Answers.Body: <p>Looking that the link you provided only 30% of the reads appear to be e. coli.  I would suggest filtering the reads for those that align to e. coli and see where that gets you.  Did the original publication describe the method they used for data cleaning? </p>
  
  <p>A: </p>
  
  <pre><code>https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?run=SRR3385327
  
  Unidentified reads: 68.82%
  Identified reads: 31.18%
  cellular organisms: 30%
  Bacteria: 30%
  Proteobacteria: 28.09%
  Gammaproteobacteria: 19.7%
  Enterobacterales: 13.3%
  Enterobacteriaceae: 12.54%
  Vibrionales: 0.01%
  Viruses: 1.17%
  </code></pre>
  
  <p>B: </p>
  
  <p>Without a good assembly canu won't be able to determine if the assembled genome is circular.  I would ignore this until you have a good example.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-12-02T22:16:57.903Z
Answers.Text: Sadly the original publication only seems to say "assembled by miniasm/minimap" (along those lines) and I'm trying to figure out the pipeline as we speak -- there's no indication on whether the reads were error corrected, etc.
Answers.UserDisplayName: thestatnoob
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-12-02T22:26:42.513Z
Answers.Text: Also, how would I filter the E. coli-identified reads?
Answers.UserDisplayName: thestatnoob
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-12-03T02:24:45.970Z
Answers.Text: I would align the reads to the plasmid reference and see how many align and what the coverage of the reference looks like.  Are there any regions of the reference that are uncovered.  how even is the coverage across the reference.  These are things that can help you understand what to expect from the assembly.
Answers.UserDisplayName: Bioathlete
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-12-03T05:40:53.023Z
Answers.Text: It seems like it's actually the opposite that @thestatnoob needs to do: exclude *E. coli* reads, instead of selecting for them
Answers.UserDisplayName: gringer
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-12-03T06:35:33.307Z
Answers.Text: Correct he has the plasmid reference he is trying to assemble.  He needs to select for those reads assuming the alignment coverage shows an assembly is even worthwhile.
Answers.UserDisplayName: Bioathlete
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-12-02T22:09:14.673Z
Answers.LastActivityDate: 2017-12-02T22:09:14.673Z
Answers.UserDisplayName: Bioathlete
Answers.Body: <p>The <em>E. coli</em> genome is closer to 4.6 Mb in size; if there are a large quantity of host reads in there, then they'll take over the reads and there won't be enough for your target plasmid. By default, it will only try to assemble from a maximum of ~40X coverage of the specified genome size (higher values do not add any more information for assembly).</p>
  
  <p>It's possible that miniasm will take all reads into account, rather than just a portion of reads. You may need to filter out the <em>E. coli</em> reads for this assembly to work properly.</p>
  
  <p>For what it's worth, Canu's default parameters <em>do</em> include a correction step.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-12-03T00:46:53.553Z
Answers.Text: Yep, but the reads aren't for an E. coli genome. It's for a plasmid within E. coli, which is definitely within the 170-180k region (in fact, I got the number from the ENA...). I'm basically just trying to reconstruct the consensus plasmid that the authors have, using canu instead of miniasm.
Answers.UserDisplayName: thestatnoob
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-12-03T00:37:18.400Z
Answers.LastActivityDate: 2017-12-03T05:39:02.240Z
Answers.LastEditDate: 2017-12-03T05:39:02.240Z
Answers.UserDisplayName: gringer
Body: <p>I've downloaded reads from <a href="https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?run=SRR3385327" rel="noreferrer">this BioProject</a>. Using canu with default parameters (no correction), I've got 4 contigs, none of which really look like the reference plasmid <a href="https://www.ebi.ac.uk/ena/data/view/KX156772" rel="noreferrer">here</a>.</p>
  
  <p>The command I used was:</p>
  
  <pre><code>canu -p ip40a -d ip40a_assembly -useGrid=false -maxThreads=6 -genomeSize=175k -nanopore-raw reads.fastq
  </code></pre>
  
  <p>Granted, the original publication for this sequence used minimap/miniasm but my questions are:</p>
  
  <p>A. Why do I get multiple contigs? Doesn't canu generate a single assembled sequence?</p>
  
  <p>B. This is a plasmid but canu says <code>suggestedCircular=false</code>. What's going on here?</p>
  
Tags: <nanopore><fastq><canu>
Title: Canu assembly not making a single consensus?
ThreadText: # Canu assembly not making a single consensus?
  
  (Asked by: username_0 on 2017-12-02)
  
  I've downloaded reads from [this BioProject](https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?run=SRR3385327). Using canu with default parameters (no correction), I've got 4 contigs, none of which really look like the reference plasmid [here](https://www.ebi.ac.uk/ena/data/view/KX156772).
  
  The command I used was:
  
  ```
  canu -p ip40a -d ip40a_assembly -useGrid=false -maxThreads=6 -genomeSize=175k -nanopore-raw reads.fastq
  
  ```
  
  Granted, the original publication for this sequence used minimap/miniasm but my questions are:
  
  A. Why do I get multiple contigs? Doesn't canu generate a single assembled sequence?
  
  B. This is a plasmid but canu says `suggestedCircular=false`. What's going on here?
  
  **Tags:** nanopore, fastq, canu
  
  
  ---
  
  ## 2 answers
  
  ### Answer by username_1 on 2017-12-02. Score: 3
  
  Looking that the link you provided only 30% of the reads appear to be e. coli. I would suggest filtering the reads for those that align to e. coli and see where that gets you. Did the original publication describe the method they used for data cleaning?
  
  A:
  
  ```
  https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?run=SRR3385327
  
  Unidentified reads: 68.82%
  Identified reads: 31.18%
  cellular organisms: 30%
  Bacteria: 30%
  Proteobacteria: 28.09%
  Gammaproteobacteria: 19.7%
  Enterobacterales: 13.3%
  Enterobacteriaceae: 12.54%
  Vibrionales: 0.01%
  Viruses: 1.17%
  
  ```
  
  B:
  
  Without a good assembly canu won't be able to determine if the assembled genome is circular. I would ignore this until you have a good example.
  
  > **Comments:**
  > **username_0**: Sadly the original publication only seems to say "assembled by miniasm/minimap" (along those lines) and I'm trying to figure out the pipeline as we speak -- there's no indication on whether the reads were error corrected, etc.  
  > **username_0**: Also, how would I filter the E. coli-identified reads?  
  > **username_1**: I would align the reads to the plasmid reference and see how many align and what the coverage of the reference looks like.  Are there any regions of the reference that are uncovered.  how even is the coverage across the reference.  These are things that can help you understand what to expect from the assembly.  
  > **username_2**: It seems like it's actually the opposite that @username_0 needs to do: exclude *E. coli* reads, instead of selecting for them  
  > **username_1**: Correct he has the plasmid reference he is trying to assemble.  He needs to select for those reads assuming the alignment coverage shows an assembly is even worthwhile.  
  
  ---
  
  ### Answer by username_2 on 2017-12-03. Score: 1
  
  The *E. coli* genome is closer to 4.6 Mb in size; if there are a large quantity of host reads in there, then they'll take over the reads and there won't be enough for your target plasmid. By default, it will only try to assemble from a maximum of ~40X coverage of the specified genome size (higher values do not add any more information for assembly).
  
  It's possible that miniasm will take all reads into account, rather than just a portion of reads. You may need to filter out the *E. coli* reads for this assembly to work properly.
  
  For what it's worth, Canu's default parameters *do* include a correction step.
  
  > **Comments:**
  > **username_0**: Yep, but the reads aren't for an E. coli genome. It's for a plasmid within E. coli, which is definitely within the 170-180k region (in fact, I got the number from the ENA...). I'm basically just trying to reconstruct the consensus plasmid that the authors have, using canu instead of miniasm.  
  
  ---
  
  