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 27 shown)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1936,1690T.36.0% 260.1 / 154.5 119intergenic (‑13/+21)ECDH1_RS04615/ECDH1_RS04620tRNA‑Met/tRNA‑Met
*DH12,169,5770TG30.0% 314.8 / 88.8 317R661R (CGA→CGCrsxCelectron transporter RsxC
*DH1936,2970AG29.7% 329.4 / 157.7 190intergenic (‑31/+3)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
*DH13,448,5260AC27.9% 778.0 / 114.3 430intergenic (‑86/+70)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
*DH1936,2930T.27.2% 388.8 / 153.0 156intergenic (‑27/+7)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
*DH12,169,5480GT26.2% 607.4 / 71.6 349P671Q (CCG→CAG) rsxCelectron transporter RsxC
*DH1613,7700CA26.1% 362.7 / 25.2 262intergenic (+631/+392)garK/yhaCglycerate 2‑kinase/pentapeptide repeats‑related protein
*DH13,184,5850CA25.0% 406.2 / 51.3 168intergenic (+8/‑30)ECDH1_RS15795/ECDH1_RS15800tRNA‑Gln/tRNA‑TERM
*DH12,326,2870GT24.9% 1031.5 / 97.1 511intergenic (+207/+200)yddM/fdnItranscriptional regulator/formate dehydrogenase nitrate‑inducible cytochrome b556(Fdn) subunit
*DH12,591,5130TG24.5% 251.5 / 37.0 226D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
*DH11,065,6050AT24.3% 400.9 / 130.7 247intergenic (+40/‑23)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
*DH11,358,8240CA23.1% 791.3 / 53.6 346intergenic (+8/‑32)ECDH1_RS06705/ECDH1_RS06710tRNA‑Ala/tRNA‑Ala
*DH13,906,8370CG22.8% 197.6 / 15.2 112intergenic (+8/‑27)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
*DH11,647,6690AC22.5% 308.7 / 31.2 233intergenic (+248/+125)dusC/ECDH1_RS08050tRNA dihydrouridine synthase DusC/hypothetical protein
*DH11,065,6180AG22.4% 464.3 / 149.3 283intergenic (+53/‑10)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
*DH13,448,5230AC22.4% 955.7 / 71.2 430intergenic (‑83/+73)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
*DH12,965,9060TG22.4% 266.8 / 28.7 205intergenic (+139/‑356)aqpZ/ybjEaquaporin/transporter
*DH13,447,4030AC21.9% 850.8 / 84.4 403G346G (GGT→GGGmalZalpha‑glycosidase
*DH11,065,6140AG21.7% 438.4 / 134.9 274intergenic (+49/‑14)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
*DH12,169,7200CG21.5% 193.4 / 24.2 157E614Q (GAA→CAA) rsxCelectron transporter RsxC

Marginal new junction evidence (lowest skew 10 of 53 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? DH1 = 2668875364 (0.690)209 (0.450) 128/268 6.8 40.3% pseudogene (20/345 nt) pinE hypothetical protein
?DH1 = 2670672 294 (0.630)coding (306/789 nt) stfP integrase
* ? DH1 2668875 =303 (0.640)166 (0.360) 126/268 7.2 35.6% pseudogene (20/345 nt) pinE hypothetical protein
?DH1 2670674 = NA (NA)noncoding (1/242 nt) repeat_region repeat region
* ? DH1 3183985 =509 (0.970)109 (0.230)
+GATTGTTTCTTC
45/276 32.9 26.1% intergenic (+379/‑1) asnB/ECDH1_RS15770 asparagine synthetase B/tRNA‑Met
?DH1 = 3184256 161 (0.310)intergenic (+2/‑33) ECDH1_RS15780/ECDH1_RS15785 tRNA‑Gln/tRNA‑Gln
* ? DH1 1355843 =127 (0.240)107 (0.210)
+C
44/298 36.0 30.0% intergenic (‑44/+3) ECDH1_RS06675/ECDH1_RS06680 tRNA‑Val/tRNA‑Val
?DH1 = 3099750 375 (0.720)noncoding (1/76 nt) ECDH1_RS15360 tRNA‑Lys
* ? DH1 = 2590941488 (0.930)50 (0.110) 28/260 40.1 15.1% intergenic (+1/‑209) ECDH1_RS12805/ECDH1_RS12810 tRNA‑Tyr/tRNA‑Tyr
?DH1 2591414 = 138 (0.300)noncoding (199/555 nt) repeat_region repeat region
* ? DH1 1355843 =127 (0.240)96 (0.180)
+C
36/298 40.1 21.6% intergenic (‑44/+3) ECDH1_RS06675/ECDH1_RS06680 tRNA‑Val/tRNA‑Val
?DH1 = 3100565 576 (1.100)noncoding (1/76 nt) ECDH1_RS15385 tRNA‑Lys
* ? DH1 1355843 =127 (0.240)222 (0.430)
+C
36/298 40.1 63.8% intergenic (‑44/+3) ECDH1_RS06675/ECDH1_RS06680 tRNA‑Val/tRNA‑Val
?DH1 = 3099542 NA (NA)noncoding (268/268 nt) repeat_region repeat region
* ? DH1 1355966 =188 (0.360)61 (0.120)
+CTTC
31/292 42.2 18.2% noncoding (76/76 nt) ECDH1_RS06685 tRNA‑Val
?DH1 = 3099750 375 (0.720)noncoding (1/76 nt) ECDH1_RS15360 tRNA‑Lys
* ? DH1 = 1355505570 (1.090)35 (0.080) 24/260 42.7 8.4% noncoding (54/348 nt) repeat_region repeat region
?DH1 3099693 = 267 (0.590)noncoding (58/76 nt) ECDH1_RS15360 tRNA‑Lys
* ? DH1 = 1519281333 (0.640)33 (0.070)
+AACCTTCTAATT
24/276 44.7 9.7% coding (788/900 nt) yfaD transposase
?DH1 1519282 = 333 (0.640)coding (787/900 nt) yfaD transposase