![]() |
Query Ensemble or USCS
Is there any automated way to get DNA sequences from UCSC genome browser or Ensemble for a list of chromosomal coordinates of interest.
I have several chromosomal coordinates sets and cannot do it manually. Thanks in advance. |
If you have a bunch of ranges, then it might be most straight forward to just download the genomes in fasta format and just write a quick parser (samtools provides an fai_fetch function in the C API that can be given a genomic range and will return the sequence).
|
At UCSC, you can upload a custom track bed file, then use the table browser and do an intersection. You could also do this with galaxy "Fetch Sequences".
On the command line, you have a number of options... bedtools getfasta comes to mind. |
All times are GMT -8. The time now is 04:35 PM. |
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.