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

Marginal read alignment evidence (highest frequency 20 of 156 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,1150CT65.9% 398.1 / 305.5 308intergenic (+29/+136)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_000913579,4570AG65.2% 442.7 / inf 426I46I (ATT→ATCybcVDUF1398 domain‑containing protein YbcV
*NC_000913580,1190AG64.3% 310.9 / inf 306intergenic (+33/+132)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_000913579,6710TC63.9% 402.2 / inf 427intergenic (‑77/‑209)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
*NC_000913580,0200GA63.8% 390.2 / inf 339R47R (AGG→AGAybcWuncharacterized protein YbcW
*NC_000913580,1200CT63.3% 356.4 / inf 311intergenic (+34/+131)ybcW/ylcIuncharacterized protein YbcW/DUF3950 domain‑containing protein YlcI
*NC_000913579,6510CT63.3% 496.3 / inf 454intergenic (‑57/‑229)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
*NC_000913579,7170AC62.8% 469.5 / inf 465intergenic (‑123/‑163)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
*NC_000913566,6510TC62.7% 242.8 / 300.1 256intergenic (‑673/‑33)intD/renDputative integrase/protein RenD
*NC_000913579,5020AG62.7% 416.3 / inf 489N31N (AAT→AACybcVDUF1398 domain‑containing protein YbcV
*NC_000913578,5710CT62.5% 338.6 / inf 332intergenic (+3/+29)rzpD/borDputative prophage endopeptidase RzpD/prophage lipoprotein BorD
*NC_000913579,8110CT62.4% 358.5 / inf 363intergenic (‑217/‑69)ybcV/ybcWDUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
*NC_000913579,5400TC62.3% 444.7 / inf 471N19D (AAC→GAC) ybcVDUF1398 domain‑containing protein YbcV
*NC_000913579,0990A.62.2% 332.7 / inf 354intergenic (‑206/+85)borD/ybcVprophage lipoprotein BorD/DUF1398 domain‑containing protein YbcV
*NC_000913578,7290AG61.9% 295.8 / inf 333A55A (GCT→GCCborDprophage lipoprotein BorD
*NC_000913578,5740AC61.9% 344.2 / inf 336intergenic (+6/+26)rzpD/borDputative prophage endopeptidase RzpD/prophage lipoprotein BorD
*NC_000913579,1270TC61.6% 274.9 / inf 336intergenic (‑234/+57)borD/ybcVprophage lipoprotein BorD/DUF1398 domain‑containing protein YbcV
*NC_000913578,6050TG61.3% 354.6 / inf 349K97Q (AAA→CAA) borDprophage lipoprotein BorD
*NC_000913579,1340TC61.3% 260.6 / inf 333intergenic (‑241/+50)borD/ybcVprophage lipoprotein BorD/DUF1398 domain‑containing protein YbcV
*NC_000913579,1460AG61.1% 281.3 / inf 353intergenic (‑253/+38)borD/ybcVprophage lipoprotein BorD/DUF1398 domain‑containing protein YbcV

Marginal new junction evidence (lowest skew 10 of 22 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3437349 =96 (0.910)28 (0.260) 22/598 3.6 22.6% coding (832/1377 nt) trkA NAD‑binding component of Trk potassium transporters
?NC_000913 4507459 = NA (NA)noncoding (1221/1221 nt) IS30 repeat region
* ? NC_000913 = 367067075 (0.710)28 (0.260) 21/598 3.8 27.2% coding (1079/1650 nt) treF cytoplasmic trehalase
?NC_000913 4507459 = NA (NA)noncoding (1221/1221 nt) IS30 repeat region
* ? NC_000913 = 343735091 (0.860)18 (0.170) 16/598 4.8 16.5% coding (833/1377 nt) trkA NAD‑binding component of Trk potassium transporters
?NC_000913 = 4508679 NA (NA)noncoding (1/1221 nt) IS30 repeat region
* ? NC_000913 = 294746793 (0.880)9 (0.120) 9/432 5.2 11.5% intergenic (+4/‑30) metZ/metW tRNA‑Met/tRNA‑Met
?NC_000913 2947688 = 71 (0.930)intergenic (+5/+69) metV/amiC tRNA‑Met/N‑acetylmuramoyl‑L‑alanine amidase C
* ? NC_000913 3670669 =79 (0.750)15 (0.140) 14/598 5.3 16.0% coding (1078/1650 nt) treF cytoplasmic trehalase
?NC_000913 = 4508679 NA (NA)noncoding (1/1221 nt) IS30 repeat region
* ? NC_000913 = 69674172 (0.680)8 (0.100) 7/442 6.1 12.1% intergenic (‑1/+15) metU/glnW tRNA‑Met/tRNA‑Gln
?NC_000913 697135 = 63 (0.810)intergenic (‑2/+378) metT/asnB tRNA‑Met/asparagine synthetase B
* ? NC_000913 = 4606166105 (0.990)5 (0.070) 5/420 6.7 6.6% intergenic (‑1/+34) leuV/leuP tRNA‑Leu/tRNA‑Leu
?NC_000913 4606403 = 68 (0.910)intergenic (‑2/+266) leuQ/rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase
* ? NC_000913 = 265354996 (0.910)8 (0.080) 8/566 7.0 8.6% intergenic (‑94/‑306) timR/ryfA putative small regulatory RNA TimR/small regulatory RNA RyfA
?NC_000913 2653872 = 80 (0.800)noncoding (18/304 nt) ryfA small regulatory RNA RyfA
* ? NC_000913 = 128690652 (0.490)8 (0.080) 7/558 7.3 14.6% coding (198/276 nt) ychS putative uncharacterized protein YchS
?NC_000913 1287557 = 45 (0.460)noncoding (66/85 nt) tyrT tRNA‑Tyr
* ? NC_000913 = 4542690173 (1.630)7 (0.070) 7/580 7.5 4.2% intergenic (+57/‑425) fimE/fimA regulator for fimA/type 1 fimbriae major subunit
?NC_000913 = 4542986 148 (1.440)intergenic (+353/‑129) fimE/fimA regulator for fimA/type 1 fimbriae major subunit