diff --git a/Makefile b/Makefile
index d02bedf93dce0bf17a9e3d3b1d0c40791a24a2ef..9b673b241dec18a0495da4bc5d8fb4db4194c686 100644
--- a/Makefile
+++ b/Makefile
@@ -5,7 +5,12 @@ test:testsuite/test_rep.sh
 	@echo "all tests passed"
 
 test_y:
-	src/phybema.py --tools mash andi dnadiff -- testdata/Yersinia/
+	src/phybema.py --tools mash andi -- testdata/Yersinia/
+
+test_e:
+	src/phybema.py --tools mash andi -- testdata/E_coli/ 
 
 clean:
 	@${RM} -r temp out
+	@${RM} -f testdata/E_coli/sequence.fasta
+	@${RM} -f testdata/Yersinia/sequence.fasta
diff --git a/src/callProg.py b/src/callProg.py
index 967c0cb77a6a04ee364f22c667b592252c9ff531..6e8c92a548ad7d0ca5113d336b328872ffac87fc 100644
--- a/src/callProg.py
+++ b/src/callProg.py
@@ -25,7 +25,7 @@ def run_prog(progname,call_list,fixed_out_file_name,progoutfilepath):
   sys.stderr.write("# starting program {}\n".format(progname))
   start_time = datetime.now()
   
-  sys.stderr.write(' '.join(call_list))
+  sys.stderr.write('# {}\n'.format(' '.join(call_list)))
   try:
     output = subprocess.run(call_list,stdout=subprocess.PIPE,
                             stderr=subprocess.PIPE, universal_newlines=True)
diff --git a/src/parser.py b/src/parser.py
index 16163eb47b3ec3255705c1f66c3fafe6f04d1587..0724f24b2b4698c1ce63ee7f932698f0910209dc 100644
--- a/src/parser.py
+++ b/src/parser.py
@@ -122,7 +122,6 @@ def create_datatuplelist(datasetpaths):
         uncompressed_file = "{}".format(mo.group(1))
         if not os.path.isfile(uncompressed_file):
           call_list = ['gzip','-k','-d',datafile]
-          sys.stderr.write('{}\n'.format(' '.join(call_list)))
           try:
             output = subprocess.run(call_list,stdout=subprocess.PIPE,
                                     stderr=subprocess.PIPE, 
@@ -143,7 +142,7 @@ def create_datatuplelist(datasetpaths):
       reftree = None
     # Check if more than one reftree was found
     elif len(reftreefilepaths) > 1:
-      sys.stderr.write("{}:There is more than one file with a {} suffix "
+      sys.stderr.write("{}: There is more than one file with a {} suffix "
                        "in directory {}\n"
                        .format(sys.argv[0], REFTREE_SUFFIXES, datasetpath))
       exit(1)
diff --git a/testdata/E_coli/.gitignore b/testdata/E_coli/.gitignore
new file mode 100644
index 0000000000000000000000000000000000000000..83836288e2db3c96f2e4892241bc2a4a262ccc1b
--- /dev/null
+++ b/testdata/E_coli/.gitignore
@@ -0,0 +1 @@
+sequence.fasta
diff --git a/testdata/E_coli/sequence.fasta.gz b/testdata/E_coli/sequence.fasta.gz
new file mode 100644
index 0000000000000000000000000000000000000000..a684adfcd2b476f812c4a1781b58bad5945016c5
Binary files /dev/null and b/testdata/E_coli/sequence.fasta.gz differ
diff --git a/testdata/E_coli/tree.nh b/testdata/E_coli/tree.nh
new file mode 100644
index 0000000000000000000000000000000000000000..7ed5c64cdc624602cebccbd8a54601d12c7ebde6
--- /dev/null
+++ b/testdata/E_coli/tree.nh
@@ -0,0 +1 @@
+(((CP000970,CU928164),(FM180568,((CP000247,AE014075),(CU928162,(CP000243,(CU928161,CP000468)))))),(CU928163,((CP000034, (BA000007, AE005174)),(((AP009048, U00096.3),(CP000802,CP000946)),((CP000266,(AE005674, AE014073)), ((CP000038,(CP000036,CP001063)),(CP000800,(CU928145,(AP009240,CU928160)))))))));