mothur

Filter.seqs (V3-V4 data)


#1

Hi all,

I’m having a difficult time trying to figure out what has gone wrong with my filter.seqs. I get filtered alignment of 30.
Here is my code:

make.file(inputdir=., type=fastq, prefix=suis)
make.contigs(file=/scratch/shill09/files/suis.files,processors=32)
summary.seqs(fasta=suis.trim.contigs.fasta)
screen.seqs(fasta=suis.trim.contigs.fasta, group=suis.contigs.groups, maxambig=0, maxlength=500)
unique.seqs(fasta=suis.trim.contigs.good.fasta)
count.seqs(name=suis.trim.contigs.good.names, group=suis.contigs.good.groups)
summary.seqs(count=suis.trim.contigs.good.count_table)
pcr.seqs(fasta=silva.bacteria2.fasta, start=6388, end=25316, keepdots=F, processors=8)
rename.file(input=silva.bacteria2.pcr.fasta, new=silva.v34.prim2.fasta)
summary.seqs(fasta=silva.v34.prim2.fasta)
align.seqs(fasta=suis.trim.contigs.good.unique.fasta, reference=silva.v34.prim2.fasta)
summary.seqs(fasta=suis.trim.contigs.good.unique.align, count=suis.trim.contigs.good.count_table)
screen.seqs(fasta=suis.trim.contigs.good.unique.align, count=suis.trim.contigs.good.count_table, summary=suis.trim.contigs.good.unique.summary, start=17571, end=18912, maxhomop=8)
summary.seqs(fasta=suis.trim.contigs.good.unique.good.align, count=suis.trim.contigs.good.good.count_table)
filter.seqs(fasta=suis.trim.contigs.good.unique.good.align, vertical=T, trump=.)

I’ve also attached the summary.seqs of the pcr.seqs, align.seqs, screen.seqs and filter.seqs!20%20PM


#2




#3

In screen.seqs you use the following…

screen.seqs(fasta=suis.trim.contigs.good.unique.align, count=suis.trim.contigs.good.count_table, summary=suis.trim.contigs.good.unique.summary, start=17571, end=18912, maxhomop=8)

You need to set start=1. The bottom output from screen.seqs that you put up showed that most of your sequences started at position 1.


#4

That worked! Thank you very much!


closed #5

This topic was automatically closed 10 days after the last reply. New replies are no longer allowed.