Biostar Beta. Not for public use.
How can I right align InDels ?
0
Entering edit mode
23 months ago
France

How can I align all indels to the right in Bam file ?

for example

 ACGT**A**AAACGTA  -->  ACGTAAA**A**CGTA

alignment • 813 views
1
Entering edit mode

Is there a particular reason that last of the A in the homopolymer is pushed in the middle of the indel and the indels are not clustered together?

I might be naive, wouldn't the moving of the indels remove all the sense from the alignment. I wouldn't change the alignments without looking into the reference, so why not use a mapper that is homopolymer aware or with that particular error profile?

0
Entering edit mode

I need it in order to count InDels afters homopolymers (since they are aligned to the left by aligners)