Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,725,643 Δ1 bp coding (104/1299 nt) kgtP ← alpha‑ketoglutarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,725,6420C.97.6% 186.6 / ‑3.7 41coding (105/1299 nt)kgtPalpha‑ketoglutarate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (22/18);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTC  >  NC_000913/2725492‑2725787
                                                                                                                                                      |                                                                                                                                                 
gCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAAc                                                                                                                                                   >  2:657462/1‑151 (MQ=255)
 cGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAAca                                                                                                                                                  <  1:231540/151‑2 (MQ=255)
     cAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCg                                                                                                                                             >  2:12673/1‑151 (MQ=255)
        aTAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGACGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGAc                                                                                                                                           >  2:582092/1‑150 (MQ=255)
                 cGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCt                                                                                                                                 >  1:464277/1‑151 (MQ=255)
                             gcATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGccccc                                                                                                                     <  1:138894/151‑1 (MQ=255)
                                 cAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACaa                                                                                                                 <  1:302623/151‑1 (MQ=255)
                                     aaTCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCg                                                                                                              >  2:51319/1‑150 (MQ=255)
                                      aTCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGcc                                                                                                            <  1:51230/151‑1 (MQ=255)
                                      aTCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGcc                                                                                                            >  2:519024/1‑151 (MQ=255)
                                       tCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGccc                                                                                                           >  1:509626/1‑151 (MQ=255)
                                         ccGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCaa                                                                                                         <  1:648845/151‑1 (MQ=255)
                                           gcagcaAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAAt                                                                                                       <  1:98647/151‑1 (MQ=255)
                                                                   agtagtTGAGTCGTCGTGTTCCCGGAAGGTAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGGGTAGACATAGAAATCGAACACGCGGCCAGACTACCCGAAGAGGCCCCCGCAAACGCCCAAATGCGGCGACGAGTTTCACTACTTg                                                                                <  2:225684/150‑1 (MQ=255)
                                                                   agtagtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTg                                                                                <  1:63773/150‑1 (MQ=255)
                                                                       gtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGt                                                                           >  2:580499/1‑151 (MQ=255)
                                                                              cgtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTg                                                                     <  1:62708/150‑1 (MQ=255)
                                                                               gtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTc                                                                   >  1:44248/1‑151 (MQ=255)
                                                                                 cgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCt                                                                  <  1:156669/150‑1 (MQ=255)
                                                                                 cgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTg                                                                 <  1:667612/151‑1 (MQ=255)
                                                                                  gtgtTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACGCTCGACCAGATTGCCTGAAGAGGCCCCCACAATCGGCCAAATGCGGCGACGAGTATCGCTACTTGTCAGTTTGCTGTCt                                                                  <  2:396379/149‑1 (MQ=255)
                                                                                           gaagggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGt                                                       <  1:305707/151‑1 (MQ=255)
                                                                                                gaagaaGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACttt                                                  >  2:9506/1‑151 (MQ=255)
                                                                                                    aagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgc                                              <  2:150251/151‑1 (MQ=255)
                                                                                                      gaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgcc                                             >  1:69496/1‑150 (MQ=255)
                                                                                                       aTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccat                                           >  1:252139/1‑151 (MQ=255)
                                                                                                             gCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCAtt                                     >  1:535017/1‑151 (MQ=255)
                                                                                                                 aGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGt                                 <  2:551932/151‑1 (MQ=255)
                                                                                                                  gTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGtc                                <  1:703046/151‑1 (MQ=255)
                                                                                                                   tAGAGTGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGtct                               >  2:386879/1‑151 (MQ=255)
                                                                                                                        tGAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGc                          >  2:377389/1‑151 (MQ=255)
                                                                                                                         gAACAGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGc                          <  2:379545/150‑1 (MQ=255)
                                                                                                                             aGAACGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGCAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTa                      >  2:439606/1‑150 (MQ=255)
                                                                                                                               aacgatgcgccATAGAAATCGACCACTCTACCAGTTTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCTCCGAACATCTCTACTTGTCAGTGTGCGGTCTGCCGTTACCGTACTTTACGGCAGGCCCTTATGTCACCTGCCGTAAt                    <  2:278668/140‑1 (MQ=255)
                                                                                                                                 cGAGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCg                 <  1:232454/151‑1 (MQ=255)
                                                                                                                                   aGTAGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGAt               >  2:752784/1‑151 (MQ=255)
                                                                                                                                      aGACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCt            >  2:62708/1‑151 (MQ=255)
                                                                                                                                       gACATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCtt           >  2:38249/1‑151 (MQ=255)
                                                                                                                                        aCATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCttt          >  1:608339/1‑151 (MQ=255)
                                                                                                                                        aCATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCttt          >  1:728996/1‑151 (MQ=255)
                                                                                                                                         cATAGAAATCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCtttt         >  1:742283/1‑151 (MQ=255)
                                                                                                                                                aTCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcggtc  >  2:156454/1‑151 (MQ=255)
                                                                                                                                                aTCGAACACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcggtc  <  1:175457/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GCGGCCAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTC  >  NC_000913/2725492‑2725787

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: