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