Currently, the way my pipeline works is that it outputs a bam for the spike-in controls, and another for the genomic contigs. Due to this, I cannot run scan_genomic_contigs() because if I run on only the spike-in bam, there are no genomic contigs present.
Can we have a function to load in spike in data and make GLM and another function to load in genomic contigs to use GLM on?
Currently, the way my pipeline works is that it outputs a bam for the spike-in controls, and another for the genomic contigs. Due to this, I cannot run scan_genomic_contigs() because if I run on only the spike-in bam, there are no genomic contigs present.
Can we have a function to load in spike in data and make GLM and another function to load in genomic contigs to use GLM on?