The filtered reads for long reads (more than 50bp) was wrongly set to 0.

yanghchen / perm

Automatically exported from code.google.com/p/perm

1 stars 0 forks source link

The filtered reads for long reads (more than 50bp) was wrongly set to 0. #20

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago

The number of filtered reads count printed at the end of mapping was
wrongly set to zero for all long reads. This will be fixed in the next
release as soon as possible.

Original issue reported on code.google.com by Chen.Yan...@gmail.com on 7 May 2010 at 2:17

GoogleCodeExporter commented 9 years ago

This has been fixed in 0.2.9.1. 
The counter is now correct for long reads. If you think there are too many reads
filtered with 'N', try --includeReadsWN. Internally, this will include reads by
encoding 'N' or '.' as one nucleotide and could be match or mismatches by 
chance. 
This will rescue some unmapped reads but may be bias the result.

Original comment by Chen.Yan...@gmail.com on 12 May 2010 at 12:14