jianhong / ATACseqQC

ATAC-seq Quality Control

Home Page:https://jianhong.github.io/ATACseqQC/articles/ATACseqQC.html

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Pls help with the error... I am using hg38

hjack123 opened this issue · comments

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

Did you checked the sequence length of gal1, genome and phastCons100way.UCSC.hg38? Are they identical to each others?