Hi.
We notice that it is possible to select only the best hit. Is there a way to select hit b?
Thank you in advance,
Michal
Hi Michal,
thanks for the message. In the world of BLAST, “hit” means a sequence with one identifier in the database. In Sequenceserver, we number the hits. So that select button applies to that hit - i.e. NBlab330C09, and should download the entire sequence.
The letters refer to different HSPs (“high-scoring pairs” - ie. chunks that align between query and the hit sequence). So for your particular example, there are 29 or more such HSP. To differentiate HSPs from hits, in SequenceServer, we label them alphabetically a, b, c, and so-forth in the picture, but also in the pairwise alignments (your screenshot shows a. and the top of b.).
If you want to independently download all of HSPs, ie. the chunks of database sequence that aligns to your query sequences, you can click on the download Alignment button. That gives you a FASTA of all HSPs. The format is a bit funny, in that it includes the query sequence, “alignment strength”, “hit sequence” (see my screenshot - where I added some whitespace to facilitate interpretation). I think you’re looking for every third sequence in there.
Currently, those sequences are not identified using letters a, b, c, in the output. Instead, we show the from-to coordinates. I think we should have a, b, c listed and have opened an issue to that extent https://github.com/wurmlab/sequenceserver/issues/554
I hope this helps
Yannick