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

GCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720552‑3720639                                             |                                            gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                       >  1:1089671/1‑50 (MQ=255)  cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      <  1:2128952/50‑1 (MQ=255)  cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      >  1:2669022/1‑50 (MQ=255)   gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      >  1:1604021/1‑49 (MQ=39)    tGAAATAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                    <  1:1845542/50‑1 (MQ=255)      aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                   <  1:904339/49‑1 (MQ=255)      aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                  <  1:2416069/50‑1 (MQ=255)      aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                  <  1:2034206/50‑1 (MQ=255)      aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                  >  1:1578650/1‑50 (MQ=255)       aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                 <  1:1606054/50‑1 (MQ=255)         cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                               >  1:563818/1‑50 (MQ=255)          aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                              >  1:643598/1‑50 (MQ=255)          aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                              >  1:188273/1‑50 (MQ=255)           aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                             <  1:910057/50‑1 (MQ=255)           aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                             <  1:1500015/50‑1 (MQ=255)            aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                            >  1:1461916/1‑50 (MQ=255)               aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         >  1:1049089/1‑50 (MQ=255)                cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:113789/1‑50 (MQ=255)                cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:2182608/1‑50 (MQ=255)                   ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                       <  1:1807315/48‑1 (MQ=255)                   ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                     >  1:859301/1‑50 (MQ=255)                   ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                     <  1:368991/50‑1 (MQ=255)                    tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                    >  1:361744/1‑50 (MQ=255)                      aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                   <  1:647515/49‑1 (MQ=255)                      aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                   <  1:966935/49‑1 (MQ=255)                      aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:2011535/1‑50 (MQ=255)                      aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:2723155/1‑50 (MQ=255)                        gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                <  1:123840/50‑1 (MQ=255)                          aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              <  1:993438/50‑1 (MQ=255)                          aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              <  1:192319/50‑1 (MQ=255)                            gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              >  1:1927268/1‑48 (MQ=255)                            gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGa            <  1:332015/50‑1 (MQ=255)                              aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg           >  1:1536568/1‑49 (MQ=255)                               aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa         <  1:986428/50‑1 (MQ=255)                                 cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        >  1:2267648/1‑49 (MQ=255)                                    tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact    <  1:1225591‑M1/50‑5 (MQ=255)                                    tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac     <  1:785277‑M1/49‑4 (MQ=255)                                    tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa      <  1:470432‑M1/48‑3 (MQ=255)                                      tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  <  1:1585772‑M1/50‑7 (MQ=255)                                       ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  >  1:2793089‑M1/1‑43 (MQ=255)                                       ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  >  1:2038540‑M1/1‑43 (MQ=255)                                             |                                            GCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720552‑3720639

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