Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,896,815 A→G 38.5% S179S (TCA→TCG argA → fused acetylglutamate kinase homolog (inactive) and amino acid N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,896,8150AG38.5% 32.0 / 45.9 70S179S (TCA→TCGargAfused acetylglutamate kinase homolog (inactive) and amino acid N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (16/27);  new base G (14/13);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATCCATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTCAACTGGCCATCAAA  >  minE/1896756‑1896875
                                                           |                                                            
gcgATCCATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGc                                                     <  1:1408747/69‑1 (MQ=255)
 cgATCCATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGc                                                     <  1:1678124/68‑1 (MQ=255)
   aTCCATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCgag                                                  <  1:2508971/69‑1 (MQ=255)
    tCCATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCgag                                                  <  1:2699109/68‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCTGTCGCTGTTTCAGTCACTGGCGagag                                                <  1:1618318/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:1579521/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:2577316/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:2845885/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:3171627/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:1789391/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:395147/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:467907/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:722996/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:955214/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:1085568/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:610045/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:518346/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGagag                                                <  1:3305150/69‑1 (MQ=255)
     ccATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGagag                                                <  1:1458751/69‑1 (MQ=255)
      cATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCGGTCACTGGCGagag                                                <  1:1955993/68‑1 (MQ=255)
      cATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGagag                                                <  1:560088/68‑1 (MQ=255)
      cATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGagag                                                <  1:174347/68‑1 (MQ=255)
          gTCTACTGGAGAGCGGTGGAATAGTGCTAATGGGGCCTGTCGCTGTTTCAGTCACTGGCGAGTGCTTTa                                           <  1:3371437/69‑1 (MQ=255)
          gTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTa                                           <  1:877240/69‑1 (MQ=255)
               cTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCt                                       <  1:327871/68‑1 (MQ=255)
               cTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCt                                       <  1:3076956/68‑1 (MQ=255)
               cTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCt                                       <  1:1140628/68‑1 (MQ=255)
                          tGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTTACCTCGGAaga                           <  1:1164392/69‑1 (MQ=255)
                          tGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAaga                           <  1:1048837/69‑1 (MQ=255)
                          tGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAaga                           <  1:2047752/69‑1 (MQ=255)
                           gCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAagag                          <  1:2606765/69‑1 (MQ=255)
                            cAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAgaga                         >  1:2091122/1‑69 (MQ=255)
                                 gTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAaga                           <  1:1580450/62‑1 (MQ=255)
                                 gTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAaga                           <  1:2826586/62‑1 (MQ=255)
                                 gTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAaga                           <  1:2783899/62‑1 (MQ=255)
                                 gTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGcc                    >  1:3138342/1‑69 (MQ=255)
                                 gTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGcc                    >  1:34078/1‑69 (MQ=255)
                                 gTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGcc                    >  1:143769/1‑69 (MQ=255)
                                 gTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGcc                    >  1:2484219/1‑69 (MQ=255)
                                 gTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAAACTGACCTCGGAAGAGATTGc                     >  1:286732/1‑68 (MQ=255)
                                     tAATGGGGCCGGTCGCTGTTTCAGTCACTTGCgag                                                  <  1:1351539/35‑1 (MQ=255)
                                          gggCCGGTCGGTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGc                     <  1:1171285/59‑1 (MQ=255)
                                          gggCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGc                     <  1:151580/59‑1 (MQ=255)
                                          gggCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGc                     <  1:1016181/59‑1 (MQ=255)
                                          gggCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGc                     <  1:1931606/59‑1 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGGGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:149216/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCGACt                 >  1:3258356/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:1985307/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:1242987/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:629937/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:1366711/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:2489919/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:2493728/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:464683/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:2076844/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:1989386/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:326352/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:216155/1‑58 (MQ=255)
                                               ggTCGCTGTTTCGGTCACTGGCGAGAGCTTTAATCTGACCTCGGAAGAGATTGCCACt                 >  1:3134856/1‑58 (MQ=255)
                                                  cGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTCAAc             >  1:652058/1‑59 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACt                 >  1:2419889/1‑54 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTc                >  1:2035645/1‑55 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTc                >  1:846883/1‑55 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTc                >  1:2181582/1‑55 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTCAACTGGCCATCaaa  >  1:2960311/1‑69 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTCAACTGGCCATCaaa  >  1:729903/1‑69 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTCAACTGGCCATCaaa  >  1:899325/1‑69 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACCCAACTGGCCATCaaa  >  1:1656649/1‑69 (MQ=255)
                                                   gCTGTTTCAGTCACTGGCGAGAGCCTTAACCTGACCTCGGAAGAGATTGCCACTc                >  1:3216522/1‑55 (MQ=255)
                                                      gTTTCAGTCACTGGCGAGAGCTTTATCCTGACCTCGGAagag                          <  1:232128/42‑1 (MQ=255)
                                                           aGTCACTGGCGAGAGCTTTAACCTGACCTCGGAaga                           <  1:336406/36‑1 (MQ=255)
                                                           |                                                            
GCGATCCATCGTCAACTGGACAGCGGTGCAATAGTGCTAATGGGGCCGGTCGCTGTTTCAGTCACTGGCGAGAGCTTTAACCTGACCTCGGAAGAGATTGCCACTCAACTGGCCATCAAA  >  minE/1896756‑1896875

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: