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)
     cAAATAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCg                                                                                                                                             >  2:12673/1‑151 (MQ=255)
        aTAGCCAACCGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGACGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGAc                                                                                                                                           >  2:582092/1‑150 (MQ=255)
                 cGCCTATTGGGCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCt                                                                                                                                 >  1:464277/1‑151 (MQ=255)
                             gcATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGccccc                                                                                                                     <  1:138894/151‑1 (MQ=255)
                                 cAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACaa                                                                                                                 <  1:302623/151‑1 (MQ=255)
                                     aaTCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCg                                                                                                              >  2:51319/1‑150 (MQ=255)
                                      aTCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGcc                                                                                                            <  1:51230/151‑1 (MQ=255)
                                      aTCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGcc                                                                                                            >  2:519024/1‑151 (MQ=255)
                                       tCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGccc                                                                                                           >  1:509626/1‑151 (MQ=255)
                                         ccGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCaa                                                                                                         <  1:648845/151‑1 (MQ=255)
                                           gcagcaAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAAt                                                                                                       <  1:98647/151‑1 (MQ=255)
                                                                   agtagtTGAGTCGTCGTGTTCCCGGAAGGTAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGGGTAGACATAGAAATCGAA‑CACGCGGCCAGACTACCCGAAGAGGCCCCCGCAAACGCCCAAATGCGGCGACGAGTTTCACTACTTg                                                                                <  2:225684/150‑1 (MQ=255)
                                                                   agtagtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTg                                                                                <  1:63773/150‑1 (MQ=255)
                                                                       gtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGt                                                                           >  2:580499/1‑151 (MQ=255)
                                                                              cgtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTg                                                                     <  1:62708/150‑1 (MQ=255)
                                                                               gtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTc                                                                   >  1:44248/1‑151 (MQ=255)
                                                                                 cgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCt                                                                  <  1:156669/150‑1 (MQ=255)
                                                                                 cgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTg                                                                 <  1:667612/151‑1 (MQ=255)
                                                                                  gtgtTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CGCTCGACCAGATTGCCTGAAGAGGCCCCCACAATCGGCCAAATGCGGCGACGAGTATCGCTACTTGTCAGTTTGCTGTCt                                                                  <  2:396379/149‑1 (MQ=255)
                                                                                           gaagggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGt                                                       <  1:305707/151‑1 (MQ=255)
                                                                                                gaagaaGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACttt                                                  >  2:9506/1‑151 (MQ=255)
                                                                                                    aagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgc                                              <  2:150251/151‑1 (MQ=255)
                                                                                                      gaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgcc                                             >  1:69496/1‑150 (MQ=255)
                                                                                                       aTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccat                                           >  1:252139/1‑151 (MQ=255)
                                                                                                             gCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCAtt                                     >  1:535017/1‑151 (MQ=255)
                                                                                                                 aGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGt                                 <  2:551932/151‑1 (MQ=255)
                                                                                                                  gTAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGtc                                <  1:703046/151‑1 (MQ=255)
                                                                                                                   tAGAGTGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGtct                               >  2:386879/1‑151 (MQ=255)
                                                                                                                        tGAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGc                          >  2:377389/1‑151 (MQ=255)
                                                                                                                         gAACAGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGc                          <  2:379545/150‑1 (MQ=255)
                                                                                                                             aGAACGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGCAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTa                      >  2:439606/1‑150 (MQ=255)
                                                                                                                               aacgatgcgccATAGAAATCG‑ACCACTCTACCAGTTTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCTCCGAACATCTCTACTTGTCAGTGTGCGGTCTGCCGTTACCGTACTTTACGGCAGGCCCTTATGTCACCTGCCGTAAt                    <  2:278668/140‑1 (MQ=255)
                                                                                                                                 cGAGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCg                 <  1:232454/151‑1 (MQ=255)
                                                                                                                                   aGTAGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGAt               >  2:752784/1‑151 (MQ=255)
                                                                                                                                      aGACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCt            >  2:62708/1‑151 (MQ=255)
                                                                                                                                       gACATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCtt           >  2:38249/1‑151 (MQ=255)
                                                                                                                                        aCATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCttt          >  1:608339/1‑151 (MQ=255)
                                                                                                                                        aCATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCttt          >  1:728996/1‑151 (MQ=255)
                                                                                                                                         cATAGAAATCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCtttt         >  1:742283/1‑151 (MQ=255)
                                                                                                                                                aTCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcggtc  >  2:156454/1‑151 (MQ=255)
                                                                                                                                                aTCGAA‑CACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcggtc  <  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: ‑