valegale / ONT_methylation

1 stars 0 forks source link

Why does the pipeline extract reads from a BAM to map them again? #4

Open hoelzer opened 2 days ago

hoelzer commented 2 days ago

Hey,

While restructuring the pipeline code, I saw that your current inputs are BAM files, and then you extract the reads to map them again with minimap2—why? There is probably some good reason for that.

We will probably anyway provide more input options in the future, such as POD5 (adding a Dorado process) or FASTQ (and then directly doing the minimap2 mapping, etc...)