I'm creating consensus sequences using a custom script but I am unsure if lowercase bases (the bases that are on the reverse strand) are already reverse complemented, or they need to be reverse complemented?
In IGV, when I see reads aligned to a position and I have a mix of A (forward reads), and 'a' (reverse reads) from single-end data, can I assume that 'a' is already reverse complemented where the reverse read actually had a 'T' at that position? Or will I have to account for this myself and really I have 'A' and 'T' matched at that position?
Thank you!
•
written
by
DNAngel
1 answer
All alignments in a BAM file are written for the + strand, so they've been reverse complemented as appropriate. IGV is just indicating their orientation with upper/lower case letters, you don't need to do anything else yourself in that regard.
Whether a reverse complemented read actually arose from the - strand will end up being library type dependent, in case you're wondering about that as well.
Log in to answer this question.
-
Calling variants with UMI-collapsed duplex consensus reads
written by blid11 •
-
Bcftools consensus generates mismatched consensus sequence
written by Duy •
-
SNV occuring 2 fold more in forward than in reverse reads should be filtered?
written by jeni •
-
ADR, ADF and dp4 counts do not count strands correct
written by g.t.los •
-
Merge paired reads for error-correction
written by Mick •
-
Illumina Overlapping R1/R2 reads, error-correction in python
written by Mick •
-
Confused about how to generate a consensus sequence after bwa
written by DNAngel
-
Batch pairwise alignment and consensus
written by lewis272727 •
-
Inline barcodes in the reverse reads
written by Picasa
-
Create Consensus Sequences For Sequence Pairs Within A Multiple Alignment?
written by Dave Lunt •
