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.1% 154.5 / ‑1.4 34coding (79/1299 nt)kgtPalpha‑ketoglutarate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (14/19);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

GCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTAAAGGAAA  >  NC_000913/2725519‑2725816
                                                                                                                                                     |                                                                                                                                                    
ccgcATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGAGGGCAAAGTAGAGTGAACAGAACGAGCAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGcc                                                                                                                                                     >  2:315235/2‑151 (MQ=255)
      cAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGccccaca                                                                                                                                                >  1:10199/1‑147 (MQ=255)
         gAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCg                                                                                                                                           <  2:79890/151‑1 (MQ=255)
                     aaaaaCACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATgcggc                                                                                                                               >  2:439854/1‑151 (MQ=255)
                       aaaCACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGcga                                                                                                                             <  2:100589/151‑1 (MQ=255)
                          cacCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGcgacg                                                                                                                           >  1:291185/1‑150 (MQ=255)
                            ccTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGt                                                                                                                        <  1:431446/151‑1 (MQ=255)
                               gctgTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTa                                                                                                                       <  1:473423/149‑1 (MQ=255)
                                      gtagtaGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCactact                                                                                                               <  2:402623/150‑1 (MQ=255)
                                      gtagtaGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACtt                                                                                                              >  1:55612/1‑151 (MQ=255)
                                        agtagtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGt                                                                                                            >  2:56950/1‑151 (MQ=255)
                                                  tcgtcgTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATAGCACAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTg                                                                                                  >  2:463558/1‑151 (MQ=255)
                                                    gtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGt                                                                                                 <  2:158556/150‑1 (MQ=255)
                                                            cccGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGtt                                                                                         <  1:550730/150‑1 (MQ=255)
                                                                 aagggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTa                                                                                   <  2:535090/151‑1 (MQ=255)
                                                                        gaagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTc                                                                              >  2:149499/1‑149 (MQ=255)
                                                                           gaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCCGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgcc                                                                          >  1:204706/1‑150 (MQ=255)
                                                                                gggCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccatgcc                                                                     <  1:556011/150‑1 (MQ=255)
                                                                                  gCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCAtt                                                                  >  1:273937/1‑151 (MQ=255)
                                                                                   cAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCAtt                                                                  >  1:105252/1‑150 (MQ=255)
                                                                                     aaGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATg                                                               >  2:394400/1‑151 (MQ=255)
                                                                                           agTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCt                                                         <  1:580571/151‑1 (MQ=255)
                                                                                            gTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTg                                                        >  2:69830/1‑151 (MQ=255)
                                                                                             tGAACAGAACGAGTAGACATAGAAATCGAACCACACGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTg                                                        <  1:95468/150‑1 (MQ=255)
                                                                                              gAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGc                                                       <  2:13657/150‑1 (MQ=255)
                                                                                                      cGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGATTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCg                                              <  2:621356/151‑1 (MQ=255)
                                                                                                       gAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCg                                              <  1:356207/150‑1 (MQ=255)
                                                                                                           aGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCt                                         >  2:365404/1‑151 (MQ=255)
                                                                                                                   aaaTCGAACCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcg                                  <  1:88682/150‑1 (MQ=255)
                                                                                                                         aaCCACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCtt                           >  1:148333/1‑151 (MQ=255)
                                                                                                                           ccACTCGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCtttt                         >  1:467510/1‑151 (MQ=255)
                                                                                                                                cGACCAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGCCGAGTATCCCTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGtt                      <  2:387850/149‑1 (MQ=255)
                                                                                                                                   ccAGATTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAtt                  <  1:540631/150‑1 (MQ=255)
                                                                                                                                       aTTACCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAtttttt              <  1:283184/150‑1 (MQ=255)
                                                                                                                                          aCCTGAAGAGG‑CCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTg           <  1:369304/150‑1 (MQ=255)
                                                                                                                                                   ggCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATTTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTAAAGGaaa  <  1:13681/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCATCAGGAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTGTAAAGGAAA  >  NC_000913/2725519‑2725816

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: ‑