how to extract specific gene coordinates by passing a list of genes using homo.sapiens library
1
0
Entering edit mode
Tharun • 0
@92c31cb5
Last seen 4 months ago
India

can we extract specific rows of data from Homo.sapiens library like
transcripts(Homo.sapiens, columns=c( "SYMBOL"=="BRCA2","ENTREZID")) or by passing specific list of gene.

missRows • 254 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 2 days ago
United States
> subset(transcripts(Homo.sapiens, columns = "SYMBOL"), SYMBOL == "BRCA2")
'select()' returned 1:1 mapping
between keys and columns
GRanges object with 4 ranges and 1 metadata column:
      seqnames            ranges
         <Rle>         <IRanges>
  [1]    chr13 32889617-32907524
  [2]    chr13 32889617-32973809
  [3]    chr13 32928998-32936830
  [4]    chr13 32945093-32953652
      strand |          SYMBOL
       <Rle> | <CharacterList>
  [1]      + |           BRCA2
  [2]      + |           BRCA2
  [3]      + |           BRCA2
  [4]      + |           BRCA2
  -------
  seqinfo: 93 sequences (1 circular) from hg19 genome

Login before adding your answer.

Traffic: 427 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