Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,108,056 G→C V169L (GTG→CTG)  pfkA → 6‑phosphofructokinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,108,0560GC96.2% 181.8 / ‑0.2 53V169L (GTG→CTG) pfkA6‑phosphofructokinase I
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (26/25);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAAGTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGAAATCAAAGCGGGTATCGCGAAAGG  >  NC_000913/4107919‑4108189
                                                                                                                                         |                                                                                                                                     
tGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTTCACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTg                                                                                                                                     <  1:196507/140‑1 (MQ=255)
tGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTg                                                                                                                                     >  1:576717/1‑140 (MQ=255)
tGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTg                                                                                                                                     >  2:667666/1‑140 (MQ=255)
 gCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGa                                                                                                                                    <  1:144631/140‑1 (MQ=255)
 gCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGa                                                                                                                                    <  1:627555/140‑1 (MQ=255)
      gCACTATCGACAACGACATCAAAGGCACTGACTACACTATCAGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGc                                                                                                                               <  2:502624/140‑1 (MQ=255)
       cACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAAGTGATGGGcc                                                                                                                              >  2:667232/1‑140 (MQ=255)
       cACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGcc                                                                                                                              >  1:700124/1‑140 (MQ=255)
        aCTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCg                                                                                                                             <  1:554718/140‑1 (MQ=255)
           aTCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGACGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTa                                                                                                                          <  2:602367/140‑1 (MQ=255)
                cAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAAGTGATGGGCCGTTATtgtg                                                                                                                     <  1:667232/140‑1 (MQ=255)
                cAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATtgtg                                                                                                                     >  1:394272/1‑140 (MQ=255)
                 aaCGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTgg                                                                                                                    <  1:561963/140‑1 (MQ=255)
                      cATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGt                                                                                                                            <  2:689127/127‑1 (MQ=255)
                      cATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGt                                                                                                                            >  1:689127/1‑127 (MQ=255)
                      cATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTgg                                                                                                                    >  2:525398/1‑135 (MQ=255)
                      cATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTgg                                                                                                                    <  1:525398/135‑1 (MQ=255)
                      cATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTgg                                                                                                                    >  2:302883/1‑135 (MQ=255)
                      cATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTgg                                                                                                                    <  1:302883/135‑1 (MQ=255)
                                            aTCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGc                                                                                         >  2:679258/1‑140 (MQ=255)
                                              cGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCg                                                                                       >  1:129709/1‑140 (MQ=255)
                                              cGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCg                                                                                       >  1:696868/1‑140 (MQ=255)
                                                gTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCa                                                                                             <  2:208738/132‑1 (MQ=255)
                                                gTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCa                                                                                             >  1:208738/1‑132 (MQ=255)
                                                gTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGt                                                                                     <  1:667666/140‑1 (MQ=255)
                                                    cttcACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCt                                                                                 >  1:257355/1‑140 (MQ=255)
                                                     ttcACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCtg                                                                                >  2:634259/1‑140 (MQ=255)
                                                     ttcACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCtg                                                                                <  2:394272/140‑1 (MQ=255)
                                                       cACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCtgtg                                                                              >  2:183044/1‑140 (MQ=255)
                                                           gcgcTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAAtt                                                                          >  2:535458/1‑140 (MQ=255)
                                                                gAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTtg                                                                     <  1:679258/140‑1 (MQ=255)
                                                                 aGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTtgt                                                                    >  1:596446/1‑140 (MQ=255)
                                                                         tGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCgg                                                            <  1:535458/140‑1 (MQ=255)
                                                                                gCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTtgtg                                                                   <  2:44855/126‑1 (MQ=255)
                                                                                gCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTtgtg                                                                   >  1:44855/1‑126 (MQ=255)
                                                                                gCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGa                                                     >  2:771808/1‑140 (MQ=255)
                                                                                gCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGa                                                     >  2:57609/1‑140 (MQ=255)
                                                                                   aTCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAAtt                                                  >  1:68269/1‑140 (MQ=255)
                                                                                     cGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCa                                                >  2:545799/1‑140 (MQ=255)
                                                                                       aCCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGc                                              <  2:68269/140‑1 (MQ=255)
                                                                                           tCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTg                                          <  2:129709/140‑1 (MQ=255)
                                                                                              gCGTGACTCCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAg                                       <  2:576717/140‑1 (MQ=255)
                                                                                                gTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCTGTGAAGAc                                     <  2:257355/140‑1 (MQ=255)
                                                                                                     acCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGt                                >  1:295582/1‑140 (MQ=255)
                                                                                                          ttcttctCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACg                           >  1:852301/1‑140 (MQ=255)
                                                                                                          ttcttctCACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACg                           <  1:634259/140‑1 (MQ=255)
                                                                                                               ctcACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGAAATc                      <  1:57609/140‑1 (MQ=255)
                                                                                                               ctcACCAGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGAAATc                      <  2:366699/140‑1 (MQ=255)
                                                                                                                     aGCGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGAAATCAAAGCg                <  2:22756/140‑1 (MQ=255)
                                                                                                                       cGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGaaa                        >  2:484269/1‑130 (MQ=255)
                                                                                                                       cGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGaaa                        <  1:484269/130‑1 (MQ=255)
                                                                                                                       cGTATTTCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGAAATCAAAGCggg              >  1:293140/1‑140 (MQ=255)
                                                                                                                            ttCCGTGGTGGAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGAAATCAAAGCGGGTATcg         >  1:781068/1‑140 (MQ=255)
                                                                                                                                   gtggAACTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGAAATCAAAGCGGGTATCGCGAAAgg  <  1:771808/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
TGCCGGGCACTATCGACAACGACATCAAAGGCACTGACTACACTATCGGTTTCTTCACTGCGCTGAGCACCGTTGTAGAAGCGATCGACCGTCTGCGTGACACCTCTTCTTCTCACCAGCGTATTTCCGTGGTGGAAGTGATGGGCCGTTATTGTGGAGATCTGACGTTGGCTGCGGCCATTGCCGGTGGCTGTGAATTCGTTGTGGTTCCGGAAGTTGAATTCAGCCGTGAAGACCTGGTAAACGAAATCAAAGCGGGTATCGCGAAAGG  >  NC_000913/4107919‑4108189

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