Logic error? The program will miss one read for each junction

dieterich-lab / FUCHS

FUCHS - FUll circle CHaracterization from rna-Seq

GNU General Public License v3.0

3 stars 6 forks source link

Open Christina-hshi opened 3 years ago

Christina-hshi commented 3 years ago

When the junction is found the first time in the ".junction" file, the corresponding supporting read is skipped or missed. https://github.com/dieterich-lab/FUCHS/blob/c2f36ad402606293e21ff12106ac40674a3b6c3a/FUCHS/get_readnames_from_DCC.py#L34

tjakobi commented 3 years ago

Dear @Christina-hshi,

Thank you for spotting this, I'll fix the issue shortly.