This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I right align InDels ?

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

for example

 ACGT**A**AAACGTA  -->  ACGTAAA**A**CGTA
alignment
• 2,800 views • link
updated by Biostar Biostar • written by Chadi Saad Chadi Saad

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

• 0 views • link
written by Chadi Saad Chadi Saad

...but should you even trust indels that are called around homopolymers? From my experience, these vanish from the VCF if you raise the MAPQ threshold to 50 or 60.

• 0 views • link
written by Kevin Blighe Kevin Blighe

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 views • link
written by Rohit Rohit

0 answers

No answers yet.

Log in to answer this question.

More posts like this