Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,072,171 C→A V161F (GTC→TTC)  yihT ← 6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,072,1710CA98.6% 258.5 / ‑3.3 73V161F (GTC→TTC) yihT6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (36/36);  total (37/36)
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.38e-01

GTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACG  >  NC_000913/4072076‑4072264
                                                                                               |                                                                                             
gTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGc                                                                                           >  1:73626/1‑100 (MQ=255)
  agagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTc                                                                                         >  2:1437010/1‑100 (MQ=255)
  agagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTc                                                                                         <  1:1115098/100‑1 (MQ=255)
   gagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCa                                                                                        >  1:1041829/1‑100 (MQ=255)
    agGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCaa                                                                                       <  1:1020891/100‑1 (MQ=255)
     ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCagt                                                                                      >  2:1400135/1‑98 (MQ=255)
     ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                      >  2:1496310/1‑100 (MQ=255)
      gTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatg                                                                                     <  2:265345/100‑1 (MQ=255)
          gcgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATCTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGc                                                                                 >  2:1622839/1‑100 (MQ=255)
                cTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCa                                                                           <  2:524351/100‑1 (MQ=255)
                    cGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGAcc                                                                       >  2:798778/1‑100 (MQ=255)
                    cGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGAcc                                                                       >  1:2299602/1‑100 (MQ=255)
                      cccAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgt                                                                     >  1:571522/1‑100 (MQ=255)
                      cccAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgt                                                                     <  1:52448/100‑1 (MQ=255)
                       ccAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgtg                                                                    <  1:1400138/100‑1 (MQ=255)
                       ccAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgtg                                                                    >  2:2026390/1‑100 (MQ=255)
                             ctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTg                                                              >  1:1222751/1‑100 (MQ=255)
                             ctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTg                                                              <  1:1314685/100‑1 (MQ=255)
                              tttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTgg                                                             <  2:1877736/100‑1 (MQ=255)
                               ttccgatcttCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                            <  1:686160/91‑1 (MQ=255)
                               ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                            <  1:1254682/100‑1 (MQ=255)
                                 ggcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcac                                                          >  1:1801523/1‑100 (MQ=255)
                                  gcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcaca                                                         >  1:922539/1‑100 (MQ=255)
                                   cggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAg                                                        >  1:2079956/1‑100 (MQ=255)
                                   cggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAg                                                        <  2:643536/100‑1 (MQ=255)
                                   cggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAg                                                        >  1:1973330/1‑100 (MQ=255)
                                    ggcATCGATTATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAAGGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAg                                                        <  2:2002531/99‑1 (MQ=255)
                                      cATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTc                                                     <  2:647016/100‑1 (MQ=255)
                                       aTCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCg                                                    <  2:1895621/100‑1 (MQ=255)
                                        tCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGcagatcggaa                                                   >  2:686158/1‑92 (MQ=255)
                                        tCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGt                                                   <  2:1973322/100‑1 (MQ=255)
                                        tCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGt                                                   >  1:1099716/1‑100 (MQ=255)
                                         cGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGtt                                                  >  1:1777718/1‑100 (MQ=255)
                                         cGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGtt                                                  <  2:1040218/100‑1 (MQ=255)
                                          gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                 <  2:922537/100‑1 (MQ=255)
                                          gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                 <  2:1801517/100‑1 (MQ=255)
                                          gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                 <  1:393812/100‑1 (MQ=255)
                                           atgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGa                                                <  2:1357450/100‑1 (MQ=255)
                                             gatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAc                                              >  1:382930/1‑100 (MQ=255)
                                              atCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAct                                             >  1:458453/1‑100 (MQ=255)
                                              atCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAct                                             >  2:1592101/1‑100 (MQ=255)
                                              atCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAct                                             <  1:765036/100‑1 (MQ=255)
                                                cGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAActct                                           >  2:696282/1‑100 (MQ=255)
                                                 gCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCtt                                          >  1:1757813/1‑100 (MQ=255)
                                                 gCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCtt                                          >  1:1174365/1‑100 (MQ=255)
                                                   ttGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTc                                        <  2:2107974/100‑1 (MQ=255)
                                                    tGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCa                                       <  2:73627/100‑1 (MQ=255)
                                                     gTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCAc                                      <  2:2299591/100‑1 (MQ=255)
                                                      ttCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCAcc                                     <  1:2012584/100‑1 (MQ=255)
                                                       tCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCa                                    <  2:1041827/100‑1 (MQ=255)
                                                       tCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCa                                    >  2:1554315/1‑100 (MQ=255)
                                                       tCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCa                                    <  2:439948/100‑1 (MQ=255)
                                                         gcgATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCata                                  >  1:714464/1‑100 (MQ=255)
                                                          cgATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCatat                                 <  2:804389/100‑1 (MQ=255)
                                                              cGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAg                             >  1:1285892/1‑100 (MQ=255)
                                                               gAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACCCTTTCACCATATCCAGa                            >  1:1181374/1‑100 (MQ=255)
                                                                aaTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGAc                           >  2:1597989/1‑100 (MQ=255)
                                                                 aTTTATCGCCACGACGCGGTGGGCGGACGAATGGCGCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACg                          <  2:458455/100‑1 (MQ=255)
                                                                       cGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctg                    <  2:2028508/100‑1 (MQ=255)
                                                                       cGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctg                    >  1:1866858/1‑100 (MQ=255)
                                                                           acgacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCa                <  2:283134/100‑1 (MQ=255)
                                                                            cgacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCAt               >  2:2221363/1‑100 (MQ=255)
                                                                             gacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATc              <  2:1400651/100‑1 (MQ=255)
                                                                             gacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATc              >  2:1682909/1‑100 (MQ=255)
                                                                               cgcgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCtt            <  2:2079948/100‑1 (MQ=255)
                                                                                gcgGTGGGCGTACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtc           >  1:14736/1‑100 (MQ=255)
                                                                                gcgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtc           <  2:571523/100‑1 (MQ=255)
                                                                                 cgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtca          <  1:1622844/100‑1 (MQ=255)
                                                                                    tGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcatca       >  1:256213/1‑100 (MQ=255)
                                                                                     gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc      <  2:61110/100‑1 (MQ=255)
                                                                                       gCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTa    >  2:740984/1‑100 (MQ=255)
                                                                                        cGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTAc   >  2:896663/1‑100 (MQ=255)
                                                                                         ggACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACg  >  1:1656073/1‑100 (MQ=255)
                                                                                               |                                                                                             
GTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACG  >  NC_000913/4072076‑4072264

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