Pls help with the error... I am using hg38
hjack123 opened this issue · comments
hjack123 commented
objs <- splitGAlignmentsByCut(gal1, txs=txs, genome=genome, outPath = outPath,
conservation=phastCons100way.UCSC.hg38)
Error in mergeNamedAtomicVectors(seqlengths(x), seqlengths(y), what = c("sequence", :
sequence chr10 has incompatible seqlengths:
- in 'x': 135534747
- in 'y': 133797422
JIANHONG OU commented
Did you checked the sequence length of gal1, genome and phastCons100way.UCSC.hg38? Are they identical to each others?