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

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641                                                  |                                              tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:2825909/50‑1 (MQ=255) tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:1771134/50‑1 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             <  1:6220953/50‑3 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgc                                              <  1:4921531/49‑2 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        <  1:2574237/50‑1 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:3097035/1‑50 (MQ=255)        gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                       <  1:6686268/50‑1 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                      <  1:4754270/50‑1 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    >  1:5610536/1‑50 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    >  1:4785865/1‑50 (MQ=255)            aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                   <  1:342687/50‑1 (MQ=255)            aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                   >  1:325235/1‑50 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  >  1:5484025/1‑50 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  <  1:3292398/50‑1 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  <  1:1393637/50‑1 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  <  1:1264081/50‑1 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  >  1:5257257/1‑50 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 <  1:1202221/50‑1 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:6114192/1‑50 (MQ=255)               aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                                <  1:1884266/50‑1 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCTGCCTCCtg                              >  1:6717645/1‑50 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              <  1:4852496/50‑1 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              <  1:2413579/50‑1 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              >  1:200659/1‑50 (MQ=255)                  tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                             >  1:5443855/1‑50 (MQ=255)                  tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                             <  1:4148601/50‑1 (MQ=255)                  tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                             <  1:5850990/50‑1 (MQ=255)                  tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                             <  1:2893249/50‑1 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                           >  1:2543598/1‑50 (MQ=255)                     cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  1:3066740/1‑50 (MQ=255)                      tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  1:3376749/1‑49 (MQ=255)                       aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                        >  1:4742628/1‑50 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                        <  1:6050508/49‑1 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                       <  1:6403032/50‑1 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                       >  1:108388/1‑50 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                       >  1:3874110/1‑50 (MQ=255)                         tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      >  1:2878416/1‑50 (MQ=255)                         tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      >  1:4446488/1‑50 (MQ=255)                          aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      <  1:541265/49‑1 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     <  1:319829/49‑1 (MQ=255)                            agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   <  1:2506759/50‑1 (MQ=255)                            agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   <  1:6460828/50‑1 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   >  1:980145/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   >  1:4631077/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:6148057/50‑1 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:4255264/50‑1 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:2689576/50‑1 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:5282280/1‑50 (MQ=255)                                agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:4060250/1‑49 (MQ=255)                                 gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGa              <  1:726169/50‑1 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:3445764/1‑50 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:3383395/50‑1 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           <  1:3456642/50‑1 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           <  1:2825869/50‑1 (MQ=255)                                     tCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:3938010/1‑49 (MQ=255)                                      cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          >  1:336018/1‑49 (MQ=255)                                       tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:2813126‑M1/1‑48 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:2042452‑M1/1‑47 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:3205418‑M1/1‑47 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:4610978‑M1/1‑47 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  1:561384/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      >  1:2642932‑M1/1‑46 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      <  1:1357447‑M1/50‑5 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      <  1:3852552‑M1/50‑5 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      >  1:4315848‑M1/1‑46 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:4979352‑M1/47‑2 (MQ=255)                                          atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  1:6663891‑M1/1‑45 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:6222298‑M1/1‑44 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:2331515‑M1/1‑44 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    <  1:3479626‑M1/50‑7 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:744497‑M1/1‑44 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:1763675‑M1/49‑6 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:2979003‑M1/1‑42 (MQ=255)                                                  |                                              TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641

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