# Using Sffinfo to parse by MID

**URL:** <https://forum.mothur.org/t/using-sffinfo-to-parse-by-mid/1819>\
**Category:** Commands in mothur\
**Created:** [April 29, 2014, 10:51pm UTC](https://forum.mothur.org/t/using-sffinfo-to-parse-by-mid/1819 "2014-04-29T22:51:39Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![Amblyomma](https://avatars.discourse-cdn.com/v4/letter/a/278dde/32.png) [@Amblyomma](https://forum.mothur.org/u/Amblyomma)\
**Post date:** [April 29, 2014, 10:51pm UTC](https://forum.mothur.org/t/using-sffinfo-to-parse-by-mid/1819/1 "2014-04-29T22:51:39Z")

</div>

I’m prepping some SFF files for upload to the SRA. Or at least I’m trying to. Files that I’ve manipulated dozens of times in mothur are suddenly giving me problems. I’ve tried this in both versions 1.31.2 and 1.33.3 on my Ubuntu 12.04 64 bit system.

sffinfo(sff=myfile.sff, fasta=F, qfile=F, trim=F, flow=F, oligos=122010.oligos)

On both versions mothur generates the individual temp sff files for each MID, but at the end of the command everything gets dumped in the scrap file. I’ve tried multiple pairs of sff and oligos files, all of which have been previously analyzed, and nothing works. I’ve grabbed copies from other drives in case some how the oligos files got corrupted. Same problem. I also tried running everything in Windows. Same problems.

Also, oddly, I’ve always named my oligos files something like 12-2010.oligos, and suddenly versions of mothur I’ve used these files with before are refusing to accept the hyphen in the name. It’s an easy enough change, but I can’t figure out why now it’s a problem.

Please help. This is already a horrendously tedious process, and I haven’t had this much trouble operating mothur is years.

---

<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:** [May 5, 2014, 4:44pm UTC](https://forum.mothur.org/t/using-sffinfo-to-parse-by-mid/1819/2 "2014-05-05T16:44:30Z")

</div>

Have you tried with pdiffs and bdiffs?

---

<div class="post-metadata">

**Author:** ![Amblyomma](https://avatars.discourse-cdn.com/v4/letter/a/278dde/32.png) [@Amblyomma](https://forum.mothur.org/u/Amblyomma)\
**Post date:** [May 5, 2014, 7:55pm UTC](https://forum.mothur.org/t/using-sffinfo-to-parse-by-mid/1819/3 "2014-05-05T19:55:51Z")

</div>

Yes, tried the standard settings from the SOP for trim.seqs. No luck.

---

<div class="post-metadata">

**Author:** ![Amblyomma](https://avatars.discourse-cdn.com/v4/letter/a/278dde/32.png) [@Amblyomma](https://forum.mothur.org/u/Amblyomma)\
**Post date:** [May 5, 2014, 9:09pm UTC](https://forum.mothur.org/t/using-sffinfo-to-parse-by-mid/1819/4 "2014-05-05T21:09:14Z")

</div>

It’s looking like we can chalk this one up to “unexpected behaviors” (it was unexpected to me anyway).

1. This produces the parsed sff files plus pooled fasta, qual, flow, and scrap files:  
sffinfo(sff=myfile.sff, oligos=122010A.oligos, bdiffs=1, pdiffs=2)

2. This puts everything in the scrap file:  
sffinfo(sff=myfile.sff, fasta=F, qfile=F, trim=F, flow=F, oligos=122010A.oligos, bdiffs=1, pdiffs=2)

3. This produces unparsed raw.fasta, raw.qual, flow, and scrap files:  
sffinfo(sff=myfile.sff, trim=F, oligos=122010A.oligos, bdiffs=1, pdiffs=2)

I am now attempting to determine whether or not the sequences are being trimmed in #1 (I want raw sequences minus the MID). For the record, I’m working with v1.33.3 on Win8 today.

[Edit] Confirmed, #1 is trimming the sequences. What is being trimmed off here? The documentation says it’s being trimmed to the “clipQualLeft and clipQualRight values”, but I cannot find where those values are described.

---

<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 7, 2014, 11:32am UTC](https://forum.mothur.org/t/using-sffinfo-to-parse-by-mid/1819/5 "2014-05-07T11:32:22Z")

</div>

> [Edit] Confirmed, #1 is trimming the sequences. What is being trimmed off here? The documentation says it’s being trimmed to the “clipQualLeft and clipQualRight values”, but I cannot find where those values are described.

If you run sffinfo(sfftxt=T) on the sff files you generated and open them, you’ll see that for each sequence the clip values are set there.
