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

TAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCA  >  NC_000913/3720547‑3720641                                                  |                                              tAAAAGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCg                                               <  1:6604045/50‑1 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             <  1:7786951/50‑3 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             >  1:2985324/1‑48 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgca                                             <  1:2924864/50‑3 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgc                                              <  1:554020/49‑2 (MQ=255)   aaaGCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGCgc                                              <  1:2676089/49‑2 (MQ=255)      gCGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTc                                         <  1:4649875/50‑1 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        <  1:2830978/50‑1 (MQ=255)       cGTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:2528173/1‑50 (MQ=255)        gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:4621916/1‑49 (MQ=39)        gTGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:8051852/1‑49 (MQ=39)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTcc                                        >  1:5598347/1‑48 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                      <  1:7562791/50‑1 (MQ=255)         tGAAACAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCa                                      <  1:2873011/50‑1 (MQ=255)           aaaCAAATCACTGTTAAAGAAAGTAATCTATATGTCTGC‑GCATTCCCAg                                     <  1:6362064/49‑1 (MQ=37)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                     <  1:1522480/49‑1 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAg                                     <  1:4675490/49‑1 (MQ=255)           aaaCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGc                                    >  1:6731757/1‑50 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  >  1:6793385/1‑50 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  <  1:4448556/50‑1 (MQ=255)             aCAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcct                                  <  1:6738012/50‑1 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 <  1:4802659/50‑1 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 <  1:7785469/50‑1 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 >  1:1799227/1‑50 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 <  1:7267086/50‑1 (MQ=255)              cAAATCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctc                                 <  1:2460787/50‑1 (MQ=255)               aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                                >  1:2871130/1‑50 (MQ=255)               aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                                >  1:3172629/1‑50 (MQ=255)               aaaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcc                                >  1:523059/1‑50 (MQ=255)                aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                               >  1:5457009/1‑50 (MQ=255)                aaTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGcctcct                               >  1:5662327/1‑50 (MQ=255)                 aTCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtg                              <  1:2700659/50‑1 (MQ=255)                  tCACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgt                             >  1:3176310/1‑50 (MQ=255)                   cACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            >  1:5386293/1‑50 (MQ=255)                   cACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            >  1:6684172/1‑50 (MQ=255)                   cACTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            >  1:2599194/1‑50 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            <  1:7705621/49‑1 (MQ=255)                    aCTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCtgtg                            <  1:2355479/49‑1 (MQ=255)                     cTATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          <  1:2253902/50‑1 (MQ=255)                      tATTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGtt                          >  1:7047465/1‑49 (MQ=255)                       aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                        >  1:3530948/1‑50 (MQ=255)                       aTTAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGa                        >  1:3419730/1‑50 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTg                         <  1:6490810/48‑1 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                       <  1:5422879/50‑1 (MQ=255)                        ttAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAt                       <  1:5028487/50‑1 (MQ=255)                         tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      >  1:8582022/1‑50 (MQ=255)                         tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      >  1:2013284/1‑50 (MQ=255)                         tAAAGAAAGTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAtt                      >  1:8226574/1‑50 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGAttt                     <  1:7572219/49‑1 (MQ=255)                           aagaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:6859/1‑50 (MQ=255)                            agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:7514162/1‑49 (MQ=255)                            agaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTc                    >  1:4193849/1‑49 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:5127369/50‑1 (MQ=255)                             gaaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCa                  <  1:7954363/50‑1 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:1691718/1‑50 (MQ=255)                              aaagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCaa                 >  1:3911192/1‑50 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:4622971/1‑50 (MQ=255)                               aagTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                <  1:6261013/50‑1 (MQ=255)                                agTAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:3236848/1‑49 (MQ=255)                                 gcaaTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAAc                >  1:4965267/3‑48 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:6924276/50‑1 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             <  1:3125863/50‑1 (MQ=255)                                  tAATCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAg             >  1:680013/1‑50 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGt            <  1:8541530/49‑1 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           <  1:1301172/50‑1 (MQ=255)                                    aTCTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTa           <  1:8030113/50‑1 (MQ=255)                                      cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          >  1:3459691/1‑49 (MQ=255)                                      cTATATTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          >  1:4429275/1‑49 (MQ=255)                                       tataTTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        >  1:4737054‑M1/1‑48 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:7855803‑M1/1‑47 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:6468385‑M1/1‑47 (MQ=255)                                        atatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       >  1:7054843‑M1/1‑47 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCGTCCTGTGTTGATTTCCAACGAGTAt          <  1:2280143/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  1:7747875/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  1:1944817/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTAt          <  1:1149709/46‑1 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaact      <  1:4955972‑M1/50‑5 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaac       <  1:1192730‑M1/49‑4 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:2601150‑M1/48‑3 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:2373674‑M1/48‑3 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaa        <  1:2213429‑M1/48‑3 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:5067334‑M1/47‑2 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:2572710‑M1/47‑2 (MQ=255)                                         tatTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATa         <  1:961739‑M1/47‑2 (MQ=255)                                          atTTCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     >  1:3301301‑M1/1‑45 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCTACGAGTATaacttt    >  1:3902652‑M1/1‑44 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    <  1:2043832‑M1/50‑7 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:1488948‑M1/1‑44 (MQ=255)                                           tttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactt     <  1:4551279‑M1/49‑6 (MQ=255)                                            ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:7384078‑M1/1‑43 (MQ=255)                                            ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:7519795‑M1/1‑43 (MQ=255)                                            ttCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaacttt    >  1:217677‑M1/1‑43 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttga  >  1:4225234‑M1/1‑42 (MQ=255)                                             tCTGC‑GCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATaactttg   <  1:6926663‑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 < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑