How To Convert Newbler Output Or Ace To Sam Format?
2
5
Entering edit mode
13.6 years ago
Casbon ★ 3.3k

Anyone converted the output from Roche's mapper (newbler) output to SAM? There are at least two ways to go: either parsing the 454PairAlign.txt or the ACE output.

glu-genetics claims to go the first way, but I can't get it to work (and I have emailed the author but received no response).

sam • 4.9k views
ADD COMMENT
1
Entering edit mode
13.6 years ago
Michael 54k

I googled this, but haven't tried it. But the amos tool looks like what you need starting with ace files.

ADD COMMENT
0
Entering edit mode

Unfortunately, with the latest release the '-s' option just gives:

/tmp/AMOS/src/Bank/bank2contig: invalid option -- 's'

Not sure if you need to compile it with samtools.

ADD REPLY
1
Entering edit mode
13.6 years ago
Casbon ★ 3.3k

I couldn't get AMOS to work since their is no SAM output option. I did get glu-genetics to work by checking out the latest code and modifying setup.py not to build the C extensions and fixing the monkey patch 1 2. You can then use

glu seq.Newbler2SAM

It seems to work ok, but the calculation of mapping qualities is definitely not optimal, being either 1 or 60.

ADD COMMENT

Login before adding your answer.

Traffic: 2979 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6