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.C99.1% 282.8 / ‑4.2 110coding (1772/1776 nt)yihKputative GTP‑binding factor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (32/77);  total (32/78)
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

Only 100 of 116 total aligned reads displayed.

AACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAA  >  NC_000913/4057683‑4057828
                                                                            |                                                                      
aaCAAGCTTTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCa                                                                         <  1:2637811/76‑1 (MQ=255)
  cAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcc                                                                      >  1:1410010/1‑75 (MQ=255)
  cAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcc                                                                      >  1:4863282/1‑75 (MQ=255)
    aGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcctg                                                                    <  2:2528606/76‑4 (MQ=255)
    aGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcctg                                                                    >  1:1804445/1‑73 (MQ=255)
    aGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcctg                                                                    <  2:4633916/76‑4 (MQ=255)
       tctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACg                                                                  >  2:3060203/1‑76 (MQ=255)
       tctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACg                                                                  >  2:1136731/1‑76 (MQ=255)
         tGGGGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGa                                                                <  2:557234/76‑1 (MQ=255)
         tGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGa                                                                >  2:42609/1‑76 (MQ=255)
          ggAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaa                                                               >  1:1555194/1‑76 (MQ=255)
           gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                              <  1:3863564/76‑1 (MQ=255)
           gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                              <  1:4154972/76‑1 (MQ=255)
           gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                              <  1:4611346/76‑1 (MQ=255)
             gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAAc                                                            <  1:2567454/76‑1 (MQ=255)
              ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACg                                                           <  2:33668/76‑1 (MQ=255)
               tCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGa                                                          <  2:3675437/76‑1 (MQ=255)
               tCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGa                                                          <  2:1483680/76‑1 (MQ=255)
                cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAACCGAt                                                         >  2:3489294/1‑76 (MQ=255)
                 ctCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtc                                                        <  2:3642464/75‑1 (MQ=255)
                 aTCGATGATGACGAACTGGTAGAAGTGACTCAGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtc                                                        <  2:1066740/76‑1 (MQ=255)
                 aTCGATGAGGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtc                                                        <  2:4453705/76‑1 (MQ=255)
                   cGATGGTGACGAACTGGGAGAAGTGACTCCGGCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgt                                                      <  2:4780169/76‑1 (MQ=255)
                   cGATGATGACGAACTGGTAGAAGCGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgc                                                      >  1:233950/1‑75 (MQ=255)
                    gatgctGCCGGACTGGTAGAAGTGGCTCCGGCCTCTATCCGTCTTGGTGAACGTCACCTGACGGAAAACGAtcttc                                                     <  2:1245299/76‑1 (MQ=255)
                    gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                     >  2:511199/1‑76 (MQ=255)
                    gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                     >  1:3936796/1‑76 (MQ=255)
                    gatgatGACGAACTGGGAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                     <  2:295558/76‑1 (MQ=255)
                     gtgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtcg                                                    <  2:684998/75‑1 (MQ=255)
                     atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtcg                                                    >  1:645269/1‑76 (MQ=255)
                     atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtcg                                                    >  1:4846948/1‑76 (MQ=255)
                      tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgc                                                   <  1:4381806/76‑1 (MQ=255)
                      tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgc                                                   <  2:3548417/76‑1 (MQ=255)
                       tatgaCGAACTGGGAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACAGAAAACGATCGTcgcc                                                  <  2:464987/75‑1 (MQ=255)
                       gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                  <  1:781049/76‑1 (MQ=255)
                       gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                  <  1:3647306/76‑1 (MQ=255)
                       gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                  >  2:3484282/1‑76 (MQ=255)
                        atgaCGAACTGGTAGAAGTGGCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccg                                                 <  2:178746/76‑1 (MQ=255)
                        atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccg                                                 >  2:3350430/1‑76 (MQ=255)
                        atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccg                                                 >  2:3937897/1‑76 (MQ=255)
                         tgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccgc                                                <  2:2695415/76‑1 (MQ=255)
                          gaCGAGCTTGTGGTAGTGACTCCGACCTCTAGCTGTATTCTTAATCGCCGCCGGTCTGAAAACTATCGTCGCcgcg                                               <  1:140288/76‑1 (MQ=255)
                            cGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGcc                                             >  1:2848523/1‑76 (MQ=255)
                             gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCa                                            <  1:4859485/76‑1 (MQ=255)
                               aCTGGTGAAAGTGACTCCGGCCTCTATCCGTGTTCGTAAACGTCACCTGACGGAAAACGATCTTCGCCGCGCCAAc                                          <  2:4774543/76‑1 (MQ=255)
                               aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAc                                          <  2:3544969/76‑1 (MQ=255)
                                cTGGTAGAAGTGACTCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAcc                                         <  1:2433180/76‑1 (MQ=255)
                                cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGCAAACGATCGGCCCCGGGCCAAcc                                         >  1:388436/1‑76 (MQ=255)
                                cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAcc                                         <  1:4136997/76‑1 (MQ=255)
                                 tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcg                                        <  1:3144946/76‑1 (MQ=255)
                                 tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcg                                        <  2:4013274/76‑1 (MQ=255)
                                  ggTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcga                                       >  1:4784429/1‑75 (MQ=255)
                                   gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgcg                                      <  2:2597521/76‑1 (MQ=255)
                                   gTAGAAGTGACTCCGACCCCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgcg                                      <  2:2609025/76‑1 (MQ=255)
                                    tAGAAGTGACTCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCgcgc                                     <  2:70299/76‑1 (MQ=255)
                                    tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCgcgc                                     <  2:1854219/76‑1 (MQ=255)
                                    tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCgcgc                                     <  2:1975108/76‑1 (MQ=255)
                                    tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCgcgc                                     <  2:4760808/76‑1 (MQ=255)
                                     aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCa                                    <  1:2497945/76‑1 (MQ=255)
                                      gAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAc                                   >  2:3336336/1‑76 (MQ=255)
                                      aaaGTGACTCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAc                                   <  1:3027815/75‑1 (MQ=255)
                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAcc                                  <  2:4283962/76‑1 (MQ=255)
                                       aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAcc                                  >  1:18426/1‑76 (MQ=255)
                                        aGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCg                                 <  1:3146068/76‑1 (MQ=255)
                                        aGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCg                                 <  2:3126109/76‑1 (MQ=255)
                                           gACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaaa                              <  1:4310556/76‑1 (MQ=255)
                                            aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCTCCAACCGCGCACCGAAAg                             <  1:1495106/76‑1 (MQ=255)
                                            aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAg                             <  1:3685474/76‑1 (MQ=255)
                                             cTCCTAACTCTATCCGTAGTCGTAAACGGCACCTAACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGa                            <  1:2787824/76‑1 (MQ=255)
                                             cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGa                            >  1:3186757/1‑76 (MQ=255)
                                              tCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGAc                           <  1:1073736/76‑1 (MQ=255)
                                               ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                          >  2:3491527/1‑76 (MQ=255)
                                               ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                          <  2:4863282/76‑1 (MQ=255)
                                                cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGACAGACGa                         >  1:3779493/1‑76 (MQ=255)
                                                cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                         <  1:3217625/76‑1 (MQ=255)
                                                cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                         >  1:3042433/1‑76 (MQ=255)
                                                cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                         <  1:1865634/76‑1 (MQ=255)
                                                cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                         >  1:3877277/1‑76 (MQ=255)
                                                 gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAt                        <  2:1518469/76‑1 (MQ=255)
                                                 gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAt                        >  2:3771564/1‑76 (MQ=255)
                                                  gccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAtt                       <  2:2159143/75‑1 (MQ=255)
                                                   ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTa                      >  2:3284988/1‑76 (MQ=255)
                                                   ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTa                      <  1:2834215/76‑1 (MQ=255)
                                                   ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTa                      <  1:2392675/76‑1 (MQ=255)
                                                    ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTaa                     <  2:540712/76‑1 (MQ=255)
                                                        aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTc                 <  1:1723032/76‑1 (MQ=255)
                                                         tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCg                <  1:41037/76‑1 (MQ=255)
                                                          gcGTATTGGTAAACGTGACCTGACGGAAACCGATCGTCGCCGCGCCAACCGCGCACCGAAAGGCGATTAATTTCGt               <  2:2114639/75‑1 (MQ=255)
                                                          ccGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGt               >  1:2119944/1‑76 (MQ=255)
                                                            gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttt             <  1:2671858/76‑1 (MQ=255)
                                                             tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttt            <  1:1368419/76‑1 (MQ=255)
                                                             tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttt            <  2:2008964/76‑1 (MQ=255)
                                                               ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGCCGATTAATTTCGtttttt          >  2:135405/1‑76 (MQ=255)
                                                                tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTa         <  2:1486141/76‑1 (MQ=255)
                                                                 cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAg        <  1:2656979/76‑1 (MQ=255)
                                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTc      <  2:2956436/76‑1 (MQ=255)
                                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTc      <  2:3795562/76‑1 (MQ=255)
                                                                   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTc      <  1:2664515/76‑1 (MQ=255)
                                                                    aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCa     <  1:2253309/76‑1 (MQ=255)
                                                                     aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAt    >  1:1535173/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 < 23 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: