[15:06:54] Command: --force --outdir /home/yquentin/work/wolbachia/prokka/wMun --addgenes --prefix wMun --locustag wMun --species Wolbachia endosymbiont --strain wMun --compliant --quiet /home/yquentin/work/wolbachia/prokka/wMun.fna
[15:06:54] Looking for 'aragorn' - found /usr/local/bioinfo/bin/aragorn
[15:06:54] Determined aragorn version is 1.2
[15:06:54] Looking for 'barrnap' - found /usr/local/bioinfo/bin/barrnap
[15:06:54] Determined barrnap version is 0.5
[15:06:54] Looking for 'blastp' - found /usr/local/bioinfo/src/NCBI_Blast+/ncbi-blast-2.6.0+/bin/blastp
[15:06:54] Determined blastp version is 2.6
[15:06:54] Looking for 'cmpress' - found /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../binaries/linux/cmpress
[15:06:54] Determined cmpress version is 1.1
[15:06:54] Looking for 'cmscan' - found /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../binaries/linux/cmscan
[15:06:54] Determined cmscan version is 1.1
[15:06:54] Looking for 'egrep' - found /bin/egrep
[15:06:54] Looking for 'find' - found /bin/find
[15:06:54] Looking for 'grep' - found /bin/grep
[15:06:54] Looking for 'hmmpress' - found /usr/local/bioinfo/bin/hmmpress
[15:06:54] Determined hmmpress version is 3.1
[15:06:54] Looking for 'hmmscan' - found /usr/local/bioinfo/bin/hmmscan
[15:06:54] Determined hmmscan version is 3.1
[15:06:54] Looking for 'less' - found /usr/bin/less
[15:06:54] Looking for 'makeblastdb' - found /usr/local/bioinfo/src/NCBI_Blast+/ncbi-blast-2.6.0+/bin/makeblastdb
[15:06:55] Determined makeblastdb version is 2.6
[15:06:55] Looking for 'minced' - found /usr/local/bioinfo/bin/minced
[15:06:55] Determined minced version is 1.6
[15:06:55] Looking for 'parallel' - found /usr/local/bioinfo/bin/parallel
[15:06:55] Determined parallel version is 20141022
[15:06:55] Looking for 'prodigal' - found /usr/local/bioinfo/bin/prodigal
[15:06:55] Determined prodigal version is 2.6
[15:06:55] Looking for 'rnammer' - found /usr/local/bioinfo/bin/rnammer
[15:06:55] Determined rnammer version is 1.2
[15:06:55] Looking for 'sed' - found /bin/sed
[15:06:55] Looking for 'signalp' - found /usr/local/bioinfo/bin/signalp
[15:06:56] Determined signalp version is 4.1
[15:06:56] Looking for 'tbl2asn' - found /home/yquentin/save/scripts/tbl2asn
[15:06:56] Determined tbl2asn version is 25.3
[15:06:56] Using genetic code table 11.
[15:06:56] Loading and checking input file: /home/yquentin/work/wolbachia/prokka/wMun.fna
[15:06:56] Wrote 256 contigs
[15:06:56] Predicting tRNAs and tmRNAs
[15:06:56] Running: aragorn -gc11  -w \/home\/yquentin\/work\/wolbachia\/prokka\/wMun\/wMun\.fna
[15:06:56] 1 tRNA-Asp c[5832,5908] 35 (gtc)
[15:06:56] 1 tRNA-Leu [7097,7180] 34 (caa)
[15:06:56] 1 tRNA-Cys c[4383,4455] 34 (gca)
[15:06:56] 1 tRNA-Glu [8676,8748] 34 (ttc)
[15:06:56] 2 tRNA-Pro [9680,9755] 36 (tgg)
[15:06:56] 3 tRNA-Val [16471,16543] 34 (tac)
[15:06:56] 1 tRNA-Met [2733,2806] 34 (cat)
[15:06:56] 1 tRNA-Leu [13019,13106] 35 (gag)
[15:06:56] 1 tRNA-Thr [5855,5928] 34 (tgt)
[15:06:56] 1 tRNA-Lys c[2725,2798] 34 (ttt)
[15:06:56] 2 tRNA-Ala [3140,3215] 34 (tgc)
[15:06:56] 1 tRNA-Ile c[2730,2806] 35 (gat)
[15:06:56] 1 tRNA-Trp c[8154,8227] 34 (cca)
[15:06:56] 1 tRNA-Ser [2802,2887] 35 (cag)
[15:06:56] 1 tRNA-Thr [599,672] 34 (ggt)
[15:06:56] 1 tRNA-Gln c[768,842] 33 (ttg)
[15:06:56] Found 16 tRNAs
[15:06:56] Predicting Ribosomal RNAs
[15:06:56] Running Barrnap with 8 threads
[15:06:57] Found 0 rRNAs
[15:06:57] Skipping ncRNA search, enable with --rfam if desired.
[15:06:57] Total of 16 tRNA + rRNA features
[15:06:57] Searching for CRISPR repeats
[15:06:58] Found 0 CRISPRs
[15:06:58] Predicting coding sequences
[15:06:58] Contigs total 867873 bp, so using single mode
[15:06:58] Running: prodigal -i \/home\/yquentin\/work\/wolbachia\/prokka\/wMun\/wMun\.fna -c -m -g 11 -p single -f sco -q
[15:07:00] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wMun_contig000034:1441..2733 on + strand
[15:07:00] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wMun_contig000053:1453..2757 on - strand
[15:07:00] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wMun_contig000186:125..706 on + strand
[15:07:00] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wMun_contig000202:833..1168 on - strand
[15:07:00] Found 840 CDS
[15:07:00] Connecting features back to sequences
[15:07:00] Option --gram not specified, will NOT check for signal peptides.
[15:07:00] Not using genus-specific database. Try --usegenus to enable it.
[15:07:00] Annotating CDS, please be patient.
[15:07:00] Will use 8 CPUs for similarity searching.
[15:07:01] There are still 840 unannotated CDS left (started with 840)
[15:07:01] Will use blast to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/kingdom/Bacteria/sprot with 8 CPUs
[15:07:01] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wMun\/proteins\.faa | parallel --gnu -j 8 --block 12919 --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\/wMun\/proteins\.bls 2> /dev/null
[15:07:13] Modify product: Probable transcriptional regulatory protein aq_1575 => putative transcriptional regulatory protein
[15:07:13] Modify product: Probable phospholipid ABC transporter-binding protein MlaD => putative phospholipid ABC transporter-binding protein MlaD
[15:07:13] Modify product: Bifunctional protein FolD => Bifunctional protein FolD protein
[15:07:13] Modify product: Actin-binding protein Smlt3054 => Actin-binding protein
[15:07:13] Modify product: Uncharacterized oxidoreductase YciK => putative oxidoreductase YciK
[15:07:13] Modify product: Transcription termination factor Rho => hypothetical protein
[15:07:13] Modify product: Transcription termination/antitermination protein NusA => hypothetical protein
[15:07:13] Modify product: SCO1 protein homolog => hypothetical protein
[15:07:13] Modify product: Probable phospholipid ABC transporter permease protein MlaE => putative phospholipid ABC transporter permease protein MlaE
[15:07:13] Modify product: Probable ribonucleotide transport ATP-binding protein mkl => putative ribonucleotide transport ATP-binding protein mkl
[15:07:13] Modify product: Probable chromosome-partitioning protein ParB => putative chromosome-partitioning protein ParB
[15:07:14] Modify product: Co-chaperone protein HscB homolog => hypothetical protein
[15:07:14] Modify product: Probable peroxiredoxin => putative peroxiredoxin
[15:07:14] Modify product: Probable protein kinase UbiB => putative protein kinase UbiB
[15:07:14] Modify product: Probable GTP-binding protein EngB => putative GTP-binding protein EngB
[15:07:14] Modify product: Probable tRNA-dihydrouridine synthase => putative tRNA-dihydrouridine synthase
[15:07:14] Modify product: Uncharacterized metal-dependent hydrolase YcfH => putative metal-dependent hydrolase YcfH
[15:07:14] Modify product: Uncharacterized HIT-like protein HP_0404 => putative HIT-like protein
[15:07:14] Modify product: Transcription termination/antitermination protein NusG => hypothetical protein
[15:07:14] Modify product: Probable acyltransferase YihG => putative acyltransferase YihG
[15:07:15] Modify product: Putative multidrug export ATP-binding/permease protein SAV1866 => Putative multidrug export ATP-binding/permease protein
[15:07:15] Modify product: GTP-binding protein TypA/BipA homolog => hypothetical protein
[15:07:15] Modify product: Actin-binding protein Smlt3054 => Actin-binding protein
[15:07:15] Cleaned 23 /product names
[15:07:15] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/proteins.faa
[15:07:15] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/proteins.bls
[15:07:15] There are still 408 unannotated CDS left (started with 840)
[15:07:15] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/HAMAP.hmm with 8 CPUs
[15:07:15] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wMun\/proteins\.faa | parallel --gnu -j 8 --block 4473 --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\/wMun\/proteins\.bls 2> /dev/null
[15:07:25] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/proteins.faa
[15:07:25] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/proteins.bls
[15:07:25] There are still 378 unannotated CDS left (started with 840)
[15:07:25] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/CLUSTERS.hmm with 8 CPUs
[15:07:25] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wMun\/proteins\.faa | parallel --gnu -j 8 --block 3908 --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\/wMun\/proteins\.bls 2> /dev/null
[15:07:48] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/proteins.faa
[15:07:48] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/proteins.bls
[15:07:48] There are still 364 unannotated CDS left (started with 840)
[15:07:48] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/Pfam.hmm with 8 CPUs
[15:07:48] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wMun\/proteins\.faa | parallel --gnu -j 8 --block 3724 --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\/wMun\/proteins\.bls 2> /dev/null
[15:08:13] Modify product: Bacterial SH3 domain => Bacterial SH3 domain protein
[15:08:13] Modify product: KAP family P-loop domain => KAP family P-loop domain protein
[15:08:13] Modify product: Tetratricopeptide repeat => Tetratricopeptide repeat protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: SPFH domain / Band 7 family => SPFH domain / Band 7 family protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Piwi/Argonaute/Zwille siRNA-binding domain => Piwi/Argonaute/Zwille siRNA-binding domain protein
[15:08:13] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein
[15:08:13] Modify product: Sodium:dicarboxylate symporter family => Sodium:dicarboxylate symporter family protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Ankyrin repeat => Ankyrin repeat protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Transposase IS116/IS110/IS902 family => Transposase IS116/IS110/IS902 family protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: RDD family => RDD family protein
[15:08:13] Modify product: RDD family => RDD family protein
[15:08:13] Modify product: RDD family => RDD family protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Tim44-like domain => Tim44-like domain protein
[15:08:13] Modify product: Uracil DNA glycosylase superfamily => Uracil DNA glycosylase superfamily protein
[15:08:13] Modify product: ETC complex I subunit conserved region => hypothetical protein
[15:08:13] Modify product: Integrase core domain => Integrase core domain protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Acetyltransferase (GNAT) family => Acetyltransferase (GNAT) family protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: MgtE intracellular N domain => MgtE intracellular N domain protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Alpha/beta hydrolase family => Alpha/beta hydrolase family protein
[15:08:13] Modify product: TENA/THI-4/PQQC family => TENA/THI-4/PQQC family protein
[15:08:13] Modify product: Transposase IS116/IS110/IS902 family => Transposase IS116/IS110/IS902 family protein
[15:08:13] Modify product: Transposase IS200 like => Transposase IS200 like protein
[15:08:13] Modify product: Integral membrane protein TerC family => Integral membrane protein TerC family protein
[15:08:13] Modify product: YGGT family => YGGT family protein
[15:08:13] Modify product: Transposase IS116/IS110/IS902 family => Transposase IS116/IS110/IS902 family protein
[15:08:13] Modify product: Transposase DDE domain => Transposase DDE domain protein
[15:08:13] Modify product: PAS domain => PAS domain protein
[15:08:13] Modify product: Zinc-finger domain => Zinc-finger domain protein
[15:08:13] Modify product: COQ9 => hypothetical protein
[15:08:13] Cleaned 60 /product names
[15:08:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/proteins.faa
[15:08:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/proteins.bls
[15:08:13] Labelling remaining 255 proteins as 'hypothetical protein'
[15:08:13] Possible /pseudo 'Multidrug resistance protein MdtL' at gnl|Prokka|wMun_contig000003 position 5869
[15:08:13] Possible /pseudo 'PD-(D/E)XK nuclease family transposase' at gnl|Prokka|wMun_contig000006 position 2708
[15:08:13] Possible /pseudo 'Transposase IS200 like protein' at gnl|Prokka|wMun_contig000007 position 4774
[15:08:13] Possible /pseudo 'Arginine--tRNA ligase' at gnl|Prokka|wMun_contig000031 position 39654
[15:08:13] Possible /pseudo 'Arginine--tRNA ligase' at gnl|Prokka|wMun_contig000031 position 40110
[15:08:13] Possible /pseudo 'RDD family protein' at gnl|Prokka|wMun_contig000034 position 8411
[15:08:13] Possible /pseudo 'Quaternary ammonium compound-resistance protein SugE' at gnl|Prokka|wMun_contig000037 position 13569
[15:08:13] Possible /pseudo 'Dephospho-CoA kinase' at gnl|Prokka|wMun_contig000044 position 1944
[15:08:13] Possible /pseudo 'Membrane protein insertase YidC' at gnl|Prokka|wMun_contig000047 position 1044
[15:08:13] Possible /pseudo 'Threonine--tRNA ligase' at gnl|Prokka|wMun_contig000053 position 684
[15:08:13] Possible /pseudo 'Transposase IS200 like protein' at gnl|Prokka|wMun_contig000055 position 5439
[15:08:13] Found 412 unique /gene codes.
[15:08:13] Fixed 18 colliding /gene names.
[15:08:13] Adding /locus_tag identifiers
[15:08:13] Assigned 856 locus_tags to CDS and RNA features.
[15:08:13] Writing outputs to /home/yquentin/work/wolbachia/prokka/wMun/
[15:08:15] Generating annotation statistics file
[15:08:15] Generating Genbank and Sequin files
[15:08:15] 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\/wMun\/wMun\.err -i \/home\/yquentin\/work\/wolbachia\/prokka\/wMun\/wMun\.fsa 2> /dev/null
[15:08:19] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/errorsummary.val
[15:08:19] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/wMun.dr
[15:08:19] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/wMun.fixedproducts
[15:08:19] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/wMun.ecn
[15:08:19] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wMun/wMun.val
[15:08:19] Output files:
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.tbl
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.gbk
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.log
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.sqn
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.fsa
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.err
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.txt
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.gff
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.ffn
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.fna
[15:08:19] /home/yquentin/work/wolbachia/prokka/wMun/wMun.faa
[15:08:19] Walltime used: 1.43 minutes
[15:08:19] If you use this result please cite the Prokka paper:
[15:08:19] Seemann T (2014) Prokka: rapid prokaryotic genome annotation. Bioinformatics. 30(14):2068-9.
[15:08:19] Type 'prokka --citation' for more details.
[15:08:19] Thank you, come again.