How to get consensus sequences in FASTA format from BAM
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?
• 5,274 views
•
link
updated
by
zx8754
•
written
by
phongphak.06 •
•
written
by
phongphak.06 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to make consensus sequence from aligned reads ?
written by Info.shi •
-
Getting consensus fasta sequence for a chromosome from a bam file
written by nagarsaggi •
-
samtools mpileup `error uninitialized value $l in numeric lt` during BAM to FASTA conversion
written by marongiu.luigi
-
bcftools: Wrong number of PL fields? nals=4 npl=6
written by richardschumaker126 •
-
How to do consensus calling from bam files
written by AW
-
Consensus from mpileup error
written by mks002
-
Error while trying to get consensus Fastq from BAM
written by shreyasibiswas88 •
-
Generating wrong consensus sequence using samtools + bcftools
written by zhuhenan0514 •
-
Problem about Getting Consensus Sequence Of Mapping Against Reference Genome Using Samtools Mpileup
written by wangyang703092
-
Samtools Mpileup Defined Regions
written by Vova Naumov
Thank you so much for your help.
So I'll try this :)