dieterich-lab / DCC

DCC uses output from the STAR read mapper to systematically detect back-splice junctions in next-generation sequencing data. DCC applies a series of filters and integrates data across replicate sets to arrive at a precise list of circRNA candidates.
https://dieterichlab.org/software/
GNU General Public License v3.0
36 stars 20 forks source link

Check for sorted BAM files fails if first chromosome has less than 100 mapped reads #38

Closed tjakobi closed 7 years ago

tjakobi commented 7 years ago

This is due to a hard-coded limit of 100 reads used for checking the sorting. If more than one chromosome is encountered during the first 100 reads the sorting check will fail.

Best solution is probably to only check for .bai files and assume that the user did not sort by read name.