Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,677,494 G→C 51.9% V19V (GTC→GTG ndk ← multifunctional nucleoside diphosphate kinase, apyrimidinic endonuclease, and 3'‑phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,677,4940GC51.9% 5.3 / 94.3 79V19V (GTC→GTGndkmultifunctional nucleoside diphosphate kinase, apyrimidinic endonuclease, and 3'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (29/9);  new base C (19/22);  total (48/31)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

AGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCATTTTTACCTCTG  >  minE/1677426‑1677561
                                                                    |                                                                   
aGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATg                                                                     <  1:2188969/69‑1 (MQ=255)
aGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATg                                                                     <  1:802414/69‑1 (MQ=255)
  gTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGAc                                                                   >  1:3430328/1‑69 (MQ=255)
  gTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGAc                                                                   >  1:2933903/1‑69 (MQ=255)
   tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                                  >  1:125211/1‑69 (MQ=255)
   tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                                  >  1:2273960/1‑69 (MQ=255)
   tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                                  >  1:2708920/1‑69 (MQ=255)
   tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                                  >  1:1328800/1‑69 (MQ=255)
   tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                                  >  1:1402522/1‑69 (MQ=255)
   tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                                  >  1:1487889/1‑69 (MQ=255)
   tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                                  >  1:3343596/1‑69 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGCCGt                                                                 >  1:1507810/1‑69 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                                  >  1:2543781/1‑68 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                                 >  1:2535205/1‑69 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                                 >  1:2403358/1‑69 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                                 >  1:311800/1‑69 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                                 >  1:1960778/1‑69 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                                 <  1:188217/69‑1 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                                 >  1:1017816/1‑69 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                                 >  1:524424/1‑69 (MQ=255)
    gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                                 >  1:1189148/1‑69 (MQ=255)
       gcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGtttt                                                              >  1:420252/1‑69 (MQ=255)
       gcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGtttt                                                              >  1:254317/1‑69 (MQ=255)
       gcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGtttt                                                              >  1:1809595/1‑69 (MQ=255)
            tttgTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:2681815/4‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGc                                                          >  1:1154584/1‑68 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:2705591/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:2785385/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:2914142/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:2486354/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:2921944/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:293515/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:2113284/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:314740/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:1577731/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:1436286/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:3195797/1‑69 (MQ=255)
            ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                         >  1:2106367/1‑69 (MQ=255)
            ttGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGc                                                          >  1:1246381/1‑68 (MQ=255)
                  ccAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACcgcg                                                   >  1:2701101/1‑69 (MQ=255)
                  ccAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACcgcg                                                   >  1:2289291/1‑69 (MQ=255)
                  ccAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACcgcg                                                   >  1:2554825/1‑69 (MQ=255)
                      caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                               >  1:1307157/1‑69 (MQ=255)
                      caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                               >  1:1622639/1‑69 (MQ=255)
                      caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                               >  1:553566/1‑69 (MQ=255)
                      caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                               >  1:3208004/1‑69 (MQ=255)
                      caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                               >  1:131195/1‑69 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:3217147/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:423649/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:809393/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:36367/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:2744907/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:936384/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:417016/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:3151051/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:3032204/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:2657581/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:1688215/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:1132268/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:1893267/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:1953831/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:1974717/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:1987740/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:2026403/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:210093/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:2447676/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:2459526/68‑1 (MQ=255)
                               aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCCCGTTCGGTTtgatg                                       <  1:2575595/68‑1 (MQ=255)
                                cccTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                       <  1:1657031/67‑1 (MQ=255)
                                 ccTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTgatgat                                     >  1:2538209/1‑68 (MQ=255)
                                    gCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTCCCGCGtt                                                 >  1:1372964/1‑53 (MQ=255)
                                     cAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGa                                  >  1:15723/1‑67 (MQ=255)
                                               gcgcgcAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTc                      <  1:2143788/69‑1 (MQ=255)
                                               gcgcgcAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTc                      <  1:901467/69‑1 (MQ=255)
                                               gcgcgcAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTc                      <  1:2284113/69‑1 (MQ=255)
                                                cgcgcAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTc                      <  1:3348606/68‑1 (MQ=255)
                                                            ttACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCAttttt         >  1:1285311/1‑69 (MQ=255)
                                                            ttACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCAttttt         <  1:1094813/69‑1 (MQ=255)
                                                            ttACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCAttttt         <  1:1735547/69‑1 (MQ=255)
                                                            ttACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCAttttt         <  1:1787871/69‑1 (MQ=255)
                                                                cAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCATTTTTAcc      <  1:2254522/68‑1 (MQ=255)
                                                                    gACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCATTTTTACCTCTg  >  1:1502718/1‑68 (MQ=255)
                                                                    gACGTTTTTTGCTAACGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCATTTTTACCTCTg  >  1:2631666/1‑68 (MQ=255)
                                                                    |                                                                   
AGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCATTTTTACCTCTG  >  minE/1677426‑1677561

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: