I went with calling Tracy from python scripts, thank you!
Hello Fellow Bioinformaticians,
I'm trying to find a tool to assemble .ab1 trace files (ABI sequencer data files), ideally in Python. However, it looks like Biopython doesn't support DNA assembly. I was wondering if you have any recommendations or suggestions?
Thank you very much!
3 answers
I'd be surprised if such a tool existed that assembles directly from ab1 files. You may be able to convert .ab1 files into .fastq format as described here. Then you can assemble the converted reads using MIRA, SPAdes, phrap or some other assembler that can use Sanger reads. You may also want to check out phred which is part of the phrap suite, as it also can read and convert .ab1 files.
Indeed, phred/phrap/consed is the standard toolkit for Sanger assembly. One can still get them, but they existed before python was even a thing.
Depending on the size of the sanger dataset and if you are fine with commercial software, you might want to consider packages like Genious.
Working with MIRA or even phrap is fun but probably very time-consuming nowadays, especially if you are not familiar with assembly file formats. Staden Package or Consed may be used as alignment editors.
I have been working with these packages years ago ..
Log in to answer this question.
If you have Mac you can use BioLabDonkey program to assemble the reads from .ab1 files directly.
About assembly of .ab1 reads in BioLabDonkey