Entering edit mode
8.2 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.