2022-06-08修改vcf文件樣本名稱，提取指定vcf樣本

1.修改vcf樣本名

使用bcftools的 rehead指令

bcftools reheader -s sample.file old.vcf -o new.sample.vcf

-s 后面的為sample.file 指定要替換的樣本名吐限，內容如下

2.提取指定樣本vcf文件

準備好樣本ID文件鲜侥，samplelist.txt，一個樣本名一行毯盈，如下

bcftools view -S?samplelist.txt? old.vcf.gz?-Ov?> new.vcf

• 人面猴

  序言：七十年代末剃毒，一起剝皮案震驚了整個濱河市，隨后出現(xiàn)的幾起案子搂赋，更是在濱河造成了極大的恐慌赘阀，老刑警劉巖，帶你破解...

  沈念sama閱讀 211,290評論 6贊 491
• 死咒

  序言：濱河連續(xù)發(fā)生了三起死亡事件脑奠，死亡現(xiàn)場離奇詭異基公，居然都是意外死亡，警方通過查閱死者的電腦和手機宋欺，發(fā)現(xiàn)死者居然都...

  沈念sama閱讀 90,107評論 2贊 385
• 救了他兩次的神仙讓他今天三更去死

  文/潘曉璐 我一進店門轰豆，熙熙樓的掌柜王于貴愁眉苦臉地迎上來，“玉大人齿诞，你說我怎么就攤上這事酸休。” “怎么了祷杈？”我有些...

  開封第一講書人閱讀 156,872評論 0贊 347
• 道士緝兇錄：失蹤的賣姜人

  文/不壞的土叔 我叫張陵斑司，是天一觀的道長。 經常有香客問我但汞，道長宿刮，這世上最難降的妖魔是什么互站？ 我笑而不...

  開封第一講書人閱讀 56,415評論 1贊 283
• ?港島之戀（遺憾婚禮）

  正文 為了忘掉前任，我火速辦了婚禮僵缺，結果婚禮上胡桃，老公的妹妹穿的比我還像新娘。我一直安慰自己磕潮，他們只是感情好翠胰，可當我...

  茶點故事閱讀 65,453評論 6贊 385
• 惡毒庶女頂嫁案：這布局不是一般人想出來的

  文/花漫 我一把揭開白布。 她就那樣靜靜地躺著揉抵，像睡著了一般亡容。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上冤今，一...

  開封第一講書人閱讀 49,784評論 1贊 290
• 城市分裂傳說

  那天闺兢，我揣著相機與錄音，去河邊找鬼戏罢。 笑死屋谭，一個胖子當著我的面吹牛，可吹牛的內容都是我干的龟糕。 我是一名探鬼主播桐磁，決...

  沈念sama閱讀 38,927評論 3贊 406
• 雙鴛鴦連環(huán)套：你想象不到人心有多黑

  文/蒼蘭香墨 我猛地睜開眼，長吁一口氣：“原來是場噩夢啊……” “哼讲岁！你這毒婦竟也來了我擂？” 一聲冷哼從身側響起，我...

  開封第一講書人閱讀 37,691評論 0贊 266
• 萬榮殺人案實錄

  序言：老撾萬榮一對情侶失蹤缓艳，失蹤者是張志新（化名）和其女友劉穎校摩，沒想到半個月后，有當地人在樹林里發(fā)現(xiàn)了一具尸體阶淘，經...

  沈念sama閱讀 44,137評論 1贊 303
• ?護林員之死

  正文 獨居荒郊野嶺守林人離奇死亡衙吩，尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...

  茶點故事閱讀 36,472評論 2贊 326
• ?白月光啟示錄

  正文 我和宋清朗相戀三年，在試婚紗的時候發(fā)現(xiàn)自己被綠了溪窒。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片坤塞。...

  茶點故事閱讀 38,622評論 1贊 340
• 活死人

  序言：一個原本活蹦亂跳的男人離奇死亡，死狀恐怖澈蚌，靈堂內的尸體忽然破棺而出摹芙，到底是詐尸還是另有隱情，我是刑警寧澤宛瞄，帶...

  沈念sama閱讀 34,289評論 4贊 329
• ?日本核電站爆炸內幕

  正文 年R本政府宣布瘫辩，位于F島的核電站，受9級特大地震影響坛悉，放射性物質發(fā)生泄漏伐厌。R本人自食惡果不足惜，卻給世界環(huán)境...

  茶點故事閱讀 39,887評論 3贊 312
• 男人毒藥：我在死后第九天來索命

  文/蒙蒙 一裸影、第九天 我趴在偏房一處隱蔽的房頂上張望挣轨。 院中可真熱鬧，春花似錦轩猩、人聲如沸卷扮。這莊子的主人今日做“春日...

  開封第一講書人閱讀 30,741評論 0贊 21
• 一樁弒父案均践，背后竟有這般陰謀

  文/蒼蘭香墨 我抬頭看了看天上的太陽晤锹。三九已至，卻和暖如春彤委，著一層夾襖步出監(jiān)牢的瞬間鞭铆，已是汗流浹背。 一陣腳步聲響...

  開封第一講書人閱讀 31,977評論 1贊 265
• 情欲美人皮

  我被黑心中介騙來泰國打工焦影， 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留车遂，地道東北人。 一個月前我還...

  沈念sama閱讀 46,316評論 2贊 360
• 代替公主和親

  正文 我出身青樓斯辰，卻偏偏與公主長得像舶担，于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子彬呻，可洞房花燭夜當晚...

  茶點故事閱讀 43,490評論 2贊 348

推薦閱讀更多精彩內容

• 使用bcftools view 拆分vcf格式文件

  首先看一下bcftools view用法 我們重點關注Subset 和output選項 我們僅使用-o -O -...

  宗肅書閱讀 6,896評論 2贊 5
• bcftools常用命令詳解

  http://samtools.github.io/bcftools/bcftools.html https://...

  生信師姐閱讀 42,369評論 0贊 39
• 如何從vcf文件中批量提取一系列基因的SNP位點？

  需求 客戶的一個簡單需求： 我有一批功能基因位點闸氮，想從重測序的群體材料中找到這些位點剪况，如何批量快速獲得？ 示例文件...

  生物信息與育種閱讀 17,310評論 1贊 19
• picard:處理HaplotypeCaller斷點的問題--gatk提取指定染色體vcf文件

  使用GATK4 call snp時總會出錯湖苞。鑒于GATK4新版的不穩(wěn)定性拯欧，即使有坑也得跳。這里是填坑的地方财骨。 重要...

  wo_monic閱讀 3,601評論 3贊 8
• bcftools 處理vcf文件，尋找多個vcf文件中突變的交集

  目標 通過varscan等找到somatic mutations之后隆箩，想看多個樣本的突變之間的交集该贾。bcftool...

  果蠅的小翅膀閱讀 9,184評論 0贊 18