vcf文件按照家系合并,批量处理-2020-01-10
作者:
爬山小虎 | 来源:发表于
2020-01-10 12:41 被阅读0次cat config | while read id
do
arr=(${id})
F=${arr[0]}
M=${arr[1]}
YS=${arr[2]}
CF=${arr[3]}
echo $F $M $Prob $YS $CF
bgzip -c $F > $F.gz
tabix -f $F.gz
bgzip -c $M > $M.gz
tabix -f $M.gz
bgzip -c $YS > $YS.gz
tabix -f $YS.gz
bgzip -c $CF > $CF.gz
tabix -f $CF.gz
/home/omics2017/hx_du/program/bcftools-1.9/bcftools merge -O v $F.gz $M.gz $YS.gz $CF.gz > ${M%%.*}.dedup.merge.vcf
done
本文标题:vcf文件按照家系合并,批量处理-2020-01-10
本文链接:https://www.haomeiwen.com/subject/wvepactx.html
网友评论