Question: how do I get a wonderful genome draft? : gap closed step
gravatar for vmicrobio
23 months ago by
vmicrobio240 wrote:

Hi all!

I'm trying to get a complete draft genome assembly of a bacteria from NGS Proton data. Here are my first steps : read cleaning using prinseq, alignment against hg19 (with bowtie2) in order to remove aligned reads, de novo assembly using Mira and scaffold building using abacas on a close reference genome. I got some gaps (NNNN sequences) on my consensus sequence. My next steps would be gap closed, iterative mapping and then genome annotation. I'm struggling with gaps, do you have any tips about removing gaps of my consensus sequences? I tried to use GMcloser, gapfiller and IMAGE but they are not especially designed for Proton technology... Any feedback on these softwares or other I missed?

Thank you very much!

ADD COMMENTlink written 23 months ago by vmicrobio240

Have you considered the possibility that the "missing" sequence is missing from your data (not sequenced or got removed in the post-processing you did) or are perhaps absent in your particular strain?

ADD REPLYlink modified 23 months ago • written 23 months ago by genomax67k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1018 users visited in the last hour