# TRIM.FLOWS Scraps All Sequences -

**URL:** https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583
**Category:** Commands in mothur
**Created:** [November 19, 2013, 3:47am UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583 "2013-11-19T03:47:21Z")
**Posts on this page:** 15
**Page:** 2

<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: [March 24, 2014, 7:08pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/21 "2014-03-24T19:08:57Z")

</div>

Are you using pdiffs and bdiffs?

---

<div class="post-metadata">

### Author: ![rmb154](https://avatars.discourse-cdn.com/v4/letter/r/df788c/32.png) [@rmb154](https://forum.mothur.org/u/rmb154)
#### Post date: [April 9, 2014, 8:12pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/22 "2014-04-09T20:12:29Z")

</div>

Yes I am using both. I have been following your SOP.

---

<div class="post-metadata">

### Author: ![rmb154](https://avatars.discourse-cdn.com/v4/letter/r/df788c/32.png) [@rmb154](https://forum.mothur.org/u/rmb154)
#### Post date: [April 22, 2014, 11:55am UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/23 "2014-04-22T11:55:48Z")

</div>

do you have any idea how to help me?

---

<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: [April 28, 2014, 3:42pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/24 "2014-04-28T15:42:07Z")

</div>

Are you running the command with bdiffs=1? You might try running the set.dir(debug=t) before the trim.seqs command. Mothur will output the aligned fragments and the barcode fragments as well as the number of diffs. This may help you determine the issue.

---

<div class="post-metadata">

### Author: ![rmb154](https://avatars.discourse-cdn.com/v4/letter/r/df788c/32.png) [@rmb154](https://forum.mothur.org/u/rmb154)
#### Post date: [May 26, 2014, 11:17am UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/25 "2014-05-26T11:17:21Z")

</div>

I have been running the trim.seqs command with bdiffs=1. When I run the set.dir(debug=t) then it is showing me that the numdiffs range from 4 to 15 on average! Is this normal? When I adjust the bdiffs=8 then I get a lot more of my sequences saved from the scrap.flows (18/24 of my samples). Is there an acceptable range to adjust bdiffs? If I go higher in bdiffs= in order to prevent the majority of sequences being scrapped, am I compromising the data analysis further down the pipeline?

Thanks for your 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: [May 27, 2014, 12:30pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/26 "2014-05-27T12:30:16Z")

</div>

No, that’s a problem. Have you tried using order=B in trim.flows and shhh.flows? Do you know what type of 454 sequencing was done on these?

pat

---

<div class="post-metadata">

### Author: ![rmb154](https://avatars.discourse-cdn.com/v4/letter/r/df788c/32.png) [@rmb154](https://forum.mothur.org/u/rmb154)
#### Post date: [May 27, 2014, 1:15pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/27 "2014-05-27T13:15:08Z")

</div>

I have tried order=B, it doesn’t improve the amount of sequences saved from the scrap.flows. It looks like there are also quite some that are scrapped based strictly on length but the majority are “lbp”. These were supposedly done with FLX… I am at a loss here. This is the first time this has happened to me that we have had such terrible sequencing…

---

<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 30, 2014, 12:32pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/28 "2014-05-30T12:32:06Z")

</div>

Can you post your sff file and oligos file somewhere (e.g. dropbox/google drive) for me to take a look?

---

<div class="post-metadata">

### Author: ![rmb154](https://avatars.discourse-cdn.com/v4/letter/r/df788c/32.png) [@rmb154](https://forum.mothur.org/u/rmb154)
#### Post date: [September 1, 2014, 7:43am UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/29 "2014-09-01T07:43:46Z")

</div>

I have send you the files via dropbox to the mothurbugs email. I included the URL for this forum. Please let me know if this email address is correct or incorrect. Thank you.

---

<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: [September 10, 2014, 1:22pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/30 "2014-09-10T13:22:37Z")

</div>

Sorry for the delay in getting back to you - we were on vacation at the end of the month and are slowly getting back to all our missed emails and forum posts…

It appears that you do have FLX and so you should be using flow order A. Like you, when I run the following…

```nohighlight
trim.flows(flow=sc.flow, oligos=sc.oligos.txt, pdiffs=1, bdiffs=2, processors=8)

```

I get a lot of sequences scrapped because of length:

```nohighlight
3 bf
856 f
108913 l
2157 lbf
7964 lf

```

The length gets flagged because you have noisy data in the first 450 flows that end up inserting ambiguous base calls (caused by reagents) or because you’re accidentally sequencing short reads (such as primer dimers).

However of the 378614 total flows, 258721 flows are good, so I’m not sure that I would be overly concerned.

---

<div class="post-metadata">

### Author: ![rmb154](https://avatars.discourse-cdn.com/v4/letter/r/df788c/32.png) [@rmb154](https://forum.mothur.org/u/rmb154)
#### Post date: [September 12, 2014, 6:58am UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/31 "2014-09-12T06:58:25Z")

</div>

Thank you!

---

<div class="post-metadata">

### Author: ![rmb154](https://avatars.discourse-cdn.com/v4/letter/r/df788c/32.png) [@rmb154](https://forum.mothur.org/u/rmb154)
#### Post date: [September 12, 2014, 7:10pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/32 "2014-09-12T19:10:07Z")

</div>

Ok, I have other similar data. From the same sequencing group. Same kind of analysis. When I run trim.flows it puts all except 25bytes of the 6.2GB of sequencing data into the scrap.flows. If this is a case of simply poor sequencing. Is there anything I can do to recover some information from this data?

---

<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: [September 15, 2014, 9:19pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/33 "2014-09-15T21:19:30Z")

</div>

Are you sure that it’s not FLX+?

---

<div class="post-metadata">

### Author: ![rmb154](https://avatars.discourse-cdn.com/v4/letter/r/df788c/32.png) [@rmb154](https://forum.mothur.org/u/rmb154)
#### Post date: [September 17, 2014, 7:28am UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/34 "2014-09-17T07:28:36Z")

</div>

I have tried both order=A and order=B. Order A got me 2 samples out of 48, order B got all files scrapped.

---

<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: [September 17, 2014, 3:29pm UTC](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583/35 "2014-09-17T15:29:10Z")

</div>

Can you try running trim.seqs with the fasta and qual file and see what happens?

See:  
[http://www.mothur.org/wiki/454\_SOP#Using\_quality\_scores](http://www.mothur.org/wiki/454_SOP#Using_quality_scores)

I suspect the data are bad and that perhaps it’s all primer dimer or there were reagent problems.

Pat

[Previous page](https://forum.mothur.org/t/trim-flows-scraps-all-sequences/1583.md?page=1)
