search

bcgsc / tigmint

⛓ Correct misassemblies using linked AND long reads
https://bcgsc.github.io/tigmint/
GNU General Public License v3.0
54 stars 13 forks source link

Add parameter for min reads per molec (for tigmint-molecule), fix missing molecule bug #20

Closed lcoombe closed 3 years ago

lcoombe commented 6 years ago

Added another parameter in tigmint-make, reads_per_molecule which can be used to set the -m parameter for tigmint-molecule.

lcoombe commented 6 years ago

Also found a bug where the last molecule from a set of alignments may be missed. If the last set of cur_reads from the BAM file does constitute a legitimate molecule, we didn't print it because we only look for molecules when we see another alignment that is a different barcode or reference -- which we don't if we hit the end of the file. Fixed by refactoring logic to produce molecule extents into a method, then calling that method when we hit the end of the file.

stale[bot] commented 3 years ago

This issue has been automatically marked as stale because it has not had any recent activity. It will be closed in 7 days if no further activity occurs. Thank you for your interest in Tigmint!