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

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641                                                  |                                              tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:1258224/50‑1 (MQ=255)  aaaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               >  1:2691545/1‑49 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             <  1:5402087/50‑3 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:1667365/1‑48 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  2:3176887/1‑48 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  2:1819570/1‑48 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:2465577/1‑48 (MQ=255)      gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                         <  2:4605776/50‑1 (MQ=255)      gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                         <  1:2852691/50‑1 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        <  1:1747411/50‑1 (MQ=255)        gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTccc                                       <  1:2717303/50‑1 (MQ=255)        gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  2:3054544/1‑49 (MQ=39)        gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  2:1089855/1‑49 (MQ=39)          gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                     >  1:4161919/1‑50 (MQ=255)          gAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                     >  1:5485665/1‑50 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                     <  2:3385237/49‑1 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    >  2:1040392/1‑50 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    <  1:1593827/50‑1 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    >  1:5293690/1‑50 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    >  1:4620395/1‑50 (MQ=255)            aaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcc                                   >  1:3468446/1‑50 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  >  1:4578384/1‑50 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:1009622/1‑50 (MQ=255)               aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                                >  2:3152568/1‑50 (MQ=255)               aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                                >  2:4666596/1‑50 (MQ=255)                aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                               <  1:2066686/50‑1 (MQ=255)                aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                               <  1:3639039/50‑1 (MQ=255)                aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                               >  1:5105276/1‑50 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              >  2:5042603/1‑50 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              >  1:3231228/1‑50 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              >  2:5371851/1‑50 (MQ=255)                  tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                             >  1:1553975/1‑50 (MQ=255)                  tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                             >  1:4930766/1‑50 (MQ=255)                   cACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            >  2:3388543/1‑50 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            <  2:4844244/49‑1 (MQ=255)                     cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  2:3098316/1‑50 (MQ=255)                     cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  2:4721219/1‑50 (MQ=255)                     cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          <  1:5023627/50‑1 (MQ=255)                      tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  2:5547486/1‑49 (MQ=255)                      tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  2:5298730/1‑49 (MQ=255)                      tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  2:724842/1‑49 (MQ=255)                      tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  2:5519301/1‑49 (MQ=255)                         tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      >  2:245188/1‑50 (MQ=255)                          caagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     >  2:1487939/2‑50 (MQ=255)                          aaagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGCGTTGAttt                     >  2:3556250/1‑50 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  2:4764452/1‑50 (MQ=255)                            agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:4585669/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTcc                   >  1:2039120/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:5343796/50‑1 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  2:1314114/1‑50 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  2:1393893/1‑50 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 <  2:438347/50‑1 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:206097/1‑50 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:4372958/50‑1 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  2:5773875/50‑1 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:457397/1‑50 (MQ=255)                                agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:2011834/1‑49 (MQ=255)                                agTAATCTATATGTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACg               <  1:5187673/50‑1 (MQ=255)                                 ggaaTTTATATTTCTGC‑GCATTCCCAGACTCCTGTGTTGGTTTCCAAc                >  2:4903129/3‑48 (MQ=25)                                 gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:1616380/1‑48 (MQ=255)                                 gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  2:4256524/1‑48 (MQ=255)                                 gTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGa              <  2:5664260/50‑1 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  2:1549734/50‑1 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:2363033/1‑50 (MQ=255)                                   aaTCTATTTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGTTTTCCAACGAg             >  2:659537/1‑49 (MQ=255)                                   aaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  2:5200346/1‑49 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt            <  2:1993239/49‑1 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  2:3039809/1‑50 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:4515814/1‑50 (MQ=255)                                     tCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  2:960232/1‑49 (MQ=255)                                     tCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           >  1:2230507/1‑49 (MQ=255)                                       tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:3242925‑M1/1‑48 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  2:3818736‑M1/1‑47 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  2:4726346/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  2:5633162/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      <  1:288561‑M1/50‑5 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      >  1:4030698‑M1/1‑46 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      <  2:4175565‑M1/50‑5 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      <  2:5351726‑M1/50‑5 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      >  2:1536551‑M1/1‑46 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:5842334‑M1/48‑3 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  2:5571449‑M1/47‑2 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  2:2750882‑M1/47‑2 (MQ=255)                                          atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  2:1140786‑M1/1‑45 (MQ=255)                                          atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  1:659442‑M1/1‑45 (MQ=255)                                          atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  2:1367168‑M1/1‑45 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:2831818‑M1/1‑44 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  2:5356687‑M1/1‑44 (MQ=255)                                            ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:1203250‑M1/1‑43 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  2:5293002‑M1/1‑42 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  2:3520113‑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 < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑