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_000913	3,720,596	0	G	.	100.0%	209.8 / NA	44	coding (5/153 nt)	hokA	small toxic polypeptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (24/20);  total (24/20)

AAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720548‑3720640                                                 |                                             aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                              >  1:2138191/1‑49 (MQ=255)  aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                            <  1:1878249/50‑3 (MQ=255)      cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:1826353/1‑50 (MQ=255)       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                      <  1:1163294/50‑1 (MQ=255)       gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:1801534/1‑49 (MQ=39)             cAAATCACTATTGAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                >  1:2318958/1‑50 (MQ=255)              aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                               >  1:770746/1‑50 (MQ=255)               aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              <  1:588146/50‑1 (MQ=255)               aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              <  1:443608/50‑1 (MQ=255)                 tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            <  1:478101/50‑1 (MQ=255)                 tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            <  1:2098859/50‑1 (MQ=255)                   aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                          >  1:16285/1‑50 (MQ=255)                   aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                          >  1:1857451/1‑50 (MQ=255)                     tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:942085/1‑49 (MQ=255)                     tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:969733/1‑49 (MQ=255)                       ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      >  1:152335/1‑50 (MQ=255)                         aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                    >  1:226723/1‑50 (MQ=255)                          aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                    <  1:1905346/49‑1 (MQ=255)                          aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:1785424/1‑50 (MQ=255)                           agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  <  1:631532/50‑1 (MQ=255)                           agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  <  1:2195512/50‑1 (MQ=255)                           agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:1456312/1‑49 (MQ=255)                            gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                 <  1:526325/50‑1 (MQ=255)                              aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:1638850/1‑50 (MQ=255)                              aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:406786/1‑50 (MQ=255)                              aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               <  1:1455371/50‑1 (MQ=255)                              aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               <  1:953837/50‑1 (MQ=255)                               agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg              <  1:2116012/50‑1 (MQ=255)                                gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGa             <  1:128459/50‑1 (MQ=255)                                 tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            <  1:1925247/50‑1 (MQ=255)                                   aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt           <  1:489863/49‑1 (MQ=255)                                   aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:613110/1‑50 (MQ=255)                                   aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:1490281/1‑50 (MQ=255)                                    tCTATATTTCTGC‑GCATTACCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:413539/1‑49 (MQ=255)                                     cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:125017/1‑49 (MQ=255)                                     cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:319408/1‑49 (MQ=255)                                     cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         >  1:1443582/1‑49 (MQ=255)                                       atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac      >  1:890584‑M1/1‑47 (MQ=255)                                       atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac      >  1:939746‑M1/1‑47 (MQ=255)                                        tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       <  1:933878‑M1/48‑3 (MQ=255)                                         atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:210203‑M1/1‑45 (MQ=255)                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  1:585513‑M1/50‑7 (MQ=255)                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  1:2223688‑M1/50‑7 (MQ=255)                                          tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  1:1491990‑M1/1‑44 (MQ=255)                                            tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg  <  1:429848‑M1/49‑8 (MQ=255)                                            tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg  <  1:787128‑M1/49‑8 (MQ=255)                                                 |                                             AAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720548‑3720640

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: ‑