Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,626,105	G→A	E331K (GAA→AAA)	ppx →	exopolyphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,626,105	0	G	A	97.5%	133.4 / ‑4.0	40	E331K (GAA→AAA)	ppx	exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (19/20);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CCGTGAACTGCGCCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCGAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAAT  >  NC_000913/2625975‑2626239                                                                                                                                    |                                                                                                                                       ccGTGAACTGCGCCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGccc                                                                                                                               >  1:314271/1‑140 (MQ=255)    tGAACTGCGCCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGAc                                                                                                                            >  1:138486/1‑140 (MQ=255)       aCTGCGCCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACggg                                                                                                                         <  2:373639/140‑1 (MQ=255)         tGCGCCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTg                                                                                                                       <  2:285846/140‑1 (MQ=255)         tGCGCCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTg                                                                                                                       <  2:150647/140‑1 (MQ=255)         cgcgcCTTTCTGGCGGGGGGTTTCGCGAAGGGGTACTTTATGAAATGGAAGGACGTTTCCGTTATCCGGATGGGGGTAGGCCCCCCCCCCGCAGCCCCGGCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTg                                                                                                                       <  2:139642/139‑1 (MQ=255)           cgcCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCt                                                                                                                     >  1:99658/1‑140 (MQ=255)           cgcCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCt                                                                                                                     >  2:335488/1‑140 (MQ=255)           cgcCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCt                                                                                                                     >  2:36326/1‑140 (MQ=255)           cgcCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGAATC‑ACTTCGACACCAAACGGCCCCGACGGGGGcg                                                                                                                     >  1:139642/1‑139 (MQ=255)           cgcCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCG‑CAACCAGGATCAACTTCGCCA‑CAAACAGGGCCCACggggggtg                                                                                                                    >  2:156352/1‑135 (MQ=255)              ttttCTGATGGGGTGTTACGCGAAGGCGTACTTTATTAAATGAAAGCGCGTTGCCGTCATGAGGATGTGCGTAGTCGCGCCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGa                                                                                                                  <  2:11526/139‑1 (MQ=255)                         ggCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGc                                                                                                       >  2:390024/1‑140 (MQ=255)                          gCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCa                                                                                                      <  1:335488/140‑1 (MQ=255)                          gCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCa                                                                                                      >  2:213747/1‑140 (MQ=255)                                    aaGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGaa                                                                                            <  2:138486/140‑1 (MQ=255)                                    aaGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGaa                                                                                            >  1:121693/1‑140 (MQ=255)                                                    aaTGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGc                                                                            >  1:85648/1‑140 (MQ=255)                                                      tGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGGAAAGCaa                                                                          >  1:333222/1‑140 (MQ=255)                                                         aaGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGGACaaaaaaaa                                                                       >  2:221780/1‑133 (MQ=255)                                                         aaGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAAAAACaaacc                                                                       >  1:160830/1‑137 (MQ=255)                                                                                  gCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAgg                                              <  1:213747/140‑1 (MQ=255)                                                                                   cGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGGACAGCAACCGAAGCTGGAGAAACCGCAAATTGAggg                                             >  2:95782/1‑139 (MQ=255)                                                                                       gTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGGACAACAAACCAAACCGGCGCGCCCGCCACCGGGGgggcgc                                         >  2:382227/1‑138 (MQ=255)                                                                                       gTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTa                                         <  1:221780/140‑1 (MQ=255)                                                                                           cACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGGAAAGCAACCCAAACCGGGGCATCCCCCACCGGGGGGGGTACttc                                     >  1:17787/1‑138 (MQ=255)                                                                                           cACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGc                                     <  2:314271/140‑1 (MQ=255)                                                                                                 cagcagCCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATggg                               <  1:36326/140‑1 (MQ=255)                                                                                                    cagcCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGgccg                            <  2:121693/140‑1 (MQ=255)                                                                                                     agcCTCGCCAACCAGTATC‑ACATCGACAGCGAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGgccgc                           <  2:382427/140‑1 (MQ=255)                                                                                                     agcCTCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGgccgc                           <  2:14744/140‑1 (MQ=255)                                                                                                         tCGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCAtg                       <  1:241657/140‑1 (MQ=255)                                                                                                          cGCCAACCAGTATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGGAAAACAAACCAAAACGGGGCCTCCGCCACCGGAGGGGGCACTGCGGTGGggcgccccgc                      >  2:60024/1‑136 (MQ=255)                                                                                                              aaCCAGTATC‑ACATCGACAGCAAACAGGCCCGTCGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCa                  <  1:358066/140‑1 (MQ=255)                                                                                                                    tATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGt            <  1:156352/140‑1 (MQ=255)                                                                                                                    tATC‑ACATCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGt            <  2:85648/140‑1 (MQ=255)                                                                                                                          caTCGACAGCAAACAGGCCCGTCGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAAGAACAAACGAAGCCGGCGCATCCCCAACCGGGGGGGCCACTTCGGTGGGGCGCCCTGCCGCCTGAGGGCgggg       >  2:358066/1‑139 (MQ=255)                                                                                                                          caTCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGt       <  1:299105/140‑1 (MQ=255)                                                                                                                          caTCGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGt       <  2:99658/140‑1 (MQ=255)                                                                                                                             cGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGa    <  1:77948/140‑1 (MQ=255)                                                                                                                             cGACAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAAAAGCAAACCAAACCGGGGCACCCCCCACCGGGGGCGGCACTGCGATGGGGCGCCCCGCTGCGTGAGGGCggggggg    >  2:75052/1‑136 (MQ=255)                                                                                                                               aCAGCAAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGGAAAGCAAACCAAAACGGGGCATCCCCAACCGGGGGGGGCACTGCGAGGGGGCGCCCCGCTGCCTGGGGGCggggggggt  >  2:77948/1‑134 (MQ=255)                                                                                                                                    |                                                                                                                                       CCGTGAACTGCGCCTTTCTGACGGGGCGTTACGCGAAGGCGTACTGTATGAAATGGAAGGACGTTTCCGTCATCAGGATGTGCGTAGTCGCACCGCCAGCAGCCTCGCCAACCAGTATC‑ACATCGACAGCGAACAGGCCCGACGGGTGCTGGATACCACTATGCAAATGTACGAACAGTGGCGGGAACAGCAACCGAAGCTGGCGCATCCGCAACTGGAGGCGCTACTGCGATGGGCCGCCATGCTGCATGAGGTCGGGTTGAAT  >  NC_000913/2625975‑2626239

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