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% 635.5 / NA 127coding (5/153 nt)hokAsmall toxic polypeptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (73/54);  total (73/54)

Only 100 of 130 total aligned reads displayed.

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641
                                                 |                                             
tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:2123982/1‑50 (MQ=255)
 aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:325154/1‑49 (MQ=255)
     gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTc                                         <  1:5821037/50‑1 (MQ=255)
     gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTc                                         <  1:374009/50‑1 (MQ=255)
     gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTc                                         >  1:6449917/1‑50 (MQ=255)
     gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTc                                         >  1:1786065/1‑50 (MQ=255)
      cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:1082704/1‑50 (MQ=255)
      cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:4705238/1‑50 (MQ=255)
      cGTGAAACAAATCACTATTAAAGAAAGAAATCTATATTTCTGCGCATTcc                                        >  1:33228/1‑50 (MQ=255)
       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTccc                                       <  1:67941/50‑1 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTcc                                        >  1:7141120/1‑48 (MQ=255)
         gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAg                                     >  1:2679118/1‑50 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAg                                     <  1:7255331/49‑1 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGc                                    <  1:7213067/50‑1 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGc                                    >  1:3271507/1‑50 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGc                                    >  1:2780367/1‑50 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGc                                    >  1:4202863/1‑50 (MQ=255)
          aaaCAAATCACTATAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcc                                   >  1:3144888/1‑50 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcct                                  >  1:6848647/1‑50 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcct                                  >  1:3893253/1‑50 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcct                                  >  1:5231177/1‑50 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcct                                  <  1:5931676/50‑1 (MQ=255)
             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctc                                 >  1:4597446/1‑50 (MQ=255)
              aaaTCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctcc                                <  1:6999301/50‑1 (MQ=255)
              aaaTCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctcc                                <  1:1223748/50‑1 (MQ=255)
              aaaTCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctcc                                >  1:5139497/1‑50 (MQ=255)
               aaTCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGcctcct                               <  1:2623901/50‑1 (MQ=255)
                aTCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCtg                              >  1:6527445/1‑50 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCtgt                             <  1:7194723/50‑1 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCtgt                             >  1:803061/1‑50 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCtgt                             <  1:1145458/50‑1 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCtgt                             >  1:5964681/1‑50 (MQ=255)
                   aCTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCtgtg                            <  1:876463/49‑1 (MQ=255)
                   aCTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCtgtg                            <  1:415239/49‑1 (MQ=255)
                    cTATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGtt                          <  1:3767469/50‑1 (MQ=255)
                     tATTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGtt                          >  1:7143525/1‑49 (MQ=255)
                      aTTAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGa                        >  1:4765905/1‑50 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAt                       <  1:2011181/50‑1 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAt                       >  1:5572182/1‑50 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAt                       >  1:6583500/1‑50 (MQ=255)
                        tAAAGAAAGTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAtt                      >  1:1446722/1‑50 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAttt                     <  1:4840164/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGAttt                     <  1:3948724/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTc                    <  1:7385339/50‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTc                    <  1:4461355/50‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTc                    <  1:3743462/50‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:2356588/1‑50 (MQ=255)
                           agaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTcc                   <  1:4872926/50‑1 (MQ=255)
                           agaaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTcc                   <  1:5625526/50‑1 (MQ=255)
                             aaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:4249125/1‑50 (MQ=255)
                             aaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:4922769/1‑50 (MQ=255)
                             aaagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:6453285/1‑50 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:6411266/1‑50 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:4554329/50‑1 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:6172259/1‑50 (MQ=255)
                              aagTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:2408853/50‑1 (MQ=255)
                               agTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACg               <  1:875841/50‑1 (MQ=255)
                                gTAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:1051780/1‑48 (MQ=255)
                                 tAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:1195658/50‑1 (MQ=255)
                                 tAATCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:2939373/50‑1 (MQ=255)
                                  aaTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:4566886/1‑49 (MQ=255)
                                  aaTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:1218100/1‑49 (MQ=255)
                                  aaTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:3548699/1‑49 (MQ=255)
                                  aaTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:3693612/1‑49 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt            <  1:2045573/49‑1 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt            <  1:1829716/49‑1 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt            <  1:3004958/49‑1 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:93850/1‑50 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:6588331/1‑50 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:5972826/1‑50 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:7595564/1‑50 (MQ=255)
                                   aTCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:5711318/1‑50 (MQ=255)
                                    tCTATATTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:42652/1‑49 (MQ=255)
                                      tataTTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:3832600‑M1/1‑48 (MQ=255)
                                      tataTTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:5622641‑M1/1‑48 (MQ=255)
                                      tataTTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:4231785‑M1/1‑48 (MQ=255)
                                       atatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:1251908‑M1/1‑47 (MQ=255)
                                       atatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:1499742‑M1/1‑47 (MQ=255)
                                       atatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:4130813‑M1/1‑47 (MQ=255)
                                        tatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  1:827085/46‑1 (MQ=255)
                                        tatTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       <  1:5278941‑M1/49‑4 (MQ=255)
                                        tatGTCTGCGCATTACCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       <  1:3400613‑M1/49‑4 (MQ=255)
                                         atTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  1:6115032‑M1/1‑45 (MQ=255)
                                         atTTCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  1:6516028‑M1/1‑45 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:1397175‑M1/1‑44 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    <  1:3539562‑M1/50‑7 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:6649273‑M1/1‑44 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:3810709‑M1/1‑44 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:4882494‑M1/1‑44 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:1644267‑M1/49‑6 (MQ=255)
                                          tttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaaatt     <  1:2590503‑M1/49‑6 (MQ=255)
                                           ttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:3929138‑M1/1‑43 (MQ=255)
                                           ttCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:5060422‑M1/1‑43 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:2740500‑M1/1‑42 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:6333392‑M1/1‑42 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:4843467‑M1/1‑42 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  1:1588065‑M1/49‑8 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  1:5719875‑M1/49‑8 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  1:3952560‑M1/49‑8 (MQ=255)
                                            tCTGCGCATTCCCAGCCTCCTGTGTTGATTACCAACGAGTATaactttg   <  1:5930307‑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 < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: