<p>Hello everyone,</p>

<p>I&#39;m working on improving the gene coverage for the chromosome of an <em>E. coli</em> strain, using existing WGS data. There are two WGS projects for this strain on NCBI (according to the assembly report, the assembly level for one of them is &quot;scaffold&quot;, and for the other is &quot;contig&quot;).</p>

<p>So far, I&#39;ve obtained the contigs from both of these projects from NCBI (which I saved as two separate multifasta files), found a decent reference genome, and re-ordered the contig sets against it using Mauve. I then ran the re-ordered contigs through GLIMMER-3, which produced annotated Genbank format files for each WGS project contig set. Basically, I&#39;ve been following the steps in this tutorial: <a href="http://www.microbialinformaticsj.com/content/3/1/2">http://www.microbialinformaticsj.com/content/3/1/2</a>.</p>

<p>After plotting the two sets of annotated contigs against each other using DNA plotter, it looks like together they can give decent gene coverage. There are gaps in the CDS in each case which are complemented by the other.</p>

<p>I have been using contigs from both WGS projects so far, but I think I&#39;m going to start using the scaffolds from the scaffold level assembly; as I&#39;ve just found out that the scaffold level project includes two plasmids (the contig level WGS project is just a chromosome), and I know which scaffolds are the plasmids, so I can exclude them.</p>

<p>How can I go about assembling a single chromosome using data from both WGS projects? Is this possible? If not, which programmes can highlight the differences between the assembled chromosomes?</p>

<p>Should I be using the scaffolds instead of the contigs for the WGS project which is at a scaffold level?</p>

<p>Are there any programmes for closing gaps between/in scaffolds using additional contigs?</p>

<p>Thanks for your time,</p>

<p>Ronan</p>

<p>If the raw sequencing reads were also released, you could try <a href="https://github.com/vice87/gam-ngs">GAM-NGS</a>. Another option is to use MIRA to perform a genome-guided assembly, use the scaffolds as reference, and the two sets of contigs as two additional strains to be assembled. A third option is to use CAP3 to merge the assemblies, but I remember reading (from more than one source, though I do not recall a single one) this is not recommended, as it introduces mis-assemblies.</p>

<p><strong>edit</strong>: I saw <a href="http://www.biostars.org/p/97652/">this</a> post listed as similar to yours, seems like a good alternative to try.</p>
