While trying to assign transcript IDs to my tags using the
TSSs <- assignTxID(TSSs, txModels = txdb, swap="thick")
I was met with the following error:
Error: seqlengths(object) not identical to seqlengths(txModels)
Unless I am mistaken, this seems to arise because
is 455 because it is created using the BSgenome.Hsapiens.UCSC.hg38 database using the
quantifyCTSSs function, while
is 595 because txdb is the TxDb.Hsapiens.UCSC.hg38.knownGene database. Unfortunately my R-Fu is not advanced enough to solve this. Is there a way to easily rectify this, or am I missing something?