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

Marginal read alignment evidence (highest frequency 20 of 251 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,7240GA60.3% ‑3.4 / 165.0 141F54F (TTC→TTTtfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,6990AC60.2% 3.9 / 180.0 162S63A (TCA→GCA) tfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,6610CG59.7% ‑1.5 / 252.4 176T75T (ACG→ACCtfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,6760AG59.4% ‑2.2 / 219.8 170D70D (GAT→GACtfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,5690AT58.1% 5.9 / 272.5 155F106Y (TTT→TAT) tfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,5710CT58.1% 0.3 / 206.6 155E105E (GAG→GAAtfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,5970CA58.0% 7.8 / 230.2 169V97F (GTT→TTT) tfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,5780TC57.7% 2.7 / 202.6 158E103G (GAA→GGA) tfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,6580AG56.6% 18.6 / 246.5 183V76V (GTT→GTCtfaQputative tail fiber assembly protein TfaQ
*NC_000913580,1150CT55.1% 82.8 / 255.0 185intergenic (+29/+136)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_0009131,634,4820TC54.0% 15.0 / 139.6 100N135S (AAC→AGC) tfaQputative tail fiber assembly protein TfaQ
*NC_0009131,634,4840GT54.0% 15.9 / 149.1 100N134K (AAC→AAAtfaQputative tail fiber assembly protein TfaQ
*NC_000913580,1190AG53.8% 30.6 / 276.9 184intergenic (+33/+132)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_000913579,2850TC53.7% 47.5 / inf 246I104V (ATT→GTT) ybcVDUF1398 domain‑containing protein YbcV
*NC_000913580,1200CT53.5% 59.0 / 267.7 185intergenic (+34/+131)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_000913578,5710CT53.1% 64.0 / inf 241intergenic (+3/+29)rzpD/borDputative prophage endopeptidase RzpD/prophage lipoprotein BorD
*NC_000913578,6050TG52.7% 55.0 / inf 241K97Q (AAA→CAA) borDprophage lipoprotein BorD
*NC_000913578,5740AC52.7% 90.1 / inf 243intergenic (+6/+26)rzpD/borDputative prophage endopeptidase RzpD/prophage lipoprotein BorD
*NC_000913579,7170AC52.6% 117.2 / inf 323intergenic (‑123/‑163)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
*NC_000913579,2710GA52.5% 44.3 / inf 240H108H (CAC→CATybcVDUF1398 domain‑containing protein YbcV

Marginal new junction evidence (lowest skew 10 of 18 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 382179 =NA (NA)29 (0.280) 22/594 3.6 NA noncoding (920/1331 nt) IS2 repeat region
?NC_000913 = 382590 NA (NA)noncoding (1331/1331 nt) IS2 repeat region
* ? NC_000913 = 4542690165 (1.600)25 (0.250) 21/580 3.7 15.3% intergenic (+57/‑425) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 = 4542986 117 (1.170)intergenic (+353/‑129) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
* ? NC_000913 4542682 =152 (1.470)25 (0.250) 20/580 3.9 15.9% intergenic (+49/‑433) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 4542996 = 117 (1.170)intergenic (+363/‑119) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
* ? NC_000913 = 566774102 (0.990)39 (0.380) 18/594 4.4 28.5% pseudogene (91/93 nt) renD protein RenD
?NC_000913 568036 = 94 (0.920)pseudogene (2/213 nt) renD protein RenD
* ? NC_000913 = 4606166121 (1.170)8 (0.110) 8/420 5.5 9.4% intergenic (‑1/+34) leuV/leuP tRNA‑Leu/tRNA‑Leu
?NC_000913 4606403 = 70 (0.970)intergenic (‑2/+266) leuQ/rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase
* ? NC_000913 = 3648146147 (1.420)8 (0.100) 8/482 6.2 6.8% intergenic (+2/‑148) agrA/agrB small RNA AgrA/small regulatory RNA AgrB
?NC_000913 3648378 = 100 (1.200)intergenic (+3/‑150) agrB/arsR small regulatory RNA AgrB/DNA‑binding transcriptional repressor ArsR
* ? NC_000913 = 258675NA (NA)9 (0.090) 7/598 7.8 9.0% noncoding (1/768 nt) IS1 repeat region
?NC_000913 = 543175 91 (0.880)intergenic (+142/+86) glxK/allE glycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase
* ? NC_000913 257908 =NA (NA)7 (0.070) 7/598 7.8 7.2% noncoding (768/768 nt) IS1 repeat region
?NC_000913 543168 = 90 (0.870)intergenic (+135/+93) glxK/allE glycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase
* ? NC_000913 = 196252397 (0.940)8 (0.080) 6/596 8.2 7.6% intergenic (+63/+57) yecT/flhE protein YecT/flagellar protein FlhE
?NC_000913 = 1979270 NA (NA)noncoding (1/768 nt) IS1 repeat region
* ? NC_000913 1096961 =83 (0.800)3 (0.040) 3/438 8.2 4.3% intergenic (+132/+604) ycdU/serX uncharacterized protein YcdU/tRNA‑Ser
?NC_000913 = 1097503 73 (0.970)intergenic (+674/+62) ycdU/serX uncharacterized protein YcdU/tRNA‑Ser