Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,725,672 (C)5→4 coding (75/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,6680C.97.7% 203.7 / ‑3.5 44coding (79/1299 nt)kgtPalpha‑ketoglutarate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (23/20);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTAA  >  NC_000913/2725524‑2725810
                                                                                                                                                |                                                                                                                                              
tCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGccccac                                                                                                                                           >  1:304832/1‑148 (MQ=255)
   ggAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATc                                                                                                                                      >  1:470227/1‑151 (MQ=255)
   ggAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATc                                                                                                                                      <  1:333683/151‑1 (MQ=255)
     aaTCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGc                                                                                                                                    >  1:491400/1‑151 (MQ=255)
        cccGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGccc                                                                                                                                  <  2:537760/150‑1 (MQ=255)
          cGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCaaa                                                                                                                               <  2:288294/151‑1 (MQ=255)
            cagcaAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATg                                                                                                                             <  1:165677/151‑1 (MQ=255)
               caAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATgcgg                                                                                                                          >  2:264994/1‑151 (MQ=255)
                  aaaCACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGcga                                                                                                                       <  1:228726/151‑1 (MQ=255)
                     cacCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGcgacg                                                                                                                     >  1:664379/1‑150 (MQ=255)
                        ctgctgTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGt                                                                                                                  >  1:594314/1‑150 (MQ=255)
                                 gtagtaGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACtt                                                                                                        >  1:335978/1‑151 (MQ=255)
                                     tagtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTc                                                                                                     >  1:26950/1‑150 (MQ=255)
                                        ttGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGtt                                                                                                 >  1:345025/1‑151 (MQ=255)
                                          gAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTg                                                                                               >  1:331022/1‑151 (MQ=255)
                                               gtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTc                                                                                          >  1:551005/1‑151 (MQ=255)
                                                   tgtTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGcc                                                                                      >  1:331642/1‑151 (MQ=255)
                                                            aagggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTa                                                                             <  2:180366/151‑1 (MQ=255)
                                                                  agaagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCa                                                                       <  1:247066/151‑1 (MQ=255)
                                                                  agaagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCa                                                                       <  1:26946/151‑1 (MQ=255)
                                                                    aagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgc                                                                     >  2:22131/1‑151 (MQ=255)
                                                                    aagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCAATCGACCAGATTACCTGAAGAGGCCACACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgc                                                                     >  2:22151/1‑151 (MQ=255)
                                                                       aTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccat                                                                  >  1:400034/1‑151 (MQ=255)
                                                                              cAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCAtt                                                            >  1:539055/1‑150 (MQ=255)
                                                                               aaaGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTAt                                                          <  2:675249/151‑1 (MQ=255)
                                                                               aaaGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTAt                                                          <  2:224613/151‑1 (MQ=255)
                                                                                 aGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGt                                                        >  1:584171/1‑151 (MQ=255)
                                                                                     gagTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGAATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTcc                                                    >  2:182648/1‑151 (MQ=255)
                                                                                       gTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTg                                                  <  2:329044/151‑1 (MQ=255)
                                                                                               aaCGAGTAGACATAGAAATCGAACCACTCGACCAGATTCCCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATc                                          <  2:664178/151‑1 (MQ=255)
                                                                                               aaCGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATc                                          <  1:707147/151‑1 (MQ=255)
                                                                                                          aTAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCtttg                                <  2:93458/150‑2 (MQ=255)
                                                                                                           tAGAAATCGACCCCCCCGACCAGATTACCTGAAGAGCCCCCACAATCGCCCAAATGCGGCGACAAGTATCCCTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCtttt                                <  2:550825/149‑1 (MQ=255)
                                                                                                             gAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcg                            >  1:141314/1‑151 (MQ=255)
                                                                                                              aaaTCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcgg                           <  2:144615/151‑1 (MQ=255)
                                                                                                              aaaTCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcgg                           >  1:6633/1‑151 (MQ=255)
                                                                                                                  cGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATgaaagtactgtaac                        <  2:399901/150‑14 (MQ=255)
                                                                                                                    aaCCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCtt                     >  2:360186/1‑151 (MQ=255)
                                                                                                                      ccACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCtttt                   <  1:261693/151‑1 (MQ=255)
                                                                                                                          tCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGttt               >  2:246168/1‑151 (MQ=255)
                                                                                                                           cGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGttt               <  2:127292/150‑1 (MQ=255)
                                                                                                                           cGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGtt                <  1:317273/149‑1 (MQ=255)
                                                                                                                               cAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAtttt          >  2:33700/1‑151 (MQ=255)
                                                                                                                               cAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAtttt          <  2:304803/151‑1 (MQ=255)
                                                                                                                                         gAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTaa  >  2:665470/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTAA  >  NC_000913/2725524‑2725810

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