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

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641                                                  |                                              tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:3405494/50‑1 (MQ=255) tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:1522797/50‑1 (MQ=255)  aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:3181180/1‑49 (MQ=255)  aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:1967643/1‑49 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             <  1:641889/50‑3 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        <  1:1062017/50‑1 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        <  1:3216334/50‑1 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:29198/1‑50 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                      <  1:673082/50‑1 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                      <  1:3119292/50‑1 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                      <  1:2124565/50‑1 (MQ=255)          gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                     >  1:2857058/1‑50 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                     <  1:3068779/49‑1 (MQ=255)            aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                   >  1:3452137/1‑50 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  <  1:1548254/50‑1 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:2509953/1‑50 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:2652369/1‑50 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              <  1:136918/50‑1 (MQ=255)                     cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  1:76872/1‑50 (MQ=255)                         tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      >  1:2945166/1‑50 (MQ=255)                          aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     >  1:2138333/1‑50 (MQ=255)                          aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     >  1:2250612/1‑50 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     <  1:1160235/49‑1 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    <  1:1911206/50‑1 (MQ=255)                            agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   <  1:3127669/50‑1 (MQ=255)                            agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:1990255/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   >  1:1736739/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:3695487/50‑1 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:1306146/1‑50 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:246700/50‑1 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:1088658/50‑1 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:3067325/1‑50 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:3156294/50‑1 (MQ=255)                                agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:3075851/1‑49 (MQ=255)                                 gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:1702175/1‑48 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:1372543/50‑1 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:2033890/50‑1 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:605453/1‑50 (MQ=255)                                   aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:8697/1‑49 (MQ=255)                                   aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:1057148/1‑49 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt            <  1:974501/49‑1 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           <  1:2888253/50‑1 (MQ=255)                                      cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          >  1:2457989/1‑49 (MQ=255)                                       tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:2547688‑M1/1‑48 (MQ=255)                                       tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:806603‑M1/1‑48 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:1530556‑M1/1‑47 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:116640‑M1/47‑2 (MQ=255)                                          atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  1:953512‑M1/1‑45 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:836130‑M1/49‑6 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:3515911‑M1/1‑42 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  1:3176127‑M1/49‑8 (MQ=255)                                                  |                                              TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641

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: ‑