[15:06:39] Command: --force --outdir /home/yquentin/work/wolbachia/prokka/wDci --addgenes --prefix wDci --locustag wDci --species Wolbachia endosymbiont --strain wDci --compliant --quiet /home/yquentin/work/wolbachia/prokka/wDci.fna [15:06:39] Looking for 'aragorn' - found /usr/local/bioinfo/bin/aragorn [15:06:39] Determined aragorn version is 1.2 [15:06:39] Looking for 'barrnap' - found /usr/local/bioinfo/bin/barrnap [15:06:39] Determined barrnap version is 0.5 [15:06:39] Looking for 'blastp' - found /usr/local/bioinfo/src/NCBI_Blast+/ncbi-blast-2.6.0+/bin/blastp [15:06:39] Determined blastp version is 2.6 [15:06:39] Looking for 'cmpress' - found /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../binaries/linux/cmpress [15:06:40] Determined cmpress version is 1.1 [15:06:40] Looking for 'cmscan' - found /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../binaries/linux/cmscan [15:06:40] Determined cmscan version is 1.1 [15:06:40] Looking for 'egrep' - found /bin/egrep [15:06:40] Looking for 'find' - found /bin/find [15:06:40] Looking for 'grep' - found /bin/grep [15:06:40] Looking for 'hmmpress' - found /usr/local/bioinfo/bin/hmmpress [15:06:40] Determined hmmpress version is 3.1 [15:06:40] Looking for 'hmmscan' - found /usr/local/bioinfo/bin/hmmscan [15:06:40] Determined hmmscan version is 3.1 [15:06:40] Looking for 'less' - found /usr/bin/less [15:06:40] Looking for 'makeblastdb' - found /usr/local/bioinfo/src/NCBI_Blast+/ncbi-blast-2.6.0+/bin/makeblastdb [15:06:40] Determined makeblastdb version is 2.6 [15:06:40] Looking for 'minced' - found /usr/local/bioinfo/bin/minced [15:06:40] Determined minced version is 1.6 [15:06:40] Looking for 'parallel' - found /usr/local/bioinfo/bin/parallel [15:06:40] Determined parallel version is 20141022 [15:06:41] Looking for 'prodigal' - found /usr/local/bioinfo/bin/prodigal [15:06:41] Determined prodigal version is 2.6 [15:06:41] Looking for 'rnammer' - found /usr/local/bioinfo/bin/rnammer [15:06:41] Determined rnammer version is 1.2 [15:06:41] Looking for 'sed' - found /bin/sed [15:06:41] Looking for 'signalp' - found /usr/local/bioinfo/bin/signalp [15:06:41] Determined signalp version is 4.1 [15:06:41] Looking for 'tbl2asn' - found /home/yquentin/save/scripts/tbl2asn [15:06:42] Determined tbl2asn version is 25.3 [15:06:42] Using genetic code table 11. [15:06:42] Loading and checking input file: /home/yquentin/work/wolbachia/prokka/wDci.fna [15:06:42] Wrote 104 contigs [15:06:42] Predicting tRNAs and tmRNAs [15:06:42] Running: aragorn -gc11 -w \/home\/yquentin\/work\/wolbachia\/prokka\/wDci\/wDci\.fna [15:06:42] 1 tRNA-Phe c[12331,12404] 34 (gaa) [15:06:42] 1 tRNA-Arg c[6947,7021] 35 (cct) [15:06:42] 1 tRNA-Thr [11633,11706] 34 (tgt) [15:06:42] 1 tRNA-Trp c[43022,43095] 34 (cca) [15:06:42] 2 tRNA-Leu [50240,50322] 35 (tag) [15:06:42] 1 tRNA-His c[18117,18192] 34 (gtg) [15:06:42] 1 tRNA-Met c[5937,6009] 33 (cat) [15:06:42] 1 tRNA-Asn [5681,5753] 33 (gtt) [15:06:42] 1 tRNA-Met c[6577,6653] 35 (cat) [15:06:42] 2 tRNA-Gly [7874,7947] 34 (gcc) [15:06:42] 3 tRNA-Ser [15066,15153] 35 (gga) [15:06:42] 1 tRNA-Thr [2318,2391] 34 (ggt) [15:06:42] 2 tRNA-Met c[20234,20310] 35 (cat) [15:06:42] 1 tRNA-Ser [2122,2207] 35 (cag) [15:06:42] 2 tRNA-Ser [49814,49906] 36 (gct) [15:06:42] 1 tRNA-Glu [14738,14810] 34 (ttc) [15:06:42] 2 tRNA-Pro [15659,15734] 36 (tgg) [15:06:42] 3 tRNA-Val [22788,22860] 34 (tac) [15:06:42] 4 tRNA-Arg [33634,33709] 36 (ccg) [15:06:42] 1 tRNA-Ser c[12235,12323] 35 (tga) [15:06:42] 2 tRNA-Lys c[18082,18155] 34 (ttt) [15:06:42] 1 tRNA-Leu c[26279,26362] 34 (caa) [15:06:42] 1 tRNA-Val [13342,13415] 34 (gac) [15:06:42] 2 tRNA-Ile [47683,47759] 35 (gat) [15:06:42] 3 tRNA-Tyr [70448,70532] 36 (gta) [15:06:43] 4 tRNA-Gly [70536,70607] 33 (tcc) [15:06:43] 5 tRNA-Asp c[124947,125023] 35 (gtc) [15:06:43] 6 tRNA-Leu [134910,134997] 35 (gag) [15:06:43] 7 tRNA-Leu c[177883,177969] 35 (taa) [15:06:43] 1 tmRNA* c[60913,61245] 195,242 ANDNFAAEDNVALAA* [15:06:43] 2 tRNA-Arg [81722,81799] 36 (acg) [15:06:43] Found 31 tRNAs [15:06:43] Predicting Ribosomal RNAs [15:06:43] Running Barrnap with 8 threads [15:06:44] 1 gnl|Prokka|wDci_contig000060 1499 16S ribosomal RNA [15:06:44] 2 gnl|Prokka|wDci_contig000062 55548 23S ribosomal RNA [15:06:44] 3 gnl|Prokka|wDci_contig000062 58378 5S ribosomal RNA [15:06:44] Found 3 rRNAs [15:06:44] Skipping ncRNA search, enable with --rfam if desired. [15:06:44] Total of 33 tRNA + rRNA features [15:06:44] Searching for CRISPR repeats [15:06:44] Found 0 CRISPRs [15:06:44] Predicting coding sequences [15:06:44] Contigs total 1240904 bp, so using single mode [15:06:44] Running: prodigal -i \/home\/yquentin\/work\/wolbachia\/prokka\/wDci\/wDci\.fna -c -m -g 11 -p single -f sco -q [15:06:46] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wDci_contig000042:5885..6013 on - strand [15:06:46] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wDci_contig000042:6000..7301 on - strand [15:06:46] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|wDci_contig000049:2288..2434 on + strand [15:06:46] Excluding CDS which overlaps existing RNA (rRNA) at gnl|Prokka|wDci_contig000060:2645..2896 on + strand [15:06:46] Excluding CDS which overlaps existing RNA (rRNA) at gnl|Prokka|wDci_contig000062:56379..56651 on + strand [15:06:46] Found 1133 CDS [15:06:46] Connecting features back to sequences [15:06:46] Option --gram not specified, will NOT check for signal peptides. [15:06:46] Not using genus-specific database. Try --usegenus to enable it. [15:06:46] Annotating CDS, please be patient. [15:06:46] Will use 8 CPUs for similarity searching. [15:06:47] There are still 1133 unannotated CDS left (started with 1133) [15:06:47] Will use blast to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/kingdom/Bacteria/sprot with 8 CPUs [15:06:47] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wDci\/proteins\.faa | parallel --gnu -j 8 --block 21501 --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\/wDci\/proteins\.bls 2> /dev/null [15:07:04] Modify product: Probable GTP-binding protein EngB => putative GTP-binding protein EngB [15:07:04] Modify product: Uncharacterized zinc protease Rv2782c => putative zinc protease [15:07:04] Modify product: Uncharacterized zinc protease Rv2782c => putative zinc protease [15:07:04] Modify product: Actin-binding protein Smlt3054 => Actin-binding protein [15:07:04] Modify product: Transcription termination factor Rho => hypothetical protein [15:07:04] Modify product: Probable tRNA-dihydrouridine synthase => putative tRNA-dihydrouridine synthase [15:07:05] Modify product: Bifunctional protein FolD => Bifunctional protein FolD protein [15:07:05] Modify product: N utilization substance protein B homolog => hypothetical protein [15:07:05] Modify product: Putative O-methyltransferase MSMEG_5073/MSMEI_4947 => Putative O-methyltransferase/MSMEI_4947 [15:07:05] Modify product: Stringent starvation protein A homolog => hypothetical protein [15:07:05] Modify product: Uncharacterized oxidoreductase YciK => putative oxidoreductase YciK [15:07:05] Modify product: Transcription termination/antitermination protein NusG => hypothetical protein [15:07:05] Modify product: Cytochrome b561 homolog 2 => hypothetical protein [15:07:05] Modify product: Probable murein peptide carboxypeptidase => putative murein peptide carboxypeptidase [15:07:05] Modify product: Uncharacterized zinc protease Rv2782c => putative zinc protease [15:07:05] Modify product: Probable phospholipid ABC transporter-binding protein MlaD => putative phospholipid ABC transporter-binding protein MlaD [15:07:05] Modify product: Uncharacterized metal-dependent hydrolase YcfH => putative metal-dependent hydrolase YcfH [15:07:05] Modify product: Uncharacterized HIT-like protein HP_0404 => putative HIT-like protein [15:07:05] Modify product: Probable phospholipid ABC transporter permease protein MlaE => putative phospholipid ABC transporter permease protein MlaE [15:07:05] Modify product: Riboflavin biosynthesis protein VVA0006 => Riboflavin biosynthesis protein [15:07:05] Modify product: Co-chaperone protein HscB homolog => hypothetical protein [15:07:05] Modify product: Actin-binding protein Smlt3054 => Actin-binding protein [15:07:06] Modify product: Probable transcriptional regulatory protein aq_1575 => putative transcriptional regulatory protein [15:07:06] Modify product: Putative multidrug export ATP-binding/permease protein SAV1866 => Putative multidrug export ATP-binding/permease protein [15:07:06] Modify product: Probable deoxyribonuclease RhsA => putative deoxyribonuclease RhsA [15:07:06] Modify product: SCO1 protein homolog => hypothetical protein [15:07:06] Modify product: Probable acyltransferase YihG => putative acyltransferase YihG [15:07:06] Modify product: Probable chromosome-partitioning protein ParB => putative chromosome-partitioning protein ParB [15:07:06] Modify product: Transcription termination/antitermination protein NusA => hypothetical protein [15:07:06] Modify product: Probable protein kinase UbiB => putative protein kinase UbiB [15:07:06] Modify product: Probable peroxiredoxin => putative peroxiredoxin [15:07:07] Modify product: Probable deoxyribonuclease RhsA => putative deoxyribonuclease RhsA [15:07:07] Cleaned 32 /product names [15:07:07] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/proteins.faa [15:07:07] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/proteins.bls [15:07:07] There are still 535 unannotated CDS left (started with 1133) [15:07:07] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/HAMAP.hmm with 8 CPUs [15:07:07] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wDci\/proteins\.faa | parallel --gnu -j 8 --block 8090 --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\/wDci\/proteins\.bls 2> /dev/null [15:07:17] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/proteins.faa [15:07:17] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/proteins.bls [15:07:18] There are still 499 unannotated CDS left (started with 1133) [15:07:18] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/CLUSTERS.hmm with 8 CPUs [15:07:18] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wDci\/proteins\.faa | parallel --gnu -j 8 --block 7242 --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\/wDci\/proteins\.bls 2> /dev/null [15:07:43] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/proteins.faa [15:07:43] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/proteins.bls [15:07:43] There are still 464 unannotated CDS left (started with 1133) [15:07:43] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/Pfam.hmm with 8 CPUs [15:07:43] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/wDci\/proteins\.faa | parallel --gnu -j 8 --block 6751 --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\/wDci\/proteins\.bls 2> /dev/null [15:08:07] Modify product: COQ9 => hypothetical protein [15:08:07] Modify product: RDD family => RDD family protein [15:08:07] Modify product: RDD family => RDD family protein [15:08:07] Modify product: Ulp1 protease family, C-terminal catalytic domain => hypothetical protein [15:08:07] Modify product: EVE domain => EVE domain protein [15:08:07] Modify product: RmlD substrate binding domain => RmlD substrate binding domain protein [15:08:07] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:07] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:07] Modify product: KAP family P-loop domain => KAP family P-loop domain protein [15:08:07] Modify product: AhpC/TSA family => AhpC/TSA family protein [15:08:07] Modify product: Na+/H+ antiporter family => Na+/H+ antiporter family protein [15:08:07] Modify product: Transposase IS116/IS110/IS902 family => Transposase IS116/IS110/IS902 family protein [15:08:07] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:07] Modify product: Transposase DDE domain => Transposase DDE domain protein [15:08:07] Modify product: Major Facilitator Superfamily => Major Facilitator Superfamily protein [15:08:07] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:07] Modify product: Transposase IS200 like => Transposase IS200 like protein [15:08:07] Modify product: Integral membrane protein TerC family => Integral membrane protein TerC family protein [15:08:07] Modify product: DSBA-like thioredoxin domain => DSBA-like thioredoxin domain protein [15:08:07] Modify product: PAS domain => PAS domain protein [15:08:07] Modify product: PAAR motif => PAAR motif protein [15:08:07] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:07] Modify product: RDD family => RDD family protein [15:08:07] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein [15:08:07] Modify product: SPFH domain / Band 7 family => SPFH domain / Band 7 family protein [15:08:07] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein [15:08:07] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:07] Modify product: Transposase IS116/IS110/IS902 family => Transposase IS116/IS110/IS902 family protein [15:08:07] Modify product: Phage capsid family => Phage capsid family protein [15:08:07] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein [15:08:07] Modify product: Uracil DNA glycosylase superfamily => Uracil DNA glycosylase superfamily protein [15:08:07] Modify product: ETC complex I subunit conserved region => hypothetical protein [15:08:08] Modify product: Transposase IS116/IS110/IS902 family => Transposase IS116/IS110/IS902 family protein [15:08:08] Modify product: Transposase DDE domain => Transposase DDE domain protein [15:08:08] Modify product: Piwi/Argonaute/Zwille siRNA-binding domain => Piwi/Argonaute/Zwille siRNA-binding domain protein [15:08:08] Modify product: Terminase-like family => Terminase-like family protein [15:08:08] Modify product: Tim44-like domain => Tim44-like domain protein [15:08:08] Modify product: Alpha/beta hydrolase family => Alpha/beta hydrolase family protein [15:08:08] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:08] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:08] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:08] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:08] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:08] Modify product: Transposase IS116/IS110/IS902 family => Transposase IS116/IS110/IS902 family protein [15:08:08] Modify product: Terminase-like family => Terminase-like family protein [15:08:08] Modify product: MgtE intracellular N domain => MgtE intracellular N domain protein [15:08:08] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:08] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:08] Modify product: Na+/H+ antiporter family => Na+/H+ antiporter family protein [15:08:08] Modify product: YGGT family => YGGT family protein [15:08:08] Modify product: SURF1 family => SURF1 family protein [15:08:08] Modify product: RDD family => RDD family protein [15:08:08] Modify product: Transposase IS116/IS110/IS902 family => Transposase IS116/IS110/IS902 family protein [15:08:08] Modify product: Ankyrin repeat => Ankyrin repeat protein [15:08:08] Cleaned 54 /product names [15:08:08] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/proteins.faa [15:08:08] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/proteins.bls [15:08:08] Labelling remaining 337 proteins as 'hypothetical protein' [15:08:08] Possible /pseudo 'Phage terminase large subunit (GpA)' at gnl|Prokka|wDci_contig000025 position 2951 [15:08:08] Possible /pseudo 'putative zinc protease' at gnl|Prokka|wDci_contig000034 position 4360 [15:08:08] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|wDci_contig000037 position 7196 [15:08:08] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|wDci_contig000037 position 9565 [15:08:08] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|wDci_contig000037 position 12538 [15:08:08] Possible /pseudo 'RDD family protein' at gnl|Prokka|wDci_contig000037 position 25592 [15:08:08] Possible /pseudo 'Ankyrin repeat protein' at gnl|Prokka|wDci_contig000041 position 24689 [15:08:08] Possible /pseudo 'Arginine/agmatine antiporter' at gnl|Prokka|wDci_contig000047 position 8827 [15:08:08] Possible /pseudo 'Ankyrin repeats (3 copies)' at gnl|Prokka|wDci_contig000047 position 19431 [15:08:08] Possible /pseudo 'Chromosome partition protein Smc' at gnl|Prokka|wDci_contig000050 position 26847 [15:08:08] Possible /pseudo 'Ankyrin repeat protein' at gnl|Prokka|wDci_contig000055 position 69279 [15:08:08] Possible /pseudo 'Ankyrin repeat protein' at gnl|Prokka|wDci_contig000056 position 577 [15:08:08] Possible /pseudo 'Amino-acid carrier protein AlsT' at gnl|Prokka|wDci_contig000062 position 18313 [15:08:08] Possible /pseudo 'Aminopyrimidine aminohydrolase' at gnl|Prokka|wDci_contig000062 position 49740 [15:08:08] Possible /pseudo 'Phage tail tube protein FII' at gnl|Prokka|wDci_contig000062 position 185494 [15:08:08] Found 550 unique /gene codes. [15:08:08] Fixed 25 colliding /gene names. [15:08:08] Adding /locus_tag identifiers [15:08:08] Assigned 1167 locus_tags to CDS and RNA features. [15:08:08] Writing outputs to /home/yquentin/work/wolbachia/prokka/wDci/ [15:08:09] Generating annotation statistics file [15:08:09] Generating Genbank and Sequin files [15:08:09] 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\/wDci\/wDci\.err -i \/home\/yquentin\/work\/wolbachia\/prokka\/wDci\/wDci\.fsa 2> /dev/null [15:08:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/errorsummary.val [15:08:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/wDci.dr [15:08:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/wDci.fixedproducts [15:08:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/wDci.ecn [15:08:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/wDci/wDci.val [15:08:13] Output files: [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.fna [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.txt [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.gff [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.fsa [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.gbk [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.tbl [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.sqn [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.faa [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.ffn [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.err [15:08:13] /home/yquentin/work/wolbachia/prokka/wDci/wDci.log [15:08:13] Walltime used: 1.57 minutes [15:08:13] If you use this result please cite the Prokka paper: [15:08:13] Seemann T (2014) Prokka: rapid prokaryotic genome annotation. Bioinformatics. 30(14):2068-9. [15:08:13] Type 'prokka --citation' for more details. [15:08:13] Thank you, come again.