How to get consensus sequences in FASTA format from BAM
0
0
Entering edit mode
2.7 years ago
phongphak.06 ▴ 10

I've searched for an answer but I only got this command:

samtools mpileup -uf reference.fa aligment.bam | bcftools view -cg - |
vcfutils.pl vcf2fq


Anyway, I cannot use it to generate consensus sequences since the changing of the version of each tool, which mean the options are also changed.

Is there any idea that can to this?

alignment sequence genome • 2.1k views
0
Entering edit mode
0
Entering edit mode

Thank you so much for your help.

So I'll try this :)