Thanks a lot for your helpful comments.
• 0 views
•
link
written
by
anikcropscience •
How to obtain 200 bp flanking sequence from both sides of BLAST search?
Hello,
I have created a local blast database and blasted some genes against the reference using
blastn -query a.fasta -db BLASTDATABASE -max_target_seqs 5 -outfmt 6 -out output
I obtain a ~98% similarity score for the best hit. Now, I want to extract the fasta sequence of that BLAST search together with the 200 bp flanking sequence from both start and stop sides for designing a primer.
Can you please give me a suggestion on how to do this for a single BLAST search?
Thank you.
• 2,741 views
•
link
updated
by
lieven.sterck
•
written
by
anikcropscience •
•
written
by
anikcropscience •
1 answer
For extracting the region + context of the hit, you're best of with using blastdbcmd. it has an option -range that you can use to extract a certain region from a hit in the blastdb .
As with the approach manaswwm proposed (for the query) you can extract the start & stop from the blast hit table and use that in the blastdbcmd (potentially recalculating them with + and - 200 bp )
• 0 views
•
link
written
by
lieven.sterck
Log in to answer this question.
More posts like this
-
Running blastn to get best hits match using (-max_target_seqs 1 -max_hsps 1 ) ?
written by sunnykevin97
-
blastn not returning all the expected results
written by selarch •
-
Working backwards - How do I retrieve transcripts from a primer list using sequence and amplicon le…
written by dacotahm •
-
Blast+ Sam outfmt parameter check
written by larvalanobium •
-
Modifying Standalone BLAST output
written by glady
-
full sequences from command line tblastn
written by peachila •
-
Local blast memory usage
written by haylogo •
-
comparing two genomes
written by dr_mm_sh •
-
Obtain NCBI Taxonomy ID from local blast output
written by ntmarshall406 •
-
Blast Database Size Influence On Number Of Significant Hits
written by PoGibas
One way to do this, given that you have start and stop coordinates of the best hit, is by using samtools.
here
a.fasta- reference file using in blastn;Chr- chromosome of the hit sequence; forstartandstopyou can add 200bp to the coordinates of the best hit and extract it from the genomehere you're extracting context for the query, not the hit.
Ok. Thank you. I will try this and see if it works.
you want to extract context region for the hit or the query ?
and spoiler alert: you can't so this with a single BLAST search! (best will be two-step approach)
So, I want to extract the context region for the hit. Because the query gene is from another species and I am trying to find the homologous sequence in my reference species.