# trim seqs

**URL:** <https://forum.mothur.org/t/trim-seqs/1620>\
**Category:** Commands in mothur\
**Created:** [December 17, 2013, 1:43pm UTC](https://forum.mothur.org/t/trim-seqs/1620 "2013-12-17T13:43:17Z")\
**Posts on this page:** 6\
**Page:** 1

<div class="post-metadata">

**Author:** ![Bjcampbell](https://avatars.discourse-cdn.com/v4/letter/b/838e76/32.png) [@Bjcampbell](https://forum.mothur.org/u/Bjcampbell)\
**Post date:** [December 17, 2013, 1:43pm UTC](https://forum.mothur.org/t/trim-seqs/1620/1 "2013-12-17T13:43:17Z")

</div>

I have MiSeq data that was successfully grouped with the make.contigs file. However, I want to use trim.seqs to get rid of the primers, but it keeps putting most of the sequences in the scrap file. I know this is because the primers are not at the start of the sequence. Is there a way I can use trim sequences to search for the primer and get rid of all the sequences prior to that? Or, do I have to put in all the barcode information (I used the double barcodes as in the MiSeq paper by Kozich et al, 2013). If so, what is the best way of doing that that includes the barcodes, spacer and primers? It is not specifically listed in the trim.seqs command information. Thanks!

---

<div class="post-metadata">

**Author:** ![pschloss](https://yyz2.discourse-cdn.com/flex036/user_avatar/forum.mothur.org/pschloss/32/4_2.png) [@pschloss](https://forum.mothur.org/u/pschloss)\
**Post date:** [December 17, 2013, 3:47pm UTC](https://forum.mothur.org/t/trim-seqs/1620/2 "2013-12-17T15:47:38Z")

</div>

If you followed the approach we used in Kozich, there won’t be any barcodes or primers on the sequences and so trim.seqs won’t be necessary. If you aren’t following the MiSeq wet lab SOP that we posted, can you give some more details of what you’ve done?

---

<div class="post-metadata">

**Author:** ![Bjcampbell](https://avatars.discourse-cdn.com/v4/letter/b/838e76/32.png) [@Bjcampbell](https://forum.mothur.org/u/Bjcampbell)\
**Post date:** [December 17, 2013, 7:06pm UTC](https://forum.mothur.org/t/trim-seqs/1620/3 "2013-12-17T19:06:46Z")

</div>

I emailed you back in November and ended up getting instructions from Jim Kosich about setting up a sample sheet. However, the sequencing provider could not make the sample sheet work last time (miseq rejected it) so they demultiplexed the run afterward. So, the sequences were already split into samples when I got the data, and I then created a stability file based on the names. They have the barcode, linker and primers still in the sequences. I would like to remove them so I can work with a smaller number of sequences- right now it is at 3.5 million, with only 22 samples (out of 96). Any suggestions would be appreciated.

---

<div class="post-metadata">

**Author:** ![westcott](https://yyz2.discourse-cdn.com/flex036/user_avatar/forum.mothur.org/westcott/32/18_2.png) [@westcott](https://forum.mothur.org/u/westcott)\
**Post date:** [December 18, 2013, 2:11pm UTC](https://forum.mothur.org/t/trim-seqs/1620/4 "2013-12-18T14:11:36Z")

</div>

Could you post a sample sequence with the barcode, linker and primer?

---

<div class="post-metadata">

**Author:** ![jayt](https://avatars.discourse-cdn.com/v4/letter/j/e5b9ba/32.png) [@jayt](https://forum.mothur.org/u/jayt)\
**Post date:** [December 30, 2013, 4:06am UTC](https://forum.mothur.org/t/trim-seqs/1620/5 "2013-12-30T04:06:22Z")

</div>

I’m also curious about this. I am following the Miseq tutorial and I am not sure what primers were used, if there are any adapters/barcodes/linkers on the sequences, and if the make.configs command does this automatically. I also cannot find the paper that was mentioned as being published in AEM. I’d like to use the tutorial to test the abilities of a cloud based tool called CloVR. It would be really nice if I could get the results for the sample test or the entire test. Please help!

---

<div class="post-metadata">

**Author:** ![pschloss](https://yyz2.discourse-cdn.com/flex036/user_avatar/forum.mothur.org/pschloss/32/4_2.png) [@pschloss](https://forum.mothur.org/u/pschloss)\
**Post date:** [January 6, 2014, 8:31pm UTC](https://forum.mothur.org/t/trim-seqs/1620/6 "2014-01-06T20:31:35Z")

</div>

> <https://www.ncbi.nlm.nih.gov/pubmed/23793624>
>
> Rapid advances in sequencing technology have changed the experimental landscape of microbial ecology. In the last 10 years, the field has moved from sequencing hundreds of 16S rRNA gene fragments per study using clone libraries to the sequencing of millions of fragments per study using next-generation sequencing technologies from 454 and Illumina. As these technologies advance, it is critical to assess the strengths, weaknesses, and overall suitability of these platforms for the interrogation of microbial communities. Here, we present an improved method for sequencing variable regions within the 16S rRNA gene using Illumina's MiSeq platform, which is currently capable of producing paired 250-nucleotide reads. We evaluated three overlapping regions of the 16S rRNA gene that vary in length (i.e., V34, V4, and V45) by resequencing a mock community and natural samples from human feces, mouse feces, and soil. By titrating the concentration of 16S rRNA gene amplicons applied to the flow cell and using a quality score-based approach to correct discrepancies between reads used to construct contigs, we were able to reduce error rates by as much as two orders of magnitude. Finally, we reprocessed samples from a previous study to demonstrate that large numbers of samples could be multiplexed and sequenced in parallel with shotgun metagenomes. These analyses demonstrate that our approach can provide data that are at least as good as that generated by the 454 platform while providing considerably higher sequencing coverage for a fraction of the cost.
