breseq  version 0.28.1  revision 52323b585d6f
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 37 shown)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1936,1690T.37.1% 239.9 / 150.6 106intergenic (‑13/+21)ECDH1_RS04615/ECDH1_RS04620tRNA‑Met/tRNA‑Met
*DH13,906,8320CA36.1% 108.7 / 18.3 92intergenic (+3/‑32)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
*DH1936,2930T.34.2% 362.5 / 210.9 163intergenic (‑27/+7)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
*DH12,169,5480GT30.6% 502.2 / 82.1 276P671Q (CCG→CAG) rsxCelectron transporter RsxC
*DH12,326,2870GT29.9% 972.9 / 171.0 544intergenic (+207/+200)yddM/fdnItranscriptional regulator/formate dehydrogenase nitrate‑inducible cytochrome b556(Fdn) subunit
*DH12,169,5770TG29.5% 263.7 / 65.6 253R661R (CGA→CGCrsxCelectron transporter RsxC
*DH1936,2970AG29.0% 391.8 / 181.0 226intergenic (‑31/+3)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
*DH13,448,5260AC27.8% 795.3 / 120.0 450intergenic (‑86/+70)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
*DH12,591,5130TG27.1% 269.1 / 59.1 297D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
*DH1613,5970GC26.4% 333.6 / 54.0 258intergenic (+458/+565)garK/yhaCglycerate 2‑kinase/pentapeptide repeats‑related protein
*DH12,169,6730AT25.9% 140.8 / 35.2 161R629R (CGT→CGArsxCelectron transporter RsxC
*DH1936,2781.T25.7% 249.0 / 135.3 102intergenic (‑12/+22)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
*DH12,169,5980GT25.6% 307.1 / 30.7 216A654A (GCC→GCArsxCelectron transporter RsxC
*DH1613,6130TG25.5% 241.7 / 35.4 225intergenic (+474/+549)garK/yhaCglycerate 2‑kinase/pentapeptide repeats‑related protein
*DH14,122,1240CA25.2% 426.7 / 18.5 291intergenic (‑15/+22)ECDH1_RS20430/ECDH1_RS20435tRNA‑Gly/tRNA‑Gly
*DH11,065,6050AT24.9% 342.5 / 112.6 218intergenic (+40/‑23)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
*DH13,906,8370CG24.2% 163.6 / 17.6 105intergenic (+8/‑27)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
*DH11,065,6180AG23.7% 377.7 / 132.5 240intergenic (+53/‑10)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
*DH12,169,5680CG23.2% 499.4 / 28.1 254A664A (GCG→GCCrsxCelectron transporter RsxC
*DH12,169,5950GT22.9% 340.5 / 13.1 236V655V (GTC→GTArsxCelectron transporter RsxC

Marginal new junction evidence (lowest skew 10 of 38 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? DH1 2668875 =353 (0.750)88 (0.190) 75/268 19.6 20.1% pseudogene (20/345 nt) pinE hypothetical protein
?DH1 2670674 = NA (NA)noncoding (1/242 nt) repeat_region repeat region
* ? DH1 = 2668875402 (0.760)75 (0.160) 65/268 23.0 16.5% pseudogene (20/345 nt) pinE hypothetical protein
?DH1 = 2670672 399 (0.850)coding (306/789 nt) stfP integrase
* ? DH1 1355966 =202 (0.380)152 (0.300)
+CTTC
58/292 28.4 37.4% noncoding (76/76 nt) ECDH1_RS06685 tRNA‑Val
?DH1 = 3099750 321 (0.610)noncoding (1/76 nt) ECDH1_RS15360 tRNA‑Lys
* ? DH1 1355966 =202 (0.380)229 (0.450)
+CTTC
58/292 28.4 53.8% noncoding (76/76 nt) ECDH1_RS06685 tRNA‑Val
?DH1 = 3099542 NA (NA)noncoding (268/268 nt) repeat_region repeat region
* ? DH1 1355966 =202 (0.380)129 (0.250)
+CTTC
56/292 29.2 25.2% noncoding (76/76 nt) ECDH1_RS06685 tRNA‑Val
?DH1 = 3100565 586 (1.110)noncoding (1/76 nt) ECDH1_RS15385 tRNA‑Lys
* ? DH1 = 1355505502 (0.950)83 (0.180) 47/260 29.3 20.2% noncoding (54/348 nt) repeat_region repeat region
?DH1 3099693 = 220 (0.480)noncoding (58/76 nt) ECDH1_RS15360 tRNA‑Lys
* ? DH1 = 3099328447 (0.850)75 (0.160) 41/260 32.1 19.8% noncoding (54/268 nt) repeat_region repeat region
?DH1 3099693 = 220 (0.480)noncoding (58/76 nt) ECDH1_RS15360 tRNA‑Lys
* ? DH1 = 2590941455 (0.860)81 (0.180) 37/260 34.1 23.2% intergenic (+1/‑209) ECDH1_RS12805/ECDH1_RS12810 tRNA‑Tyr/tRNA‑Tyr
?DH1 2591414 = 142 (0.310)noncoding (199/555 nt) repeat_region repeat region
* ? DH1 3183985 =509 (0.970)74 (0.150)
+GATTGTTTCTTC
36/276 36.5 19.9% intergenic (+379/‑1) asnB/ECDH1_RS15770 asparagine synthetase B/tRNA‑Met
?DH1 = 3184256 141 (0.270)intergenic (+2/‑33) ECDH1_RS15780/ECDH1_RS15785 tRNA‑Gln/tRNA‑Gln
* ? DH1 = 1519281367 (0.700)52 (0.110)
+AACCTTCTAATT
32/276 38.7 13.4% coding (788/900 nt) yfaD transposase
?DH1 1519282 = 367 (0.700)coding (787/900 nt) yfaD transposase