TSS sites of GRCM38
0
0
Entering edit mode
5.8 years ago
kanwarjag ★ 1.2k

How can I get TSS file (bed/ GTF) for GRCM38? I can see mm9 file but not GRCm38

Thanks

ChIP-Seq TSS • 2.9k views
ADD COMMENT
1
Entering edit mode
ADD REPLY
0
Entering edit mode

I tried to use the above TSS file but it is not compatible with GRCM38 and giver error that ist column annotations-NC.... can not be recognized.

ADD REPLY
0
Entering edit mode

Use it with what? What was the error?

ADD REPLY
0
Entering edit mode

Did you try the following in R?

library(GenomicFeatures)
chromInfoDF <- read.table("chrNameLength.txt")
chromInfoDF <- data.frame(chrom = chromInfoDF[,1], length = chromInfoDF[,2], is_circular = chromInfoDF[,1] == "MT")
txdb <- makeTxDbFromGFF(file = "Mus_musculus.GRCm38.83.gtf", format = "gtf", dataSource = "Ensembl 83", organism = "Mus musculus", chrominfo = chromInfoDF)
tss = promoters(txdb, upstream = 1500, downstream = 1500, columns = c("tx_name", "gene_id"))

NOTE: chrNameLength.txt file has two columns: chromosome number and length of chromosome, you can export the tss object as bed

ADD REPLY

Login before adding your answer.

Traffic: 2105 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6