Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,029,310 T→G 62.5% A221A (GCT→GCG malY → bifunctional beta‑cystathionase, PLP‑dependent and regulator of maltose regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,029,3100TG62.5% 68.7 / 56.6 77A221A (GCT→GCGmalYbifunctional beta‑cystathionase, PLP‑dependent and regulator of maltose regulon
Reads supporting (aligned to +/- strand):  ref base T (19/10);  new base G (21/27);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 9.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTG  >  minE/1029245‑1029374
                                                                 |                                                                
tCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCt                                                                  >  1:2880159/1‑66 (MQ=255)
tCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTcg                                                                >  1:1153538/1‑68 (MQ=255)
tCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTcg                                                                >  1:377521/1‑68 (MQ=255)
  cGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTcgcg                                                              <  1:3350597/68‑1 (MQ=255)
            cATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGga                                                            <  1:788278/60‑1 (MQ=255)
              tatggatatggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGga                                                            <  1:438896/58‑1 (MQ=255)
                  gatatggTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTa                                             >  1:1823211/1‑69 (MQ=255)
                                         gCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCaa                      <  1:2308592/69‑1 (MQ=255)
                                         gCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCaa                      <  1:870298/69‑1 (MQ=255)
                                         gCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCaa                      <  1:332972/69‑1 (MQ=255)
                                         gCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGGCGGGCTCGAAAAGTTTCaa                      <  1:595008/69‑1 (MQ=255)
                                              ttCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTc                 >  1:1796706/1‑69 (MQ=255)
                                                   tGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGcc            <  1:322748/69‑1 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:2230893/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:226395/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:3404353/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:2133167/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:2126928/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:2370614/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:2125556/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:2408726/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:3405953/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:463618/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:495727/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:562885/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:745870/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:843544/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:969830/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:196247/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:1082083/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:1413137/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:160232/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:1860617/1‑54 (MQ=255)
                                                   tGGAGTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTCGGGATCAAAAAGt                           >  1:2051095/1‑54 (MQ=255)
                                                       gTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTg         <  1:558033/68‑1 (MQ=255)
                                                       gTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTg         <  1:1598490/68‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTACCGTc                                        <  1:3207790/37‑1 (MQ=38)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:3302178/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1336242/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:2295784/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:50113/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1485225/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1491348/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1269107/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1544940/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:35247/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1589752/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:181419/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:2179984/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1742449/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:3324522/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:2697089/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1810274/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1016440/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1125127/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:2464614/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:2529130/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:798101/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1754112/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:2785373/37‑1 (MQ=255)
                                                       gTAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:755697/37‑1 (MQ=255)
                                                       gTAATGGGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:2934902/37‑1 (MQ=38)
                                                        tAATGTGGCGCGCGGAGACTGGGCGTTGCTAACGTc                                        <  1:1919162/36‑1 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:1194399/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:1183579/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:70096/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:3330679/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:89905/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:1936235/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:541465/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:1542236/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:2923790/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:1800516/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:2627285/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:231430/1‑69 (MQ=255)
                                                          aTGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCg     >  1:2270326/1‑69 (MQ=255)
                                                             tgGCTCGCGGAGACTGGGCGTTGCTATCGTCGGGCTCGAAAAGTTTCTATATTCCCGCCCTGACCGGTg  >  1:2460350/1‑69 (MQ=255)
                                                             tgGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTg  >  1:2675221/1‑69 (MQ=255)
                                                                 |                                                                
TCCGATGAAATCCATATGGATATGGTTTGGGGCGAGCAGCCGCATATTCCCTGGAGTAATGTGGCTCGCGGAGACTGGGCGTTGCTAACGTCGGGCTCGAAAAGTTTCAATATTCCCGCCCTGACCGGTG  >  minE/1029245‑1029374

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: