Blastn / Tblastn : Mapping The Features Of The Query To The Hit.

Question

[Pierre Lindenbaum]

I'm blasting+ (blastn+ or tblastn) an annotated sequence (a Genbank.xml sequence (nucleotide or protein) or an Uniprot.xml entry) against a DNA database.

Is there a standard tool to map the features of the query to the hit ? The best output would be a BED file.

Thanks.

Answer 1

[Pierre Lindenbaum]

I've answered my own question on my blog: http://plindenbaum.blogspot.fr/2013/07/mapping-annotations-of-query-sequence.html

$ java -jar dist/blastmapannots.jar I=P04514.xml B=blast_P04514_vs_AY065842.xml

AY065842    25    961    Non-structural_protein_3    943    +    25961    255,255,255    1    936    25
AY065842    34    469    RNA-binding    970    +    34    469    255,255,255    1    435    34
AY065842    472    640    Dimerization    947    +    472    640    255,255,255    1    168    472
AY065842    532    724    Interaction_with_ZC3H7B    917    +    532    724    255,255,255    1    192    532
AY065842    646    961    Interaction_with_EIF4G1    905    +    646    961    255,255,255    1    315    646
AY065842    520    733    coiled-coil_region    916    +    520    733    255,255,255    1    213    520