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 38 shown)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,993,5330AG69.1% 747.4 / inf 492L55S (TTG→TCG) rcdADNA‑binding transcriptional regulator
*DH1936,1690T.40.9% 227.1 / 176.8 112intergenic (‑13/+21)ECDH1_RS04615/ECDH1_RS04620tRNA‑Met/tRNA‑Met
*DH1936,2970AG40.0% 197.6 / 189.2 167intergenic (‑31/+3)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
*DH1936,2781.T38.8% 117.3 / 139.9 67intergenic (‑12/+22)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
*DH1936,2930T.38.3% 260.8 / 192.5 131intergenic (‑27/+7)ECDH1_RS04620/ECDH1_RS04625tRNA‑Met/tRNA‑Met
*DH13,906,8320CA31.7% 180.4 / 18.8 117intergenic (+3/‑32)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
*DH11,065,6050AT31.1% 308.6 / 148.5 226intergenic (+40/‑23)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
*DH13,448,5260AC29.5% 732.7 / 119.4 414intergenic (‑86/+70)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
*DH12,169,5480GT28.6% 509.4 / 62.1 274P671Q (CCG→CAG) rsxCelectron transporter RsxC
*DH12,591,5130TG28.5% 216.8 / 54.4 261D86A (GAC→GCC) ECDH1_RS12815hypothetical protein
*DH12,169,5770TG26.7% 311.3 / 75.1 272R661R (CGA→CGCrsxCelectron transporter RsxC
*DH1613,5970GC26.3% 410.3 / 60.1 280intergenic (+458/+565)garK/yhaCglycerate 2‑kinase/pentapeptide repeats‑related protein
*DH12,326,2870GT26.3% 1021.1 / 113.7 522intergenic (+207/+200)yddM/fdnItranscriptional regulator/formate dehydrogenase nitrate‑inducible cytochrome b556(Fdn) subunit
*DH11,065,6180AG25.8% 395.6 / 157.3 263intergenic (+53/‑10)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
*DH1613,7700CA25.4% 327.8 / 24.1 211intergenic (+631/+392)garK/yhaCglycerate 2‑kinase/pentapeptide repeats‑related protein
*DH13,906,8370CG25.4% 228.8 / 25.1 129intergenic (+8/‑27)ECDH1_RS19315/ECDH1_RS19320tRNA‑Leu/tRNA‑Leu
*DH13,448,5230AC24.8% 888.7 / 81.4 414intergenic (‑83/+73)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
*DH11,065,6140AG23.8% 389.2 / 143.2 250intergenic (+49/‑14)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
*DH12,169,7200CG23.5% 159.6 / 28.0 142E614Q (GAA→CAA) rsxCelectron transporter RsxC
*DH1102,6180GC23.2% 889.8 / 56.4 328intergenic (‑70/+128)trmL/lldDtRNA (cytosine(34)‑2'‑O)‑methyltransferase TrmL/alpha‑hydroxy‑acid oxidizing enzyme

Marginal new junction evidence (lowest skew 10 of 63 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? DH1 2668875 =228 (0.510)209 (0.470) 147/268 3.1 48.0% pseudogene (20/345 nt) pinE hypothetical protein
?DH1 2670674 = NA (NA)noncoding (1/242 nt) repeat_region repeat region
* ? DH1 2732072 =370 (0.760)160 (0.330) 124/296 8.7 30.2% coding (17/1041 nt) plsX phosphate acyltransferase
?DH1 3271770 = NA (NA)noncoding (3/1345 nt) repeat_region repeat region
* ? DH1 = 161689466 (0.940)165 (0.340) 123/296 8.9 29.6% intergenic (+254/+25) glyS/insK glycine‑‑tRNA ligase subunit beta/transposase
?DH1 258177 = 326 (0.670)intergenic (+93/+72) yhhJ/yhhI inner membrane transport permease YhhJ/H repeat‑associated protein yhhI
* ? DH1 = 2732077361 (0.750)167 (0.350) 118/292 9.5 31.6% coding (12/1041 nt) plsX phosphate acyltransferase
?DH1 = 3855990 NA (NA)noncoding (1346/1346 nt) repeat_region repeat region
* ? DH1 3970988 =NA (NA)170 (0.340) 113/300 11.2 24.7% noncoding (1/770 nt) repeat_region repeat region
?DH1 = 4170051 519 (1.050)coding (1162/1422 nt) melB melibiose carrier protein
* ? DH1 4079177 =NA (NA)147 (0.300) 110/300 11.8 21.1% noncoding (2/771 nt) repeat_region repeat region
?DH1 4104586 = 551 (1.110)coding (721/732 nt) rlmB 23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
* ? DH1 2830573 =NA (NA)151 (0.310) 106/298 12.5 22.7% noncoding (4/773 nt) repeat_region repeat region
?DH1 = 4170050 513 (1.040)coding (1163/1422 nt) melB melibiose carrier protein
* ? DH1 3592224 =NA (NA)135 (0.270) 96/298 14.9 19.8% noncoding (183/951 nt) repeat_region repeat region
?DH1 4104586 = 551 (1.110)coding (721/732 nt) rlmB 23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
* ? DH1 = 3581520NA (NA)152 (0.310) 58/300 26.7 22.7% noncoding (770/976 nt) repeat_region repeat region
?DH1 = 4170051 519 (1.050)coding (1162/1422 nt) melB melibiose carrier protein
* ? DH1 3183985 =437 (0.880)97 (0.210)
+GATTGTTTCTTC
46/276 29.0 25.6% intergenic (+379/‑1) asnB/ECDH1_RS15770 asparagine synthetase B/tRNA‑Met
?DH1 = 3184256 175 (0.350)intergenic (+2/‑33) ECDH1_RS15780/ECDH1_RS15785 tRNA‑Gln/tRNA‑Gln