This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Can Ncbi Blast Only Search A Special Region In Database For Hits?

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
• 3,282 views • link
updated by vkkodali_ncbi vkkodali_ncbi • written by litiancheng.gansu • litiancheng.gansu

2 answers

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
• 0 views • link
written by vkkodali_ncbi vkkodali_ncbi

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.

• 0 views • link
written by vaskin90 vaskin90

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.

• 0 views • link
written by litiancheng.gansu • litiancheng.gansu

Log in to answer this question.

More posts like this