Convert blast tabular output to bed format
1
1
Entering edit mode
3.0 years ago
Denis ▴ 200

I need to convert blast output in tabular format (format 6) to bed format. I'm interested only in qstart and qend coordinates (the 7th & 8th columns respectively). As far i know bed format is 0-based, but blast probably uses 1-based coordinates. Am i right?

So to do so with AWK i need to execute the code below:

cat balst_tabular_output.txt | awk '{print($2"\t"$7-1"\t"$8)}' > blast_otput.bed
genome • 2.2k views
ADD COMMENT
0
Entering edit mode

I've manually checked blast output and corresponding query sequences. So, blast definitely uses 1-based coordinates.

ADD REPLY
3
Entering edit mode
3.0 years ago
cschu181 ★ 2.6k

Yes, BLAST is 1-based and your conversion should be right (keeping qend 1-based makes it compatible with the open interval style that bed uses.)

ADD COMMENT

Login before adding your answer.

Traffic: 1623 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6