Answers.Body: <p><a href="https://github.com/lh3/miniasm" rel="nofollow noreferrer">Miniasm</a> will probably work for this; it's designed to be an ultra-fast, greedy assembler, which has been used for people doing real-time assemblies from MinION data. Because of its greedy nature, it can misassemble regions where the overlapping fragments are ambiguous, so it's worth using some other method in addition to miniasm to confirm that the assembly is correct.</p>
  <p>By default it looks for overlap regions that have a coverage of 3 reads when joining (option <code>-c</code>), but I assume that could be reduced down to 2 reads, as in your example (which has a minimum overlapping coverage of 2).</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2020-10-19T04:25:06.710Z
Answers.Text: Thanks for the suggestion. This actually wound up working in a roundabout way!
Answers.UserDisplayName: RD Ward
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2019-12-17T01:38:23.023Z
Answers.LastActivityDate: 2020-10-19T09:35:30.610Z
Answers.LastEditDate: 2020-10-19T09:35:30.610Z
Answers.UserDisplayName: gringer
Body: <p>I wrote a &quot;script&quot; to pull out reads from a huge fastq file in an iterative manner, by finding homology to the previous sequence.</p>
  <p>It should be relatively easy to overlap them and assemble the gene, however traditional assemblers like SPAdes fail due to <em>obviously</em> low coverage.</p>
  <p>I am hesitant to just feed it more faked reads, when I know I have all the data necessary to just overlap them into a single sequence.</p>
  <p>The reads are set up like this:</p>
  <pre><code>GATCAAACATC
       AACATCAGTTAG
           TCAGTTAG
                TAGAGGATAGC
  </code></pre>
  <p>What kind of assembler is out there where I can feed it this kind of terribly-formatted data?</p>
  
Tags: <phylogenetics><assembly>
Title: Finding simple sequence from reads with significant overlap
ThreadText: # Finding simple sequence from reads with significant overlap
  
  (Asked by: username_0 on 2019-12-15)
  
  I wrote a "script" to pull out reads from a huge fastq file in an iterative manner, by finding homology to the previous sequence.
  
  It should be relatively easy to overlap them and assemble the gene, however traditional assemblers like SPAdes fail due to *obviously* low coverage.
  
  I am hesitant to just feed it more faked reads, when I know I have all the data necessary to just overlap them into a single sequence.
  
  The reads are set up like this:
  
  ```
  GATCAAACATC
       AACATCAGTTAG
           TCAGTTAG
                TAGAGGATAGC
  
  ```
  
  What kind of assembler is out there where I can feed it this kind of terribly-formatted data?
  
  **Tags:** phylogenetics, assembly
  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_1 on 2019-12-17. Score: 1 | Accepted answer
  
  [Miniasm](https://github.com/lh3/miniasm) will probably work for this; it's designed to be an ultra-fast, greedy assembler, which has been used for people doing real-time assemblies from MinION data. Because of its greedy nature, it can misassemble regions where the overlapping fragments are ambiguous, so it's worth using some other method in addition to miniasm to confirm that the assembly is correct.
  
  By default it looks for overlap regions that have a coverage of 3 reads when joining (option `-c`), but I assume that could be reduced down to 2 reads, as in your example (which has a minimum overlapping coverage of 2).
  
  > **Comments:**
  > **username_0**: Thanks for the suggestion. This actually wound up working in a roundabout way!  
  
  ---
  
  