Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,863	A→G	L204L (CTT→CTC)	dgoK ←	2‑oxo‑3‑deoxygalactonate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,863	0	A	G	85.7%	52.3 / 1.7	21	L204L (CTT→CTC)	dgoK	2‑oxo‑3‑deoxygalactonate kinase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (10/8);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

GGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCA  >  NC_000913/3873730‑3873992                                                                                                                                      |                                                                                                                                  gggTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTg                                                                                                                              <  2:663119/139‑1 (MQ=255)            tgtTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATgg                                                                                                                   >  2:655124/1‑139 (MQ=255)            tgtTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATgg                                                                                                                   >  2:421763/1‑139 (MQ=255)                gATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATggcggg                                                                                                               <  1:655124/139‑1 (MQ=255)                                    aTACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCaa                                                                                           >  2:703505/1‑139 (MQ=255)                                         ggCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCa                                                                                      <  1:421763/139‑1 (MQ=255)                                          gCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAg                                                                                     >  2:257749/1‑139 (MQ=255)                                           cGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGc                                                                                    <  1:252156/139‑1 (MQ=255)                                                                   cAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATggcg                                                                                                                 <  1:254336/86‑1 (MQ=255)                                                                   cAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATggcg                                                                                                                 >  2:254336/1‑86 (MQ=255)                                                                    agagagAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTcc                                                           >  1:248116/1‑139 (MQ=255)                                                                      agagaAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTAtt                                                                                                         <  1:193462/91‑1 (MQ=255)                                                                      agagaAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTAtt                                                                                                         >  2:193462/1‑91 (MQ=255)                                                                                                       ttCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCt                                                                                   <  2:18501/80‑1 (MQ=255)                                                                                                       ttCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCt                                                                                   >  1:18501/1‑80 (MQ=255)                                                                                                            aGCACATGCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGAGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACa                   >  2:413240/1‑139 (MQ=255)                                                                                                               acaTGCGAGGCGCGAACTTC‑AAGAGCTGCGGCAATATTGCGGGCGTATTAAGACCACGCTCAAAGCCAACTGTGAAGGCATCGGCAGAGGTTTCCCGCGGCGGCAAACCTTCGCCCATCAATGAGTGATTTAACAATaa               >  1:196285/1‑139 (MQ=255)                                                                                                                   gagAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTgcgc                                        >  2:429172/3‑111 (MQ=255)                                                                                                                   gCGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAatg            >  1:646989/1‑139 (MQ=255)                                                                                                                   ccGAGGCGCGAACTTCAAAGAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTgcgc                                        <  1:429172/110‑1 (MQ=255)                                                                                                                             aaCTTCAAAGAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCa  <  1:700568/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                  GGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGGCGGCAAACCTGCGCCAATCAATGAGTGATTTAACAGTAAATGATGTAATTCA  >  NC_000913/3873730‑3873992

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