[BioC] converting position from '-' strand to '+' strand

michael watson (IAH-C) michael.watson at bbsrc.ac.uk
Mon Apr 5 22:20:55 CEST 2010


Begin from the 5' end of what? Chromosomes don't have a 5' "end" 
________________________________________
From: bioconductor-bounces at stat.math.ethz.ch [bioconductor-bounces at stat.math.ethz.ch] On Behalf Of Tim Smith [tim_smith_666 at yahoo.com]
Sent: 05 April 2010 20:07
To: bioc
Subject: [BioC] converting position from '-' strand to '+' strand

Apologies if this seems like a trivial question.

I wanted to have a consistent set of locations and wanted to all the locations to begin from the 5' end. How can I convert locations that are given for the '-' strand?
For example:
-----------------------------
library(biomaRt)

mart.obj <- useMart(biomart = 'ensembl', dataset = 'hsapiens_gene_ensembl')

atb <- c('ensembl_gene_id', 'chromosome_name', 'start_position', 'end_position', 'strand')

mir.locs <- getBM(attributes=atb, filters="biotype", values="miRNA", mart=mart.obj)
mir.locs[1:5,]
>  ensembl_gene_id chromosome_name start_position end_position strand
1 ENSG00000222732               5      171706206    171706319      1
2 ENSG00000207864               9       97847727     97847823      1
3 ENSG00000221173               9      129338809    129338909     -1
4 ENSG00000222961               5       32379501     32379581     -1
5 ENSG00000221058              18       51612956     51613026     -1


----------------------------
Is there a quick way that I can convert the last 3 rows so that they reflect positions from the 5' strand?
many thanks!



        [[alternative HTML version deleted]]

_______________________________________________
Bioconductor mailing list
Bioconductor at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list