[11:48:51] Command: --force --outdir /home/yquentin/work/wolbachia/prokka/wLsig --addgenes --prefix wLsig --locustag wLsig --species Wolbachia endosymbiont --strain wLsig --compliant /home/formation/public_html/M2_Phylogenomique/data/Wolbachia/Additional/Litomosoides_sigmodontis_wolbachia_2.0.fna
[11:48:51] Looking for 'aragorn' - found /usr/local/bioinfo/bin/aragorn
[11:48:51] Determined aragorn version is 1.2
[11:48:51] Looking for 'barrnap' - found /usr/local/bioinfo/bin/barrnap
[11:48:51] Determined barrnap version is 0.5
[11:48:51] Looking for 'blastp' - found /usr/local/bioinfo/src/NCBI_Blast+/ncbi-blast-2.6.0+/bin/blastp
[11:48:51] Determined blastp version is 2.6
[11:48:51] Looking for 'cmpress' - found /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../binaries/linux/cmpress
[11:48:51] Determined cmpress version is 1.1
[11:48:51] Looking for 'cmscan' - found /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../binaries/linux/cmscan
[11:48:51] Determined cmscan version is 1.1
[11:48:51] Looking for 'egrep' - found /bin/egrep
[11:48:51] Looking for 'find' - found /bin/find
[11:48:51] Looking for 'grep' - found /bin/grep
[11:48:51] Looking for 'hmmpress' - found /usr/local/bioinfo/bin/hmmpress
[11:48:51] Determined hmmpress version is 3.1
[11:48:51] Looking for 'hmmscan' - found /usr/local/bioinfo/bin/hmmscan
[11:48:51] Determined hmmscan version is 3.1
[11:48:51] Looking for 'less' - found /usr/bin/less
[11:48:51] Looking for 'makeblastdb' - found /usr/local/bioinfo/src/NCBI_Blast+/ncbi-blast-2.6.0+/bin/makeblastdb
[11:48:51] Determined makeblastdb version is 2.6
[11:48:51] Looking for 'minced' - found /usr/local/bioinfo/bin/minced
[11:48:51] Determined minced version is 1.6
[11:48:51] Looking for 'parallel' - found /usr/local/bioinfo/bin/parallel
[11:49:00] Determined parallel version is 20141022
[11:49:00] Looking for 'prodigal' - found /usr/local/bioinfo/bin/prodigal
[11:49:00] Determined prodigal version is 2.6
[11:49:00] Looking for 'rnammer' - found /usr/local/bioinfo/bin/rnammer
[11:49:00] Determined rnammer version is 1.2
[11:49:00] Looking for 'sed' - found /bin/sed
[11:49:00] Looking for 'signalp' - found /usr/local/bioinfo/bin/signalp
[11:49:00] Determined signalp version is 4.1
[11:49:00] Looking for 'tbl2asn' - found /home/yquentin/save/scripts/tbl2asn
[11:49:00] Determined tbl2asn version is 25.3
[11:49:00] Using genetic code table 11.
[11:49:00] Loading and checking input file: /home/formation/public_html/M2_Phylogenomique/data/Wolbachia/Additional/Litomosoides_sigmodontis_wolbachia_2.0.fna
[11:49:01] Wrote 10 contigs
[11:49:01] Predicting tRNAs and tmRNAs
[11:49:01] Running: aragorn -gc11  -w \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/wLsig\.fna
[11:49:01] 1 tRNA-Ile [12815,12889] 35 (gat)
[11:49:01] 2 tRNA-Arg c[40168,40244] 36 (cct)
[11:49:01] 3 tRNA-Glu c[49883,49956] 35 (ttc)
[11:49:01] 4 tRNA-Leu [63875,63962] 35 (gag)
[11:49:01] 5 tRNA-Phe c[140023,140097] 35 (gaa)
[11:49:01] 6 tRNA-Asp [145565,145639] 35 (gtc)
[11:49:01] 7 tRNA-Leu c[146732,146819] 35 (taa)
[11:49:01] 8 tRNA-Met [201741,201816] 34 (cat)
[11:49:01] 9 tRNA-Arg c[219945,220020] 36 (acg)
[11:49:01] 10 tRNA-Gln c[245244,245318] 33 (ttg)
[11:49:01] 11 tRNA-Cys c[317150,317222] 34 (gca)
[11:49:01] 12 tRNA-Ser c[326800,326888] 36 (gga)
[11:49:01] 13 tmRNA* c[348076,348411] 198,245 ANDNFAAENDVAVAA*
[11:49:01] 14 tRNA-Trp c[361054,361127] 34 (cca)
[11:49:01] 15 tRNA-Gly c[362401,362473] 34 (tcc)
[11:49:01] 16 tRNA-Tyr c[362477,362560] 36 (gta)
[11:49:01] 17 tRNA-Ser c[384212,384300] 35 (tga)
[11:49:01] 18 tRNA-Met c[415942,416017] 35 (cat)
[11:49:01] 19 tRNA-His [422796,422871] 34 (gtg)
[11:49:01] 20 tRNA-Gly c[502590,502663] 34 (gcc)
[11:49:01] 21 tRNA-Arg [538382,538456] 35 (ccg)
[11:49:01] 22 tRNA-Lys c[577422,577495] 34 (ttt)
[11:49:01] 1 tRNA-Ser [48439,48524] 35 (cag)
[11:49:01] 2 tRNA-Pro c[58308,58383] 36 (tgg)
[11:49:01] 3 tRNA-Val [98888,98960] 34 (tac)
[11:49:01] 4 tRNA-Arg [137740,137814] 35 (tct)
[11:49:01] 5 tRNA-Met [236616,236690] 35 (cat)
[11:49:01] 1 tRNA-Leu [5224,5307] 34 (caa)
[11:49:01] 2 tRNA-Leu [79501,79583] 35 (tag)
[11:49:01] 3 tRNA-Ser c[109906,109997] 35 (gct)
[11:49:01] 4 tRNA-Thr c[115858,115931] 34 (ggt)
[11:49:01] 5 tRNA-Val [128347,128420] 34 (gac)
[11:49:01] 6 tRNA-Asn c[131454,131526] 33 (gtt)
[11:49:01] 7 tRNA-Thr c[133391,133463] 34 (tgt)
[11:49:01] 1 tRNA-Ala [21593,21666] 34 (tgc)
[11:49:01] Found 35 tRNAs
[11:49:01] Predicting Ribosomal RNAs
[11:49:01] Running Barrnap with 8 threads
[11:49:02] 1 gnl|Prokka|wLsig_contig000001 22292 23S ribosomal RNA
[11:49:02] 2 gnl|Prokka|wLsig_contig000001 25150 5S ribosomal RNA (partial)
[11:49:02] 3 gnl|Prokka|wLsig_contig000001 107779 16S ribosomal RNA
[11:49:02] Found 3 rRNAs
[11:49:02] Skipping ncRNA search, enable with --rfam if desired.
[11:49:02] Total of 37 tRNA + rRNA features
[11:49:02] Searching for CRISPR repeats
[11:49:02] Found 0 CRISPRs
[11:49:02] Predicting coding sequences
[11:49:02] Contigs total 1048936 bp, so using single mode
[11:49:02] Running: prodigal -i \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/wLsig\.fna -c -m -g 11 -p single -f sco -q
[11:49:03] Excluding CDS which overlaps existing RNA (rRNA) at gnl|Prokka|wLsig_contig000001:23150..23278 on + strand
[11:49:03] Excluding CDS which overlaps existing RNA (rRNA) at gnl|Prokka|wLsig_contig000001:108925..109176 on + strand
[11:49:03] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wLsig_contig000001:145385..145630 on - strand
[11:49:03] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wLsig_contig000001:200446..201741 on + strand
[11:49:03] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wLsig_contig000001:245309..245644 on - strand
[11:49:03] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wLsig_contig000003:79486..79644 on + strand
[11:49:03] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wLsig_contig000003:115821..115940 on - strand
[11:49:03] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wLsig_contig000003:128383..129342 on + strand
[11:49:03] Found 880 CDS
[11:49:03] Connecting features back to sequences
[11:49:03] Option --gram not specified, will NOT check for signal peptides.
[11:49:03] Not using genus-specific database. Try --usegenus to enable it.
[11:49:03] Annotating CDS, please be patient.
[11:49:03] Will use 8 CPUs for similarity searching.
[11:49:04] There are still 880 unannotated CDS left (started with 880)
[11:49:04] Will use blast to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/kingdom/Bacteria/sprot with 8 CPUs
[11:49:04] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/proteins\.faa | parallel --gnu -j 8 --block 15090 --recstart '>' --pipe blastp -query - -db /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/kingdom/Bacteria/sprot -evalue 1e-06 -num_threads 1 -num_descriptions 1 -num_alignments 1 -seg no > \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/proteins\.bls 2> /dev/null
[11:49:26] Modify product: GTP-binding protein TypA/BipA homolog => hypothetical protein
[11:49:26] Modify product: Uncharacterized zinc protease Rv2782c => putative zinc protease
[11:49:26] Modify product: Transcription termination/antitermination protein NusG => hypothetical protein
[11:49:27] Modify product: Uncharacterized zinc protease Rv2782c => putative zinc protease
[11:49:27] Modify product: Probable dipeptidase PepE => putative dipeptidase PepE
[11:49:27] Modify product: N utilization substance protein B homolog => hypothetical protein
[11:49:27] Modify product: Uncharacterized ABC transporter ATP-binding protein TM_0288 => putative ABC transporter ATP-binding protein
[11:49:27] Modify product: Uncharacterized oxidoreductase YciK => putative oxidoreductase YciK
[11:49:27] Modify product: Probable transcriptional regulatory protein aq_1575 => putative transcriptional regulatory protein
[11:49:27] Modify product: SCO1 protein homolog => hypothetical protein
[11:49:27] Modify product: Cytochrome b561 homolog 2 => hypothetical protein
[11:49:27] Modify product: Transcription termination factor Rho => hypothetical protein
[11:49:27] Modify product: Bifunctional protein FolD => Bifunctional protein FolD protein
[11:49:27] Modify product: Probable UbiX-like flavin prenyltransferase => putative UbiX-like flavin prenyltransferase
[11:49:27] Modify product: Stringent starvation protein A homolog => hypothetical protein
[11:49:27] Modify product: Putative O-methyltransferase MSMEG_5073/MSMEI_4947 => Putative O-methyltransferase/MSMEI_4947
[11:49:27] Modify product: Probable tRNA-dihydrouridine synthase => putative tRNA-dihydrouridine synthase
[11:49:27] Modify product: Transcription termination/antitermination protein NusA => hypothetical protein
[11:49:27] Modify product: Probable murein peptide carboxypeptidase => putative murein peptide carboxypeptidase
[11:49:27] Modify product: Co-chaperone protein HscB homolog => hypothetical protein
[11:49:27] Modify product: Uncharacterized HIT-like protein HP_0404 => putative HIT-like protein
[11:49:27] Modify product: Probable GTP-binding protein EngB => putative GTP-binding protein EngB
[11:49:27] Modify product: Probable peroxiredoxin => putative peroxiredoxin
[11:49:27] Modify product: Probable protein kinase UbiB => putative protein kinase UbiB
[11:49:27] Cleaned 24 /product names
[11:49:27] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/proteins.faa
[11:49:27] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/proteins.bls
[11:49:27] There are still 365 unannotated CDS left (started with 880)
[11:49:27] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/HAMAP.hmm with 8 CPUs
[11:49:27] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/proteins\.faa | parallel --gnu -j 8 --block 3592 --recstart '>' --pipe hmmscan --noali --notextw --acc -E 1e-06 --cpu 1 /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/HAMAP.hmm /dev/stdin > \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/proteins\.bls 2> /dev/null
[11:49:39] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/proteins.faa
[11:49:39] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/proteins.bls
[11:49:39] There are still 340 unannotated CDS left (started with 880)
[11:49:39] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/CLUSTERS.hmm with 8 CPUs
[11:49:39] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/proteins\.faa | parallel --gnu -j 8 --block 3211 --recstart '>' --pipe hmmscan --noali --notextw --acc -E 1e-06 --cpu 1 /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/CLUSTERS.hmm /dev/stdin > \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/proteins\.bls 2> /dev/null
[11:50:03] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/proteins.faa
[11:50:03] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/proteins.bls
[11:50:03] There are still 321 unannotated CDS left (started with 880)
[11:50:03] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/Pfam.hmm with 8 CPUs
[11:50:03] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/proteins\.faa | parallel --gnu -j 8 --block 2915 --recstart '>' --pipe hmmscan --noali --notextw --acc -E 1e-06 --cpu 1 /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/Pfam.hmm /dev/stdin > \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/proteins\.bls 2> /dev/null
[11:50:19] Modify product: ETC complex I subunit conserved region => hypothetical protein
[11:50:19] Modify product: Tim44-like domain => Tim44-like domain protein
[11:50:19] Modify product: Peptidase M16 inactive domain => Peptidase M16 inactive domain protein
[11:50:20] Modify product: SPFH domain / Band 7 family => SPFH domain / Band 7 family protein
[11:50:20] Modify product: Zinc-finger domain => Zinc-finger domain protein
[11:50:20] Modify product: PD-(D/E)XK nuclease superfamily => PD-(D/E)XK nuclease superfamily protein
[11:50:20] Modify product: YGGT family => YGGT family protein
[11:50:20] Modify product: PAS domain => PAS domain protein
[11:50:20] Modify product: Uracil DNA glycosylase superfamily => Uracil DNA glycosylase superfamily protein
[11:50:20] Modify product: BioY family => BioY family protein
[11:50:20] Modify product: COQ9 => hypothetical protein
[11:50:20] Modify product: AhpC/TSA family => AhpC/TSA family protein
[11:50:20] Modify product: SURF1 family => SURF1 family protein
[11:50:20] Modify product: RDD family => RDD family protein
[11:50:20] Modify product: Integral membrane protein TerC family => Integral membrane protein TerC family protein
[11:50:20] Modify product: Sodium:alanine symporter family => Sodium:alanine symporter family protein
[11:50:20] Modify product: EVE domain => EVE domain protein
[11:50:20] Modify product: Alpha/beta hydrolase family => Alpha/beta hydrolase family protein
[11:50:20] Modify product: RDD family => RDD family protein
[11:50:20] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein
[11:50:20] Cleaned 20 /product names
[11:50:20] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/proteins.faa
[11:50:20] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/proteins.bls
[11:50:20] Labelling remaining 278 proteins as 'hypothetical protein'
[11:50:20] Possible /pseudo 'tRNA pseudouridine synthase B' at gnl|Prokka|wLsig_contig000001 position 324988
[11:50:20] Possible /pseudo 'ATP-dependent DNA helicase RecG' at gnl|Prokka|wLsig_contig000001 position 335592
[11:50:20] Possible /pseudo 'ATP-dependent DNA helicase RecG' at gnl|Prokka|wLsig_contig000001 position 335916
[11:50:20] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|wLsig_contig000001 position 444592
[11:50:20] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|wLsig_contig000001 position 447681
[11:50:20] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|wLsig_contig000001 position 450152
[11:50:20] Possible /pseudo 'Enolase' at gnl|Prokka|wLsig_contig000001 position 461494
[11:50:20] Possible /pseudo 'Enolase' at gnl|Prokka|wLsig_contig000001 position 461850
[11:50:20] Possible /pseudo 'GTPase Obg' at gnl|Prokka|wLsig_contig000001 position 462808
[11:50:20] Possible /pseudo 'ATP synthase subunit b' at gnl|Prokka|wLsig_contig000001 position 515291
[11:50:20] Possible /pseudo 'Arginine/agmatine antiporter' at gnl|Prokka|wLsig_contig000002 position 154093
[11:50:20] Possible /pseudo 'Elongation factor Tu' at gnl|Prokka|wLsig_contig000003 position 38721
[11:50:20] Possible /pseudo 'Chaperone protein HscA' at gnl|Prokka|wLsig_contig000007 position 943
[11:50:20] Found 487 unique /gene codes.
[11:50:20] Fixed 19 colliding /gene names.
[11:50:20] Adding /locus_tag identifiers
[11:50:20] Assigned 918 locus_tags to CDS and RNA features.
[11:50:20] Writing outputs to /home/yquentin/work/wolbachia/prokka/wLsig/
[11:50:21] Generating annotation statistics file
[11:50:21] Generating Genbank and Sequin files
[11:50:21] Running: tbl2asn -V b -a r10k -l paired-ends -M n -N 1 -y 'Annotated using prokka 1.10 from http://www.vicbioinformatics.com' -Z \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/wLsig\.err -i \/home\/yquentin\/work\/wolbachia\/prokka\/wLsig\/wLsig\.fsa 2> /dev/null
[11:50:23] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/errorsummary.val
[11:50:23] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.dr
[11:50:23] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.fixedproducts
[11:50:23] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.ecn
[11:50:23] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.val
[11:50:23] Output files:
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.gff
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.faa
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.tbl
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.fsa
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.err
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.log
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.ffn
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.gbk
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.sqn
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.fna
[11:50:23] /home/yquentin/work/wolbachia/prokka/wLsig/wLsig.txt
[11:50:23] Walltime used: 1.53 minutes
[11:50:23] If you use this result please cite the Prokka paper:
[11:50:23] Seemann T (2014) Prokka: rapid prokaryotic genome annotation. Bioinformatics. 30(14):2068-9.
[11:50:23] Type 'prokka --citation' for more details.
[11:50:23] Thank you, come again.