Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,678,816 T→C E184E (GAA→GAG yphF ← putative sugar ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,8160TC98.1% 187.0 / ‑3.7 53E184E (GAA→GAGyphFputative sugar ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/24);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CGCACCGCCCGATTCCCCCATAATGGCGTTGAGATCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACTTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCAAATTCCAGCGGATCGCCGACC  >  NC_000913/2678675‑2678963
                                                                                                                                             |                                                                                                                                                   
cGCACCGCCCGATTCCCCCATAATGGCGTTGAGATCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAAt                                                                                                                                            <  2:688388/151‑1 (MQ=255)
  cACCGCCCGATTCCCCCATAATGGCGTTGAGATCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATcc                                                                                                                                          <  2:559518/151‑1 (MQ=255)
  cACCGCCCGATTCCCCCATAATGGCGTTGAGATCCGGCGTGGAGATAATTAGGTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATcc                                                                                                                                          >  2:600818/1‑151 (MQ=255)
          gATTCCCCCATAATGGCGTTGAGATCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACTTCTTCAAATCCTTTACGTc                                                                                                                                  >  1:147737/1‑151 (MQ=255)
          gATTCCCCCATAATGGCGTTGAGATCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTc                                                                                                                                  <  1:539074/151‑1 (MQ=255)
                     aaTGGCGTTGAGATCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGcacaca                                                                                                                        <  2:441751/150‑1 (MQ=255)
                     aaTGGCGTTGAGATCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGcacaca                                                                                                                        >  1:273361/1‑150 (MQ=255)
                              gagaTCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCaa                                                                                                                >  2:223355/1‑149 (MQ=255)
                                      gCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGc                                                                                                       >  2:698231/1‑150 (MQ=255)
                                       cGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCa                                                                                                      >  2:376280/1‑150 (MQ=255)
                                           gagaTAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATtg                                                                                                  >  2:513104/1‑150 (MQ=255)
                                              aTAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATtgatg                                                                                               >  1:143730/1‑150 (MQ=255)
                                                aaTTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGAc                                                                                             <  1:761666/150‑1 (MQ=255)
                                                aaTTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACg                                                                                            <  1:629043/151‑1 (MQ=255)
                                                 aTTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACgg                                                                                           >  2:190940/1‑151 (MQ=255)
                                                 aTTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACgg                                                                                           >  2:648234/1‑151 (MQ=255)
                                                  ttAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGc                                                                                          >  2:662143/1‑151 (MQ=255)
                                                    aGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCa                                                                                         <  1:169903/150‑1 (MQ=255)
                                                    aGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCa                                                                                         >  2:212715/1‑150 (MQ=255)
                                                    aGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCa                                                                                         >  2:240583/1‑150 (MQ=255)
                                                        tttCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAAtttt                                                                                    >  2:28355/1‑151 (MQ=255)
                                                          tCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTc                                                                                   <  1:204020/150‑1 (MQ=255)
                                                           cACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCgg                                                                                 >  2:627429/1‑151 (MQ=255)
                                                              cAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTg                                                                              <  1:153739/151‑1 (MQ=255)
                                                              cAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTg                                                                              >  1:96948/1‑151 (MQ=255)
                                                              cAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTg                                                                              <  2:680644/151‑1 (MQ=255)
                                                                     aaTCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAAtt                                                                       >  1:110248/1‑151 (MQ=255)
                                                                         gCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTAt                                                                   >  2:393097/1‑151 (MQ=255)
                                                                         gCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTAt                                                                   >  2:125357/1‑151 (MQ=255)
                                                                         gCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTAt                                                                   >  2:758299/1‑151 (MQ=255)
                                                                           tttATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTg                                                                 >  2:593800/1‑151 (MQ=255)
                                                                            ttATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTAtt                                                                  <  1:775842/149‑1 (MQ=255)
                                                                             tATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTgg                                                                <  2:253329/150‑1 (MQ=255)
                                                                                  aaaaCAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCaata                                                           <  1:648176/150‑1 (MQ=255)
                                                                                    aaCAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCaataaaa                                                        <  1:642126/151‑1 (MQ=255)
                                                                                     aCAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCaataaaat                                                       <  1:167890/151‑1 (MQ=255)
                                                                                      cAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCaataaaata                                                      <  1:112310/151‑1 (MQ=255)
                                                                                       aGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCaataaaata                                                      <  2:552870/150‑1 (MQ=255)
                                                                                              tCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCg                                               <  2:99367/150‑1 (MQ=255)
                                                                                              tCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGCCAATTTTATTGGCAATAAAATAATCAGCg                                               >  1:155094/1‑150 (MQ=255)
                                                                                               ccTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCgg                                              <  1:278217/150‑1 (MQ=255)
                                                                                               ccTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGc                                             <  2:204849/151‑1 (MQ=255)
                                                                                                  gATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCa                                            <  1:571539/149‑1 (MQ=255)
                                                                                                             aTTTGCGCGCCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGtt                               <  1:660092/151‑1 (MQ=255)
                                                                                                                     gcCGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTAcc                        >  1:159062/1‑150 (MQ=255)
                                                                                                                      ccGGGAACGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCa                       >  2:679026/1‑150 (MQ=255)
                                                                                                                            aCGCGGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGGTACCCAGTTTTTTACCAAATTcc                 >  1:731651/1‑150 (MQ=255)
                                                                                                                               cgGGATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCAAATTCCAGc              <  1:764884/150‑1 (MQ=255)
                                                                                                                                 ggATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCAAATTCCAGCGGa           <  1:60739/151‑1 (MQ=255)
                                                                                                                                  gATTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCAAATTCCAGCGGAt          >  2:98579/1‑151 (MQ=255)
                                                                                                                                   aTTTTAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCAAATTCCAGCGGATc         >  2:719365/1‑151 (MQ=255)
                                                                                                                                      ttAAAACCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCAAATTCCAGCGGATCGcc      >  2:371250/1‑151 (MQ=255)
                                                                                                                                          aaCCTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCAAATTCCAGCGGATCGCCGAcc  >  2:258241/1‑151 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CGCACCGCCCGATTCCCCCATAATGGCGTTGAGATCCGGCGTGGAGATAATTAGTTTTTCACCAACGGAAATCGCTTTATCTAAAACAGTCCCTTCCTGATTAGCGACAATTTGCGCGCCGGGAACGCGGGATTTTAAAACTTCTTCAAATCCTTTACGTCGCTGCACACAAACTTCAAAGGCTTCGCAATTGATGACGGCAATTTTCGGCTGGTCAATTTTATTGGCAATAAAATAATCAGCGGCAGCGTTACCCAGTTTTTTACCAAATTCCAGCGGATCGCCGACC  >  NC_000913/2678675‑2678963

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