5、提取snp indel 位点

le final.snp.list | perl -lane '{$a+=1;print "$a $F[0] $F[1] $F[1]"}' | less >snp_site

le final.indel.vcf |grep -v '^#' | less -S|perl -lane '{$a+=1;$b=$F[1]+length($F[3]);print "$a $F[0] $F[1] $b"}' | less -S >indel_site

原文地址:https://www.cnblogs.com/renping/p/7356586.html