search

biomedicalinformaticsgroup / Sargasso

Sargasso disambiguates mixed-species high-throughput sequencing data.
http://biomedicalinformaticsgroup.github.io/Sargasso/
Other
8 stars 4 forks source link

Change minmatch threshold option to max unmatched? #31

Closed lweasel closed 8 years ago

lweasel commented 8 years ago

I think we might change the minmatch threshold option to the maximum number of bases that are allowed not to match? That would mean we can specify a default value (e.g. 0) that is independent of the length of the reads. Obviously there will have to be a little bit of extra processing in filter_sample_reads.py to turn the combination of the length of a particular read and the option value into a minmatch threshold for that read.

Sound reasonable?

s-heron commented 8 years ago

Sounds reasonable indeed