hi, dear guys
I performance blastn (-m 8) using a query file of many sequences, and for each query sequence, the output contains many fragmental hits of significance.
however, these hits have no overlap, and what is interesting is that most gaps < 300bp (much shorter than full-length of the query sequence).
so, how can i concatenate those closely related hits into one via setting a value (e.g 300bp) when these hits match the same subject (different regions), ——also to reduce the number of output hits per query.
for example:
http://www.freeimagehosting.net/bohcp
are there any scripts or tools for this purpose?
all your replies are welcome!
•
written
by
xiongtl2013 •
2 answers
You can use Biopython to parse the blast output and then concatenate the sequences that match your criteria. I would not recommend parsing the tabular output though, instead re-run blast and get the results in xml format since the that is easier to parse using a script.
In the biopython tutorial the chapters you would be interested in are 3, 4, 5, and 7.
Log in to answer this question.
-
Compare genotype genome sequences at basepair level
written by kashiff007
-
Joining blast hits together into single fasta sequence
written by Wilber0x •
-
Merge blastn hits sharing same subject
written by vinicius •
-
sorting blast results by sseq coverage
written by nheb1977 •
-
comparing two genomes
written by dr_mm_sh •
-
Blastall gives several hits to the same subject sequence
written by gislevestergaard •
-
Blast - Formatting Output
written by timjoncooper
-
Errors In A Bioperl Script, Please Help Me To Have A Check, Thanks In Advance!
written by xiongtl2013 •
-
How To Merge Contiguous Blast Hsps! (-M 8 Tab)
written by xiongtl2013 •
-
Blast Short Sequences Containing Repeats
written by grafabian •
