Generate consensus from a BAM file to each of the variants
0
0
Entering edit mode
3 months ago
fhsantanna ▴ 550

I found a putative SARS-CoV-2 co-infection case. After examining the bam file, I found that some positions present two variants. Is it possible to build two consensus, one to the major (most frequent) and one to the minor form (less frequent) of the viral sample?

bam coinfection consensus ivar • 278 views
2
Entering edit mode

How about using GATK FastaAlternateReferenceMaker (LINK).

0
Entering edit mode

Seems promising, I'll give you a feedback.

0
Entering edit mode

Dear GenoMax , let me see if I understood. I managed to generate the alternative reference, now do I have to map the reads of my sample to this reference in order to generate a new consensus?

Since the alternative reference was generated over the mn908947 (wuhan), there is no gaps in the sequence, which are originally found in my sequence of interest.

Or should I map the reads of my sequence of interest to its own consensus, then generate the alternative sequence using FastaAlternateReferenceMaker?