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 28 shown) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | DH1 | 936,169 | 0 | T | . | 40.5% | 157.3 / 116.3 | 80 | intergenic (‑13/+21) | ECDH1_RS04615/ECDH1_RS04620 | tRNA‑Met/tRNA‑Met |
* | DH1 | 3,906,837 | 0 | C | G | 33.3% | 104.2 / 19.6 | 76 | intergenic (+8/‑27) | ECDH1_RS19315/ECDH1_RS19320 | tRNA‑Leu/tRNA‑Leu |
* | DH1 | 936,293 | 0 | T | . | 32.0% | 229.6 / 120.2 | 99 | intergenic (‑27/+7) | ECDH1_RS04620/ECDH1_RS04625 | tRNA‑Met/tRNA‑Met |
* | DH1 | 613,770 | 0 | C | A | 30.1% | 246.9 / 21.1 | 177 | intergenic (+631/+392) | garK/yhaC | glycerate 2‑kinase/pentapeptide repeats‑related protein |
* | DH1 | 2,169,548 | 0 | G | T | 29.6% | 351.1 / 71.8 | 198 | P671Q (CCG→CAG) | rsxC | electron transporter RsxC |
* | DH1 | 2,169,577 | 0 | T | G | 29.3% | 197.7 / 53.8 | 188 | R661R (CGA→CGC) | rsxC | electron transporter RsxC |
* | DH1 | 3,448,526 | 0 | A | C | 28.6% | 578.4 / 84.8 | 345 | intergenic (‑86/+70) | malZ/proY | alpha‑glycosidase/proline‑specific permease ProY |
* | DH1 | 613,597 | 0 | G | C | 28.3% | 281.0 / 38.7 | 206 | intergenic (+458/+565) | garK/yhaC | glycerate 2‑kinase/pentapeptide repeats‑related protein |
* | DH1 | 936,297 | 0 | A | G | 27.6% | 257.0 / 107.3 | 138 | intergenic (‑31/+3) | ECDH1_RS04620/ECDH1_RS04625 | tRNA‑Met/tRNA‑Met |
* | DH1 | 3,184,585 | 0 | C | A | 26.9% | 245.1 / 34.8 | 135 | intergenic (+8/‑30) | ECDH1_RS15795/ECDH1_RS15800 | tRNA‑Gln/tRNA‑TERM |
* | DH1 | 2,326,287 | 0 | G | T | 26.5% | 746.6 / 82.3 | 376 | intergenic (+207/+200) | yddM/fdnI | transcriptional regulator/formate dehydrogenase nitrate‑inducible cytochrome b556(Fdn) subunit |
* | DH1 | 613,613 | 0 | T | G | 24.8% | 177.4 / 23.6 | 165 | intergenic (+474/+549) | garK/yhaC | glycerate 2‑kinase/pentapeptide repeats‑related protein |
* | DH1 | 2,591,513 | 0 | T | G | 24.1% | 242.2 / 38.7 | 241 | D86A (GAC→GCC) | ECDH1_RS12815 | hypothetical protein |
* | DH1 | 613,601 | 0 | G | T | 24.1% | 288.1 / 16.9 | 194 | intergenic (+462/+561) | garK/yhaC | glycerate 2‑kinase/pentapeptide repeats‑related protein |
* | DH1 | 936,278 | 1 | . | T | 24.1% | 138.4 / 65.7 | 54 | intergenic (‑12/+22) | ECDH1_RS04620/ECDH1_RS04625 | tRNA‑Met/tRNA‑Met |
* | DH1 | 4,122,124 | 0 | C | G | 23.4% | 309.9 / 38.8 | 245 | intergenic (‑15/+22) | ECDH1_RS20430/ECDH1_RS20435 | tRNA‑Gly/tRNA‑Gly |
* | DH1 | 3,447,403 | 0 | A | C | 22.2% | 612.5 / 53.1 | 293 | G346G (GGT→GGG) | malZ | alpha‑glycosidase |
* | DH1 | 613,589 | 0 | G | T | 22.2% | 402.8 / 13.8 | 223 | intergenic (+450/+573) | garK/yhaC | glycerate 2‑kinase/pentapeptide repeats‑related protein |
* | DH1 | 2,169,598 | 0 | G | T | 22.0% | 254.5 / 10.3 | 160 | A654A (GCC→GCA) | rsxC | electron transporter RsxC |
* | DH1 | 2,965,906 | 0 | T | G | 21.8% | 242.5 / 20.2 | 192 | intergenic (+139/‑356) | aqpZ/ybjE | aquaporin/transporter |
Marginal new junction evidence (lowest skew 10 of 36 shown) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | DH1 | 1355966 = | 152 (0.390) | 116 (0.310) +CTTC |
51/292 | 24.9 | 37.0% | noncoding (76/76 nt) | ECDH1_RS06685 | tRNA‑Val |
? | DH1 | = 3099750 | 253 (0.650) | noncoding (1/76 nt) | ECDH1_RS15360 | tRNA‑Lys | |||||
* | ? | DH1 | 1355966 = | 152 (0.390) | 90 (0.240) +CTTC |
49/292 | 25.7 | 23.8% | noncoding (76/76 nt) | ECDH1_RS06685 | tRNA‑Val |
? | DH1 | = 3100565 | 441 (1.130) | noncoding (1/76 nt) | ECDH1_RS15385 | tRNA‑Lys | |||||
* | ? | DH1 | 1355966 = | 152 (0.390) | 187 (0.490) +CTTC |
49/292 | 25.7 | 55.8% | noncoding (76/76 nt) | ECDH1_RS06685 | tRNA‑Val |
? | DH1 | = 3099542 | NA (NA) | noncoding (268/268 nt) | repeat_region | repeat region | |||||
* | ? | DH1 | 1355844 = | 118 (0.300) | 78 (0.200) +C |
46/298 | 27.6 | 29.7% | intergenic (‑45/+2) | ECDH1_RS06675/ECDH1_RS06680 | tRNA‑Val/tRNA‑Val |
? | DH1 | = 3099750 | 254 (0.650) | noncoding (1/76 nt) | ECDH1_RS15360 | tRNA‑Lys | |||||
* | ? | DH1 | = 1355505 | 365 (0.940) | 62 (0.180) | 37/260 | 27.6 | 20.0% | noncoding (54/348 nt) | repeat_region | repeat region |
? | DH1 | 3099693 = | 181 (0.540) | noncoding (58/76 nt) | ECDH1_RS15360 | tRNA‑Lys | |||||
* | ? | DH1 | = 3970042 | 448 (1.150) | 42 (0.110) | 38/282 | 29.6 | 9.4% | intergenic (‑128/‑86) | fimA/ECDH1_RS19620 | type‑1 fimbrial protein, A chain/hypothetical protein |
? | DH1 | = 3970338 | 385 (1.050) | intergenic (+22/+58) | ECDH1_RS19620/fimE | hypothetical protein/tyrosine recombinase | |||||
* | ? | DH1 | 3183985 = | 360 (0.920) | 68 (0.190) +GATTGTTTCTTC |
35/276 | 30.4 | 22.4% | intergenic (+379/‑1) | asnB/ECDH1_RS15770 | asparagine synthetase B/tRNA‑Met |
? | DH1 | = 3184256 | 151 (0.390) | intergenic (+2/‑33) | ECDH1_RS15780/ECDH1_RS15785 | tRNA‑Gln/tRNA‑Gln | |||||
* | ? | DH1 | = 1519086 | 248 (0.640) | 46 (0.120) | 35/300 | 33.0 | 16.1% | coding (71/186 nt) | ECDH1_RS07485 | hypothetical protein |
? | DH1 | 1519099 = | 232 (0.600) | coding (58/186 nt) | ECDH1_RS07485 | hypothetical protein | |||||
* | ? | DH1 | = 1519281 | 297 (0.760) | 41 (0.110) +AACCTTCTAATT |
30/276 | 33.1 | 13.1% | coding (788/900 nt) | yfaD | transposase |
? | DH1 | 1519282 = | 297 (0.760) | coding (787/900 nt) | yfaD | transposase | |||||
* | ? | DH1 | 1355966 = | 152 (0.390) | 92 (0.240) +CTC |
30/294 | 35.0 | 31.6% | noncoding (76/76 nt) | ECDH1_RS06685 | tRNA‑Val |
? | DH1 | = 3099750 | 254 (0.650) | noncoding (1/76 nt) | ECDH1_RS15360 | tRNA‑Lys |