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

AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720549‑3720640                                                |                                             aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                            <  1:5284550/50‑3 (MQ=255)    gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                        <  1:2717465/50‑1 (MQ=255)      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                      <  1:3106766/50‑1 (MQ=255)      gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:2590189/1‑49 (MQ=39)       tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                       >  1:5850175/1‑48 (MQ=255)         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   >  1:1479685/1‑50 (MQ=255)         aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                   >  1:5694685/1‑50 (MQ=255)              aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                              <  1:4780792/50‑1 (MQ=255)                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            <  1:2797215/50‑1 (MQ=255)                tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                            >  1:706051/1‑50 (MQ=255)                  aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                           <  1:2874028/49‑1 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         >  1:996500/1‑50 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         <  1:4667440/50‑1 (MQ=255)                   cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                         <  1:1403811/50‑1 (MQ=255)                     aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                       >  1:1331657/1‑50 (MQ=255)                      ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                      <  1:5099771/50‑1 (MQ=255)                        aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                    >  1:4527948/1‑50 (MQ=255)                        aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                     <  1:3761532/49‑1 (MQ=255)                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                    <  1:2377904/49‑1 (MQ=255)                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:1255708/1‑50 (MQ=255)                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:1742813/1‑50 (MQ=255)                         aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   <  1:2603455/50‑1 (MQ=255)                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  <  1:6218451/50‑1 (MQ=255)                          agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                   >  1:2489420/1‑49 (MQ=255)                           gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                  >  1:1469168/1‑49 (MQ=255)                           gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                 <  1:2859708/50‑1 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                <  1:1002328/50‑1 (MQ=255)                            aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                <  1:3480243/50‑1 (MQ=255)                             aagTAATCTATATTTCTGC‑GCATTCTCAGCCTCCTGTGTTGATTTCCAAc               >  1:3826688/1‑50 (MQ=38)                             aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               <  1:4707807/50‑1 (MQ=255)                               gTAATCTATATTTCTGC‑TCATTCCCAGCCTCCTGTGTTGATTTCCAACGa             <  1:3791572/50‑1 (MQ=38)                               gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc               >  1:5818238/1‑48 (MQ=255)                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            <  1:3139282/50‑1 (MQ=255)                                tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg            <  1:3636601/50‑1 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          <  1:2646501/50‑1 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          >  1:5172616/1‑50 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          <  1:4359131/50‑1 (MQ=255)                                  aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa          <  1:35115/50‑1 (MQ=255)                                     tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa       >  1:4808155‑M1/1‑48 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt         <  1:1740441/46‑1 (MQ=255)                                       tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact     >  1:3001412‑M1/1‑46 (MQ=255)                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:1589565‑M1/1‑45 (MQ=255)                                        atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt    >  1:569821‑M1/1‑45 (MQ=255)                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  1:5173799‑M1/1‑44 (MQ=255)                                         tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   <  1:6126004‑M1/50‑7 (MQ=255)                                          ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt   >  1:1516606‑M1/1‑43 (MQ=255)                                           tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg  <  1:5438168‑M1/49‑8 (MQ=255)                                                |                                             AAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGC  >  NC_000913/3720549‑3720640

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