diff --git a/_modules/functions_PhageTerm.py b/_modules/functions_PhageTerm.py index ea44ef4835d0552b57488ef30d3e74f7defdd6b4..3610616415a62ed5a7d419bfb8c6ecb2ab8dbd9f 100644 --- a/_modules/functions_PhageTerm.py +++ b/_modules/functions_PhageTerm.py @@ -57,7 +57,8 @@ from reportlab.lib.units import inch from reportlab.lib import colors from reportlab.lib.utils import ImageReader -from utilities import reverseComplement,changeCase,hybridCoverage,applyCoverage,correctEdge +from utilities import reverseComplement,hybridCoverage,applyCoverage,correctEdge +from common_readsCoverage_processing import picMax diff --git a/_modules/main_utils.py b/_modules/main_utils.py index 7eba44e34f78e51a5c1401591cb638e4c7d112ce..c3542e9b787e5c6fbc1dd111bce171cca1c157ed 100644 --- a/_modules/main_utils.py +++ b/_modules/main_utils.py @@ -255,9 +255,19 @@ class InputDerivedDataArgs: if fparms.virome: self.refseq_liste, self.refseq_name, refseq_rejected = ["N" * self.mean_virome], ["Test_virome"], 0 # HOST sequence recovery + # # HOST sequence recovery + # if len(refseq_liste) == 1 and host != "": + # hostseq = genomeFastaRecovery(host, limit_reference, edge, 1) + # if len(hostseq[0]) != 0 and len(hostseq[0]) < len(refseq_liste[0]): + # print "\nHost length < Phage length : removing host sequence." + # hostseq = "" + # else: + # hostseq = "" + # if len(refseq_liste) > 1: + # print "\nWARNING: Host analysis impossible with multiple fasta input\n" if len(self.refseq_liste) == 1 and inputRaw.host != "": self.hostseq = genomeFastaRecovery(inputRaw.host, fparms.limit_reference, fparms.edge, 1) - if len(self.hostseq[0]) != 0 and len(self.hostseq[0]) < len(self.refseq_liste[0]): + if len(self.hostseq[0]) != 0 and len(self.hostseq[0]) > len(self.refseq_liste[0]): print "\nHost length < Phage length : removing host sequence." self.hostseq = "" else: