Filter question

[pamelacamejo]

Hi,

I run the following command to classify some 16S amplicon sequences:

hmmufotu gg_97_otus_GTR unmapped_otus_rev.fasta -s 0 -t 0 > hmmufotu_gg_97_unmapped.out

I noticed that some reads have very poor alignment coverage with the database, and even that, they have a taxonomic assignation. Do you recommend to use any filter for taxonomic assignation after running hmmufotu?

[e00011027]

Hi, Looks like you are trying to map some unmapped reads. First make sure you have the latest HmmUFOtu installed, and second make sure HmmUFOtu guessed your 16S read orientation correctly. Sometimes 16S primers are miss labeled some forward reads are reverse reads and vice versa. If everything is OK you could try to filer reads in the hmmufotu-summ step. There are a few identity options and try something like 0.5 is reasonable. Best Qi

…

On Tue, Jun 25, 2019, 2:31 PM pamelacamejo ***@***.***> wrote: Hi, I run the following command to classify some 16S amplicon sequences: hmmufotu gg_97_otus_GTR unmapped_otus_rev.fasta -s 0 -t 0 > hmmufotu_gg_97_unmapped.out I noticed that some reads have very poor alignment coverage with the database, and even that, they have a taxonomic assignation. Do you recommend to use any filter for taxonomic assignation after running hmmufotu? — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#4?email_source=notifications&email_token=ABXZ44OO6CE75D3NARFACMDP4JQAHA5CNFSM4H3LCKDKYY3PNVWWK3TUL52HS4DFUVEXG43VMWVGG33NNVSW45C7NFSM4G3T2TXQ>, or mute the thread <https://github.com/notifications/unsubscribe-auth/ABXZ44KPVALI4PLXSCHWBADP4JQAHANCNFSM4H3LCKDA> .