hg19转hg38
my_liftover <- function(peaks,from_to = 'hg19_hg38'){
if(from_to == 'hg19_hg38') {
path <- '/mdshare/node8/txmdata/scATAC_seq/AtacAnnoR/data/genome_files/hg19ToHg38.over.chain'
ch <- rtracklayer::import.chain(path)
one_liftover <- function(peak){Signac::GRangesToString(rtracklayer::liftOver(Signac::StringToGRanges(peak), ch)[[1]][1])}
# one_liftover('chr19-44821468-44827058')
if(length(peaks) > 1000){
out_peaks <- pbmcapply::pbmclapply(peaks,one_liftover,mc.cores = 10) %>% unlist()
} else{
out_peaks <- lapply(peaks,one_liftover) %>% unlist()
}
}
return(out_peaks)
# which(hg38_peaks == '--')
}
peaks <- c('chr19-44821468-44827058','chr1-44821468-44827058')
my_liftover(peaks)
