WGLab / LongReadSum

MIT License
12 stars 2 forks source link

Add RRMS support #47

Closed jonperdomo closed 1 year ago

jonperdomo commented 1 year ago

This adds support for RRMS data filtering via a provided CSV with the appropriate columns (read_id and decision).

longreadsum rrms --csv rrms_results.csv --input input.bam --output output_directory --threads 12

The update simply calls the BAM module once for each type of read (accepted and rejected), and filters the CSV read ids based on this flag. Thus two QC reports are created with suffix _rejected.html (left) or _accepted.html (right):

image

image

jonperdomo commented 1 year ago

I also increase plot font sizes and set the max. read and base quality values to 100 in histograms.