Hi everyone,
I am trying to process some CAGE data from Drosophila using CAGEfightR but I am having the following error. Has anybody seen this before? I cannot find anything related on internet.
This is the code I am using:
```
library(CAGEfightR)
bsg <- BSgenome.Dmelanogaster.UCSC.dm6
txdb <- TxDb.Dmelanogaster.UCSC.dm6.ensGene
odb <- org.Dm.eg.db::org.Dm.eg
bw_plus <- "/n/projects/sga/analysis/CAGE_data/cage_pos.bw"
bw_minus <- "/n/projects/sga/analysis/CAGE_data/cage_neg.bw"
bw_plus <- BigWigFileList(bw_plus)
bw_minus <- BigWigFileList(bw_minus)
names(bw_plus) <- names(bw_minus) <- "kc167_CAGE-seq_rep1"
genomeInfo <- SeqinfoForUCSCGenome("dm6")
CTSSs <- quantifyCTSSs(plusStrand=bw_plus, minusStrand=bw_minus, genome=genomeInfo)
Checking supplied genome compatibility...
Iterating over 2 genomic tiles in 1 samples using 38 worker(s)...
Importing CTSSs from plus strand...
Importing CTSSs from minus strand...
Merging strands... Error in .check_GenomicRanges_for_coercion_to_GPos(from, "UnstitchedGPos") :
all the ranges in the object to coerce to UnstitchedGPos must have a width of 1
Thanks!