Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,177,171 T→A 9.3% intergenic (+43/‑187) tufB → / → secE translation elongation factor Tu 2/Sec translocon subunit SecE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,177,1710TA 9.3% 135.8 / 8.1 54intergenic (+43/‑187)tufB/secEtranslation elongation factor Tu 2/Sec translocon subunit SecE
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base A (2/3);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-02

ATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCG  >  NC_000913/4177033‑4177311
                                                                                                                                          |                                                                                                                                             
atgcccccGAGCGCGATGGACGACGGTCTGCGGTTCGCAATCCGTGAAGGCGGCCGTACCGTAGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGa                                                                                                                                           <  1:411174/138‑1 (MQ=255)
aTCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCACCATTTGa                                                                                                                                           >  2:68667/1‑143 (MQ=255)
    aCCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGcc                                                                                                                                       >  2:741115/1‑143 (MQ=255)
      ccGAGCGCGAGGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCt                                                                                                                                     <  1:725034/143‑1 (MQ=255)
         aTCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCtttt                                                                                                                                  <  2:560837/143‑1 (MQ=255)
          tCGCGATGGACGACGGTCTGCGCTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCttttt                                                                                                                                 <  2:219844/143‑1 (MQ=255)
          tCGCGAGGTACGACGGTCTGCGTTTCGCAATCCGTGAAGCCGGCCGCACCGTTGGCGCAAGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttta                                                                                                                                 <  1:302998/143‑2 (MQ=25)
          gcgcgACGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCAAATGATGCCCtttta                                                                                                                                 <  2:511853/142‑2 (MQ=255)
             gggtGGACGGCGGCGTGCGTGCGGCACTTCGGGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAGGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                                                              <  2:730201/140‑1 (MQ=37)
             cgATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCa                                                                                                                              >  2:813105/1‑143 (MQ=255)
                gttccggcgGTCTGCGTTTCGCAAGCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGc                                                                                                                           <  2:741973/136‑1 (MQ=255)
                     gacgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAACGGGCATCATTGGATGCCCTTTTTGCACGCTTTCg                                                                                                                      >  1:241401/1‑143 (MQ=255)
                     gacgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGACGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTATTTGCACGCTTTCg                                                                                                                      >  2:658326/1‑143 (MQ=255)
                      acgGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTCGCGCAGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGt                                                                                                                     <  1:112025/143‑1 (MQ=255)
                        gggcTGCGTTTCGTAGTCCTTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAc                                                                                                                   <  2:605973/140‑1 (MQ=255)
                         gTCTGCGTTTCGTAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGGTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                                                                  <  1:469001/143‑1 (MQ=255)
                         gTCTGCGTTTCGCAATCCGTGAGGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTAcc                                                                                                                  <  2:491803/143‑1 (MQ=255)
                          tCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAGCAGTGGACGCAATGCGCACGAAAAGGGCATCGTTCGATGCCATTTTTGCACGCAATCGtgcgg                                                                                                                 >  1:214806/1‑139 (MQ=21)
                             gCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCGCGCTTTCGTACCAGaa                                                                                                             >  2:342130/1‑143 (MQ=255)
                             gCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGaa                                                                                                             >  1:554261/1‑143 (MQ=255)
                                tttCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCCCGCTTTCGTACCAGAACCt                                                                                                          >  1:309444/1‑143 (MQ=255)
                                 ttCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTTCCGATAACATTTGCAGGCAGACGCACTAAAAGGGCATCATTCGATACCCTTTTTGCACGCTTTCGTACCAGAACCTg                                                                                                         >  1:375639/1‑143 (MQ=25)
                                  ccGCCATCCGTGAAGGCGTCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTCACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCCTACCAGACCCTgg                                                                                                        <  2:684204/142‑1 (MQ=38)
                                    gCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCCGAACCTGGCt                                                                                                      >  2:139126/1‑143 (MQ=255)
                                          cGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCtcatccg                                                                                                >  2:471766/1‑141 (MQ=255)
                                             ggtggcggcCGTGCCGTAGGCGCGGGCGTTGTAGCAAAAGTTCTTAGCTAATTGCCGATAGCATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCTCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGa                                                                                             <  2:143576/140‑1 (MQ=255)
                                                      ggTACCGTGGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCGATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTg                                                                                    <  1:474044/142‑1 (MQ=255)
                                                                   cgGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGc                                                                       <  1:123304/143‑1 (MQ=255)
                                                                      gccggggaGCAAAAGTTCGGAGCTAAGTGCCGATAACATTTGACGCAATGCGCACTAAAAGGCCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTga                                                                    <  1:96088/136‑1 (MQ=255)
                                                                         ttGTAGCAAAAGTTCTTAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATTTATTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGAc                                                                 <  1:512196/143‑1 (MQ=255)
                                                                         ttGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACATGGCTCATCAGTGCTTTTCTTTGTCATAATCATTGCTGAGAc                                                                 >  2:649782/1‑143 (MQ=255)
                                                                           gTAGCCAAAGTTCTGAGGTCAGCGCCGATAACATTTGACGCACTGGACAACCAAGGGGCAACATTTGCTGCTCTTTTTGCACGCTTTCGTACCAGACAATAGGCAATTCAACATTTTCTTTGTCATACTCATTGCTGAGACAg                                                               <  2:540196/143‑1 (MQ=255)
                                                                            tAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGACTCATCAGTGATATTCTTTGTCATAATCACTGCTGGGACaag                                                              >  2:477889/1‑141 (MQ=255)
                                                                               cAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCCTCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGCTCTTCTTTCTCATACTCATTGCTGAGACaaggtc                                                           >  2:608511/1‑138 (MQ=255)
                                                                                     ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGa                                                     <  2:495992/143‑1 (MQ=255)
                                                                                     ttCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCGCTGTTGa                                                     >  2:51178/1‑143 (MQ=255)
                                                                                        tGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAggg                                                  <  1:813105/143‑1 (MQ=255)
                                                                                           gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt                                               <  1:667239/143‑1 (MQ=255)
                                                                                           gCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGt                                               >  1:690307/1‑143 (MQ=255)
                                                                                                  gacGATAACATTCGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATcc                                        <  1:642693/141‑1 (MQ=255)
                                                                                                   ccGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGACACAGGCTCTGTTGAGGGCGTATAATCCg                                       >  2:30000/1‑143 (MQ=255)
                                                                                                           aTTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTa                               >  1:5341/1‑143 (MQ=255)
                                                                                                            cttGACGCGAGGCGCACTAAAGGGGCATCATTTGATGCCCTTTTGGTACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTaa                              <  2:639638/142‑1 (MQ=255)
                                                                                                                acacAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTCGCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAGGCTaagccg                          >  2:13882/4‑139 (MQ=255)
                                                                                                                     aTGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGttt                     <  1:214641/143‑1 (MQ=255)
                                                                                                                      tgagcACTCTAAGGGCATCATTTGATGCCCGTTTTGCTCGCTTTCGTACCAGAACCTGGCTCATCAGGGATTTTCTTTGTCATATTCATTGCTGAGACAGGCTCTGTTGAGGGCGTATACTCCGAAAAGCTAATACGCGTTTc                    <  2:717028/140‑1 (MQ=255)
                                                                                                                          cACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTATTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGAtt                >  2:621132/1‑143 (MQ=255)
                                                                                                                             tAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCATAACCTGGCTCATCAGTGATTTTCTTGGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGAGTATAATCCGAAACGCTAATACGCGTTTCGATGTgg             >  1:174034/1‑143 (MQ=255)
                                                                                                                             tAAAAGGGCATCATTTGATGCCCTTTTTGCAAGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCGGAGACAGGCTCTGTCGAGGGCGTATACTCCGAAAAGGTAATGCGAGATTCAATGTgg             >  2:490625/1‑143 (MQ=255)
                                                                                                                             tAAAAGGGCATCATTTGATGCCCTTTGTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTGCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTgg             <  1:132040/143‑1 (MQ=255)
                                                                                                                              aaaaGGGCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATGTGGt            >  2:507566/1‑143 (MQ=255)
                                                                                                                                 aGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTg         >  1:33485/1‑143 (MQ=255)
                                                                                                                                  gggCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGc        <  2:645226/143‑1 (MQ=255)
                                                                                                                                  gggCATCAAATGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATAAGCGTTTCGATTTGGTTTgg        >  1:790607/1‑142 (MQ=255)
                                                                                                                                     catcatTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTc     >  1:745239/1‑143 (MQ=255)
                                                                                                                                      atcatTTGATGCCCTTGTTGCACGCTTTCGTACCAGACCCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTcg    <  1:173981/143‑1 (MQ=255)
                                                                                                                                       tcagTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGACTcgc   <  2:544469/143‑1 (MQ=255)
                                                                                                                                        catTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGACGGCGTATAATCCGCAAAGCTAATACGCCTTACGATTCGGTTTGCCtgggg  >  2:641596/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                             
ATCCACCCGATCGCGATGGACGACGGTCTGCGTTTCGCAATCCGTGAAGGCGGCCGTACCGTTGGCGCGGGCGTTGTAGCAAAAGTTCTGAGCTAATTGCCGATAACATTTGACGCAATGCGCACTAAAAGGGCATCATTTGATGCCCTTTTTGCACGCTTTCGTACCAGAACCTGGCTCATCAGTGATTTTCTTTGTCATAATCATTGCTGAGACAGGCTCTGTTGAGGGCGTATAATCCGAAAAGCTAATACGCGTTTCGATTTGGTTTGCCTCGCG  >  NC_000913/4177033‑4177311

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: