Biostar Beta. Not for public use.
Can Ncbi Blast Only Search A Special Region In Database For Hits?
1
Entering edit mode
6.2 years ago

I want to blast my 36SE reads to a special region of HG19, does NCBI BLAST can do this special region search? If so, how can I do this?

for example, I only want my reads to be blasted to chr1:5000-10000 of hg19, how can I do this without creating a special small region?

ncbi blast • 1.6k views
ADD COMMENTlink
1
Entering edit mode
12 months ago
vkkodali ♦ 1.1k
United States

If you have the fasta sequence of chr1 on disk, you can use the -subject and -subject_loc options. From blastn -help:

-subject <File_In>
Subject sequence(s) to search
* Incompatible with:  db, gilist, seqidlist, negative_gilist,
db_soft_mask, db_hard_mask
-subject_loc <String>
Location on the subject sequence in 1-based offsets (Format: start-stop)
* Incompatible with:  db, gilist, seqidlist, negative_gilist,
db_soft_mask, db_hard_mask, remote
ADD COMMENTlink
0
Entering edit mode
5.7 years ago
vaskin90 • 280
Milan, Italy

Why creating a small region is a problem? If you have chr1 locally on your hard drive you can open it with a sequence editor and copy the region.

ADD COMMENTlink
0
Entering edit mode

Yes, I can create a small region sequence database, but I don't like to because I have to convert the relative coordinate to hg19 coordinate. You may ask me why not covert relative coordinate to hg19 coordinate by add a offset, but I insist on blasting regionly because I have soon many targets regions and I have to uniform all coordinate to detect something.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1