Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • RDW
    replied
    You could get the hg18 (build 36) start and end co-ordinates of accessions like AY293855 from the relevant table at NCBI:

    ftp://ftp.ncbi.nih.gov/genomes/MapView/Homo_sapiens/

    e.g.:

    ftp://ftp.ncbi.nih.gov/genomes/MapVi...eq_extra.md.gz

    (other accessions may be in other tables). If you have internal co-ordinates relative to the AY293855 sequence, you can then get the equivalent hg18 co-ordinates by simple arithmetic. However, some of the numbers quoted above are too big for the 139947 bp AY293855 sequence, so perhaps they refer to something else, like a clone or contig.

    Leave a comment:


  • Bert
    replied
    Hello Mathew,

    Ok, so as far as I understand now you want to know where the region of GenBank / EMBL / DDBJ sequence AY293855 with coordinates start 140618 and stop 140769 is located on the NCBI36 (hg18) genome assembly.

    I would suggest the following steps:

    1 - Retrieve the sequence for AY293855
    2 - Get from this sequence the sequence corresponding to bp 140618 - 140769
    3 - BLAST this sequence against GRCh37 (hg19) at e.g. NCBI or Ensembl
    4 - Lift the resulting coordinates over from GRCh37 (hg19) to NCBI36 (hg18)

    It would of course be easier to directly BLAST to the NCBI36 (hg18) assembly, but I don't see this option in NCBI BLAST. In Ensembl it certainly isn't possible to BLAST to older assemblies.

    As we have annotated several Illumina arrays in Ensembl, another option would be to do step 3 using Ensembl BLAST and then directly look in Ensembl if there are any probes mapped to your regions of interest, instead of having to use the NCBI36 (hg18) based Illumina files.

    Hope this helps.

    Leave a comment:


  • Pankaj kumar
    replied
    Yes bret was right.
    I think you need to extract specific sequence from each clone. For example
    PEG3 AC006115 19 163172 163493
    PEG3 AC006115 19 161667 161938
    Same AC number but two different probe. You need to extract the particular region from each clone and need to do the blast search with hg19.

    thanks
    Pankaj

    Leave a comment:


  • mathew
    replied
    Bert,
    You are right I have locations which have been published in paper with the corresponding gene bank accession # (Ay...,AC.. etc) I only have handful of these. I have sent to you partial list via email. I was thinking perhaps I can blast theses positions (of gene bank accession #) but still the location which I get is not corresponding to the actual coordinates cooresponding to particular gene. Any help will be appreciated.

    Leave a comment:


  • Bert
    replied
    Pankaj,

    Please read my reply carefully. Mathew has a GenBank / EMBL / DDBJ ID and some coordinates. If you would look up the ID (AY293855), you would see that this is the sequence for the human IGF2R gene (http://www.ebi.ac.uk/ena/data/view/AY293855). When you then would look up where this gene is located in GRCh37 (hg19) you would see that it is on chromosome 6 at position 160390131-160534539 (http://www.ensembl.org/Homo_sapiens/...0131-160534539). Now, that are quite different numbers from 140600-1407890, right? So, the coordinates Mathew has cannot be genomic coordinates. I suspect they are relative to e.g. a clone or so.

    Cheers,
    Bert

    Leave a comment:


  • Pankaj kumar
    replied
    Hi Bert,
    He has genomic co-ordinate.
    see "Chr no (6 ) and start (140600) and stop (1407890)"

    Leave a comment:


  • Bert
    replied
    What Matthew has are not genomic coordinates as far as I can see .... AY293855 is a sequence from the IGF2R gene, and that is located in GRCh37 (hg19) at position 160390131-160534539 (http://www.ensembl.org/Homo_sapiens/...0131-160534539).

    Cheers,
    Bert

    Leave a comment:


  • Pankaj kumar
    replied
    Hi,
    check the ucsc genome browser.

    You need to give input as chr# start_position end_position
    select your assembly and you are done.
    Pankaj

    Leave a comment:


  • Bert
    replied
    Mathew,

    It's still not clear to me what your coordinates represent. Relative to which sequence are they?

    And how many of these ID's do you have?

    Cheers,
    Bert

    Leave a comment:


  • mathew
    replied
    genomic coordinates

    What I have is Genebank accession Id (like AY293855) - Chr no (6 ) and start (140600) and stop (1407890) so I want to find out corresponding genomic coordinates in Hg18.
    Thanks Bert

    Leave a comment:


  • Bert
    replied
    Hello Mathew,

    From the title and content of your post it's not clear to me what you want to accomplish .... What data do you exactly have and what data do you want to retrieve?

    Cheers,
    Bert

    Leave a comment:


  • getting genomic coordinates from gene accesion information

    I have list of gene accession numbers (e.g AY.... - Chr no start --- and stop ---).
    Is there any way that I can get genomic location of these regions in Hg18.
    I have tried USCS geneom browser Table that will not work.
    Thanks

Latest Articles

Collapse

  • seqadmin
    Exploring the Dynamics of the Tumor Microenvironment
    by seqadmin




    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
    07-08-2024, 03:19 PM
  • seqadmin
    Exploring Human Diversity Through Large-Scale Omics
    by seqadmin


    In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
    06-25-2024, 06:43 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 07-10-2024, 07:30 AM
0 responses
23 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 09:45 AM
0 responses
198 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 08:54 AM
0 responses
209 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-02-2024, 03:00 PM
0 responses
191 views
0 likes
Last Post seqadmin  
Working...
X