Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	832,557	T→C	S251P (TCC→CCC)	hcxB →	hydroxycarboxylate dehydrogenase HcXB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	832,557	0	T	C	97.8%	148.2 / NA	46	S251P (TCC→CCC)	hcxB	hydroxycarboxylate dehydrogenase HcXB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/22);  minor base A (1/0);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGACGTTAACGGCGTGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGA  >  NZ_CP009273/832425‑832681                                                                                                                                     |                                                                                                                             aTTGACGTTAACGGCGTGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCcgg                                                                                                                          <  2:281774/137‑1 (MQ=255)                 tGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCAt                                                                                                          >  2:754932/1‑137 (MQ=255)                 tGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCa                                                                                                        >  2:382557/1‑139 (MQ=255)                 tGCCGACGACCAATCCGGCGGGAATGCAGGAGGCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTCCGCCCTTGCAGCGATGCGTGTAATTTTTTGCGGGGCGCTTACCGGCGGGACAACGACGcct                                                                                                          >  2:328773/1‑135 (MQ=255)                      acgacCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   <  1:382557/139‑1 (MQ=255)                           cAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTa                                                                                              <  2:843890/139‑1 (MQ=255)                                  gcggTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAAc                                                                                         >  1:145932/1‑137 (MQ=255)                                           cAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATg                                                                              <  1:620119/139‑1 (MQ=255)                                                gTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCAtt                                                                         >  2:102566/1‑139 (MQ=255)                                                   gCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   <  1:109066/110‑1 (MQ=255)                                                   gCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   >  2:109066/1‑110 (MQ=255)                                                   gCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCt                                                                       >  1:756563/1‑138 (MQ=255)                                                           gTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCAt                                                              <  1:328773/139‑1 (MQ=255)                                                           gTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCAt                                                              <  2:613044/139‑1 (MQ=255)                                                              cGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAACCCAGTCCCGATGCCATTCTTCACTGCATgcc                                                           >  2:836052/1‑137 (MQ=255)                                                               gCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGAc                                                           >  1:5324/1‑138 (MQ=255)                                                               gCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGAc                                                           <  2:5324/138‑1 (MQ=255)                                                                      aCCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatca                                                   >  1:68435/1‑139 (MQ=255)                                                                           tGCTGATCATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatcat                                               <  2:243308/138‑1 (MQ=255)                                                                           tGCCGTACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCCTCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatcat                                               >  1:137068/1‑138 (MQ=255)                                                                           tGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAacgac                                                                                                              >  1:542545/1‑75 (MQ=255)                                                                           tGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAacgac                                                                                                              <  2:542545/75‑1 (MQ=255)                                                                           tGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatcat                                               >  1:822552/1‑138 (MQ=255)                                                                                            aCGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCgg                             >  1:215032/1‑139 (MQ=255)                                                                                             cGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCa                                                         <  2:422221/110‑1 (MQ=255)                                                                                             cGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCa                                                         >  1:422221/1‑110 (MQ=255)                                                                                              gCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTgcggcg                           <  2:137068/139‑1 (MQ=255)                                                                                                ccTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAACCCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGgcgcg                         <  1:102566/139‑1 (MQ=255)                                                                                                   tGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCa                                                               <  1:468587/98‑1 (MQ=255)                                                                                                   tGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCa                                                               >  2:468587/1‑98 (MQ=255)                                                                                                          aTGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGt                >  1:23081/1‑138 (MQ=255)                                                                                                          aTGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGt                >  1:533328/1‑138 (MQ=255)                                                                                                          aTGTGTGAAATTCTTGGCGGGGCGCCTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGt                >  2:124308/1‑138 (MQ=255)                                                                                                             tgtgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAacga                                                                                                               >  1:29267/1‑40 (MQ=255)                                                                                                             tgtgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAacga                                                                                                               <  2:29267/40‑1 (MQ=255)                                                                                                             tgtgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCAt                                                              >  1:465322/1‑89 (MQ=255)                                                                                                             tgtgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCAt                                                              <  2:465322/89‑1 (MQ=255)                                                                                                               tgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgcg          <  2:68435/139‑1 (MQ=255)                                                                                                               tgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgcg          >  2:833300/1‑139 (MQ=255)                                                                                                                 aaaTTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACgcgc         <  1:833300/138‑1 (MQ=255)                                                                                                                         ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   <  2:482826/40‑1 (MQ=255)                                                                                                                         ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   >  1:482826/1‑40 (MQ=255)                                                                                                                         ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGa                                                                                <  2:4779/59‑1 (MQ=255)                                                                                                                         ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGa                                                                                >  1:4779/1‑59 (MQ=255)                                                                                                                         ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCg   <  1:434211/136‑1 (MQ=255)                                                                                                                         ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGa  <  1:754932/137‑1 (MQ=255)                                                                                                                                     |                                                                                                                             ATTGACGTTAACGGCGTGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGA  >  NZ_CP009273/832425‑832681

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