Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	387,798	G→C	10.8%	G313A (GGC→GCC)	glyA‑I →	serine hydroxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	387,798	0	G	C	10.8%	107.6 / 4.0	37	G313A (GGC→GCC)	glyA‑I	serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (18/15);  new base C (2/2);  total (20/17)
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.80e-01

ATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCTG  >  NC_002947/387656‑387946                                                                                                                                               |                                                                                                                                                     aTGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGa                                                                                                                                                >  4:119710/1‑149 (MQ=255)     aCGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCACGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGGTGCCGGTGGCACCgcgagc                                                                                                                                            >  8:230037/1‑144 (MQ=255)             gCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGtt                                                                                                                                    <  8:183179/149‑1 (MQ=255)                gCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacat                                                                                                                                       >  8:72888/1‑141 (MQ=37)                gCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacat                                                                                                                                       <  7:72888/143‑3 (MQ=255)                 ccAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGCTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTg                                                                                                                                <  1:264513/149‑1 (MQ=255)                     ggCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCa                                                                                                                            <  4:85713/149‑1 (MQ=255)                     ggCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCa                                                                                                                            <  6:79331/149‑1 (MQ=255)                        aGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGcc                                                                                                                         <  3:119710/149‑1 (MQ=255)                               ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCg                                                                                                                  >  3:236182/1‑149 (MQ=255)                               ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacat                                                                                                                                       <  5:136301/128‑3 (MQ=255)                               ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacat                                                                                                                                       >  6:136301/1‑126 (MQ=37)                                     gAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAggg                                                                                                            >  8:196041/1‑149 (MQ=255)                                        gcgcTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTg                                                                                                                       <  2:50164/135‑1 (MQ=255)                                        gcgcTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTg                                                                                                                       >  1:50164/1‑135 (MQ=255)                                                  cTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCaaa                                                                                               <  5:276430/149‑1 (MQ=255)                                                  cTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCa                                                                                                 >  8:1697/1‑147 (MQ=255)                                                   tGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGCCCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGcc                                                                                                          >  6:258985/1‑137 (MQ=255)                                                   tGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGcc                                                                                                          <  5:258985/137‑1 (MQ=255)                                                       ttCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgc                                                                                          >  3:116405/1‑149 (MQ=255)                                                                        gcaAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTc                                                                         <  7:196041/149‑1 (MQ=255)                                                                        gcaAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTc                                                                         >  5:58627/1‑149 (MQ=255)                                                                            gTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgc                                                                     <  6:58627/149‑1 (MQ=255)                                                                               aTCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgcgc                                                                   <  6:177984/148‑1 (MQ=255)                                                                                  gAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGa                                                                                                       >  1:170882/1‑109 (MQ=255)                                                                                  gAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGa                                                                                                       <  2:170882/109‑1 (MQ=255)                                                                                   aaaaCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGAt                                                                                                                     <  2:170379/94‑1 (MQ=255)                                                                                   aaaaCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGAt                                                                                                                     >  1:170379/1‑94 (MQ=255)                                                                                           cAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaa                                                      <  4:116405/149‑1 (MQ=255)                                                                                                 aTGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGaa                                                >  8:14408/1‑149 (MQ=255)                                                                                                      ccAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCg                                           >  8:55410/1‑149 (MQ=255)                                                                                                      ccAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTCGGTCGTGCGCATATCACCGTCAACAAGAACGCCg                                           >  4:213028/1‑149 (MQ=255)                                                                                                               cATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAAc                                   >  2:35504/1‑148 (MQ=255)                                                                                                                                    ggTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCa             >  4:150083/1‑149 (MQ=255)                                                                                                                                        tCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTc         >  5:106617/1‑149 (MQ=255)                                                                                                                                          cGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCg                                      <  4:230916/118‑1 (MQ=255)                                                                                                                                          cGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCg                                      >  3:230916/1‑118 (MQ=255)                                                                                                                                               gCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCTg  <  7:1697/149‑1 (MQ=255)                                                                                                                                               gCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCTg  <  3:150083/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     ATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGGGCCTG  >  NC_002947/387656‑387946

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