Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,677,486 T→C 54.0% N22S (AAT→AGT)  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,4860TC54.0% 11.1 / 61.7 76N22S (AAT→AGT) ndkmultifunctional nucleoside diphosphate kinase, apyrimidinic endonuclease, and 3'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (28/7);  new base C (19/22);  total (47/29)
Fisher's exact test for biased strand distribution p-value = 4.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

CAACGGTCAGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCATTTTT  >  minE/1677418‑1677554
                                                                    |                                                                    
cAACGGTCAGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGatat                                                                      >  1:3367315/1‑69 (MQ=255)
cAACGGTCAGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGatat                                                                      >  1:3118515/1‑69 (MQ=255)
        aGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATg                                                              <  1:2188969/69‑1 (MQ=255)
        aGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATg                                                              <  1:802414/69‑1 (MQ=255)
          gTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGAc                                                            >  1:2933903/1‑69 (MQ=255)
          gTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGAc                                                            >  1:3430328/1‑69 (MQ=255)
           tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                           >  1:2708920/1‑69 (MQ=255)
           tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                           >  1:125211/1‑69 (MQ=255)
           tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                           >  1:2273960/1‑69 (MQ=255)
           tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                           >  1:1328800/1‑69 (MQ=255)
           tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                           >  1:3343596/1‑69 (MQ=255)
           tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                           >  1:1402522/1‑69 (MQ=255)
           tGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                           >  1:1487889/1‑69 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGCCGt                                                          >  1:1507810/1‑69 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACg                                                           >  1:2543781/1‑68 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                          <  1:188217/69‑1 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                          >  1:1960778/1‑69 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                          >  1:2403358/1‑69 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                          >  1:524424/1‑69 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                          >  1:1017816/1‑69 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                          >  1:2535205/1‑69 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                          >  1:311800/1‑69 (MQ=255)
            gcagcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGt                                                          >  1:1189148/1‑69 (MQ=255)
               gcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGtttt                                                       >  1:254317/1‑69 (MQ=255)
               gcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGtttt                                                       >  1:1809595/1‑69 (MQ=255)
               gcaTTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGtttt                                                       >  1:420252/1‑69 (MQ=255)
                    tttgTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:2681815/4‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGc                                                   >  1:1154584/1‑68 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:314740/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:2486354/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:293515/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:2113284/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:3195797/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:2921944/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:2914142/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:1436286/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:1577731/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:2785385/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:2705591/1‑69 (MQ=255)
                    ttGGTGCCGACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCt                                                  >  1:2106367/1‑69 (MQ=255)
                    ttGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGc                                                   >  1:1246381/1‑68 (MQ=255)
                          ccAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACcgcg                                            >  1:2289291/1‑69 (MQ=255)
                          ccAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACcgcg                                            >  1:2701101/1‑69 (MQ=255)
                          ccAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACcgcg                                            >  1:2554825/1‑69 (MQ=255)
                              caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                        >  1:1622639/1‑69 (MQ=255)
                              caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                        >  1:3208004/1‑69 (MQ=255)
                              caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                        >  1:553566/1‑69 (MQ=255)
                              caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                        >  1:131195/1‑69 (MQ=255)
                              caaTTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCg                                        >  1:1307157/1‑69 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:3217147/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:3032204/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:3151051/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:2459526/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:36367/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:417016/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:423649/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:809393/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:936384/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:2657581/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:1132268/68‑1 (MQ=255)
                                       aCCCTGCAGCTTCAAAGCGCGCAAAGATACTACCAATCACGTTTTTTGCTACCGCGTTCGGTTtgatg                                <  1:1688215/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:2744907/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:2284113/69‑1 (MQ=255)
                                                       gcgcgcAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTc               <  1:901467/69‑1 (MQ=255)
                                                        cgcgcAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTc               <  1:3348606/68‑1 (MQ=255)
                                                                    ttACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCAttttt  <  1:1787871/69‑1 (MQ=255)
                                                                    ttACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCAttttt  <  1:1735547/69‑1 (MQ=255)
                                                                    ttACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCAttttt  >  1:1285311/1‑69 (MQ=255)
                                                                    ttACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCAttttt  <  1:1094813/69‑1 (MQ=255)
                                                                    |                                                                    
CAACGGTCAGGTGCAGCATTTTGGTGCCAACAATTTTGAACCCTGCAGCTTCAAAGCGCGCAAAGATATTACCAATGACGTTTTTTGCTACCGCGTTCGGTTTGATGATGGAAAAAGTACGTTCAATAGCCATTTTT  >  minE/1677418‑1677554

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: