Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,776,458 A→G 9.8% Y410H (TAT→CAT)  arpA ← regulator of acetyl CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,776,4580AG 9.8% 170.6 / 10.2 61Y410H (TAT→CAT) arpAregulator of acetyl CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base G (3/3);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATA  >  minE/2776389‑2776527
                                                                     |                                                                     
cAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCatat                                                                      <  1:2915484/71‑1 (MQ=255)
 aGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATc                                                                     <  1:954080/71‑1 (MQ=255)
 aGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATc                                                                     <  1:1778629/71‑1 (MQ=255)
 aGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATc                                                                     <  1:3216379/71‑1 (MQ=255)
  gAAAAACCATTTTTATCTTTTGCCTCCAGAATACGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATcc                                                                    <  1:1325821/71‑1 (MQ=255)
         cATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCt                                                             <  1:3158400/71‑1 (MQ=255)
         cATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCt                                                             <  1:2939884/71‑1 (MQ=255)
          aTTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCt                                                             >  1:2343815/1‑70 (MQ=255)
          aTTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCt                                                             >  1:2752133/1‑70 (MQ=255)
          aTTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTg                                                            <  1:1683292/71‑1 (MQ=255)
             tttATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGAc                                                          <  1:273176/70‑1 (MQ=255)
                aTCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAAt                                                       >  1:2753545/1‑70 (MQ=255)
                 tCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                                      <  1:2787268/70‑1 (MQ=255)
                   ttttGCCTCCAGAATATGCATGAGATTTTTTTTCAAGAGTAATCCTTCATATCCTGTTTCTGACAATGAAt                                                   >  1:1820162/1‑71 (MQ=255)
                     ttGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGa                                                     <  1:1634919/67‑1 (MQ=255)
                     ttGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGa                                                     <  1:304186/67‑1 (MQ=255)
                        ccTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGaa                                               <  1:2074221/70‑1 (MQ=255)
                         cTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGaaaa                                             <  1:494387/71‑1 (MQ=255)
                          tCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGt                                                                 <  1:2387397/50‑1 (MQ=255)
                          tCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCt                                                             <  1:1501633/54‑1 (MQ=255)
                          tCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGa                                                <  1:1134987/67‑1 (MQ=255)
                           ccAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGaaaa                                             <  1:2041092/69‑1 (MQ=255)
                             aGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTg                                                            <  1:1298871/52‑1 (MQ=255)
                             aGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATa                                           >  1:2880293/1‑69 (MQ=255)
                              gAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGa                                                <  1:108573/63‑1 (MQ=255)
                               aaTATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATGTc                                                                     <  1:192364/41‑1 (MQ=255)
                               aaTATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATGTc                                                                     <  1:308458/41‑1 (MQ=255)
                               aaTATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATGTc                                                                     <  1:2327027/41‑1 (MQ=255)
                                 tatGCATGAGATTTTTTTTCGAGAGTAATCCTTCATGTCCt                                                                   >  1:1277888/1‑41 (MQ=255)
                                 tatGCATGAGATTTTTTTTCGAGAGTAATCCTTCATGTCCt                                                                   >  1:1777589/1‑41 (MQ=255)
                                 tatGCATGAGATTTTTTTTCGAGAGTAATCCTTCATGTCCt                                                                   >  1:1552078/1‑41 (MQ=255)
                                 tatGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCa                                     <  1:296540/71‑1 (MQ=255)
                                   tGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGtt                                                                >  1:16576/1‑42 (MQ=255)
                                       tGAGATTTTTTTTCGAGAGTAATCCTTCATATCCt                                                                   <  1:2007593/35‑1 (MQ=255)
                                       tGAGATTTTTTTTCGAGAGTAATCCTTCATATCCt                                                                   <  1:2043692/35‑1 (MQ=255)
                                       tGAGATTTTTTTTCGAGAGTAATCCTTCATATCCt                                                                   <  1:279954/35‑1 (MQ=255)
                                       tGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGa                                                <  1:1908813/54‑1 (MQ=255)
                                       tGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                               <  1:3182976/71‑1 (MQ=255)
                                         agaTTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATc                             >  1:3091958/1‑71 (MQ=255)
                                         agaTTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATc                             >  1:3028894/1‑71 (MQ=255)
                                          gaTTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATcc                            >  1:401410/1‑71 (MQ=255)
                                          gaTTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATcc                            >  1:2880246/1‑71 (MQ=255)
                                           aTTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCg                           >  1:953057/1‑71 (MQ=255)
                                            ttttttttCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGc                          <  1:2380177/71‑1 (MQ=255)
                                             tttttttCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGaaaa                                             >  1:3316219/1‑51 (MQ=255)
                                             tttttttCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGatat                              >  1:2507993/1‑66 (MQ=255)
                                             tttttttCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGc                          >  1:712287/1‑70 (MQ=255)
                                             tttttttCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCa                         >  1:1098866/1‑71 (MQ=255)
                                               tttttCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAt                        <  1:928967/70‑1 (MQ=255)
                                                        agTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAAtt              <  1:2069110/71‑1 (MQ=255)
                                                        agTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAAt               >  1:1504253/1‑70 (MQ=255)
                                                          tAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTta            <  1:2675203/71‑1 (MQ=255)
                                                            aTCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTta            >  1:2845919/1‑69 (MQ=255)
                                                            aTCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTta            >  1:1524194/1‑69 (MQ=255)
                                                             tCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAg         >  1:1900039/1‑71 (MQ=255)
                                                               cTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAc                                   <  1:565857/43‑1 (MQ=255)
                                                               cTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGa                                 >  1:2002387/1‑45 (MQ=255)
                                                               cTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGc        >  1:137986/1‑70 (MQ=255)
                                                               cTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGc        >  1:2175161/1‑70 (MQ=255)
                                                                 tCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAt                        >  1:1435312/1‑52 (MQ=255)
                                                                 tCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAga        >  1:402731/1‑67 (MQ=255)
                                                                  cATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAt                        >  1:3236600/1‑51 (MQ=255)
                                                                   atatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCaaa    <  1:1661006/70‑1 (MQ=255)
                                                                   atatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCaaa    <  1:613269/70‑1 (MQ=255)
                                                                   atatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCaaa    <  1:2874635/70‑1 (MQ=255)
                                                                    tatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAc                                   <  1:2155390/38‑1 (MQ=255)
                                                                    tatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAc                                   <  1:2518012/38‑1 (MQ=255)
                                                                    tatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAc                                   <  1:2573526/38‑1 (MQ=255)
                                                                    tatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAc                                   <  1:1134726/38‑1 (MQ=255)
                                                                    tatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCa                         <  1:950493/48‑1 (MQ=255)
                                                                    tatCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAata  <  1:900584/71‑1 (MQ=255)
                                                                     |                                                                     
CAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATA  >  minE/2776389‑2776527

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: