[10:41:12] Command: --force --outdir /home/yquentin/work/wolbachia/prokka/Ob_Wba --addgenes --prefix Ob_Wba --locustag Ob_Wba --species Wolbachia endosymbiont --strain Ob_Wba --compliant --quiet /home/yquentin/work/wolbachia/prokka/Ob_Wba.fna
[10:41:12] Looking for 'aragorn' - found /usr/local/bioinfo/bin/aragorn
[10:41:12] Determined aragorn version is 1.2
[10:41:12] Looking for 'barrnap' - found /usr/local/bioinfo/bin/barrnap
[10:41:12] Determined barrnap version is 0.5
[10:41:12] Looking for 'blastp' - found /usr/local/bioinfo/src/NCBI_Blast+/ncbi-blast-2.6.0+/bin/blastp
[10:41:12] Determined blastp version is 2.6
[10:41:12] Looking for 'cmpress' - found /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../binaries/linux/cmpress
[10:41:12] Determined cmpress version is 1.1
[10:41:12] Looking for 'cmscan' - found /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../binaries/linux/cmscan
[10:41:12] Determined cmscan version is 1.1
[10:41:12] Looking for 'egrep' - found /bin/egrep
[10:41:12] Looking for 'find' - found /bin/find
[10:41:12] Looking for 'grep' - found /bin/grep
[10:41:13] Looking for 'hmmpress' - found /usr/local/bioinfo/bin/hmmpress
[10:41:13] Determined hmmpress version is 3.1
[10:41:13] Looking for 'hmmscan' - found /usr/local/bioinfo/bin/hmmscan
[10:41:13] Determined hmmscan version is 3.1
[10:41:13] Looking for 'less' - found /usr/bin/less
[10:41:13] Looking for 'makeblastdb' - found /usr/local/bioinfo/src/NCBI_Blast+/ncbi-blast-2.6.0+/bin/makeblastdb
[10:41:13] Determined makeblastdb version is 2.6
[10:41:13] Looking for 'minced' - found /usr/local/bioinfo/bin/minced
[10:41:13] Determined minced version is 1.6
[10:41:13] Looking for 'parallel' - found /usr/local/bioinfo/bin/parallel
[10:41:13] Determined parallel version is 20141022
[10:41:13] Looking for 'prodigal' - found /usr/local/bioinfo/bin/prodigal
[10:41:13] Determined prodigal version is 2.6
[10:41:13] Looking for 'rnammer' - found /usr/local/bioinfo/bin/rnammer
[10:41:13] Determined rnammer version is 1.2
[10:41:13] Looking for 'sed' - found /bin/sed
[10:41:13] Looking for 'signalp' - found /usr/local/bioinfo/bin/signalp
[10:41:13] Determined signalp version is 4.1
[10:41:13] Looking for 'tbl2asn' - found /home/yquentin/save/scripts/tbl2asn
[10:41:13] Determined tbl2asn version is 25.3
[10:41:13] Using genetic code table 11.
[10:41:13] Loading and checking input file: /home/yquentin/work/wolbachia/prokka/Ob_Wba.fna
[10:41:14] Wrote 120 contigs
[10:41:14] Predicting tRNAs and tmRNAs
[10:41:14] Running: aragorn -gc11  -w \/home\/yquentin\/work\/wolbachia\/prokka\/Ob_Wba\/Ob_Wba\.fna
[10:41:14] 1 tRNA-Thr c[3387,3460] 34 (tgt)
[10:41:14] 1 tRNA-Arg [1411,1486] 36 (tct)
[10:41:14] 1 tRNA-Trp [1198,1271] 34 (cca)
[10:41:14] 1 tRNA-Val [1705,1778] 34 (gac)
[10:41:14] 1 tRNA-Arg [2641,2718] 36 (acg)
[10:41:14] 2 tRNA-Leu [14144,14230] 35 (taa)
[10:41:14] 1 tRNA-Cys [1942,2014] 34 (gca)
[10:41:14] 1 tRNA-Leu [3842,3925] 34 (caa)
[10:41:14] 1 tRNA-Ile [572,648] 35 (gat)
[10:41:14] 1 tRNA-Gly c[64,137] 34 (gcc)
[10:41:14] 2 tRNA-Met [1354,1430] 35 (cat)
[10:41:14] 1 tRNA-Asn [1232,1304] 33 (gtt)
[10:41:14] 1 tRNA-His c[1025,1100] 34 (gtg)
[10:41:14] 1 tRNA-Gln c[989,1061] 33 (ttg)
[10:41:14] 2 tRNA-Met c[2021,2093] 33 (cat)
[10:41:14] 1 tRNA-Ser c[18511,18598] 35 (gga)
[10:41:14] 1 tRNA-Met c[24168,24244] 35 (cat)
[10:41:14] 1 tRNA-Ser c[7840,6] 35 (gct)
[10:41:14] tRNA c[7840,6] has start(7840) > end (6) - skipping.
[10:41:14] 1 tRNA-Ser c[4010,4098] 35 (tga)
[10:41:14] 2 tRNA-Lys c[10065,10138] 34 (ttt)
[10:41:14] 3 tRNA-Ala [10439,10514] 34 (tgc)
[10:41:14] 1 tRNA-Thr [153,227] 35 (ggt)
[10:41:14] 1 tRNA-Leu [6991,7073] 35 (tag)
[10:41:14] 1 tRNA-Val c[5980,6052] 34 (tac)
[10:41:14] 2 tRNA-Pro c[13112,13187] 36 (tgg)
[10:41:14] 3 tRNA-Glu c[14036,14108] 34 (ttc)
[10:41:14] 1 tRNA-Phe c[2159,2232] 34 (gaa)
[10:41:14] 1 tRNA-Arg c[3903,3977] 35 (cct)
[10:41:14] 1 tRNA-Leu c[7569,7656] 35 (gag)
[10:41:14] 2 tRNA-Asp [17858,17934] 35 (gtc)
[10:41:14] 3 tRNA-Gly c[77867,77938] 33 (tcc)
[10:41:14] 4 tRNA-Tyr c[77942,78025] 36 (gta)
[10:41:14] Found 31 tRNAs
[10:41:14] Predicting Ribosomal RNAs
[10:41:14] Running Barrnap with 8 threads
[10:41:14] 1 gnl|Prokka|Ob_Wba_contig000003 23546 16S ribosomal RNA
[10:41:14] 2 gnl|Prokka|Ob_Wba_contig000021 8483 23S ribosomal RNA (partial)
[10:41:14] 3 gnl|Prokka|Ob_Wba_contig000052 4962 5S ribosomal RNA
[10:41:14] Found 3 rRNAs
[10:41:14] Skipping ncRNA search, enable with --rfam if desired.
[10:41:14] Total of 34 tRNA + rRNA features
[10:41:14] Searching for CRISPR repeats
[10:41:15] Found 0 CRISPRs
[10:41:15] Predicting coding sequences
[10:41:15] Contigs total 1121046 bp, so using single mode
[10:41:15] Running: prodigal -i \/home\/yquentin\/work\/wolbachia\/prokka\/Ob_Wba\/Ob_Wba\.fna -c -m -g 11 -p single -f sco -q
[10:41:16] Excluding CDS which overlaps existing RNA (rRNA) at gnl|Prokka|Ob_Wba_contig000003:24692..24943 on + strand
[10:41:16] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|Ob_Wba_contig000032:1052..1387 on - strand
[10:41:16] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|Ob_Wba_contig000032:1969..2097 on - strand
[10:41:16] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|Ob_Wba_contig000032:2084..3385 on - strand
[10:41:16] Excluding CDS which overlaps existing RNA (tRNA) at gnl|Prokka|Ob_Wba_contig000100:124..273 on + strand
[10:41:16] Found 1027 CDS
[10:41:16] Connecting features back to sequences
[10:41:16] Option --gram not specified, will NOT check for signal peptides.
[10:41:16] Not using genus-specific database. Try --usegenus to enable it.
[10:41:16] Annotating CDS, please be patient.
[10:41:16] Will use 8 CPUs for similarity searching.
[10:41:17] There are still 1027 unannotated CDS left (started with 1027)
[10:41:17] Will use blast to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/kingdom/Bacteria/sprot with 8 CPUs
[10:41:17] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/Ob_Wba\/proteins\.faa | parallel --gnu -j 8 --block 19139 --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\/Ob_Wba\/proteins\.bls 2> /dev/null
[10:41:27] Modify product: Transcription termination factor Rho => hypothetical protein
[10:41:27] Modify product: Actin-binding protein Smlt3054 => Actin-binding protein
[10:41:27] Modify product: Uncharacterized zinc protease Rv2782c => putative zinc protease
[10:41:27] Modify product: Uncharacterized zinc protease Rv2782c => putative zinc protease
[10:41:27] Modify product: Actin-binding protein Smlt3054 => Actin-binding protein
[10:41:27] Modify product: Uncharacterized metal-dependent hydrolase YcfH => putative metal-dependent hydrolase YcfH
[10:41:27] Modify product: Actin-binding protein Smlt3054 => Actin-binding protein
[10:41:27] Modify product: Transcription termination/antitermination protein NusG => hypothetical protein
[10:41:27] Modify product: Bifunctional protein FolD => Bifunctional protein FolD protein
[10:41:27] Modify product: Actin-binding protein Smlt3054 => Actin-binding protein
[10:41:27] Modify product: Transcription termination/antitermination protein NusA => hypothetical protein
[10:41:27] Modify product: SCO1 protein homolog => hypothetical protein
[10:41:28] Modify product: Uncharacterized zinc protease Rv2782c => putative zinc protease
[10:41:28] Modify product: Probable murein peptide carboxypeptidase => putative murein peptide carboxypeptidase
[10:41:28] Modify product: Cytochrome b561 homolog 2 => hypothetical protein
[10:41:28] Modify product: Probable tRNA-dihydrouridine synthase => putative tRNA-dihydrouridine synthase
[10:41:28] Modify product: Probable GTP-binding protein EngB => putative GTP-binding protein EngB
[10:41:28] Modify product: Co-chaperone protein HscB homolog => hypothetical protein
[10:41:28] Modify product: Probable phospholipid ABC transporter-binding protein MlaD => putative phospholipid ABC transporter-binding protein MlaD
[10:41:28] Modify product: N utilization substance protein B homolog => hypothetical protein
[10:41:28] Modify product: Probable protein kinase UbiB => putative protein kinase UbiB
[10:41:28] Modify product: Probable peroxiredoxin => putative peroxiredoxin
[10:41:28] Modify product: Probable ribonucleotide transport ATP-binding protein mkl => putative ribonucleotide transport ATP-binding protein mkl
[10:41:28] Modify product: Probable phospholipid ABC transporter permease protein MlaE => putative phospholipid ABC transporter permease protein MlaE
[10:41:28] Modify product: Stringent starvation protein A homolog => hypothetical protein
[10:41:28] Modify product: Probable acyltransferase YihG => putative acyltransferase YihG
[10:41:28] Modify product: Riboflavin biosynthesis protein VVA0006 => Riboflavin biosynthesis protein
[10:41:28] Modify product: Uncharacterized oxidoreductase YciK => putative oxidoreductase YciK
[10:41:28] Modify product: Probable chromosome-partitioning protein ParB => putative chromosome-partitioning protein ParB
[10:41:28] Modify product: Putative multidrug export ATP-binding/permease protein SAV1866 => Putative multidrug export ATP-binding/permease protein
[10:41:28] Modify product: Probable transcriptional regulatory protein aq_1575 => putative transcriptional regulatory protein
[10:41:28] Modify product: Uncharacterized oxidoreductase Rv0945 => putative oxidoreductase
[10:41:28] Cleaned 32 /product names
[10:41:28] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/proteins.faa
[10:41:28] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/proteins.bls
[10:41:28] There are still 450 unannotated CDS left (started with 1027)
[10:41:28] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/HAMAP.hmm with 8 CPUs
[10:41:28] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/Ob_Wba\/proteins\.faa | parallel --gnu -j 8 --block 6406 --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\/Ob_Wba\/proteins\.bls 2> /dev/null
[10:41:35] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/proteins.faa
[10:41:35] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/proteins.bls
[10:41:35] There are still 424 unannotated CDS left (started with 1027)
[10:41:35] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/CLUSTERS.hmm with 8 CPUs
[10:41:35] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/Ob_Wba\/proteins\.faa | parallel --gnu -j 8 --block 5906 --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\/Ob_Wba\/proteins\.bls 2> /dev/null
[10:41:51] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/proteins.faa
[10:41:51] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/proteins.bls
[10:41:52] There are still 406 unannotated CDS left (started with 1027)
[10:41:52] Will use hmmer3 to search against /usr/local/bioinfo/src/PROKKA/prokka-1.10/bin/../db/hmm/Pfam.hmm with 8 CPUs
[10:41:52] Running: cat \/home\/yquentin\/work\/wolbachia\/prokka\/Ob_Wba\/proteins\.faa | parallel --gnu -j 8 --block 5645 --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\/Ob_Wba\/proteins\.bls 2> /dev/null
[10:42:09] Modify product: COQ9 => hypothetical protein
[10:42:09] Modify product: Transposase DDE domain => Transposase DDE domain protein
[10:42:09] Modify product: DSBA-like thioredoxin domain => DSBA-like thioredoxin domain protein
[10:42:09] Modify product: Ankyrin repeat => Ankyrin repeat protein
[10:42:09] Modify product: Alpha/beta hydrolase family => Alpha/beta hydrolase family protein
[10:42:09] Modify product: Na+/H+ antiporter family => Na+/H+ antiporter family protein
[10:42:09] Modify product: Alpha/beta hydrolase family => Alpha/beta hydrolase family protein
[10:42:09] Modify product: Piwi/Argonaute/Zwille siRNA-binding domain => Piwi/Argonaute/Zwille siRNA-binding domain protein
[10:42:09] Modify product: Major Facilitator Superfamily => Major Facilitator Superfamily protein
[10:42:09] Modify product: Ankyrin repeat => Ankyrin repeat protein
[10:42:09] Modify product: Sodium:dicarboxylate symporter family => Sodium:dicarboxylate symporter family protein
[10:42:09] Modify product: Na+/H+ antiporter family => Na+/H+ antiporter family protein
[10:42:09] Modify product: RmlD substrate binding domain => RmlD substrate binding domain protein
[10:42:09] Modify product: Transposase IS200 like => Transposase IS200 like protein
[10:42:09] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein
[10:42:09] Modify product: DnaJ domain => DnaJ domain protein
[10:42:09] Modify product: Transposase IS200 like => Transposase IS200 like protein
[10:42:09] Modify product: Transposase DDE domain => Transposase DDE domain protein
[10:42:09] Modify product: PD-(D/E)XK nuclease superfamily => PD-(D/E)XK nuclease superfamily protein
[10:42:09] Modify product: Tim44-like domain => Tim44-like domain protein
[10:42:09] Modify product: RDD family => RDD family protein
[10:42:09] Modify product: YGGT family => YGGT family protein
[10:42:09] Modify product: EVE domain => EVE domain protein
[10:42:09] Modify product: AhpC/TSA family => AhpC/TSA family protein
[10:42:09] Modify product: Ankyrin repeat => Ankyrin repeat protein
[10:42:09] Modify product: Transposase IS66 family => Transposase IS66 family protein
[10:42:09] Modify product: Phage capsid family => Phage capsid family protein
[10:42:09] Modify product: PAS domain => PAS domain protein
[10:42:09] Modify product: Phosphotransferase enzyme family => Phosphotransferase enzyme family protein
[10:42:09] Modify product: KAP family P-loop domain => KAP family P-loop domain protein
[10:42:09] Modify product: Transposase IS200 like => Transposase IS200 like protein
[10:42:09] Modify product: RDD family => RDD family protein
[10:42:09] Modify product: SURF1 family => SURF1 family protein
[10:42:09] Modify product: Terminase-like family => Terminase-like family protein
[10:42:09] Modify product: SPFH domain / Band 7 family => SPFH domain / Band 7 family protein
[10:42:09] Modify product: TrbC/VIRB2 family => TrbC/VIRB2 family protein
[10:42:09] Modify product: Alpha/beta hydrolase family => Alpha/beta hydrolase family protein
[10:42:09] Modify product: Uracil DNA glycosylase superfamily => Uracil DNA glycosylase superfamily protein
[10:42:09] Modify product: ETC complex I subunit conserved region => hypothetical protein
[10:42:09] Modify product: Integral membrane protein TerC family => Integral membrane protein TerC family protein
[10:42:09] Cleaned 40 /product names
[10:42:09] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/proteins.faa
[10:42:09] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/proteins.bls
[10:42:09] Labelling remaining 297 proteins as 'hypothetical protein'
[10:42:09] Possible /pseudo 'putative zinc protease' at gnl|Prokka|Ob_Wba_contig000001 position 33616
[10:42:09] Possible /pseudo 'Amino-acid carrier protein AlsT' at gnl|Prokka|Ob_Wba_contig000008 position 6624
[10:42:09] Possible /pseudo 'Aminopyrimidine aminohydrolase' at gnl|Prokka|Ob_Wba_contig000021 position 2669
[10:42:09] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|Ob_Wba_contig000030 position 4644
[10:42:09] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|Ob_Wba_contig000030 position 7541
[10:42:09] Possible /pseudo 'TrbL/VirB6 plasmid conjugal transfer protein' at gnl|Prokka|Ob_Wba_contig000030 position 9931
[10:42:09] Possible /pseudo 'Ankyrin repeats (3 copies)' at gnl|Prokka|Ob_Wba_contig000040 position 2271
[10:42:09] Possible /pseudo 'Arginine/agmatine antiporter' at gnl|Prokka|Ob_Wba_contig000040 position 9503
[10:42:09] Possible /pseudo 'Exopolysaccharide synthesis, ExoD' at gnl|Prokka|Ob_Wba_contig000041 position 14933
[10:42:09] Possible /pseudo 'Transposase, Mutator family' at gnl|Prokka|Ob_Wba_contig000118 position 3291
[10:42:09] Possible /pseudo 'Glutamine--fructose-6-phosphate aminotransferase [isomerizing]' at gnl|Prokka|Ob_Wba_contig000120 position 19520
[10:42:09] Found 525 unique /gene codes.
[10:42:09] Fixed 24 colliding /gene names.
[10:42:09] Adding /locus_tag identifiers
[10:42:09] Assigned 1061 locus_tags to CDS and RNA features.
[10:42:09] Writing outputs to /home/yquentin/work/wolbachia/prokka/Ob_Wba/
[10:42:10] Generating annotation statistics file
[10:42:10] Generating Genbank and Sequin files
[10:42:10] 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\/Ob_Wba\/Ob_Wba\.err -i \/home\/yquentin\/work\/wolbachia\/prokka\/Ob_Wba\/Ob_Wba\.fsa 2> /dev/null
[10:42:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/errorsummary.val
[10:42:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.dr
[10:42:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.fixedproducts
[10:42:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.ecn
[10:42:13] Deleting unwanted file: /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.val
[10:42:13] Output files:
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.gbk
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.tbl
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.faa
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.txt
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.fsa
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.fna
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.log
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.err
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.sqn
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.ffn
[10:42:13] /home/yquentin/work/wolbachia/prokka/Ob_Wba/Ob_Wba.gff
[10:42:13] Walltime used: 1.02 minutes
[10:42:13] If you use this result please cite the Prokka paper:
[10:42:13] Seemann T (2014) Prokka: rapid prokaryotic genome annotation. Bioinformatics. 30(14):2068-9.
[10:42:13] Type 'prokka --citation' for more details.
[10:42:13] Share and enjoy!