Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,720,597 Δ1 bp coding (4/153 nt) hokA ← small toxic polypeptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,720,5960G.100.0% 842.1 / NA 170coding (5/153 nt)hokAsmall toxic polypeptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (88/82);  total (88/82)

Only 100 of 180 total aligned reads displayed.

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641
                                                 |                                             
tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:5627635/50‑1 (MQ=255)
tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:7583799/1‑50 (MQ=255)
tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:6800367/1‑50 (MQ=255)
 aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:1103008/1‑49 (MQ=255)
  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:3419454/1‑48 (MQ=255)
  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:1912324/1‑48 (MQ=255)
  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:7805478/1‑48 (MQ=255)
      cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:7248023/1‑50 (MQ=255)
      cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:5275510/1‑50 (MQ=255)
      cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        <  1:4200934/50‑1 (MQ=255)
      cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:7728747/1‑50 (MQ=255)
      cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        <  1:4457362/50‑1 (MQ=255)
       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTccc                                       <  1:5858266/50‑1 (MQ=255)
       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:2276129/1‑49 (MQ=39)
       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:6401108/1‑49 (MQ=39)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:3657126/1‑48 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:7710972/1‑48 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCa                                      <  1:7942925/50‑1 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCa                                      <  1:3586394/50‑1 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCa                                      <  1:2864130/50‑1 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCa                                      <  1:5805103/50‑1 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCa                                      <  1:5192238/50‑1 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCa                                      <  1:4378785/50‑1 (MQ=255)
         gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAg                                     >  1:5948997/1‑50 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGc                                    >  1:2007804/1‑50 (MQ=255)
           aaCAAATCGCTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcc                                   <  1:249611/50‑1 (MQ=255)
           aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcc                                   <  1:7549433/50‑1 (MQ=255)
           aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcc                                   >  1:3532693/1‑50 (MQ=255)
           aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcc                                   >  1:7659346/1‑50 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcct                                  <  1:2879909/50‑1 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcct                                  >  1:5338019/1‑50 (MQ=255)
             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctc                                 >  1:4728187/1‑50 (MQ=255)
              aaaTCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctcc                                >  1:1463500/1‑50 (MQ=255)
              aaaTCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctcc                                >  1:8149210/1‑50 (MQ=255)
               aaTCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctcct                               <  1:6730466/50‑1 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCtgt                             <  1:3970324/50‑1 (MQ=255)
                   aCTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTgtgt                           >  1:4710600/1‑50 (MQ=255)
                   aCTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTgtgt                           >  1:2951862/1‑50 (MQ=255)
                    cTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGtt                          >  1:2687251/1‑50 (MQ=255)
                    cTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGtt                          <  1:2124389/50‑1 (MQ=255)
                     tATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGtt                          >  1:288757/1‑49 (MQ=255)
                     tATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGtt                          >  1:7780288/1‑49 (MQ=255)
                     tATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGtt                          >  1:152306/1‑49 (MQ=255)
                      aTTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGa                        >  1:8244977/1‑50 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTg                         <  1:7213039/48‑1 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAt                       >  1:4829683/1‑50 (MQ=255)
                         aaagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAtt                      <  1:2214025/49‑1 (MQ=255)
                         aaagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAtt                      <  1:2114545/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAttt                     <  1:3525677/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAttt                     <  1:5095383/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAttt                     <  1:4805616/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAttt                     <  1:837143/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAttt                     <  1:6492035/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTc                    <  1:1026280/50‑1 (MQ=255)
                            gaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTcc                   >  1:4072734/1‑49 (MQ=255)
                            gaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTcc                   >  1:3337065/1‑49 (MQ=255)
                            gaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:4397285/50‑1 (MQ=255)
                            gaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:3087303/50‑1 (MQ=255)
                             aaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:5236296/1‑50 (MQ=255)
                             aaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCaa                 <  1:5877821/50‑1 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:7162755/50‑1 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:3130072/50‑1 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:6245659/50‑1 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:301249/1‑50 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:809097/50‑1 (MQ=255)
                               agTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:1058109/1‑49 (MQ=255)
                               agTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:1184941/1‑49 (MQ=255)
                               agTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACg               <  1:4939227/50‑1 (MQ=255)
                               agTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACg               <  1:701070/50‑1 (MQ=255)
                                gTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:845508/1‑48 (MQ=255)
                                gTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:4628218/1‑48 (MQ=255)
                                 tAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:2243485/50‑1 (MQ=255)
                                  aaTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:6192417/1‑49 (MQ=255)
                                  aaTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:1467783/1‑49 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:4667705/1‑50 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:7816794/1‑50 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:4533565/1‑50 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:3928176/1‑50 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           <  1:2843480/50‑1 (MQ=255)
                                    tCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:7286151/1‑49 (MQ=255)
                                    tCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:2225150/1‑49 (MQ=255)
                                     cTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          >  1:5593918/1‑49 (MQ=255)
                                       atatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:7881539‑M1/1‑47 (MQ=255)
                                       atatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:2777589‑M1/1‑47 (MQ=255)
                                        tatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      >  1:7179842‑M1/1‑46 (MQ=255)
                                        tatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      >  1:2266312‑M1/1‑46 (MQ=255)
                                        tatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       <  1:7650532‑M1/49‑4 (MQ=255)
                                        tatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:2633174‑M1/48‑3 (MQ=255)
                                        tatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:5695138‑M1/47‑2 (MQ=255)
                                         atTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  1:7837203‑M1/1‑45 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:2392885‑M1/1‑44 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:5427513‑M1/1‑44 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:1844204‑M1/49‑6 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:50126‑M1/49‑6 (MQ=255)
                                           ttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:2700757‑M1/1‑43 (MQ=255)
                                           ttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:6867340‑M1/1‑43 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:6316424‑M1/1‑42 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:7227223‑M1/1‑42 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  1:3136127‑M1/49‑8 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  1:6893107‑M1/49‑8 (MQ=255)
                                                 |                                             
TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641

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