This is a test version of Biostars. For the public version, visit https://www.biostars.org.
spliced alignment algorithm

The picture is a view of an rna-seq bam file from tophat alignment. On the left side, in red, are some alignments shifted away from the leftmost exon. These "shifted-exons" are sometimes 1 base long, a base that can be found in many other locations including the real exon.
Why would the aligner pick that specific shifted location?
Image at: https://ibb.co/gxv9Cv

spliced alignment

rna-seq alignment splice

Perhaps someone can comment on this possible explanation: Might this be a result of a 2-pass alignment.The first pass determined a "shifted exon" based on one read that was mapped unambiguously with 5-6 bases to that region. The second pass then mapped ambiguous 1-base exons to that location. .

0 answers

No answers yet.

Log in to answer this question.