ATACseqQC icon indicating copy to clipboard operation
ATACseqQC copied to clipboard

Adjust the read start sites of mouse data

Open vibook opened this issue 6 years ago • 1 comments

Hi, I installed ATACseqQC to offset ATAC-seq aligning. The example code in the Guide is for human data. I downloaded BSgenome.Mmusculus.UCSC.mm9 but I'm struggling with another part.

seqlev <- "chr1" ## subsample data for quick run which <- as(seqinfo(Hsapiens)[seqlev], "GRanges")

I want to offset all sample data and I cannot figure out how to change seqinfo() instead of Hsapiens.

Thank you.

vibook avatar Nov 07 '18 19:11 vibook

library(BSgenome.Mmusculus.UCSC.mm9) which <- as(seqinfo(Mmusculus)[seqlev], "GRanges")

jianhong avatar Nov 12 '18 21:11 jianhong