breseq  version 0.27.1  revision 87c22d663cc3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 1418 shown)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8620GA66.1% 92.2 / 112.4 121intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
*NC_000913357,8900CA61.2% 4.0 / 48.1 50L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
*NC_0009132,332,8230AG56.7% 54.3 / 107.6 93pseudogene (1580/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
*NC_000913354,2880TG56.2% 69.4 / 133.8 91V528G (GTG→GGG) prpEpropionate‑‑CoA ligase
*NC_0009133,133,3020TG54.4% 129.8 / 191.6 125H219P (CAC→CCC) yghSputative ATP‑binding protein
*NC_0009134,604,3440CT54.4% 23.5 / 147.6 125S62F (TCT→TTT) bglJbgl operon transcriptional activator
*NC_0009133,665,7020CA54.4% 67.8 / 201.2 149A37S (GCG→TCG) gadXacid resistance regulon transcriptional activator; autoactivator
*NC_0009133,665,6980TC52.4% 72.6 / 199.7 145D38G (GAT→GGT) gadXacid resistance regulon transcriptional activator; autoactivator
*NC_000913993,9170AC51.9% 54.3 / 81.1 52V43G (GTG→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
*NC_0009131,222,2330GA51.5% 38.5 / 91.0 68pseudogene (60/75 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
*NC_0009131,554,0000TG51.2% 81.9 / 105.5 87Y557S (TAC→TCC) maeAmalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme
*NC_0009131,356,2010TG50.7% 99.2 / 123.4 70T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
*NC_00091329,4960TA50.0% 74.2 / 72.2 86intergenic (+301/‑155)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
*NC_0009132,386,8680TA49.6% 78.3 / 86.9 120intergenic (+37/‑66)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
*NC_0009132,571,0360TG49.2% 75.0 / 22.1 67N239T (AAC→ACC) eutEaldehyde oxidoreductase, ethanolamine utilization protein
*NC_0009132,386,8630CA49.1% 105.2 / 65.9 121intergenic (+32/‑71)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
*NC_0009132,386,8830CT48.9% 74.5 / 120.0 131intergenic (+52/‑51)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
*NC_0009134,622,2500AG48.0% 105.6 / 57.1 103E161G (GAG→GGG) yjjJputative protein kinase
*NC_000913650,4580AG47.7% 114.1 / 24.0 90T11T (ACT→ACCcitEcitrate lyase, citryl‑ACP lyase (beta) subunit
*NC_0009132,903,7150GA47.4% 36.9 / 84.8 78intergenic (+341/+1032)ygcE/queEputative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis

Marginal new junction evidence (lowest skew 10 of 217 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780579 (0.700)17 (0.170) 14/266 5.0 19.8% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 67 (0.670)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
* ? NC_000913 = 381689NA (NA)18 (0.170) 15/280 5.0 NA noncoding (430/1331 nt) IS2 repeat region
?NC_000913 = 381743 NA (NA)noncoding (484/1331 nt) IS2 repeat region
* ? NC_000913 4409174 =NA (NA)10 (0.100) 10/268 6.4 6.9% noncoding (22/98 nt) RIP331 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site RIP331 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site
?NC_000913 4409235 = 135 (1.340)noncoding (83/98 nt) RIP331 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site RIP331 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site
* ? NC_000913 1207790 =73 (0.650)11 (0.110) 10/266 6.4 14.1% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 69 (0.690)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
* ? NC_000913 = 3880655124 (1.110)13 (0.120) 10/278 6.6 10.0% coding (567/1074 nt) recF gap repair protein
?NC_000913 = 3880717 118 (1.130)coding (505/1074 nt) recF gap repair protein
* ? NC_000913 4005815 =122 (1.090)9 (0.080) 9/282 7.1 7.3% intergenic (+84/‑49) pldA/recQ outer membrane phospholipase A/ATP‑dependent DNA helicase
?NC_000913 4005837 = 114 (1.070)intergenic (+106/‑27) pldA/recQ outer membrane phospholipase A/ATP‑dependent DNA helicase
* ? NC_000913 = 3529405133 (1.190)11 (0.110) 9/278 7.1 7.6% coding (58/402 nt) hslR ribosome‑associated heat shock protein Hsp15
?NC_000913 = 3529461 142 (1.360)coding (114/402 nt) hslR ribosome‑associated heat shock protein Hsp15
* ? NC_000913 = 373065120 (1.070)10 (0.100) 8/262 7.1 8.6% coding (145/951 nt) mhpF acetaldehyde‑CoA dehydrogenase II, NAD‑binding
?NC_000913 = 373047 107 (1.090)coding (127/951 nt) mhpF acetaldehyde‑CoA dehydrogenase II, NAD‑binding
* ? NC_000913 326549 =115 (1.030)10 (0.100) 9/278 7.1 8.3% coding (699/1671 nt) betA choline dehydrogenase, a flavoprotein
?NC_000913 326652 = 115 (1.100)coding (596/1671 nt) betA choline dehydrogenase, a flavoprotein
* ? NC_000913 = 4062810119 (1.060)9 (0.090) 8/280 7.5 7.4% coding (564/1266 nt) yihN MFS transporter family protein
?NC_000913 = 4062816 114 (1.080)coding (570/1266 nt) yihN MFS transporter family protein