This is a test version of Biostars. For the public version, visit https://www.biostars.org.
A Way To Get Dna Sequence From Specific Coordinates Using An Online Service

Is there an easy way, using some online service, to retrieve a bit of sequence from the C. elegans genome. As an example, I could enter chromsome 4 and position 5,000,000 to 5,002,000 and get the DNA sequence in that location.

I know how to do this from the command line, but I'd like a much easier way for people who don't use the command line.

sequence extraction
• 10,162 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by KCC KCC

3 answers

Well the first way that comes to mind is using Ensembl or Ensembl Genomes:

  1. Go to the C. elegans Ensembl: http://www.ensembl.org/Caenorhabditis_elegans/, or the C. elegans Ensembl Genomes: http://metazoa.ensembl.org/Caenorhabditis_elegans/, these are using the same genome build so the choice is up to you.
  2. Search for the required genome region by typing it into the search box, in your case the region is "IV:5000000-5002000", and click the "Go" button
  3. Click the "Export data" item in the left-hand menu
  4. Select the required format, etc. and click "Next"
  5. Select how you want to view/download the sequence data for the region.

You could of course also do this by constructing the URL used to access the result, for example:

http://www.ensembl.org/Caenorhabditis_elegans/Export/Output/Location?db=core;flank3_display=0;flank5_display=0;output=fasta;r=IV:5000000-5002000;strand=feature;coding=yes;cdna=yes;peptide=yes;utr3=yes;exon=yes;intron=yes;genomic=unmasked;utr5=yes;_format=Text

or:

http://metazoa.ensembl.org/Caenorhabditis_elegans/Export/Output/Location?db=core;flank3_display=0;flank5_display=0;output=fasta;r=IV:5000000-5002000;strand=feature;coding=yes;cdna=yes;peptide=yes;utr3=yes;exon=yes;intron=yes;genomic=unmasked;utr5=yes;_format=Text

But that is only really worth it if you are going to fetch a lot of sections from these genomes.

I suspect that the same thing can be achieved using WormBase but I haven't used it much so I'm not sure how to go about it there.

• 164 views • link
updated by Ram Ram • written by Hamish Hamish

A really easy way to get this data would be using the Ensembl REST API sequence_region module. Details are here:

http://beta.rest.ensembl.org/documentation/info/sequence_region

In essence you can make a URL like this:

http://beta.rest.ensembl.org/sequence/region/caenorhabditis_elegans/IV:5000000..5002000:1?content-type=text/plain

Which will just take you to a page with the sequence.

Change the species and coordinates to suit your purpose.

• 24 views • link
written by Emily Emily

Log in to answer this question.

More posts like this