classify.seqs

Hi,

I am doing the pipeline of my Illumina paired end samples, but I got some unclear issues.

First, I don’t know why, but when I try to search for some topic in the forum and I go to Google, it shows me some results. However, I can’t access and it has been showing me the following message:

“Information
Sorry but you are not permitted to use the search system.”
So, basically I have to check page by page in order to check if I find a topic related to mine, but this is though.

Moving to my main issue… I am using the version v.1.39.4 and I amfollowing Myseq SOP tutorial. I ran align.seqs and in the end I got a lot of unclassified sequences.
This was the command I ran:

classify.seqs(fasta=stability.trim.contigs.good.unique.good.filter.unique.precluster.pick.fasta, count=stability.trim.contigs.good.unique.good.filter.unique.precluster.denovo.vsearch.pick.count_table, reference=trainset14_032015.pds.fasta, taxonomy=trainset14_032015.pds.tax, cutoff=80)

When I openned the tax summary file, I got it:


taxlevel rankID taxon daughterlevels total Day_01 Day_12 Day_17 Day_29 Day_30 Day_38 Day_42 Day_56 Day_68 0 0 Root 1 392548 43739 44545 43862 43225 54282 43321 34222 39640 45712 1 0.1 Bacteria 3 392548 43739 44545 43862 43225 54282 43321 34222 39640 45712 2 0.1.1 Actinobacteria 1 487 19 188 132 40 64 21 7 12 4 3 0.1.1.1 Actinobacteria 2 487 19 188 132 40 64 21 7 12 4 4 0.1.1.1.1 Actinobacteria_unclassified 1 1 0 0 0 0 1 0 0 0 0 5 0.1.1.1.1.1 Actinobacteria_unclassified 1 1 0 0 0 0 1 0 0 0 0 6 0.1.1.1.1.1.1 Actinobacteria_unclassified 0 1 0 0 0 0 1 0 0 0 0 4 0.1.1.1.2 Actinomycetales 1 486 19 188 132 40 63 21 7 12 4 5 0.1.1.1.2.1 Actinomycetales_unclassified 1 486 19 188 132 40 63 21 7 12 4 6 0.1.1.1.2.1.1 Actinomycetales_unclassified 0 486 19 188 132 40 63 21 7 12 4 2 0.1.2 Bacteria_unclassified 1 391292 43632 44314 43666 43114 54118 43224 34143 39541 45540 3 0.1.2.1 Bacteria_unclassified 1 391292 43632 44314 43666 43114 54118 43224 34143 39541 45540 4 0.1.2.1.1 Bacteria_unclassified 1 391292 43632 44314 43666 43114 54118 43224 34143 39541 45540 5 0.1.2.1.1.1 Bacteria_unclassified 1 391292 43632 44314 43666 43114 54118 43224 34143 39541 45540 6 0.1.2.1.1.1.1 Bacteria_unclassified 0 391292 43632 44314 43666 43114 54118 43224 34143 39541 45540 2 0.1.3 Proteobacteria 1 769 88 43 64 71 100 76 72 87 168 3 0.1.3.1 Proteobacteria_unclassified 1 769 88 43 64 71 100 76 72 87 168 4 0.1.3.1.1 Proteobacteria_unclassified 1 769 88 43 64 71 100 76 72 87 168 5 0.1.3.1.1.1 Proteobacteria_unclassified 1 769 88 43 64 71 100 76 72 87 168 6 0.1.3.1.1.1.1 Proteobacteria_unclassified 0 769 88 43 64 71 100 76 72 87 168

I guess this is not normal :frowning:

What could I try to do? Is this a command error or a bug in this version?

Thank you!

What region? What mothur version?

Do you have an example sequence that you think should be getting classified but isn’t?

Pat