FLX+ SOP

Hi,

I have a recent sff file generated through V1-4 FLX+ amplicon chemistry from my sequence provider.

Now if I fed it into the initial SOP commands:
sffinfo(sff=sample.sff, flow=T)
trim.flows(flow=sample.flow, oligos=sample.oligos, pdiffs=2, bdiffs=1, minflows=1050, maxflows=1050, order=B, processors=5)
shhh.flows(file=sample.flow.files) <— approximately 11-hours run (lookup file: GS FLX Titanium)
trim.seqs(fasta=sample.shhh.fasta, name=sample.shhh.names, oligo=sample.oligos, pdiffs=2, bdiffs=1, maxhomop=8, minlength=200, flip=T)
summary.seqs(fasta=sample.shhh.trim.fasta, name=sample.shhh.trim.names)

Plus if I try to proceed further into summary.seqs() here, I get an error prompt mentioning blank input files.

So this left me no choice but to find an alternative quality score approach:
trim.seqs(fasta=sample.fasta, oligos=sample.oligos, qfile=sample.qual, maxambig=0, maxhomop=8, flip=T, bdiffs=1, pdiffs=2, qwindowaverage=35, qwindowsize=50) <— “no namefile provided, fasta closely match with name=sample.shhh.trim.names”
summary.seqs(fasta=sample.trim.fasta, name=sample.trim.names)
unique.seqs(fasta=sample.trim.fasta, name=sample.trim.names)
summary.seqs(fasta=sample.trim.unique.fasta, name=sample.trim.unique.names)
align.seqs(fasta=sample.trim.unique.fasta, reference=silva.bacteria.fasta)
summary.seqs(fasta=sample.trim.unique.align, name=sample.trim.unique.names)
screen.seqs(fasta=sample.trim.unique.align, name=sample.shhh.trim.unique.names, group=sample.groups, end=27659, optimize=start, criteria=95)
summary.seqs(fasta=sample.trim.unique.good.align, name=sample.trim.unique.good.names) <— “blank input files, aborting…command not complete”

I just can’t find any solutions here onward (literally buried in a pile of technical problems I know I feel completely unaware of).

Anybody around for a desperate help, please?? :cry: :cry: :cry:

Daniel

You must use the order=b parameter in the shhh.flows command or the trim.seqs command will dump all your sequences (blank input file error)
I had the same problem and it work now !

Thanks, appreciate your help! Worked out all well!