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

AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTG  >  NC_000913/3720549‑3720639                                                |                                            aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                           >  1:688138/1‑48 (MQ=255)     cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      <  1:3360216/50‑1 (MQ=255)     cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      >  1:3751052/1‑50 (MQ=255)     cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      >  1:3818690/1‑50 (MQ=255)      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                      >  1:1388181/1‑49 (MQ=39)         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                  >  1:975803/1‑50 (MQ=255)         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                  <  1:3047180/50‑1 (MQ=255)           aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                <  1:244356/50‑1 (MQ=255)            cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                               <  1:3085591/50‑1 (MQ=255)              aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                             <  1:3024793/50‑1 (MQ=255)               aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                            <  1:173123/50‑1 (MQ=255)               aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                            >  1:2520450/1‑50 (MQ=255)                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                           >  1:1872402/1‑50 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         >  1:4090280/1‑50 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                         <  1:1238930/50‑1 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:1343268/1‑50 (MQ=255)                    tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                        >  1:901406/1‑49 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                      >  1:628511/1‑50 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                       <  1:3694440/48‑1 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                       <  1:3676938/48‑1 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                       <  1:1136904/48‑1 (MQ=255)                       tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                    >  1:2739251/1‑50 (MQ=255)                        aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                   >  1:3674788/1‑50 (MQ=255)                        aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                    <  1:2259145/49‑1 (MQ=255)                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:3171312/1‑50 (MQ=255)                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:844341/1‑50 (MQ=255)                          agacagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:4450366/1‑49 (MQ=255)                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                 <  1:3447922/50‑1 (MQ=255)                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                 <  1:4203388/50‑1 (MQ=255)                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                  >  1:2055438/1‑49 (MQ=255)                           gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                <  1:3594624/50‑1 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:1452638/1‑50 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:2599886/1‑50 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               <  1:2381369/50‑1 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa               >  1:4192529/1‑50 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              >  1:1838266/1‑50 (MQ=255)                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc              >  1:4237393/1‑49 (MQ=255)                              agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg             <  1:3400172/50‑1 (MQ=255)                                 aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg           >  1:110786/1‑49 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt          <  1:3460520/49‑1 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa         >  1:3747225/1‑50 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact    >  1:2251513‑M1/1‑46 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact    <  1:3482985‑M1/50‑5 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac     <  1:3505412‑M1/49‑4 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac     <  1:1445235‑M1/49‑4 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa       <  1:2169187‑M1/47‑2 (MQ=255)                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  >  1:3726077‑M1/1‑44 (MQ=255)                                          ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt  >  1:3820325‑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 < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑