vcf-query:
通过 vcf-query 提取DP (reads depth).
~/zengs/Tools/vcftools/perl/vcf-query -f ‘%CHROM\t%POS\t%REF\t%ALT\t%INFO/DP\n‘ my.vcf | awk ‘{if($1==1) print $1"\t"$2"\t"$3"\t"$4"\t"$5 }‘ > my.vcf.DP
原文：http://www.cnblogs.com/jluzhsai/p/6041373.html