Screen.seqs command criteria

Hello! Im new in mothur and I’m wondering which is the best parameters (start and end) to use the screen.seqs command, since I lost around 80% of the sequences with the following criteria:

  • screen.seqs(fasta=stability.trim.contigs.good.unique.align, count=stability.trim.contigs.good.count_table, summary=stability.trim.contigs.good.unique.summary, start=10368, end=25319, maxhomop=8)

It took 1192 secs to screen 1053264 sequences, removed 797534.

Its really bad to use this parameters? Im mean its valid to use the sequences as I did it or I will loose useful information?


Thanks !!

Why " start=10368, end=25319"? Half of your reads finish before that, so they get dumped. An then those starting after 10368 too… You should look at your alignment, and decide from which base to which base is the right alignment.

Pick the first 20 lines of your alignment (you do not need to check the whole dataset), check the alignment in MEGA or similar, and compare it to the region you are trying to get. I have to say I am surprised to see different starting points… Usually those are very consistent once oligos are trimmed. You might have leftover bases from your oligos?

1 Like