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)
         gAATCCCGCAGCAAAAACACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCg                                                                                                                                           <  2:79890/151‑1 (MQ=255)
                     aaaaaCACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATgcggc                                                                                                                               >  2:439854/1‑151 (MQ=255)
                       aaaCACCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGcga                                                                                                                             <  2:100589/151‑1 (MQ=255)
                          cacCTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGcgacg                                                                                                                           >  1:291185/1‑150 (MQ=255)
                            ccTGCTGTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGt                                                                                                                        <  1:431446/151‑1 (MQ=255)
                               gctgTTTGTAGTAGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTa                                                                                                                       <  1:473423/149‑1 (MQ=255)
                                      gtagtaGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCactact                                                                                                               <  2:402623/150‑1 (MQ=255)
                                      gtagtaGTTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACtt                                                                                                              >  1:55612/1‑151 (MQ=255)
                                        agtagtTGAGTCGTCGTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGt                                                                                                            >  2:56950/1‑151 (MQ=255)
                                                  tcgtcgTGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATAGCACAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTg                                                                                                  >  2:463558/1‑151 (MQ=255)
                                                    gtcgtGTTCCCGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGt                                                                                                 <  2:158556/150‑1 (MQ=255)
                                                            cccGGAAGGGAAGAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGtt                                                                                         <  1:550730/150‑1 (MQ=255)
                                                                 aagggaagAAGATGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTa                                                                                   <  2:535090/151‑1 (MQ=255)
                                                                        gaagaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTc                                                                              >  2:149499/1‑149 (MQ=255)
                                                                           gaTGTGGGCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCCGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgcc                                                                          >  1:204706/1‑150 (MQ=255)
                                                                                gggCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAgccatgcc                                                                     <  1:556011/150‑1 (MQ=255)
                                                                                  gCAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCAtt                                                                  >  1:273937/1‑151 (MQ=255)
                                                                                   cAAAGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCAtt                                                                  >  1:105252/1‑150 (MQ=255)
                                                                                     aaGTAGAGTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATg                                                               >  2:394400/1‑151 (MQ=255)
                                                                                           agTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCt                                                         <  1:580571/151‑1 (MQ=255)
                                                                                            gTGAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTg                                                        >  2:69830/1‑151 (MQ=255)
                                                                                             tGAACAGAACGAGTAGACATAGAAATCGAACCACACGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTg                                                        <  1:95468/150‑1 (MQ=255)
                                                                                              gAACAGAACGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGc                                                       <  2:13657/150‑1 (MQ=255)
                                                                                                      cGAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGATTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCg                                              <  2:621356/151‑1 (MQ=255)
                                                                                                       gAGTAGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCg                                              <  1:356207/150‑1 (MQ=255)
                                                                                                           aGACATAGAAATCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCt                                         >  2:365404/1‑151 (MQ=255)
                                                                                                                   aaaTCGAACCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTgtcg                                  <  1:88682/150‑1 (MQ=255)
                                                                                                                         aaCCACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCtt                           >  1:148333/1‑151 (MQ=255)
                                                                                                                           ccACTCGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCtttt                         >  1:467510/1‑151 (MQ=255)
                                                                                                                                cGACCAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGCCGAGTATCCCTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGtt                      <  2:387850/149‑1 (MQ=255)
                                                                                                                                   ccAGATTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAtt                  <  1:540631/150‑1 (MQ=255)
                                                                                                                                       aTTACCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTAtttttt              <  1:283184/150‑1 (MQ=255)
                                                                                                                                          aCCTGAAGAGGCCCCACAATCGCCCAAATGCGGCGACGAGTATCACTACTTGTCAGTTTGCTGTCTGCCGTTACAGTACTTTCAGCCATGCCATTATGTCTCCTGCCGTAATCCGATGCTTTTGTCGGTCGCTTTTGTTTATTTTTTTTg           <  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: