Produce fasta file from genome coordinates
0
0
Entering edit mode
2.9 years ago
claudio • 0

I wish to query blastn using genome coordinates rather than a specific gene but it seems that I would need a fasta sequence in order to do so. Is there a way to generate a fasta sequence from genomic coordinates such as the following: chr11:5,300,000-5,600,000

fasta coordinates • 925 views
ADD COMMENT
1
Entering edit mode

if you have reference genome as fasta file, use samtools faidx

ADD REPLY
1
Entering edit mode

you can use the built-in functionality of blast to do this.

check out the blastdbcmd command (entry + range )

ADD REPLY
0
Entering edit mode

Is there a way to do this from the ncbi website or another website because I don't have the capacity on my pc to download large database files.

ADD REPLY
0
Entering edit mode

OK, valid remark indeed.

You can download a fasta version of your hits from the NCBI blast output webpage. Not sure if you can just 'randomly' get a range .

ADD REPLY

Login before adding your answer.

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