Test of Blast 2.2.31 -outfmt 15

test.sam

@HD     VN:1.2  GO:query
@SQ     SN:sp|P08100|OPSD_HUMAN LN:348
@SQ     SN:sp|P08100|OPSD_HUMAN LN:348
@SQ     SN:sp|P08100|OPSD_HUMAN LN:348
@SQ     SN:sp|P08100|OPSD_HUMAN LN:348
@SQ     SN:sp|P08100|OPSD_HUMAN LN:348
@SQ     SN:sp|P08100|OPSD_HUMAN LN:348
ref|NP_001009242.1|     0       sp|P08100|OPSD_HUMAN    1       255     348M    *       0       0       *       *       AS:i:1808       EV:f:0  NM:i:0  PI:f:96.55      BS:f:701.049
sp|P56514.1|OPSD_BUFBU  0       sp|P08100|OPSD_HUMAN    1       255     333M1I8M5I7M    *       0       0       *       *       AS:i:1560       EV:f:0  NM:i:6  PI:f:84.77      BS:f:605.52
gb|ADB45242.1|  0       sp|P08100|OPSD_HUMAN    11      255     328M    *       0       0       *       *       AS:i:1625       EV:f:0  NM:i:0  PI:f:94.82      BS:f:630.558
gb|ADB45229.1|  0       sp|P08100|OPSD_HUMAN    11      255     328M    *       0       0       *       *       AS:i:1625       EV:f:0  NM:i:0  PI:f:94.82      BS:f:630.558
prf||0811197A   0       sp|P08100|OPSD_HUMAN    1       255     190M1D157M      *       0       0       *       *       AS:i:1680       EV:f:0  NM:i:1  PI:f:93.37      BS:f:651.744
dbj|BAB21486.1| 0       sp|P08100|OPSD_HUMAN    1       255     328M1I20M5H     *       0       0       *       *       AS:i:1512       EV:f:0  NM:i:1  PI:f:81.32      BS:f:587.03

test.sh

# In referecne to: http://blastedbio.blogspot.com/2015/07/ncbi-working-on-sam-output-from-blast.html
# Data from the following page:
#https://github.com/peterjc/galaxy_blast/tree/master/test-data

makeblastdb -dbtype prot -in four_human_proteins.fasta -parse_seqids -out testdb
blastp -query rhodopsin_proteins.fasta -parse_deflines -db testdb -evalue 0.0001 -outfmt 15 -out test.sam