Why is `send` lower than `sstart` from blast format 6 output?
3
2
Entering edit mode
14 months ago
O.rka ▴ 290

Is this behavior a result of strand specificity? For example, on the first line would "FUN_003800-T1" be on the negative strand while "FUN_003804-T1" on the positive strand so the ssend would be a smaller value than the sstart?

enter image description here

blast genomics • 591 views
ADD COMMENT
3
Entering edit mode
14 months ago
Asaf 8.6k

Yes, the match is on the negative strand.

ADD COMMENT
0
Entering edit mode

Ok, that's what I thought but I couldn't find it in the documentation and just wanted to be sure. Can you put this as an answer so I can mark it correct?

ADD REPLY
0
Entering edit mode
14 months ago
Fatima ▴ 930

Yes the query is aligned with the reversed complement of the target(subject). Blast has an option -strand that can be set to plus, minus, both. If you only want matches from the positive strand you can use -strand plus.

Please see the link below for more information on matching with negative strand:

https://biology.stackexchange.com/a/8162

ADD COMMENT
0
Entering edit mode
7 days ago
David • 0

I just wanted to add for the record that the sstart and send coordinates are all relative to the same reference strand you used for BLAST, regardless of which strand your query sequence aligned to. In other words, for all hits where (send-sstart) < 0, you can simply flip the two coordinates (i.e. assign send to sstart and vice versa) if you need consistency (e.g. when creating a GRanges object which requires send > sstart).

ADD COMMENT

Login before adding your answer.

Traffic: 2363 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