Hi, I am not sure if this is the right place to ask but here is my question.
When I take trancription start positions from the site below the total tss count is around 37000. http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/refGene.txt.gz When I query at ensembl the total tss count is around 137000. library(biomaRt) mart <- useMart( "ensembl", dataset = "hsapiens_gene_ensembl" ) tss <- getBM( attributes=c( "transcript_start", "transcript_end", "strand" ), filters="chromosome_name", values=c(1:23), mart=mart ) So is this the set of TSSs that includes alternative versions of genes or I am missing something ? Thank you, Best regards, Kirti Prakash _______________________________________________ Bioc-sig-sequencing mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
