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

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.6k views
1
Entering edit mode
12.7 years ago

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

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.

1
Entering edit mode
12.7 years ago
Casbon ★ 3.2k

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.