Distinguish no-hits and hits that don't pass the post-BLAST filtering - Githubissues

ohnosequences / mg7

Configurable and scalable 16S metagenomics data analysis

https://goo.gl/y3rZFD

GNU Affero General Public License v3.0

3 stars 3 forks source link

Distinguish no-hits and hits that don't pass the post-BLAST filtering #130

Open laughedelic opened 7 years ago

laughedelic commented 7 years ago

Currently both the sequences that don't have BLAST hits and those that do, but are filtered afterwards (qcovs & pident thresholds) are saved in no-hits FASTA files. So for the latter ones we just through away the BLAST results which in some cases could be useful to rerun the filter with different threshold.

we should distinguish sequences with no hits and the hits that don't pass the filter
there should be an easy way to rerun the filter on the latter ones
same about rerunning BLAST with different parameters on the former ones

/cc @rtobes

eparejatobes commented 6 years ago

I like https://github.com/ohnosequences/mg7/issues/130#issue-242943320; this issue should focus on simply producing separate output for each case, though.

rtobes commented 6 years ago

I agree.