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.2% 192.5 / ‑4.7 55T74A (ACG→GCG) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/28);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCC  >  NC_000913/4218684‑4218949
                                                                                                                                   |                                                                                                                                      
gTGCGCAGAAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                                 <  2:646111/139‑1 (MQ=255)
     cagcagCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                               <  1:889215/136‑1 (MQ=255)
     cagcagCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                               >  2:889215/1‑136 (MQ=255)
           cGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtac                                                                                                                      >  2:685170/1‑139 (MQ=255)
            gTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                     <  1:455414/139‑1 (MQ=255)
                 gAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                                <  1:306393/139‑1 (MQ=255)
                    cAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTa                                                                                                             <  2:822914/139‑1 (MQ=255)
                      ggCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTaaa                                                                                                           <  1:685170/139‑1 (MQ=255)
                       gCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTa                                                                                                             >  1:788363/1‑136 (MQ=255)
                       gCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTa                                                                                                             <  2:788363/136‑1 (MQ=255)
                                            cagGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGc                                                                                     >  2:532213/1‑139 (MQ=255)
                                              ggCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACAcc                                                                                               <  1:179936/127‑1 (MQ=255)
                                              ggCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACAcc                                                                                               >  2:179936/1‑127 (MQ=255)
                                                 gATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTa                                                                                                   <  1:280461/120‑1 (MQ=255)
                                                 gATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTa                                                                                                   >  2:280461/1‑120 (MQ=255)
                                                             tATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                     <  1:362923/90‑1 (MQ=255)
                                                             tATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                     >  2:362923/1‑90 (MQ=255)
                                                                 cATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCgaga                                                                <  2:183853/139‑1 (MQ=255)
                                                                      ttACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCg                                                           >  2:716563/1‑139 (MQ=255)
                                                                             catcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                                 >  1:436657/1‑62 (MQ=255)
                                                                             catcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                                 <  2:436657/62‑1 (MQ=255)
                                                                             catcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCt                                                    >  1:951967/1‑139 (MQ=255)
                                                                              atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCa                                                                                                                             <  1:766627/65‑1 (MQ=255)
                                                                              atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCa                                                                                                                             >  2:766627/1‑65 (MQ=255)
                                                                              atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                                  <  2:309337/76‑1 (MQ=255)
                                                                              atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                                  >  1:309337/1‑76 (MQ=255)
                                                                              atcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCtt                                                   >  2:593823/1‑139 (MQ=255)
                                                                               tcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCttt                                                  >  2:731398/1‑139 (MQ=255)
                                                                                   gTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTa                                              <  1:716563/139‑1 (MQ=255)
                                                                                      ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                                 <  2:188223/53‑1 (MQ=255)
                                                                                      ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGAcgc                                                                                                                                 >  1:188223/1‑53 (MQ=255)
                                                                                      ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACAcc                                                                                               <  2:688977/87‑1 (MQ=255)
                                                                                      ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACAcc                                                                                               >  1:688977/1‑87 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCa                                                                                                      >  1:96573/1‑78 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCa                                                                                                      <  2:96573/78‑1 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                           >  1:470392/1‑137 (MQ=255)
                                                                                        tctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                           <  2:470392/137‑1 (MQ=255)
                                                                                          tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTa                                              >  2:30344/1‑132 (MQ=255)
                                                                                          tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTa                                              <  1:30344/132‑1 (MQ=255)
                                                                                             tcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTa                                    >  1:315334/1‑139 (MQ=255)
                                                                                                 ggAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGt                                <  1:731398/139‑1 (MQ=255)
                                                                                                           cgcTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTaaa                                                                                                           <  2:173679/54‑1 (MQ=255)
                                                                                                           cgcTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTaaa                                                                                                           >  1:173679/1‑54 (MQ=255)
                                                                                                             cTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGac                    >  1:874646/1‑139 (MQ=255)
                                                                                                                   tACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGgaga                                                       >  2:382111/1‑98 (MQ=255)
                                                                                                                   tACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGgaga                                                       <  1:382111/98‑1 (MQ=255)
                                                                                                                   tACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgc              <  2:173982/139‑1 (MQ=255)
                                                                                                                     cTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGc                                                     <  2:196476/98‑1 (MQ=255)
                                                                                                                     cTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGc                                                     >  1:196476/1‑98 (MQ=255)
                                                                                                                     cTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctc            <  1:593823/139‑1 (MQ=255)
                                                                                                                        aCGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACcgctcgct         <  2:951967/139‑1 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGgaga                                                       >  1:1000120/1‑91 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGgaga                                                       <  2:1000120/91‑1 (MQ=255)
                                                                                                                          ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTa       >  2:674626/1‑139 (MQ=255)
                                                                                                                               aaGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTccc  >  2:598556/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                      
GTGCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTCCC  >  NC_000913/4218684‑4218949

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