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

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)