Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

USEARCH chimera checking output seems wrong #69

Open
spacocha opened this issue Oct 6, 2014 · 2 comments
Open

USEARCH chimera checking output seems wrong #69

spacocha opened this issue Oct 6, 2014 · 2 comments

Comments

@spacocha
Copy link

spacocha commented Oct 6, 2014

Currently using UCHIME to do dbotu_chimera checking seems to be wrong. The printed output from uchime on the mock community samples says there are 36 nonchimeric sequences, but the output is only 34 sequences. I'm not sure what happened to the other 2 sequences, but usearch doesn't throw an error and it's not clear why there is this discrepancy. Obviously not about SmileTrain, but something to try to figure out.

@swo
Copy link
Contributor

swo commented Oct 6, 2014

can't do anything about this without any example data...

@spacocha
Copy link
Author

spacocha commented Oct 6, 2014

/net/radiodurans/alm/spacocha/raw_data/mock_mixes/dbOTU_dir

python ~/lib/SmileTrain/otu_caller.py -f all_R1_2file.fastq -r all_R2_2file.fastq -b vanja_59sam_mapping_rc.txt2 --split -n 1 --merge --demultiplex --qfilter --dereplicate --index --dbotu --maxee 0.5 --alignref /home/spacocha/tmpdir/silva.bacteria.m22534.m6428.filter.fasta --k_fold 10 --dbotu_chimeras --dbotu_split --JS_cutoff 0.02 --dbotu_id 0.08 --local

Although that command has the --JS_cutoff in it...

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants