# Trim.seqs without .qual files

**URL:** https://forum.mothur.org/t/trim-seqs-without-qual-files/888
**Category:** Commands in mothur
**Created:** [May 7, 2012, 3:29pm UTC](https://forum.mothur.org/t/trim-seqs-without-qual-files/888 "2012-05-07T15:29:38Z")
**Posts on this page:** 4
**Page:** 1

<div class="post-metadata">

### Author: ![dbaptista](https://avatars.discourse-cdn.com/v4/letter/d/c68b51/32.png) [@dbaptista](https://forum.mothur.org/u/dbaptista)
#### Post date: [May 7, 2012, 3:29pm UTC](https://forum.mothur.org/t/trim-seqs-without-qual-files/888/1 "2012-05-07T15:29:38Z")

</div>

Hi All,

I have a bit of a problem. So I have some old data that never received the .qual files. I followed the Schloss SOP recommending to use the keepfirst command to trim my files to 200 bases. Unfortunately, when I do this, I loose all my samples when I filter using “trump=.”. When I increase the the keepfirst to 300 bases, I’m able to trim and keep the majority of the sequences in my alignment. The Schloss protocol mentions that sequence quality drops after 250 bases which is why they recommend the 200 cutoff. I’m wondering if anyone has any suggestions. Is a 300 base cutoff acceptable in this situation? Anything else I can do to get around this problem? Any advise will help. I did not try a 250 base cutoff, but I could. Thanks.

-Damon

---

<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: [May 8, 2012, 12:32pm UTC](https://forum.mothur.org/t/trim-seqs-without-qual-files/888/2 "2012-05-08T12:32:44Z")

</div>

Damon,

Are you running screen.seqs prior to running filter.seqs?

Pat

---

<div class="post-metadata">

### Author: ![dbaptista](https://avatars.discourse-cdn.com/v4/letter/d/c68b51/32.png) [@dbaptista](https://forum.mothur.org/u/dbaptista)
#### Post date: [May 8, 2012, 8:53pm UTC](https://forum.mothur.org/t/trim-seqs-without-qual-files/888/3 "2012-05-08T20:53:20Z")

</div>

Hi Pat,

Yes, I run screen.seqs prior. I used the recommendations in your SOP screen.seqs(optimize=start, criteria=95, fasta=current, name=current). I used an “end” value just below my 2.5%-tile. It was strange since it seemed like the summary of the filtered alignment had a negative value for the sequence length.

These are archaea samples and I used the silva.archaea.fasta file as a reference. I haven’t tried the trim with my data set that is for bacteria yet. Maybe that will trim and align better? Any suggestions?

Thanks,

Damon

---

<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: [June 4, 2012, 6:05pm UTC](https://forum.mothur.org/t/trim-seqs-without-qual-files/888/4 "2012-06-04T18:05:40Z")

</div>

Could you send your filtered alignment file to [mothur.bugs@gmail.com](mailto:mothur.bugs@gmail.com), so I can track down the negative values issue for you?
