This is a test version of Biostars. For the public version, visit https://www.biostars.org.
TSS sites of GRCM38

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

Thanks

chip-seq tss
• 3,561 views • link
written by kanwarjag kanwarjag

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.

• 1 views • link
written by kanwarjag kanwarjag

Use it with what? What was the error?

• 1 views • link
written by GenoMax GenoMax

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

• 1 views • link
written by anshupa.vssut • anshupa.vssut

0 answers

No answers yet.

Log in to answer this question.

More posts like this