美文网首页
vcf文件按照家系合并,批量处理-2020-01-10

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