Hi All,

Could you please let me know how to get rsnumbers using using IntronVariants()?

intron_regions     <- IntronVariants()

head(locateVariants(vcf_mod, enstxdb, intron_regions))
'select()' returned many:1 mapping between keys and columns
GRanges object with 6 ranges and 9 metadata columns:
     seqnames    ranges strand | LOCATION  LOCSTART    LOCEND   QUERYID        TXID         CDSID      GENEID       PRECEDEID
         <Rle> <IRanges>  <Rle> | <factor> <integer> <integer> <integer> <character> <IntegerList> <character> <CharacterList>
[1]    chr19    572764      + |   intron        63        63        82           1          <NA>         BSG            <NA>
[2]    chr19    572770      + |   intron        69        69        83           1          <NA>         BSG            <NA>
[3]    chr19    572786      + |   intron        85        85        84           1          <NA>         BSG            <NA>Best,

I would expect the rsnumbers to appear as rownames in this output. I can provide details of my input to locateVariants if it would help. The rsnumbers ARE in the vcf.

Jon

Hi,

If the rsnumbers are the ID column of the VCF file they should be the rownames in the readVcf() output. This (which you may already know) is described in the ?readVcf man page:

     rowRanges: The CHROM, POS, ID and REF fields are used to create a
          ‘GRanges’ object. Ranges are created using POS as the start
          value and width of the reference allele (REF). By default,
          the IDs become the rownames ('row.names = FALSE' to turn this
          off). If IDs are missing (i.e., ‘.’) a string of
          CHROM:POS_REF/ALT is used instead.  The ‘genome’ argument is
          stored in the seqinfo of the ‘GRanges’ and can be accessed
          with ‘genome(<VCF>)’.

Are the rsnumbers in the 'vcf_mod' object? Is the question why didn't the rownames from 'vcf_mod' propagate to the output of locateVariants()?

Valerie

I edited my answer after re-reading the question.