Hi,
I've been using the blastn (version 2.2.28+) standalone tool against a custom formatted genome via:
blastn -db BLASTDB -word_size 7 -query input.fa -out filename -perc_identity 100 -outfmt 6 -max_target_seqs 2
To discard non-perfect hits and show only the 2 top hits.
The output file has a great format however is there a way to add an extra column that contains the actual target-seq (sequence of the matched hit)? Such that the fields are:
query id, subject id, % identity, alignment length, mismatches, gap opens, q. start, q. end, s. start, s. end, evalue, bit score, sequence
Thanks!
- TJC