How to concatenate scaffolds in a plant mitochondrial genome assembly
1
0
Entering edit mode
5.9 years ago

Dear All,

I have an Illumina PE data (2X 2,3GB, reads of 150bp) and I am trying to assemble a mitochondrial genome of ~500 Kb expected size. I used ABySS for contig generation and SSPACE for scaffolding. Finally, I got 65 scaffolds and i used MAUVE to order them against a reference genome.

My question is : How to concatenate scaffolds in order to obtain the complete genome sequence? Do I need to look for overlapping between them or just to merge scaffolds into a single fasta file ?

I tried CAP3 and BLAST for overlapping, but i did't find any overlaps.

Thank you in advance!

assembly mitochondrial genome • 1.7k views
ADD COMMENT
0
Entering edit mode

Hi, you must have solved the problem? Can you share your experience here? Thats would be very helpful!

Recently, I am trying assembling a plant mitochondrial genome. To date, I already got ~17 'mito' contigs, ~650kb and hope to merge them further if possible!

Cheers!

ADD REPLY
0
Entering edit mode

Please open a new question, describing in depth what you've done so far and what the problems are. The OP of this question has not been logged in in a year, you won't get an answer from them.

ADD REPLY
1
Entering edit mode
5.9 years ago
lagartija ▴ 160

Seams difficult if you don't have overlapping sequences. Maybe you could try scaffolding your assemblage with your PE library with NpScarf to see if it's more complete than with SSPACE;

ADD COMMENT

Login before adding your answer.

Traffic: 2807 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