I was wondering if anyone would be able to help me with a problem I am having. I am using the example data on mothur to work through the MiSeq SOP. Everything has ran smoothly until I have gotten to here:
mothur > screen.seqs(fasta=stability.trim.contigs.unique.align, count=stability.trim.contigs.count_table, start=1969, end=11551)
It gives me an error that a sequence isn’t in my count table, but when I look in the file that specific sequence is there. I also noticed that I am missing a stability.trim.contigs.good.count_table file which I need for the following line and the sequential ones after it:
mothur > unique.seqs(fasta=stability.trim.contigs.unique.good.filter.fasta, count=stability.trim.contigs.good.count_table)
I’m not sure why this file can’t be generated. If anyone could help me, I would really appreciate it!