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

AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720549‑3720639                                                |                                            aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                             >  1:568918/1‑48 (MQ=255)  aaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                           >  1:1834528/1‑47 (MQ=255)     cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                       >  1:241613/1‑49 (MQ=39)      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      <  1:511576/49‑1 (MQ=39)           aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                 >  1:690584/1‑49 (MQ=255)               aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                             <  1:1419555/49‑1 (MQ=255)               aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                             >  1:1730453/1‑49 (MQ=255)                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                            <  1:118014/49‑1 (MQ=255)                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                            >  1:608491/1‑49 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         >  1:1230639/1‑49 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         >  1:1703106/1‑49 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:404557/1‑48 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:477000/1‑48 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:661023/1‑48 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:1085377/1‑48 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                      >  1:208653/1‑49 (MQ=255)                       tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                     >  1:1953440/1‑49 (MQ=255)                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                   >  1:1256298/1‑49 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               <  1:729257/49‑1 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:1956465/1‑49 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:1370566/1‑49 (MQ=255)                               gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg             <  1:285309/49‑1 (MQ=255)                                 aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg           >  1:1107079/1‑49 (MQ=255)                                     tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt        >  1:54696/1‑48 (MQ=255)                                      atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa      >  1:432446‑M1/1‑47 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa         <  1:1661522/45‑1 (MQ=255)                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact    >  1:1855694‑M1/1‑45 (MQ=255)                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt   <  1:261961‑M1/49‑6 (MQ=255)                                          ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  >  1:311985‑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 < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑