Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,101,728 (G)7→8 coding (39/984 nt) kdgT → 2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,101,7211.G95.8% 75.8 / ‑3.1 24coding (32/984 nt)kdgT2‑dehydro‑3‑deoxy‑D‑gluconate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (11/12);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCC  >  NC_000913/4101589‑4101835
                                                                                                                                     |                                                                                                                   
tttATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATccc                                                                                                                      <  1:84160/133‑1 (MQ=255)
  tATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCgg                                                                                                                   <  1:519013/133‑1 (MQ=255)
   aTTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATCCAGGTAAAACGCTCGACTGAGAAAATCCCggg                                                                                                                  >  2:256475/1‑133 (MQ=255)
   aTTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTCATGGAAATCCAGATAAAACGCTCGAGTGAGCAAATCCCggg                                                                                                                  >  2:228549/1‑133 (MQ=255)
         tAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAAGGCTAGATTGAGAAAATCCCgggggggga                                                                                                           >  2:265476/1‑131 (MQ=255)
          aaaTGAATGTCTATCTTCGTTTCCATAAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCCCGATTGAGAAAATCCCggggggggat                                                                                                          <  2:549209/133‑4 (MQ=255)
                 tgTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTc                                                                                                    >  2:434623/1‑133 (MQ=255)
                 tgTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTc                                                                                                    >  1:376849/1‑133 (MQ=255)
                 tgTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGAAAAAACGCTCGATGGAGAAAATCCCGGGGGGCGATGATGCTc                                                                                                    >  2:252125/1‑133 (MQ=255)
                       tCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCg                                                                                              <  1:434623/133‑1 (MQ=255)
                        cTTCGTTGCCATCAACACTTATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGc                                                                                             >  2:359056/1‑133 (MQ=255)
                             tttCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTc                                                                                        >  2:269096/1‑133 (MQ=255)
                             tttCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTc                                                                                        >  1:346629/1‑133 (MQ=255)
                                      acacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCa                                                                               <  1:359056/133‑1 (MQ=255)
                                      acacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCa                                                                               <  1:100543/133‑1 (MQ=255)
                                       cacTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCAc                                                                              <  2:60747/133‑1 (MQ=255)
                                              gCTCCATTGAGGAATTACACATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCca                                                                       <  1:269096/133‑1 (MQ=255)
                                               cTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGc                                                                         <  2:99223/130‑1 (MQ=255)
                                               cTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGc                                                                         >  1:99223/1‑130 (MQ=255)
                                                 ccATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCAc                                                                              >  1:460631/1‑123 (MQ=255)
                                                 ccATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCAc                                                                              <  2:460631/123‑1 (MQ=255)
                                                  cATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACAcc                                                                   >  2:392059/1‑133 (MQ=255)
                                                       aGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTtctc                                                              >  2:50441/1‑133 (MQ=255)
                                                                                atatGCCGCCAGGTGCTGGAAAAGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGa                                     <  1:265476/133‑1 (MQ=255)
                                                                                    gCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGTGCGCGGCAATATTCTGGATCAt                                 >  1:289118/1‑133 (MQ=255)
                                                                                       gACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCa                              <  1:15066/133‑1 (MQ=255)
                                                                                            gTGATGGAAATGCAGGTAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAAc                         <  1:252125/133‑1 (MQ=255)
                                                                                                           aTAAAACGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGAATCATTCACCAAGGGGATGATTACCGGt          >  2:222071/1‑133 (MQ=255)
                                                                                                             aaaaCGCTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTAc        <  2:346629/133‑1 (MQ=255)
                                                                                                                   cTCGATTGAGAAAATCCCGGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGcc  <  2:376849/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                   
TTTATTGATTAAATGAATGTCTATCTTCGTTTCCATCAACACTGATGCTCCATTGAGGAATTACGCATCAGCCCTTAAAAATATGCCGACAGGTGATGGAAATGCAGATAAAACGCTCGATTGAGAAAATCCCGGGGGGGATGATGCTCGTCCCGCTATTCCTTGGCGCACTGTGCCACACCTTCTCGCCGGGGGCGGGGAAATATTTTGGATCATTCACCAACGGGATGATTACCGGTACGGTGCC  >  NC_000913/4101589‑4101835

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