Skip to content

Instantly share code, notes, and snippets.

@zachcp

zachcp/test.sam

Last active July 3, 2018 13:57
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save zachcp/3781489afcdf26d3eeda to your computer and use it in GitHub Desktop.
Save zachcp/3781489afcdf26d3eeda to your computer and use it in GitHub Desktop.
Download ZIP
Test of Blast 2.2.31 -outfmt 15
Raw
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
Raw
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
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment