Skip to content

Latest commit

 

History

History
15 lines (11 loc) · 822 Bytes

README.md

File metadata and controls

15 lines (11 loc) · 822 Bytes
Raw

alt text

Simple python script to run BLAST on a local database over many sequences and sum up the results. Also includes proper output to analyse how different statistics such as % identity, alignment length, missmatches...etc. covariate.

To run the script:

python3 blastsumer.py -i <input file or directory in .fa (fasta) format> -o <output directory (does not need to exist)>

Optional parameters include:

  • -v or --verbose add this argument for more verbose output in the terminal while running
  • --pid this is a threshold for minimum % ID of alignment of BLAST hits. Default set to 80.0 (80%)
  • --eval this is a threshold for maximum e-value of BLAST hits. Default set to 1e-50.