vcf-tools 笔记

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