breseq  version 0.33.1  revision 8505477f25b3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 209 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913417,5330TG54.0% 8.0 / 30.7 38V131G (GTG→GGG) phoBresponse regulator in two‑component regulatory system with PhoR
*NC_0009132,386,8620GA38.6% 71.4 / 52.6 70intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
*NC_000913811,7210TG37.2% 227.2 / 26.8 122V67G (GTG→GGG) bioCmalonyl‑ACP O‑methyltransferase, SAM‑dependent
*NC_0009133,862,6640CG37.0% 23.0 / 13.2 50pseudogene (940/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
*NC_0009132,684,7790AC34.7% 65.4 / 11.1 50G243G (GGT→GGGglyAserine hydroxymethyltransferase
*NC_00091375,0230AC34.2% 127.6 / 17.2 73V153G (GTG→GGG) thiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
*NC_0009133,745,0190CG33.3% 38.0 / 11.5 42G72G (GGC→GGGyiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
*NC_000913540,5580TG33.0% 192.0 / 37.7 100intergenic (+50/‑7)allB/ybbYallantoinase/putative uracil/xanthine transporter
*NC_0009131,265,4730TG33.0% 209.3 / 20.6 94G154G (GGT→GGGprfApeptide chain release factor RF‑1
*NC_000913433,6490TG32.5% 72.1 / 19.3 83G65G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
*NC_000913354,2880TG31.6% 132.3 / 16.6 57V528G (GTG→GGG) prpEpropionate‑‑CoA ligase
*NC_0009133,782,9080TG31.2% 80.5 / 11.9 74T252P (ACC→CCC) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
*NC_0009132,509,8800AC31.2% 302.5 / 29.0 128D84A (GAC→GCC) yfeOputative ion channel protein
*NC_0009132,483,2610CA30.7% 169.4 / 16.2 76A27S (GCG→TCG) emrKmultidrug resistance efflux pump membrane fusion protein
*NC_0009131,383,7220GA30.6% 124.0 / 19.2 63F80F (TTC→TTTycjWLacI family putative transcriptional repressor
*NC_000913753,1250TA30.6% 213.2 / 46.2 98intergenic (‑330/+60)ybgD/gltAputative fimbrial‑like adhesin protein/citrate synthase
*NC_0009133,760,3840AC30.3% 120.3 / 19.8 76V289G (GTG→GGG) selAselenocysteine synthase
*NC_0009132,878,9470CT30.2% 212.3 / 25.6 106E276K (GAA→AAA) ygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
*NC_0009131,346,8530GT30.1% 192.4 / 25.5 93V20L (GTA→TTA) ymiBuncharacterized protein
*NC_00091329,4960TA30.0% 107.3 / 20.4 50intergenic (+301/‑155)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase

Marginal new junction evidence (lowest skew 10 of 30 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2713829120 (0.890)10 (0.080) 8/326 8.1 8.2% coding (934/1335 nt) srmB ATP‑dependent RNA helicase
?NC_000913 = 2713908 110 (0.870)coding (1013/1335 nt) srmB ATP‑dependent RNA helicase
* ? NC_000913 = 3698975131 (0.970)8 (0.060) 6/330 9.0 6.2% coding (762/1680 nt) bcsG DUF3260 family cellulose production inner membrane protein
?NC_000913 = 3699049 117 (0.910)coding (836/1680 nt) bcsG DUF3260 family cellulose production inner membrane protein
* ? NC_000913 4542682 =72 (0.530)5 (0.040) 5/330 9.6 4.9% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 4542996 = 124 (0.960)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
* ? NC_000913 = 449838110 (0.070)16 (0.120) 5/342 9.8 62.0% noncoding (201/1331 nt) IS2 repeat region
?NC_000913 = 4498773 NA (NA)noncoding (593/1331 nt) IS2 repeat region
* ? NC_000913 = 4068629183 (1.350)4 (0.030) 4/324 10.0 2.6% coding (648/2037 nt) yihQ putative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase
?NC_000913 = 4068723 132 (1.050)coding (554/2037 nt) yihQ putative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase
* ? NC_000913 = 4095482187 (1.380)4 (0.030) 4/328 10.1 2.5% coding (501/1260 nt) rhaA L‑rhamnose isomerase
?NC_000913 = 4095710 132 (1.030)coding (273/1260 nt) rhaA L‑rhamnose isomerase
* ? NC_000913 = 323503095 (0.700)5 (0.040) 4/338 10.2 4.6% coding (846/1137 nt) rlmG 23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
?NC_000913 = 3235119 117 (0.890)coding (757/1137 nt) rlmG 23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
* ? NC_000913 = 2283632112 (0.830)4 (0.030) 4/336 10.2 3.5% coding (316/1008 nt) yejK spermidine nucleoid‑associated protein
?NC_000913 = 2283824 110 (0.840)coding (124/1008 nt) yejK spermidine nucleoid‑associated protein
* ? NC_000913 1060875 =130 (0.960)4 (0.030) 4/336 10.2 2.8% coding (1620/2547 nt) torA trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
?NC_000913 1061210 = 149 (1.140)coding (1955/2547 nt) torA trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
* ? NC_000913 992281 =129 (0.950)4 (0.030) 4/336 10.2 3.3% coding (1660/2613 nt) pepN aminopeptidase N
?NC_000913 992374 = 109 (0.830)coding (1753/2613 nt) pepN aminopeptidase N