This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Dealing With Newbler'S Misalignments?

I have been using newbler for a while now, and the rate of misalignments is very high. I have found that these are correlated with the homopolymer problem, but should be aligned correctly.

For example, this is produced in the 454PairAlign.txt (top is read, bottom is reference): 

TAA-TTTT--GTATTTTTGTA
TAATTTTTGTA-TTTTTTGTA

I have omitted the full reads for brevity.

The read has undercalled the two homopolymers. However, the best alignment suggested is clearly wrong since the GTA in the middle should be aligned, with two homopolymer gaps.

Has anyone found a way of getting newbler to correctly align such cases (eg different parameters)?

It seems like this must be a problem with the aligner, which appears to prefer three gaps plus two substitutions to two gaps (in homopolymers).

• 3,353 views • link
updated by Biostar Biostar • written by Casbon • Casbon

Roche guys please don't read further (and this is not an answer btw) : Maybe the best way to deal with the errors would be to avoid newbler if possible as an alignment tool. I tried it about 1.5 years ago when it even had even more bugs, I would try and check for an open source tool instead and possibly ignore homopolymers.

• 0 views • link
written by Michael Dondrup Michael Dondrup

not an answer but there are these mysterious parameters mentioned in the manual:- 

[-rst] "Repeat score threshold parameter"
[-sl] "Seed length"
[-ss] "Seed step" parameter

but I am not that clear how newbler works. I wish they would explain it more clearly in the manual.

• 0 views • link
updated by Ram Ram • written by Jake • Jake

Michael, I tried BWA-sw but found even more problems in the alignments. Any suggestions on a good long read aligner?

• 0 views • link
written by Casbon • Casbon

Didn't see your comment until now, tried Lastz or blat?

• 0 views • link
written by Michael Dondrup Michael Dondrup

1 answer

In conversations with Roche, they claim the 454PairAlign.txt alignments are not the final alignments used by Newbler for variant calling. This means, until Roche support a decent output format for their alignments (like SAM), you cannot get the underlying alignments from newbler.

• 0 views • link
written by Casbon • Casbon

Log in to answer this question.

More posts like this