jbloomlab / dms_variants

Analyze deep mutational scanning of barcoded variants.
Other
6 stars 9 forks source link

add outer flank regions for `IlluminaBarcodeParser` #88

Closed jbloom closed 7 months ago

jbloom commented 7 months ago

Added ability to parse second upstream / downstream region in IlluminaBarcodeParser by adding upstream2 and downstream2 parameters. Also modified IlluminaBarcodeParser so that reads will only be parsed if they are long enough to fully cover the region containing the barcodes and specified upstream / downstream sequences. Based on docs, this is how it was supposed to function before but did not. Additionally, this adds another row ("reads too short") to the fates from the barcode parser, as well as the outer_flank_fates option to report just failing the additional upstream and downstream regions.

Related to allowing the IlluminaBarcodeParser to be used to address this issue: https://github.com/jbloomlab/seqneut-pipeline/issues/40