vcf-tools 筆記

vcf-query:spa

經過 vcf-query 提取DP (reads depth).blog

~/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
相關文章
相關標籤/搜索