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,1710CA96.2% 274.7 / ‑1.1 78V161F (GTC→TTC) yihT6‑deoxy‑6‑sulphofructose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (42/33);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGCC  >  NC_000913/4072075‑4072266
                                                                                                |                                                                                               
tGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATgg                                                                                              <  2:872926/100‑1 (MQ=255)
tGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATgg                                                                                              <  2:1546843/100‑1 (MQ=255)
 gTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGc                                                                                             <  1:132301/100‑1 (MQ=255)
  tAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCt                                                                                            <  2:2224505/100‑1 (MQ=255)
   agagGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTc                                                                                           >  2:406845/1‑100 (MQ=255)
      ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                        <  1:262848/100‑1 (MQ=255)
      ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAat                                                                                        >  2:40131/1‑100 (MQ=255)
      ggTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCAGACGAATGGCTCAat                                                                                        >  1:2364607/1‑100 (MQ=255)
         cAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATTGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAatgat                                                                                     <  1:2344764/100‑1 (MQ=255)
            cgcCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCt                                                                                  <  2:2183076/100‑1 (MQ=255)
              ccACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTa                                                                                <  1:2028074/100‑1 (MQ=255)
               cACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAc                                                                               <  1:1543250/100‑1 (MQ=255)
                aCTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAcc                                                                              >  2:107652/1‑100 (MQ=255)
                aCTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTAcc                                                                              <  1:1206167/100‑1 (MQ=255)
                  tGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGCCGGACGAATGGCTCAATGATGCTTACCAg                                                                            >  2:829186/1‑100 (MQ=255)
                   gTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGa                                                                           >  1:1739604/1‑100 (MQ=255)
                      gCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCg                                                                        >  1:1884339/1‑100 (MQ=255)
                      gCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCg                                                                        <  2:2252662/100‑1 (MQ=255)
                       cccAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCgt                                                                       >  1:2237731/1‑100 (MQ=255)
                             tctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGt                                                                 <  2:1514859/100‑1 (MQ=255)
                              ctTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTg                                                                <  2:1178560/100‑1 (MQ=255)
                               gatctcggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCGATGATGCTTACCAGACCGTGTGAGTgg                                                               <  1:182819/95‑1 (MQ=255)
                                ttGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGc                                                              <  2:2187592/100‑1 (MQ=255)
                                ttGGCGGCATCGATGATCGCTTGTCCGCGATCGAATTTATCGCCACGACGCGGTGGGTGGACGACTGGCTCAATGAGGCTTACCAGTCCGTGTGAGTg                                                                <  2:2065906/98‑1 (MQ=255)
                                    cggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCGATGATGCTTACCAGACCGTGTGAGTggagatc                                                          >  2:182818/1‑95 (MQ=255)
                                     ggcATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGt                                                         <  2:861781/100‑1 (MQ=255)
                                       cATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTc                                                       >  1:547315/1‑100 (MQ=255)
                                        aTCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGAACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCg                                                      >  1:2300808/1‑100 (MQ=255)
                                          cGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGtt                                                    >  2:2213305/1‑100 (MQ=255)
                                          cGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGtt                                                    <  2:977978/100‑1 (MQ=255)
                                          cGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGtt                                                    <  1:10143/100‑1 (MQ=255)
                                          cGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGtt                                                    >  1:279275/1‑100 (MQ=255)
                                           gatgatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTg                                                   <  1:618411/100‑1 (MQ=255)
                                              gatCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAc                                                >  2:1597133/1‑100 (MQ=255)
                                               atCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAAct                                               <  1:2233740/100‑1 (MQ=255)
                                                   cTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCttt                                           >  1:1007070/1‑100 (MQ=255)
                                                        tCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCa                                      >  2:517031/1‑100 (MQ=255)
                                                           cgATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCatat                                   >  1:1084550/1‑100 (MQ=255)
                                                            naTCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTAACCAGACCGTGTGAGAGGCACAGTACGTTGAACTCTTTCACCATATc                                  >  2:1618871/2‑100 (MQ=255)
                                                            gATCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATc                                  <  2:1022956/100‑1 (MQ=255)
                                                              tCGAATTTATCGCCACGACGCGGTGGGCGGACGTATGGCTTAATGATGCTTACCAGGCCGTGTGAGTGGCACAGTTCGTTGAACTCTCTCAACATATcc                                 <  2:1265274/99‑1 (MQ=255)
                                                              tCGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCa                                <  1:574048/100‑1 (MQ=255)
                                                               cGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAg                               >  1:1702006/1‑100 (MQ=255)
                                                               cGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAg                               >  2:1999616/1‑100 (MQ=255)
                                                               cGAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAg                               >  2:1565677/1‑100 (MQ=255)
                                                                gAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGa                              >  1:2086089/1‑100 (MQ=255)
                                                                gAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGa                              <  2:1667058/100‑1 (MQ=255)
                                                                gAATTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGa                              >  2:2086606/1‑100 (MQ=255)
                                                                  aTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACg                            >  2:1566550/1‑100 (MQ=255)
                                                                  aTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACg                            >  2:1949913/1‑100 (MQ=255)
                                                                  aTTTATCGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACg                            <  1:829185/100‑1 (MQ=255)
                                                                        cGCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctg                      >  2:896912/1‑100 (MQ=255)
                                                                         gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctgc                     >  1:1651877/1‑100 (MQ=255)
                                                                         gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctgc                     >  1:1502747/1‑100 (MQ=255)
                                                                         gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctgc                     >  1:1967675/1‑100 (MQ=255)
                                                                         gCCACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTtgctgc                     >  2:1049480/1‑100 (MQ=255)
                                                                          ccACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcg                    >  2:485329/1‑100 (MQ=255)
                                                                          ccACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcg                    >  2:75468/1‑100 (MQ=255)
                                                                          ccACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcg                    >  1:1353585/1‑100 (MQ=255)
                                                                          ccACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcg                    <  1:278938/100‑1 (MQ=255)
                                                                           cACGACGCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTgcgc                   >  1:1119828/1‑100 (MQ=255)
                                                                            acgacgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCa                  <  2:547316/100‑1 (MQ=255)
                                                                             cgacgCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCAt                 >  2:2403251/1‑100 (MQ=255)
                                                                               acgCGGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCt               <  1:40131/100‑1 (MQ=255)
                                                                                  cgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtca            <  2:704075/100‑1 (MQ=255)
                                                                                  cgGTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtca            <  1:2339326/100‑1 (MQ=255)
                                                                                  cgGTGGGCGAACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtca            <  2:2300806/100‑1 (MQ=255)
                                                                                   ggTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcat           >  2:1033029/1‑100 (MQ=255)
                                                                                   ggTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTAAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcat           >  2:2005164/1‑100 (MQ=255)
                                                                                    gTGGGCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcatc          <  2:2249107/100‑1 (MQ=255)
                                                                                      gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTtcatca         <  2:1967674/99‑1 (MQ=255)
                                                                                      gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc        >  1:361114/1‑100 (MQ=255)
                                                                                      gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc        >  1:149577/1‑100 (MQ=255)
                                                                                      gggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCAc        >  2:2099531/1‑100 (MQ=255)
                                                                                       ggCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACt       >  2:715581/1‑100 (MQ=255)
                                                                                        gCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCCTCTTCATCACTa      >  1:696436/1‑100 (MQ=255)
                                                                                        gCGGACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTa      >  1:1759752/1‑100 (MQ=255)
                                                                                          ggACGAATGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACg    >  1:1439112/1‑100 (MQ=255)
                                                                                            acgaaTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGcc  <  1:1565676/100‑1 (MQ=255)
                                                                                                |                                                                                               
TGTAGAGGTCAGCGCCACTGTCGCCCAGCTCTTTGGCGGCATCGATGATCGCTTGTTCGCGATCGAATTTATCGCCACGACGCGGTGGGCGGACGACTGGCTCAATGATGCTTACCAGACCGTGTGAGTGGCACAGTTCGTTGAACTCTTTCACCATATCCAGACGTTGCTGCGCATCTTCATCACTACGCC  >  NC_000913/4072075‑4072266

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