Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,218,815 A→G T74A (ACG→GCG)  aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,218,8150AG98.3% 199.6 / NA 58T74A (ACG→GCG) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (27/30);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCT  >  NC_000913/4218684‑4218942
                                                                                                                                   |                                                                                                                               
gTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                          <  1:449682/139‑1 (MQ=255)
  gcgcAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAGGTGCGGACGCgg                                                                                                                        <  1:1034257/139‑1 (MQ=255)
  gcgcAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                        <  2:536297/139‑1 (MQ=255)
     cagcagCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                          <  2:932575/134‑1 (MQ=255)
     cagcagCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                          >  1:932575/1‑134 (MQ=255)
      agcagcGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGTGGACGCGGCAtt                                                                                                                    <  1:440945/139‑1 (MQ=255)
               tCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                           <  2:412448/139‑1 (MQ=255)
                cGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgt                                                                                                          >  2:831294/1‑139 (MQ=255)
                                         cagcagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgcg                                                                                                                         >  1:1058956/1‑99 (MQ=255)
                                         cagcagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGt                                                                                 >  2:106757/1‑139 (MQ=255)
                                         cagcagGCGATGAAAAACCGTATCCATCTTTACGAGCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgcg                                                                                                                         <  2:1058956/99‑1 (MQ=255)
                                                                              atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtt                                            >  1:710812/1‑139 (MQ=255)
                                                                              atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtt                                            >  2:546638/1‑139 (MQ=255)
                                                                              atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtt                                            >  2:730401/1‑139 (MQ=255)
                                                                                caCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                        >  2:104275/1‑61 (MQ=255)
                                                                                caCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                        <  1:104275/61‑1 (MQ=255)
                                                                                      ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCa                                                                                                                      >  2:599268/1‑57 (MQ=255)
                                                                                      ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCa                                                                                                                      <  1:599268/57‑1 (MQ=255)
                                                                                      ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                    <  2:515249/139‑1 (MQ=255)
                                                                                      ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                    <  1:730401/139‑1 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtactac                                                                                                            <  2:278262/65‑1 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtactac                                                                                                            >  1:278262/1‑65 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTAcc                                                                     <  2:544724/104‑1 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTAcc                                                                     >  1:544724/1‑104 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGc                                              <  1:568533/127‑1 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGc                                              >  2:568533/1‑127 (MQ=255)
                                                                                          tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                          <  1:699926/49‑1 (MQ=255)
                                                                                          tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                          <  1:557318/49‑1 (MQ=255)
                                                                                          tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                          >  2:557318/1‑49 (MQ=255)
                                                                                          tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                          >  2:699926/1‑49 (MQ=255)
                                                                                          tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGTCATTTTTACTACTTGTTTATGTACATTACACCCGTCTGTTGCCGGATTACCCCCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCgt                                <  2:381143/139‑1 (MQ=255)
                                                                                            gtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGc                                              >  2:1079297/1‑123 (MQ=255)
                                                                                            gtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGc                                              <  1:1079297/123‑1 (MQ=255)
                                                                                            gtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCgtgt                              >  2:826175/1‑139 (MQ=255)
                                                                                             tcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                   <  1:672417/101‑1 (MQ=255)
                                                                                             tcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACACGGCTGTTGCCGGATTACCcg                                                                   >  2:672417/1‑101 (MQ=255)
                                                                                                 ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                    >  2:54258/1‑128 (MQ=255)
                                                                                                 ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                    <  1:54258/128‑1 (MQ=255)
                                                                                                 ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGt                         >  1:690722/1‑139 (MQ=255)
                                                                                                 ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGt                         >  2:477858/1‑139 (MQ=255)
                                                                                                 ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGt                         >  1:1104544/1‑139 (MQ=255)
                                                                                                     cAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGt                     <  2:531900/139‑1 (MQ=255)
                                                                                                       aCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTa                             >  1:404408/1‑129 (MQ=255)
                                                                                                       aCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTa                             <  2:404408/129‑1 (MQ=255)
                                                                                                           cgcgGCATTACTAACGGCAAAAGTACGGACGTGGTATTTTTAGTACGTGATAAAGAGCATTACATGCGGCCGTTGCCGGTTTACTCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTg               >  2:777774/1‑139 (MQ=255)
                                                                                                           cgcTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTaaa                                                                                                    >  2:540011/1‑54 (MQ=255)
                                                                                                           cgcTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTaaa                                                                                                    <  1:540011/54‑1 (MQ=255)
                                                                                                            gcTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGa              <  2:690722/139‑1 (MQ=255)
                                                                                                               gCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGgaga                                                >  1:817377/1‑102 (MQ=255)
                                                                                                               gCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGgaga                                                <  2:817377/102‑1 (MQ=255)
                                                                                                                  ttACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcg        <  2:1104544/139‑1 (MQ=255)
                                                                                                                        aCGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctcgct  <  1:477858/139‑1 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCa                                                                                               >  2:9910/1‑44 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCa                                                                                               <  1:9910/44‑1 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGagag                                               <  1:111128/92‑1 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGagag                                               >  2:111128/1‑92 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                    >  1:235870/1‑103 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                    <  2:235870/103‑1 (MQ=255)
                                                                                                                                   |                                                                                                                               
GTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCT  >  NC_000913/4218684‑4218942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑