bcftools:根據(jù)染色體提取vcf中的信息

根據(jù)染色體提取vcf中的信息:

bcftools view -r chr20,chr21,chr22 NA12878.vcf.gz -o NA12878.vcf.chr20_22.gz -O z
解釋:

提取金標(biāo)準(zhǔn)變異NA12878.vcf.gz中20烘浦、21和22染色體的變異刑顺,保存到NA12878.vcf.chr20_22.gz文件中
-O z:表示使用壓縮形式保存文件

提取之后的文件作為其它軟件輸入時(shí)經(jīng)常需要索引,可按如下建立索引:

bcftools index -t NA12878.vcf.chr20_22.gz
解釋:

-t:產(chǎn)生TBI格式的索引(如果沒有-t選項(xiàng)做个,默認(rèn)會(huì)產(chǎn)生csi的索引)
命令結(jié)束后會(huì)產(chǎn)生NA12878.vcf.chr20_22.gz.tbi的索引文件
————————————————
版權(quán)聲明:本文為CSDN博主「bigcindy」的原創(chuàng)文章括勺,遵循CC 4.0 BY-SA版權(quán)協(xié)議结序,轉(zhuǎn)載請(qǐng)附上原文出處鏈接及本聲明宴猾。
原文鏈接:https://blog.csdn.net/Jwenxue/article/details/109330562

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末酌住,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子刃麸,更是在濱河造成了極大的恐慌醒叁,老刑警劉巖,帶你破解...
    沈念sama閱讀 207,248評(píng)論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異把沼,居然都是意外死亡啊易,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,681評(píng)論 2 381
  • 文/潘曉璐 我一進(jìn)店門饮睬,熙熙樓的掌柜王于貴愁眉苦臉地迎上來租谈,“玉大人,你說我怎么就攤上這事捆愁「钊ィ” “怎么了?”我有些...
    開封第一講書人閱讀 153,443評(píng)論 0 344
  • 文/不壞的土叔 我叫張陵昼丑,是天一觀的道長劫拗。 經(jīng)常有香客問我,道長矾克,這世上最難降的妖魔是什么页慷? 我笑而不...
    開封第一講書人閱讀 55,475評(píng)論 1 279
  • 正文 為了忘掉前任,我火速辦了婚禮胁附,結(jié)果婚禮上酒繁,老公的妹妹穿的比我還像新娘。我一直安慰自己控妻,他們只是感情好州袒,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,458評(píng)論 5 374
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著弓候,像睡著了一般郎哭。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上菇存,一...
    開封第一講書人閱讀 49,185評(píng)論 1 284
  • 那天夸研,我揣著相機(jī)與錄音,去河邊找鬼依鸥。 笑死亥至,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的贱迟。 我是一名探鬼主播姐扮,決...
    沈念sama閱讀 38,451評(píng)論 3 401
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼衣吠!你這毒婦竟也來了茶敏?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 37,112評(píng)論 0 261
  • 序言:老撾萬榮一對(duì)情侶失蹤缚俏,失蹤者是張志新(化名)和其女友劉穎惊搏,沒想到半個(gè)月后贮乳,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,609評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡胀屿,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,083評(píng)論 2 325
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了包雀。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片宿崭。...
    茶點(diǎn)故事閱讀 38,163評(píng)論 1 334
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖才写,靈堂內(nèi)的尸體忽然破棺而出葡兑,到底是詐尸還是另有隱情,我是刑警寧澤赞草,帶...
    沈念sama閱讀 33,803評(píng)論 4 323
  • 正文 年R本政府宣布讹堤,位于F島的核電站,受9級(jí)特大地震影響厨疙,放射性物質(zhì)發(fā)生泄漏洲守。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,357評(píng)論 3 307
  • 文/蒙蒙 一沾凄、第九天 我趴在偏房一處隱蔽的房頂上張望梗醇。 院中可真熱鬧,春花似錦撒蟀、人聲如沸叙谨。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,357評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽手负。三九已至,卻和暖如春姑尺,著一層夾襖步出監(jiān)牢的瞬間竟终,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,590評(píng)論 1 261
  • 我被黑心中介騙來泰國打工切蟋, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留衡楞,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 45,636評(píng)論 2 355
  • 正文 我出身青樓敦姻,卻偏偏與公主長得像瘾境,于是被迫代替她去往敵國和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子镰惦,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,925評(píng)論 2 344

推薦閱讀更多精彩內(nèi)容