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%	233.9 / NA	48	coding (5/153 nt)	hokA	small toxic polypeptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (26/22);  total (26/22)

AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720549‑3720639                                                |                                            aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                             >  1:120954/1‑48 (MQ=255) aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                             >  1:2374801/1‑48 (MQ=255) aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                             >  1:533464/1‑48 (MQ=255)     cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                       >  1:746090/1‑49 (MQ=39)        gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      >  1:1551102/1‑47 (MQ=255)         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                   >  1:300613/1‑49 (MQ=255)               aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                             >  1:801586/1‑49 (MQ=255)                 cACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                           >  1:674099/1‑49 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                          >  1:2347440/1‑49 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                          >  1:1876370/1‑49 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                           <  1:2565884/48‑1 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         >  1:2083111/1‑49 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         <  1:1059395/49‑1 (MQ=255)                    tATTAAAGACAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:2605230/1‑49 (MQ=255)                    tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:371034/1‑49 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:33218/1‑48 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:2251151/1‑48 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:2095139/1‑48 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        <  1:1916027/47‑1 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                       <  1:1890751/48‑1 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                       <  1:326928/48‑1 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                      <  1:600841/49‑1 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                      >  1:1537700/1‑49 (MQ=255)                        aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                     <  1:1213356/48‑1 (MQ=255)                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  <  1:159195/49‑1 (MQ=255)                           gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:1932411/1‑48 (MQ=255)                           gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:1732139/1‑48 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                <  1:2336803/49‑1 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:2232663/1‑49 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:1562769/1‑49 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               <  1:2542169/49‑1 (MQ=255)                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              >  1:460831/1‑49 (MQ=255)                               gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg             <  1:1396666/49‑1 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg           >  1:1613271/1‑48 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg           >  1:1716707/1‑48 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt          <  1:841019/49‑1 (MQ=255)                                     tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        >  1:1895003/1‑48 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa         <  1:1351178/45‑1 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        <  1:2500732/46‑1 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        <  1:1779726/46‑1 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        <  1:1621347/46‑1 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac     <  1:2651105‑M1/49‑4 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa      <  1:270745‑M1/48‑3 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa      <  1:2359070‑M1/48‑3 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa       <  1:2435992‑M1/47‑2 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa       <  1:829440‑M1/47‑2 (MQ=255)                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact    >  1:2320739‑M1/1‑45 (MQ=255)                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt   <  1:479162‑M1/49‑6 (MQ=255)                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt   >  1:820723‑M1/1‑44 (MQ=255)                                          ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  >  1:1671213‑M1/1‑43 (MQ=255)                                          ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  >  1:1481047‑M1/1‑43 (MQ=255)                                                |                                            AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720549‑3720639

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