Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,949,029 T→C 25.8% A40A (GCA→GCG yjjP ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,949,0290TC25.8% 88.9 / 30.4 81A40A (GCA→GCGyjjPpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (19/41);  new base C (19/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 2.53e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

TTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGAAATAA  >  minE/2948966‑2949097
                                                               |                                                                    
tttAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCg                                                                 >  1:2441216/1‑69 (MQ=255)
tttAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCg                                                                 >  1:1474122/1‑69 (MQ=255)
tttAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCg                                                                 >  1:1807594/1‑69 (MQ=255)
    aTAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGAc                                                               <  1:3207823/67‑1 (MQ=255)
    aTAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGAc                                                               <  1:2705771/67‑1 (MQ=255)
    aTAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGAc                                                               <  1:435751/67‑1 (MQ=255)
          gTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTc                                                         >  1:2406012/1‑67 (MQ=255)
             aGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAg                                                    <  1:1019658/69‑1 (MQ=255)
             aGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAg                                                    <  1:3073478/69‑1 (MQ=255)
             aGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAg                                                    <  1:2677414/69‑1 (MQ=255)
             aGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAg                                                    <  1:2411492/69‑1 (MQ=255)
              gCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAg                                                    <  1:636654/68‑1 (MQ=255)
                aCTATGGCGTTCGAATAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGaa                                                 <  1:1533563/69‑1 (MQ=255)
                aCTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGaa                                                 <  1:892442/69‑1 (MQ=255)
                aCTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGaa                                                 <  1:356668/69‑1 (MQ=255)
                aCTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGaa                                                 <  1:1346807/69‑1 (MQ=255)
                aCTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGaa                                                 <  1:1528515/69‑1 (MQ=255)
                aCTATGGCGTGCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGaa                                                 <  1:3257170/69‑1 (MQ=255)
                   aTGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAg                                               >  1:1547165/1‑68 (MQ=255)
                         ttCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGCCCCAGTCGTg                                                      >  1:2607127/1‑55 (MQ=255)
                         ttCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAgg                                                   >  1:135601/1‑58 (MQ=255)
                         ttCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGaa                                                 >  1:1113495/1‑60 (MQ=255)
                         ttCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCtcatca                                        >  1:3044870/1‑69 (MQ=255)
                         ttCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCtcatca                                        >  1:3082763/1‑69 (MQ=255)
                         ttCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCtcatca                                        >  1:398744/1‑69 (MQ=255)
                         ttCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCtcat                                          >  1:2546524/1‑67 (MQ=255)
                         ttCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCtaatca                                        >  1:2696841/1‑69 (MQ=255)
                               gaggTAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCt                                             <  1:1573179/58‑1 (MQ=255)
                               gagaTAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCt                                             <  1:597368/58‑1 (MQ=255)
                               gagaTAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCt                                             <  1:75755/58‑1 (MQ=255)
                               gagaTAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCt                                             <  1:1528759/58‑1 (MQ=255)
                               gagaTAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCt                                             <  1:896593/58‑1 (MQ=255)
                               gagaTAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCt                                             <  1:162228/58‑1 (MQ=255)
                               gagaTAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACACAGTCGTGAGGAAAGCt                                             <  1:1571492/58‑1 (MQ=255)
                                    aGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACgcgc                              >  1:2361965/1‑68 (MQ=255)
                                    aGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACgcgc                              >  1:1921621/1‑68 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:2833823/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:3027435/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:35430/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:2687286/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:264278/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:392147/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:41752/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:1471305/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:2364348/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:1825114/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:1701014/69‑1 (MQ=255)
                                       aCTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTc                          <  1:713739/69‑1 (MQ=255)
                                             gACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCa                      <  1:1016262/67‑1 (MQ=255)
                                             gACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCa                      <  1:99487/67‑1 (MQ=255)
                                             gACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCa                      <  1:2788916/67‑1 (MQ=255)
                                               cGCTGTCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCAt                  <  1:2532768/69‑1 (MQ=255)
                                                gCTGTCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCAt                  <  1:3284522/68‑1 (MQ=255)
                                                    tCCATTCCCAACGCCCGCCCCAGTCGTGAGGAAAGCtcatc                                         >  1:2756318/1‑41 (MQ=37)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATcaaccaa                                   >  1:713827/1‑47 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCAt                  >  1:2318045/1‑64 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:2898252/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:99128/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:1514105/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:1608525/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:169335/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:600206/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:526813/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:2433384/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:2533083/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:3406589/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:337170/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:2768694/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:2929252/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:3167471/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:314488/1‑68 (MQ=255)
                                                    tCCATTCCCAACGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTg              >  1:2849486/1‑68 (MQ=255)
                                                          cccAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGa       >  1:91220/1‑69 (MQ=255)
                                                          cccAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGa       >  1:3176089/1‑69 (MQ=255)
                                                            cAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCttt                           >  1:87721/1‑47 (MQ=255)
                                                            cAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGaa      >  1:1134434/1‑68 (MQ=255)
                                                            cAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGaa      <  1:1288917/68‑1 (MQ=255)
                                                            cAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGaa      <  1:2896322/68‑1 (MQ=255)
                                                            cAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGaa      <  1:1663295/68‑1 (MQ=255)
                                                             aGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGAaat    <  1:472303/69‑1 (MQ=255)
                                                             aGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGAaat    <  1:2735486/69‑1 (MQ=255)
                                                               tGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGAaataa  >  1:2204290/1‑69 (MQ=255)
                                                               tGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGAaataa  >  1:3218172/1‑69 (MQ=255)
                                                               |                                                                    
TTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCGTGAGGAAAGCTCATCAACCAACGCGCTTTCCGCACCATGTTGCAAAAGAAATAA  >  minE/2948966‑2949097

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