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

AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720549‑3720639                                                |                                            aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgc                                            <  1:1483351/49‑2 (MQ=255) aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                             <  1:1349985/48‑1 (MQ=255)  aaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                           >  1:1124663/1‑47 (MQ=255)  aaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                           <  1:350463/49‑3 (MQ=255)  aaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                           >  1:1701821/1‑47 (MQ=255)       tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                     <  1:1943329/49‑1 (MQ=255)        gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      >  1:1416879/1‑47 (MQ=255)         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                   >  1:1940310/1‑49 (MQ=255)           aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                 >  1:130872/1‑49 (MQ=255)               aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                             <  1:1983895/49‑1 (MQ=255)                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                            >  1:1066050/1‑49 (MQ=255)                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                            >  1:468033/1‑49 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                          <  1:986319/49‑1 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                          >  1:1113740/1‑49 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                          >  1:1529893/1‑49 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         >  1:1624860/1‑49 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        <  1:830918/47‑1 (MQ=255)                       tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                     >  1:23861/1‑49 (MQ=255)                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                   >  1:1936911/1‑49 (MQ=255)                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:1630603/1‑49 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                 >  1:32032/1‑48 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                 >  1:1347947/1‑48 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:1117787/1‑49 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:873665/1‑49 (MQ=255)                              agTAATCTATATTTCTGC‑GCATTCTCAGCCTCCTGTGTTGATTTCCAAc              <  1:218317/49‑1 (MQ=38)                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              <  1:278253/49‑1 (MQ=255)                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              <  1:1575128/49‑1 (MQ=255)                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              >  1:88980/1‑47 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt          <  1:729905/49‑1 (MQ=255)                                   tCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa         >  1:60002/1‑49 (MQ=255)                                    cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa         >  1:2058494/1‑48 (MQ=255)                                      atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa      >  1:756839‑M1/1‑47 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        <  1:658021/46‑1 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac     <  1:762532‑M1/49‑4 (MQ=255)                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact    >  1:617817‑M1/1‑45 (MQ=255)                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact    >  1:1231385‑M1/1‑45 (MQ=255)                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  >  1:33249‑M1/1‑42 (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: ‑