diff --git a/Makefile b/Makefile
index 2d830af..53d2070 100644
--- a/Makefile
+++ b/Makefile
@@ -3,6 +3,10 @@
 # Even if you never use these commands here, they provide a good example of ways to use the python
 #  scripts included in this project.
 
+# --- put test first, so `make` runs this and exits
+test:
+	python bin/test_requirements.py
+
 #---------------Bio 496, Project ---------------
 
 prot_and_dna: diamond_db bowtie2_index
@@ -119,9 +123,6 @@ protein_hits:
 
 
 
-test:
-	python bin/test_requirements.py
-
 demo: clean
 	#Download 100 sample reads
 	fastq-dump --outdir Input/SRA_datasets/ -X 100000 --skip-technical --readids  --dumpbase --clip SRR3403834
diff --git a/bin/bam_stats.py b/bin/bam_stats.py
old mode 100644
new mode 100755
index 0b08b4f..077b8ca
--- a/bin/bam_stats.py
+++ b/bin/bam_stats.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import argparse
 import pysam
 import os
diff --git a/bin/conserved_regions_csv.py b/bin/conserved_regions_csv.py
old mode 100644
new mode 100755
index e531c12..50615aa
--- a/bin/conserved_regions_csv.py
+++ b/bin/conserved_regions_csv.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3'
+
 import os
 import sys
 import argparse
diff --git a/bin/coverage_depth_csv.py b/bin/coverage_depth_csv.py
old mode 100644
new mode 100755
index 2f6e12d..0ec5a08
--- a/bin/coverage_depth_csv.py
+++ b/bin/coverage_depth_csv.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import argparse
 import os
@@ -92,4 +94,4 @@ def condense_list(big_list, condense_size=100):
             outfile.write(','.join(x) + '\n')
         with open(args.output + '.row_labels.txt', 'a') as outfile:
             outfile.write(args.output + '\n')
-    sys.exit(0)
\ No newline at end of file
+    sys.exit(0)
diff --git a/bin/csv_to_heatmap.py b/bin/csv_to_heatmap.py
old mode 100644
new mode 100755
index 1cc2c02..67aad41
--- a/bin/csv_to_heatmap.py
+++ b/bin/csv_to_heatmap.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import numpy as np
 from numpy import vstack  # VSTACK((a,b)) stacks B on TOP of A
 import seaborn as sns; sns.set()
diff --git a/bin/genome_coverage.py b/bin/genome_coverage.py
old mode 100644
new mode 100755
index 6326fa5..b564ef3
--- a/bin/genome_coverage.py
+++ b/bin/genome_coverage.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import pysam
 import argparse
diff --git a/bin/hits_by_organism.py b/bin/hits_by_organism.py
old mode 100644
new mode 100755
index dff4aeb..ba14ade
--- a/bin/hits_by_organism.py
+++ b/bin/hits_by_organism.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import os
 import argparse
 import sys
@@ -30,4 +32,4 @@
             organisms[line[1]] = organisms.get(line[1], 0) + int(line[2])
     with open(args.output, 'w') as outfile:
         for k in sorted(organisms, key=organisms.get, reverse=True):
-            outfile.write(delim.join([str(k),str(organisms[k])]) + '\n')
\ No newline at end of file
+            outfile.write(delim.join([str(k),str(organisms[k])]) + '\n')
diff --git a/bin/m8_hits.py b/bin/m8_hits.py
old mode 100644
new mode 100755
index 06ba7ca..9abdb5c
--- a/bin/m8_hits.py
+++ b/bin/m8_hits.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import os
 import argparse
diff --git a/bin/m8_reader.py b/bin/m8_reader.py
old mode 100644
new mode 100755
index d6d6071..e0cc515
--- a/bin/m8_reader.py
+++ b/bin/m8_reader.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 """
 M8_Reader! This package is designed to read .m8 files produced by diamond blastx, specifically with a command
 to include the actual DNA sequence matched to the protein. This information is not provided by default when using
diff --git a/bin/m8_stats.py b/bin/m8_stats.py
old mode 100644
new mode 100755
index a7be806..3f47fb7
--- a/bin/m8_stats.py
+++ b/bin/m8_stats.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import os
 import argparse
@@ -37,4 +39,4 @@
             line = [str(x) for x in line]
             outfile.write('\t'.join(line))
             outfile.write('\n')
-        outfile.close()
\ No newline at end of file
+        outfile.close()
diff --git a/bin/pbs_aggregator.py b/bin/pbs_aggregator.py
old mode 100644
new mode 100755
index 37ab1c8..813dc18
--- a/bin/pbs_aggregator.py
+++ b/bin/pbs_aggregator.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import os
 import argparse
diff --git a/bin/rapsearch_fastq_fixer.py b/bin/rapsearch_fastq_fixer.py
old mode 100644
new mode 100755
index 0350663..8f6dfbd
--- a/bin/rapsearch_fastq_fixer.py
+++ b/bin/rapsearch_fastq_fixer.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import os
 import sys
 import argparse
@@ -45,4 +47,4 @@
         ret_code = subprocess.call(mv_command)
         if ret_code != 0:
             sys.stderr.write("Error running command \"{}\"\n".format(' '.join(mv_command)))
-            sys.exit(1)
\ No newline at end of file
+            sys.exit(1)
diff --git a/bin/remove_short_reads.py b/bin/remove_short_reads.py
old mode 100644
new mode 100755
index 7eb3eaf..80f278f
--- a/bin/remove_short_reads.py
+++ b/bin/remove_short_reads.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import os
 import argparse
diff --git a/bin/sam_hits.py b/bin/sam_hits.py
old mode 100644
new mode 100755
index 1084fe3..408bbb0
--- a/bin/sam_hits.py
+++ b/bin/sam_hits.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import os
 import argparse
@@ -16,4 +18,4 @@
         sys.exit(1)
 
     samfiles = Sam_Reader(args.input)
-    samfiles.hits(write_file=args.output)
\ No newline at end of file
+    samfiles.hits(write_file=args.output)
diff --git a/bin/sam_stats.py b/bin/sam_stats.py
old mode 100644
new mode 100755
index b01619e..3cdef51
--- a/bin/sam_stats.py
+++ b/bin/sam_stats.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import os
 import argparse
diff --git a/bin/single_pbs.py b/bin/single_pbs.py
old mode 100644
new mode 100755
index a9468b7..75b8c36
--- a/bin/single_pbs.py
+++ b/bin/single_pbs.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import os
 import argparse
@@ -17,4 +19,4 @@
         sys.exit(1)
 
     data = Sam_Reader(args.input)
-    data.per_base_stats(write_file=args.output)
\ No newline at end of file
+    data.per_base_stats(write_file=args.output)
diff --git a/bin/split_bam.py b/bin/split_bam.py
old mode 100644
new mode 100755
index c3f3237..3bb1cd8
--- a/bin/split_bam.py
+++ b/bin/split_bam.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import argparse
 import pysam
@@ -87,4 +89,4 @@ def split_sam(infile, minimum_match_length, out_directory):
     if not args.length:
         args.length = 0
 
-    split_sam(args.input, args.length, args.output)
\ No newline at end of file
+    split_sam(args.input, args.length, args.output)
diff --git a/bin/target_sites.py b/bin/target_sites.py
old mode 100644
new mode 100755
index d0cd469..081ae9a
--- a/bin/target_sites.py
+++ b/bin/target_sites.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import argparse
 import sys
 import operator
diff --git a/bin/test_requirements.py b/bin/test_requirements.py
old mode 100644
new mode 100755
index 9562c4c..93a3a0e
--- a/bin/test_requirements.py
+++ b/bin/test_requirements.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import sys
 import subprocess
 
diff --git a/bin/xml_scrape.py b/bin/xml_scrape.py
old mode 100644
new mode 100755
index b6d4e08..d381785
--- a/bin/xml_scrape.py
+++ b/bin/xml_scrape.py
@@ -1,3 +1,5 @@
+#!/usr/bin/env python3
+
 import os
 import sys
 import argparse