diff --git a/_modules/functions_PhageTerm.py b/_modules/functions_PhageTerm.py
index 63db8b31c46153491da95b91c6a9698b8e8c6536..a20be630d08f0d9faf60941029b21af75adaf017 100644
--- a/_modules/functions_PhageTerm.py
+++ b/_modules/functions_PhageTerm.py
@@ -49,6 +49,7 @@ import pandas as pd
 # String
 #import cStringIO
 import io
+import gzip
 
 # PDF report building
 import time
@@ -127,11 +128,13 @@ def readsCoverage(inRawDArgs,refseq,inDArgs,fParms,return_dict, core_id,line_sta
         sys.stdout.flush()
     
     # Mapping
-    filin            = open(inRawDArgs.fastq)
+    #filin            = open(inRawDArgs.fastq)
+    filin = gzip.open(inRawDArgs.fastq, "rt") if inRawDArgs.fastq.endswith(".gz") else open(inRawDArgs.fastq)
     line             = filin.readline()
     
     if inRawDArgs.paired != "":
-        filin_paired = open(inRawDArgs.paired)
+        #filin_paired = open(inRawDArgs.paired)
+        filin_paired = gzip.open(inRawDArgs.paired, "rt") if inRawDArgs.paired.endswith(".gz") else open(inRawDArgs.paired)
         line_paired  = filin_paired.readline()
     count_line_tmp=0
     while line and count_line!=count_line_tmp: