Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_pgi 2,263,081 G→A 8.6% L139L (CTG→TTG)  fruB ← fused fructose‑specific PTS enzymes: IIA component/HPr component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi2,263,0810GA 8.6% 125.4 / 4.0 50L139L (CTG→TTG) fruBfused fructose‑specific PTS enzymes: IIA component/HPr component
Reads supporting (aligned to +/- strand):  ref base G (17/26);  new base A (2/2);  total (21/29)
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.47e-01

TTGGTGACGAAAGTGGCGTCAACTGCCCCCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCG  >  NC_000913_3_pae_pgi/2262933‑2263229
                                                                                                                                                     |                                                                                                                                                     
ttGGTGACGAAAGTGGCTTCAACTGCCCCCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAg                                                                                                                                                       <  2:391231/149‑1 (MQ=255)
  ggTGGCGAAACGCGCGTCAACTGCCCGCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCTGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTa                                                                                                                                                     <  2:849239/149‑1 (MQ=255)
     gccGAACGTGGCCTCGACGCCCCCCGCGTCGCTCAGACGCGCAGTTGTCATCCCCTGAAGTGTCAGCACATCGCTGGCGACGATATTCAGGGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATg                                                                                                                                                 <  2:70303/147‑1 (MQ=255)
             tGGCGTCAACTGCCCCCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTACTCACTGTGCTTTTCGCCCATCATTAATGCCCGGAGt                                                                                                                                         >  2:508994/1‑149 (MQ=255)
               gCGTCAACTGCCCCCGCCTCTTTCAGACCCACAGCGTTCGCGATTTCTTTTACGAGCACATCGCTGGTGATGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGttc                                                                                                                                       <  1:731897/149‑1 (MQ=255)
                cGTCAACTGCCCCCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGttct                                                                                                                                      <  1:90350/149‑1 (MQ=255)
                 gTCAACTGCCCCCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGttctt                                                                                                                                     <  2:604838/149‑1 (MQ=255)
                  cccagtGCTCCCGCCTTTTTCAGACGCGCAGCGTTCAGTTTTTGATGTGTTAGCATATCGATGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGttcttc                                                                                                                                    <  1:31439/144‑1 (MQ=255)
                    aaCTGCCCCCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGATGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTctg                                                                                                                                  <  1:559845/149‑1 (MQ=255)
                       tGCCCCCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCACCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACCCTGCTTTTCGCCCATCAGTAACGCGCGCAGTTCTTctgctc                                                                                                                               >  2:563287/1‑148 (MQ=255)
                                ctctTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCCGATCGCTGGCGCCGATATCCAGTGTCAGCATTTCGTTGGCGAGCTTCAGCTGCTCACTCGGCTTTTCGCCCATCAGTAGTGCGCGAAGTTCTACTACTGGGGTTGCTg                                                                                                                      >  1:712440/1‑149 (MQ=255)
                                 tctTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGGTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGa                                                                                                                     <  2:509634/149‑1 (MQ=255)
                                 tctTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGa                                                                                                                     >  1:886383/1‑149 (MQ=255)
                                 tctTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGCCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGGTGCTGa                                                                                                                     >  1:875438/1‑149 (MQ=255)
                                  ctgtCAGACGCGCAGCCTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCCCACTCTGCTTCTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGAc                                                                                                                    <  2:420037/146‑1 (MQ=255)
                                   tttCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGGTCCCTGGCGACGATAAAGTTTTCCATCATTTCGTTGTCGAGCTTCAGCTTCTCACGCTGCGTTTCGCCCATCACGACTGCGCGAAGTTCTTCTGGTGTTGGTGCGGACt                                                                                                                   >  1:396376/1‑149 (MQ=255)
                                          cgcgCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAgttg                                                                                                            <  1:40890/149‑1 (MQ=255)
                                            cgcAGGGTTCGCCTCCTGAAGAGTCAGCATATGGCTGGCGACGAGATCCAGTGGCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTATGATTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAgttgtt                                                                                                          <  2:165284/149‑1 (MQ=255)
                                              cAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCa                                                                                                        <  1:650402/149‑1 (MQ=255)
                                              cAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTGTCGCCCATCAGTAATGCGCGATGTTCTTCTGCTGTTGTTGCTGCCTCCAGCTGTTCa                                                                                                        >  1:99889/1‑149 (MQ=255)
                                                               agagTCAGCAGATCGCTGGTGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGTTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCt                                                                                       <  1:372543/149‑1 (MQ=255)
                                                               agagCCAGCAGCTTGCTGGGCGCGTGATCGATTGTCCGCCTGTCGTTCTCGAGTTTGAACTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGcc                                                                                       <  2:68383/149‑2 (MQ=255)
                                                                 gtggcggcagTTCGCTGGCGACGATATCCAGTTTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATCCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTc                                                                                      <  2:118846/143‑1 (MQ=255)
                                                                   tCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAAGGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGt                                                                                   >  1:896543/1‑149 (MQ=255)
                                                                          gATCGCTGGCGAGCATAGCCCGTGTCCGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTgg                                                                             <  2:598984/149‑1 (MQ=255)
                                                                                 tgctactatatCCAATTCCAGCATTTTGTTGTCGAGCTTCAGCTTATCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGCCTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCt                                                                      <  2:301788/142‑1 (MQ=255)
                                                                                     acgaTGTCCAGTTGCAGCATTTCGTGGTCGGTCTTCAGCTGCTCTCTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCg                                                                  <  2:138383/149‑1 (MQ=255)
                                                                                          atCCAGGTTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCa                                                             <  1:602613/149‑1 (MQ=255)
                                                                                            ccAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGttctt                                                                                                                                     <  1:797987/75‑1 (MQ=255)
                                                                                            ccAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGttctt                                                                                                                                     >  2:797987/1‑75 (MQ=255)
                                                                                             cAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCGCAGTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGc                                                          <  2:937974/149‑1 (MQ=255)
                                                                                                  tCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAgg                                                     >  2:203758/1‑148 (MQ=255)
                                                                                                  tCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAgc                                                     >  1:667743/1‑148 (MQ=255)
                                                                                                  ccAGCATTTCGCTGTCGAGCTTCATCTTCTCACTGTGCTTTTTGCGCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTTCTGACTGCACTTGGTCAGCAACGGAATCATCGCTCACTACGTGGGTCAGCTGGAGTAGCAGGCCCAGa                                                     <  2:71282/148‑1 (MQ=255)
                                                                                                     gCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGt                                                  >  1:308356/1‑149 (MQ=255)
                                                                                                                       tCAGCTGCTCACTCTGCTTTTCGCCCATCAGGAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGTCTGCAGTTGTTCAGCAACGGAATCCTCGCTCAGTACGTGGGTCAGCTGGCGTTGCAGGCCCATATGTTCCTCCGAGCTGgcagca                                >  1:266004/1‑149 (MQ=255)
                                                                                                                          gctgctCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATa                             >  2:600759/1‑149 (MQ=255)
                                                                                                                          gctgctCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATa                             >  1:472101/1‑149 (MQ=255)
                                                                                                                              ctCCCTCTGCTTTTCGCCCATCCGTACTGCGCGAAATTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCACCGGAATCATCGCTGAATACGTGGGCCAGCGTGCGTAGCCGGCCCAGATGTTCATCCGAGCTGCCACCAATACCgc                         >  2:116775/1‑148 (MQ=255)
                                                                                                                              ctCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAgtt                                                                                                             >  1:306729/1‑65 (MQ=255)
                                                                                                                              ctCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAgtt                                                                                                             <  2:306729/65‑1 (MQ=255)
                                                                                                                                taCTCTGCCTTTCACCCATCACTAACGCGCGAACTTCTCTTGCTGTTGTTGCTGACTTCAGTTTTTCAGCAACGGAATCATCGCGCAGTACGGGGGTTAGCTGGCGTAGCAGGCCCAGATTTTCATCCGAGCGGGCAGCAATACCAAtt                       <  2:643231/148‑1 (MQ=255)
                                                                                                                                cACTCTGCTTTTCGCCCATCAATAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTgg                                                                    <  1:413022/104‑1 (MQ=255)
                                                                                                                                cACTCTGCTTTTCGCCCATCAATAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTgg                                                                    <  1:412997/104‑1 (MQ=255)
                                                                                                                                cACTCTGCTTTTCGCCCATCAATAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTgg                                                                    >  2:413022/1‑104 (MQ=255)
                                                                                                                                cACTCTGCTTTTCGCCCATCAATAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTgg                                                                    >  2:412997/1‑104 (MQ=255)
                                                                                                                                  ctctGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGc                     <  1:361925/149‑1 (MQ=255)
                                                                                                                                         tttAGCCCATAAGTATATGGCGATGTTCTGCTGGTGTGTTGGCGCATTTCAGTTGTTCAGCAAGCGAAGCCTCGCTCAGTTCGTGGGGTAGCTGGCGTCGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAACACTGACTGCCACGTAc              <  2:270739/149‑1 (MQ=255)
                                                                                                                                            cGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCt           <  2:886383/149‑1 (MQ=255)
                                                                                                                                                 aTCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTg      >  1:770677/1‑149 (MQ=255)
                                                                                                                                                 aTCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACCCTACCTg      >  2:918295/1‑149 (MQ=255)
                                                                                                                                                     gTAATGCGCGAAGCTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAAGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACATGACCg  >  1:270358/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGGTGACGAAAGTGGCGTCAACTGCCCCCGCCTCTTTCAGACGCGCAGCGTTCAGCGCCTGAAGAGTCAGCAGATCGCTGGCGACGATATCCAGTGTCAGCATTTCGTTGTCGAGCTTCAGCTGCTCACTCTGCTTTTCGCCCATCAGTAATGCGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCG  >  NC_000913_3_pae_pgi/2262933‑2263229

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