Entering edit mode
7.3 years ago
Bumble
▴
10
If I understand the Blast manual correctly, the 'qseqid qseq' command will write the parts of the query sequence that match to a sequence in the database into the output. I would however like to get the whole query sequence. Is this possible, maybe with one of the length commands?
I don't think this is possible with BLAST itself. I do recall doing something similar with the
blastdbcmd
. Can't help much more, but you might be able to track something down from there.The parts of the sequences you will see, in the alignments are defined by the query start/end positions, and the subject start/end positions. The length command specifies the min alignment length to report.