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

Only 100 of 136 total aligned reads displayed.

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720547‑3720640
                                                 |                                            
tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                              <  1:2781313/50‑1 (MQ=255)
 aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                              >  1:1285917/1‑49 (MQ=255)
  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                            >  2:3203002/1‑48 (MQ=255)
  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                            >  2:5320539/1‑48 (MQ=255)
  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgc                                             <  1:1310986/49‑2 (MQ=255)
       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                      <  1:4596501/50‑1 (MQ=255)
       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                      <  1:3309691/50‑1 (MQ=255)
       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                      <  2:2306594/50‑1 (MQ=255)
       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:790867/1‑49 (MQ=39)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  2:418642/1‑48 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                     <  1:5215321/50‑1 (MQ=255)
        tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                     <  1:297881/50‑1 (MQ=255)
         gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                    >  2:1421449/1‑50 (MQ=255)
         gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                    >  2:1620810/1‑50 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   <  1:220995/50‑1 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   >  2:3001223/1‑50 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   >  2:1663012/1‑50 (MQ=255)
          aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   <  1:4046298/50‑1 (MQ=255)
           aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                  <  1:2062842/50‑1 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                 <  2:5426348/50‑1 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                 >  1:3034939/1‑50 (MQ=255)
            aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                 <  1:3474324/50‑1 (MQ=255)
             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                <  1:5660722/50‑1 (MQ=255)
             cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                >  1:1086802/1‑50 (MQ=255)
              aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                               >  1:1482217/1‑50 (MQ=255)
               aaTCACTATTAAAGAAAGAAATCTATATTTCTGC‑GCATTCCCAGcctcct                              >  2:3627792/1‑50 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcctct                            <  2:4588223/50‑3 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            >  2:611130/1‑50 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            <  1:965949/50‑1 (MQ=255)
                 tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            >  2:4630629/1‑50 (MQ=255)
                  cACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                           >  2:3274963/1‑50 (MQ=255)
                  cACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                           >  2:3741393/1‑50 (MQ=255)
                   aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                           <  2:3974675/49‑1 (MQ=255)
                   aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                          <  1:4906858/50‑1 (MQ=255)
                   aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                          >  1:2832102/1‑50 (MQ=255)
                    cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         <  1:1769920/50‑1 (MQ=255)
                     tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                           <  2:278938/47‑1 (MQ=255)
                     tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  2:1688593/1‑49 (MQ=255)
                     tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:5152181/1‑49 (MQ=255)
                      aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                       >  2:4287933/1‑50 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      >  2:3286037/1‑50 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      >  1:1468112/1‑50 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      <  2:1923776/50‑1 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      >  2:2799995/1‑50 (MQ=255)
                       ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      >  1:5444196/1‑50 (MQ=255)
                        tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     >  1:2677946/1‑50 (MQ=255)
                        tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     >  1:218638/1‑50 (MQ=255)
                         aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     <  1:37790/49‑1 (MQ=255)
                         aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     <  1:2310947/49‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:2912855/1‑50 (MQ=255)
                          aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   <  1:509306/50‑1 (MQ=255)
                          aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   <  2:832067/50‑1 (MQ=255)
                           agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  2:4628453/1‑49 (MQ=255)
                           agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  2:5045963/1‑49 (MQ=255)
                            gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  >  2:774808/1‑49 (MQ=255)
                            gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  >  2:1314823/1‑49 (MQ=255)
                            gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  >  1:29283/1‑49 (MQ=255)
                             aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                <  2:1891962/50‑1 (MQ=255)
                             aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                >  2:227932/1‑50 (MQ=255)
                             aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                >  2:5257833/1‑50 (MQ=255)
                              aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  2:561727/1‑50 (MQ=255)
                              aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  2:2016509/1‑50 (MQ=255)
                              aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  2:1061358/1‑50 (MQ=255)
                               agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  2:1821695/1‑49 (MQ=255)
                               agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  2:4103014/1‑49 (MQ=255)
                               agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  2:3307384/1‑49 (MQ=255)
                                 tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  2:415509/1‑50 (MQ=255)
                                 tAATCTATATTACTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  2:1093680/1‑50 (MQ=255)
                                  aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  1:3771185/1‑49 (MQ=255)
                                  aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            >  2:4373458/1‑49 (MQ=255)
                                   aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt           <  1:1430623/49‑1 (MQ=255)
                                   aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:3045737/1‑50 (MQ=255)
                                   aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          <  2:1522873/50‑1 (MQ=255)
                                   aTCTATATTTCTGC‑ACATTCCCAGCCTCCTGTTTTGATTTCCATCGAGt           <  2:4408776/49‑1 (MQ=255)
                                    tCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:5319051/1‑49 (MQ=255)
                                     cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  2:1088868/1‑49 (MQ=255)
                                     cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:2307553/1‑49 (MQ=255)
                                      tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       >  1:4466843‑M1/1‑48 (MQ=255)
                                      tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       >  2:1039145‑M1/1‑48 (MQ=255)
                                      tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       >  2:4252346‑M1/1‑48 (MQ=255)
                                       atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac      >  2:5088171‑M1/1‑47 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         <  2:2760549/46‑1 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         <  1:4692283/46‑1 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         <  2:1033055/46‑1 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact     >  2:5630764‑M1/1‑46 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact     <  2:623694‑M1/50‑5 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact     <  1:3570148‑M1/50‑5 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact     <  1:3069898‑M1/50‑5 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac      <  1:5067777‑M1/49‑4 (MQ=255)
                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa        <  1:1557764‑M1/47‑2 (MQ=255)
                                         atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:745932‑M1/1‑45 (MQ=255)
                                         atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:3940851‑M1/1‑45 (MQ=255)
                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  2:5192174‑M1/50‑7 (MQ=255)
                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  1:3939471‑M1/50‑7 (MQ=255)
                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  1:2215512‑M1/50‑7 (MQ=255)
                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  2:4012691‑M1/50‑7 (MQ=255)
                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  2:86925‑M1/1‑44 (MQ=255)
                                           ttCTGC‑GCATTCCCCGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  2:3621965‑M1/1‑43 (MQ=255)
                                           ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  2:3009884‑M1/1‑43 (MQ=255)
                                            tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg  <  2:4121632‑M1/49‑8 (MQ=255)
                                                 |                                            
TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720547‑3720640

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