vcf-query:
通过 vcf-query 提取DP (reads depth).
~/zengs/Tools/vcftools/perl/vcf-query -f '%CHROM %POS %REF %ALT %INFO/DP
' my.vcf | awk '{if($1==1) print $1" "$2" "$3" "$4" "$5 }' > my.vcf.DP
vcf-query:
通过 vcf-query 提取DP (reads depth).
~/zengs/Tools/vcftools/perl/vcf-query -f '%CHROM %POS %REF %ALT %INFO/DP
' my.vcf | awk '{if($1==1) print $1" "$2" "$3" "$4" "$5 }' > my.vcf.DP