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

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641                                                  |                                              tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:5078543/50‑1 (MQ=255) tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:1700259/1‑50 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:2229544/1‑48 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             <  1:1537322/50‑3 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgc                                              <  1:283357/49‑2 (MQ=255)      gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                         <  1:768875/50‑1 (MQ=255)      gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                         >  1:1648653/1‑50 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:1909721/1‑50 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:893027/1‑50 (MQ=255)        gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                       <  1:3093366/50‑1 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:3077138/1‑48 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:3526295/1‑48 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                      <  1:758152/50‑1 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                      <  1:4054221/50‑1 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    <  1:1094812/50‑1 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    <  1:2820163/50‑1 (MQ=255)            aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                   >  1:666396/1‑50 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  >  1:3597193/1‑50 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  <  1:2848096/50‑1 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  >  1:1444924/1‑50 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:805408/1‑50 (MQ=255)               aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                                <  1:285627/50‑1 (MQ=255)               aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                                >  1:4179650/1‑50 (MQ=255)                aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                               >  1:2035593/1‑50 (MQ=255)                aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                               <  1:759093/50‑1 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              >  1:3579234/1‑50 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              <  1:4109285/50‑1 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              <  1:407311/50‑1 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            <  1:3706505/49‑1 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                           >  1:2421496/1‑50 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                           >  1:2752956/1‑50 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                           <  1:314444/50‑1 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTgtgt                           >  1:3985486/1‑50 (MQ=255)                      tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  1:213668/1‑49 (MQ=255)                       aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                        >  1:1466822/1‑50 (MQ=255)                       aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                        >  1:1849113/1‑50 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                         <  1:1984318/48‑1 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                       <  1:3018417/50‑1 (MQ=255)                         tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      >  1:2886235/1‑50 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     <  1:3083577/49‑1 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:377344/1‑50 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    <  1:2809361/50‑1 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:3521777/1‑50 (MQ=255)                            agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:4880224/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   >  1:2674869/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   >  1:870812/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:4360353/50‑1 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:2048382/50‑1 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:3454016/50‑1 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:1204044/1‑50 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 <  1:2989007/50‑1 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:690470/50‑1 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:2965122/1‑50 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:4032523/1‑50 (MQ=255)                                agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:1322997/1‑49 (MQ=255)                                agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:423931/1‑49 (MQ=255)                                agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg               <  1:1256200/50‑1 (MQ=255)                                agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg               <  1:3151779/50‑1 (MQ=255)                                 gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:339252/1‑48 (MQ=255)                                 gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:3576036/1‑48 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:1531165/50‑1 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:1720173/1‑50 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:3173286/1‑50 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt            <  1:878145/49‑1 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt            <  1:4967164/49‑1 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:3457120/1‑50 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           <  1:4087507/50‑1 (MQ=255)                                     tCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:3317376/1‑49 (MQ=255)                                       tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:4628976‑M1/1‑48 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:365314‑M1/1‑47 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:2291996‑M1/1‑47 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:4788152‑M1/1‑47 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  1:2906571/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  1:3538478/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:1988187‑M1/48‑3 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:3232669‑M1/48‑3 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:4229047‑M1/48‑3 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:3310565‑M1/48‑3 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:2149935‑M1/47‑2 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:2225657‑M1/47‑2 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:3873883‑M1/47‑2 (MQ=255)                                          atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  1:3123751‑M1/1‑45 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    <  1:4345428‑M1/50‑7 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:142952‑M1/49‑6 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:1260863‑M1/49‑6 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:2540778‑M1/1‑42 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  1:3593621‑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 < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑