Thanks for your quick answer. Before receiving it, I tried using the Maximilian Press 's version compatible with python3. The process went until the end, i.e without any error messages but I reckon some of the outputs were incorrect : the corrected fasta was empty for all the tested phages (~70). However, I got some seemingly good reports : I compared one of them with one obtained on the galaxy platform of the institute and they were identical. However I am still very interested in the corrected sequence.
After reading you reply I tried, as you recommended me, and I got the same error message. Nevertheless, I got for some phages a corrected fasta sequence without reports. (Seems congruent with the following error message) :
Calculating coverage values, please wait (may take a while)...
100.0 %
Finished calculating coverage values, the remainder should be completed rapidly
Bonjour !
I tried running the script on a recently sequenced phage. The input are paired reads, the command line is as follow :
After calculating coverage values until 100%, I get the following error message :
I tried running the script with the test data as follow :
/home/user/software/phageterm-1.0.12/PhageTerm.py -f Mu-like_R1.fastq -p Mu-like_R2.fastq -r Mu-like.fasta
but I get the same error.
Is there a way this can be fixed ?
Best,
Robby
Bonsoir :)
When you run the installation test, could try without your own reads, as follow :
$ /home/user/software/phageterm-1.0.12/PhageTerm.py -t C5
Tell me if this command works properly on your system. If not, maybe you should force launching phageterm with python2 :
$ python2 PhageTerm.py -t C5
And with your own reads :
$ python2 PhageTerm.py -f K10PH82C1_S15_R1_001.fastq -p K10PH82C1_S15_R2_001.fastq -r K10PH82C1_sequence.fasta -n K10PH82C1 -c 4
Keep me posted if it worked with one of those ways to run the software.
Good luck
P.S. Do you expect your phage to be a Mu-like phage?
Hi again,
Thanks for your quick answer. Before receiving it, I tried using the Maximilian Press 's version compatible with python3. The process went until the end, i.e without any error messages but I reckon some of the outputs were incorrect : the corrected fasta was empty for all the tested phages (~70). However, I got some seemingly good reports : I compared one of them with one obtained on the galaxy platform of the institute and they were identical. However I am still very interested in the corrected sequence.
After reading you reply I tried, as you recommended me, and I got the same error message. Nevertheless, I got for some phages a corrected fasta sequence without reports. (Seems congruent with the following error message) :
Calculating coverage values, please wait (may take a while)...
100.0 %
Finished calculating coverage values, the remainder should be completed rapidly
I guess I could combine the output from those 2 different sources.
Thanks for your help and congrats for the work done.
Ps : I actually had no expectation regarding the phage packaging mechanism.
Hi again Concha,
Could you send me an e-mail at julian.garneau@pasteur.fr
WIll be easier to exchange on there,
Thanks!