samtools view -bS 1.sam|samtools sort|samtools mpileup -uf 1.bin.29.fa >1.bcf
1.sam:输入的sam文件
2.1.bin.29.fa:参考基因组的fa文件
查看bcf命令:
bcftools view 1.bcf
samtools view -bS 1.sam|samtools sort|samtools mpileup -uf 1.bin.29.fa >1.bcf
1.sam:输入的sam文件
2.1.bin.29.fa:参考基因组的fa文件
查看bcf命令:
bcftools view 1.bcf