Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,057,759 +C coding (1773/1776 nt) yihK → putative GTP‑binding factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,057,7581.C98.6% 367.2 / ‑3.3 141H591P (CAC→CCC) yihKputative GTP‑binding factor
Reads supporting (aligned to +/- strand):  ref base . (0/2);  major base C (48/89);  minor base T (0/2);  total (48/93)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 147 total aligned reads displayed.

AACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAA  >  NC_000913/4057683‑4057828
                                                                              |                                                                      
aaCAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCa                                                                         <  1:3985124/76‑1 (MQ=255)
   aaGCTCTGGAGTTCATTGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcct                                                                     >  2:5017298/1‑74 (MQ=255)
   aaGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcct                                                                     >  1:5819057/1‑74 (MQ=255)
   aaGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcct                                                                     <  2:2127725/76‑3 (MQ=255)
       tctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACg                                                                  >  1:4377241/1‑76 (MQ=255)
       tctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACg                                                                  <  2:3909724/76‑1 (MQ=255)
       tctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACg                                                                  <  2:1547316/76‑1 (MQ=255)
       tctGGAGTTCATCGATGATGACGAACTGGGAAAAGTGCCTCCGACCTCTATCCGTAGTCGTAAACCTCATCTGACg                                                                  <  1:4060714/76‑1 (MQ=255)
        ctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACgg                                                                 <  2:3515128/76‑1 (MQ=255)
         tGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGa                                                                <  1:2349041/76‑1 (MQ=255)
          ggAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaa                                                               <  2:2437943/76‑1 (MQ=255)
          ggAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGCCACCTGACGGaa                                                               >  2:752179/1‑76 (MQ=255)
           gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                              <  2:5373952/76‑1 (MQ=255)
           gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                              >  1:3893628/1‑76 (MQ=255)
           gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                              <  1:3980503/76‑1 (MQ=255)
           gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                              <  1:2766573/76‑1 (MQ=255)
             gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAAc                                                            <  1:356487/76‑1 (MQ=255)
              ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACg                                                           <  1:4867609/76‑1 (MQ=255)
               tCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGa                                                          <  2:1932684/76‑1 (MQ=255)
                cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAt                                                         <  2:2560373/76‑1 (MQ=255)
                 aTCGATGATTACGAACTGGTAGAAGTTCCTCCGCCCTCTATCCGTCTTCGTAAGCGTCCCCAGACGGAATACGacc                                                        <  2:4280070/76‑3 (MQ=255)
                 aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtc                                                        <  2:4955971/76‑1 (MQ=255)
                   cGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgt                                                      <  1:5595328/76‑1 (MQ=255)
                    gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                     <  1:3775154/76‑1 (MQ=255)
                    gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                     >  2:2625990/1‑76 (MQ=255)
                    gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                     >  1:5054885/1‑76 (MQ=255)
                     atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtcg                                                    <  1:5149320/76‑1 (MQ=255)
                     atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtcg                                                    >  1:1928447/1‑76 (MQ=255)
                      tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgc                                                   >  1:396849/1‑76 (MQ=255)
                      tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgc                                                   <  2:1257617/76‑1 (MQ=255)
                       gatgaCGCACTGGTAGAAGTGACTCCGACCTCTATCCGTCTTCGTATACGTCCCCTGACGGAAAACTATCGTggcc                                                  <  1:4924098/76‑1 (MQ=255)
                       gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                  <  2:4997827/76‑1 (MQ=255)
                       gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                  <  1:4072512/76‑1 (MQ=255)
                       gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                  <  2:2326104/76‑1 (MQ=255)
                        atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccg                                                 <  1:1654017/76‑1 (MQ=255)
                         tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGCTCGTcgccgc                                                >  2:5692018/1‑76 (MQ=255)
                         tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccgc                                                <  1:4829000/76‑1 (MQ=255)
                         tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccgc                                                <  2:5603445/76‑1 (MQ=255)
                         tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccgc                                                <  1:4359404/76‑1 (MQ=255)
                          gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCcgcg                                               >  1:2500897/1‑76 (MQ=255)
                            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGcc                                             >  1:5382849/1‑76 (MQ=255)
                            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGcc                                             <  1:5280463/76‑1 (MQ=255)
                             gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCa                                            <  2:2861281/76‑1 (MQ=255)
                              aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCaa                                           >  2:3796905/1‑76 (MQ=255)
                                cTGGTAGAATGGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAcc                                         <  1:5052594/76‑1 (MQ=255)
                                cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAcc                                         <  1:1017254/76‑1 (MQ=255)
                                 tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcg                                        >  1:3232475/1‑76 (MQ=255)
                                  ggTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgc                                       <  2:3816436/76‑1 (MQ=255)
                                  ggTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgc                                       >  1:5744445/1‑76 (MQ=255)
                                  ggTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgc                                       <  1:5812524/76‑1 (MQ=255)
                                   gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgcg                                      <  1:1715964/76‑1 (MQ=255)
                                    tAGAAGTGACTCCGACCTCTAGCCGTATTCGTAAACGGCACCTGACGGAAAAAGATCGTCGCCGCGCCCACCgcgc                                     >  2:559463/1‑76 (MQ=255)
                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCTAACCGCGCAcc                                  <  1:1559795/76‑1 (MQ=255)
                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAcc                                  <  1:5176619/76‑1 (MQ=255)
                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAcc                                  <  2:2572317/76‑1 (MQ=255)
                                        aGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCg                                 >  1:657010/1‑76 (MQ=255)
                                         gTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGa                                <  2:3535120/76‑1 (MQ=255)
                                         aagACTCCGACCTCTATCCGTATTCGCAAACGTCACCTGACGGAAAACGATCGTCGCAGTGCCAACCGCGCACCGa                                >  1:1226260/3‑76 (MQ=255)
                                          tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaa                               <  1:1929404/76‑1 (MQ=255)
                                          ggACTCCTACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaa                               <  1:277505/75‑1 (MQ=255)
                                          gcaCTCGGACCTCTATCCGTATTCGTAAAGGTCACCTGTCGGAAAGCGGTCGTCGCCGCGCCAACCGCGCGCCGaa                               <  2:4418862/74‑1 (MQ=255)
                                           gACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaaa                              >  2:1192959/1‑76 (MQ=255)
                                           gACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaaa                              <  2:2582963/76‑1 (MQ=255)
                                            aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAg                             <  2:4664253/76‑1 (MQ=255)
                                            aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAg                             <  1:5335412/76‑1 (MQ=255)
                                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGa                            <  2:1744735/76‑1 (MQ=255)
                                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGa                            <  1:2996363/76‑1 (MQ=255)
                                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGa                            <  1:4446731/76‑1 (MQ=255)
                                               ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                          <  1:2427373/76‑1 (MQ=255)
                                               ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                          >  1:4630325/1‑76 (MQ=255)
                                               ccGACCTCTAGCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                          <  2:5143727/76‑1 (MQ=255)
                                                cgcccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                         <  1:2878188/73‑1 (MQ=255)
                                                cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCTCCGCGCCAACCGCGCACCGAAAGACGa                         <  2:5462193/76‑1 (MQ=255)
                                                 gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAt                        <  1:1513943/76‑1 (MQ=255)
                                                 gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAt                        >  1:4413266/1‑76 (MQ=255)
                                                  aCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAtt                       <  1:1076084/76‑1 (MQ=255)
                                                  aCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAtt                       <  1:3623000/76‑1 (MQ=255)
                                                   ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTa                      <  2:3273862/76‑1 (MQ=255)
                                                   ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAGCGTCGCCGCGCCAACCCTGCAACGAAAGGCGATTa                      >  1:993510/1‑76 (MQ=255)
                                                    ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTaa                     >  1:5327558/1‑76 (MQ=255)
                                                     tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAt                    <  1:839428/76‑1 (MQ=255)
                                                      ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAtt                   <  2:4241940/76‑1 (MQ=255)
                                                       tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAttt                  <  2:1873938/76‑1 (MQ=255)
                                                          ccGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGt               >  2:2970497/1‑76 (MQ=255)
                                                          ccGTATTCGGAAACGGCGCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGt               <  2:5200985/76‑1 (MQ=255)
                                                            gTATTCGTAAACGTCACCTGACGGCAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttt             >  2:4219688/1‑76 (MQ=255)
                                                            gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttt             >  1:2791583/1‑76 (MQ=255)
                                                            gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttt             >  2:2771129/1‑76 (MQ=255)
                                                             tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCAGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttt            <  1:2326769/76‑1 (MQ=255)
                                                              aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttttt           >  2:2569892/1‑76 (MQ=255)
                                                              aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttttt           >  1:207399/1‑76 (MQ=255)
                                                               ttCGGAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttttt          <  2:3878642/76‑1 (MQ=255)
                                                                tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTa         >  1:2055014/1‑76 (MQ=255)
                                                                  gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACTGAAAGACGATTAATTTCGTTTTTTAGt       <  1:380208/76‑1 (MQ=255)
                                                                  gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGt       <  2:494498/76‑1 (MQ=255)
                                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgc                                       <  1:4525087/43‑1 (MQ=37)
                                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgc                                       >  2:4525086/1‑43 (MQ=37)
                                                                      aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATa   <  1:5604165/76‑1 (MQ=255)
                                                                      aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATa   >  2:4773873/1‑76 (MQ=255)
                                                                       cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaa  >  1:3181501/1‑76 (MQ=255)
                                                                              |                                                                      
AACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAA  >  NC_000913/4057683‑4057828

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