breseq  version 0.39.0  
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 240 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,208,1320GA79.4% 217.9 / 67.7 128intergenic (‑600/‑15)purH/rrsEbifunctional AICAR transformylase/IMP cyclohydrolase/16S ribosomal RNA
*NC_0009131,634,4040AG60.0% ‑3.2 / 47.6 45I161T (ATC→ACC) tfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,3880TG57.5% ‑3.7 / 48.4 40A166A (GCA→GCCtfaQputative tail fiber assembly protein TfaQ
*NC_000913580,1150CT56.5% 137.6 / inf 246intergenic (+29/+136)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_000913580,0200GA54.9% 96.7 / inf 276R47R (AGG→AGAybcWuncharacterized protein YbcW
*NC_000913580,1190AG54.8% 86.2 / inf 248intergenic (+33/+132)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_000913580,1200CT54.6% 91.2 / inf 249intergenic (+34/+131)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_000913579,8110CT54.0% 72.7 / inf 311intergenic (‑217/‑69)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
*NC_0009131,634,4820TC54.0% 23.3 / 173.7 126N135S (AAC→AGC) tfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,4840GT54.0% 24.0 / 189.1 126N134K (AAC→AAAtfaQputative tail fiber assembly protein TfaQ
*NC_000913581,1850CT52.8% 10.0 / 213.6 144L118L (CTG→TTG) nohDputative DNA‑packaging protein NohD
*NC_000913579,4570AG51.9% 52.3 / inf 342I46I (ATT→ATCybcVDUF1398 domain‑containing protein YbcV
*NC_000913578,8110TG51.7% 53.7 / inf 230Q28P (CAG→CCG) borDprophage lipoprotein BorD
*NC_0009131,634,3820TC51.4% 5.7 / 50.0 36K168K (AAA→AAGtfaQputative tail fiber assembly protein TfaQ
*NC_000913581,1330AG51.1% 61.1 / 218.9 140K100K (AAA→AAGnohDputative DNA‑packaging protein NohD
*NC_000913581,2270AC51.1% ‑5.9 / 203.9 139I132L (ATC→CTC) nohDputative DNA‑packaging protein NohD
*NC_0009131,634,5690AT51.1% 69.7 / inf 187F106Y (TTT→TAT) tfaQputative tail fiber assembly protein TfaQ
*NC_000913578,8090TC50.4% 0.2 / inf 232T29A (ACA→GCA) borDprophage lipoprotein BorD
*NC_000913581,0420TG50.3% 100.9 / 237.2 155L70R (CTG→CGG) nohDputative DNA‑packaging protein NohD
*NC_000913579,6510CT50.3% ‑2.2 / inf 352intergenic (‑57/‑229)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW

Marginal new junction evidence (lowest skew 10 of 52 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 1399593109 (1.070)31 (0.300) 25/596 3.3 22.1% intergenic (‑67/+128) fnr/ogt DNA‑binding transcriptional dual regulator FNR/methylated‑DNA‑‑[protein]‑cysteine S‑methyltransferase
?NC_000913 = 2290113 NA (NA)noncoding (1/1195 nt) IS5 repeat region
* ? NC_000913 = 103413589 (0.870)26 (0.250) 24/598 3.4 22.6% coding (882/1794 nt) hyaB hydrogenase 1 large subunit
?NC_000913 = 1470516 NA (NA)noncoding (1221/1221 nt) IS30 repeat region
* ? NC_000913 = 566774110 (1.080)58 (0.570) 19/594 4.4 34.0% pseudogene (91/93 nt) renD protein RenD
?NC_000913 568036 = 116 (1.140)pseudogene (2/213 nt) renD protein RenD
* ? NC_000913 = 460616689 (0.870)12 (0.170) 11/420 4.6 16.6% intergenic (‑1/+34) leuV/leuP tRNA‑Leu/tRNA‑Leu
?NC_000913 4606403 = 58 (0.810)intergenic (‑2/+266) leuQ/rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase
* ? NC_000913 2068941 =NA (NA)22 (0.220) 18/598 4.6 16.9% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 = 2319850 108 (1.060)coding (27/2850 nt) rcsC histidine kinase RcsC
* ? NC_000913 1063705 =96 (0.940)25 (0.240) 15/598 5.4 20.7% coding (71/921 nt) cbpA curved DNA‑binding protein
?NC_000913 = 1469240 NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 280111NA (NA)15 (0.150) 15/598 5.4 14.4% noncoding (181/181 nt) IS30 repeat region
?NC_000913 = 1034135 89 (0.870)coding (882/1794 nt) hyaB hydrogenase 1 large subunit
* ? NC_000913 = 2290113NA (NA)15 (0.150) 14/598 5.6 16.5% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 4603969 = 76 (0.740)coding (493/726 nt) yjjQ DNA‑binding transcriptional repressor YjjQ
* ? NC_000913 = 271761NA (NA)14 (0.140) 14/596 5.6 13.6% noncoding (1221/1221 nt) IS30 repeat region
?NC_000913 = 1034135 89 (0.870)coding (882/1794 nt) hyaB hydrogenase 1 large subunit
* ? NC_000913 = 265341987 (0.850)12 (0.130) 12/548 5.7 12.1% noncoding (37/100 nt) timR putative small regulatory RNA TimR
?NC_000913 2653746 = 95 (1.010)intergenic (‑291/‑109) timR/ryfA putative small regulatory RNA TimR/small regulatory RNA RyfA