I was trying to extract the 5' UTR and 3'UTR coordinates for a particular isoform. My code is:
mart = useMart("ensembl")
getBM(attributes = c("5_utr_start","5_utr_end","3_utr_start","3_utr_end"),filters = "hgnc_symbol",value = "AAGAB",mart = ensembl)
Is there a way that I can specify that for the gene (AAGAB), I want information for isoform 001 (ie. AAGAB-001)? Else, how do I process the data given by my query above so that I can get this information?
Or should I be doing something else altogether?