Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,767,014	C→T	G558S (GGT→AGT)	ydiJ ←	putative FAD‑linked oxidoreductase YdiJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,767,014	0	C	T	96.3%	86.7 / ‑6.2	27	G558S (GGT→AGT)	ydiJ	putative FAD‑linked oxidoreductase YdiJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/12);  minor base A (1/0);  total (15/12)
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.97e-01

AGGCAGTTCTTGTTCCAGTTTGAGTGGATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACCGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTT  >  NC_000913/1766778‑1767251                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              aGGCAGTTCTTGTTCCAGTTTGAGTGGATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGt                                                                                                                                                                                                                                 >  2:147509/1‑251 (MQ=255)    cAGTTCTTGTTCCAGTTTGAGTGGATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTccc                                                                                                                                                                                                                              >  2:131854/1‑250 (MQ=255)                    ttGAGTGGATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGgcg                                                                                                                                                                                                               <  1:94969/250‑1 (MQ=255)                                    cgcgTTCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCGCTTTCATCATCGGCGCATCGCTCCCTTCTGGc                                                                                                                                                                                              <  1:186396/251‑1 (MQ=255)                                    cgcgATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATCGGGATCTGCCGATCGAATGTACTGCGCTTCACGGCGTCCACTTTCATCATCGGCGCCTCTAGACCTTCTGGc                                                                                                                                                                                              >  2:71082/1‑251 (MQ=255)                                           cGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGCCCTTCTGGCGGGCaa                                                                                                                                                                                        >  2:229434/1‑250 (MQ=255)                                                      gCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTggg                                                                                                                                                                            >  2:104807/1‑251 (MQ=255)                                                            aTTCACGCACCAGCGTTGCTCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGt                                                                                                                                                                      <  2:27286/251‑1 (MQ=255)                                                            aTTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGt                                                                                                                                                                      <  1:147509/251‑1 (MQ=255)                                                              tCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACGCCAGCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGTCCTGCTGTCGGGCAAATCGTCCCTGGGGTGAGTc                                                                                                                                                                     <  2:146188/250‑1 (MQ=255)                                                                                 gcCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCTGGTCAAATGcc                                                                                                                                                 <  1:182918/251‑1 (MQ=255)                                                                                                      ggTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACttttcccccttct                                                                                                                            <  1:104807/251‑11 (MQ=255)                                                                                                              gTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCTCAGGTCTGCaaaa                                                                                                                      >  2:54178/1‑249 (MQ=255)                                                                                                               tGATCTTCATCGACGGACACATCGGACTAAGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAgg                                                                                                                   >  1:202551/1‑250 (MQ=255)                                                                                                                                      gtaCTACCGGCATCAAAGTTGAAGCATAATCCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGTATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTGTCATCATCGGCGCATCGAGATCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCGCTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCt                                                                                            <  2:63602/249‑1 (MQ=255)                                                                                                                                          tACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCCGGCTGTAt                                                                                        >  1:33123/1‑251 (MQ=255)                                                                                                                                                                 aaaCCGTTGCCGTTACACTCCCTCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGGTGAGTCGGGTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGGCAAAAGCTCCCCCCCGCAAAACCCCGGGCTGTAGTCCGCGCGAAAACCTTTGACGTg                                                                 >  2:197277/1‑251 (MQ=255)                                                                                                                                                                      gTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGc                                                            <  1:106823/251‑1 (MQ=255)                                                                                                                                                                        ggCCGTTACACTCCATCGCACCGCGCCACTGCCGGCGTACCGCAATGGGGATCTGCTGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCATATCTTCCCTGGGTTGAGTCGGTTGTGCGGGGCAAATGCCGCTTTCCCTTTGCGCAGTTCTGCAAAAAGTTCCCCACCGACAAACGCCCGGCTGTATTCAGCGCGAACACCTTTGCCGTGCTCGccc                                                          >  2:167430/2‑251 (MQ=255)                                                                                                                                                                          ccGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGCTCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATCCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGccccc                                                        >  1:114090/1‑250 (MQ=255)                                                                                                                                                                           cGTTACACTACATCGCATTCCGCCAATGCTCGCGTACCGCAATGGGGATCCATCGATCGAATGTACAGCGCTTCACCGCGTCCACTTTCATCATAGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGCTGCGTCGGTTATGCGGGTCACAGCCCGCCTTCCCATTTAGCAGTCCAGCCAAAATGTCCACACCGAAAAACCACGGGCTGTATTTACCGCGAAAACCTTTGCCggtcccgccacac                                                       >  2:155527/1‑239 (MQ=255)                                                                                                                                                                                          cACCGCGCCACTGCTGGCGTACCGCAATGGGGTTCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGGTATGCGGGTCAAATGCCGCTTTCCCTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAACACGCCGCGCTGTATTCAGCGCGAAAACCTTTTCCGTGCCCGCCCCCCAACAAACCACCGTAt                                        >  2:14730/1‑251 (MQ=255)                                                                                                                                                                                          cACCGCGCCACTGCTGGCGAACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGCTGAGTCGGTTATGCGGGTCAAATTCCGCTTTCCCTTTGCGCAGTTCTTCAAAAAACTTCTCACCGAAAAACGCCGGGATGTATTCAGCCCGCACACCTTTGCCGTGCTCGCCCCACAACAAACCAACGTAt                                        >  2:78483/1‑251 (MQ=255)                                                                                                                                                                                             cgcgCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTc                                     <  2:76019/251‑1 (MQ=255)                                                                                                                                                                                                                              gCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGAtt       <  1:219162/248‑1 (MQ=255)                                                                                                                                                                                                                              gCCGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTGCGCAGTCAGCGCCACCACGTCATCAGAGATTTGc    <  2:33123/251‑1 (MQ=255)                                                                                                                                                                                                                                cGATCGAATGTACTGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCtt  >  1:27318/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              AGGCAGTTCTTGTTCCAGTTTGAGTGGATCAACGCCGCGATCCGCCAACAAACGCAGCCATTCACGCACCAGCGTTGCGCGCCCTTTCGGTGAATGAATCCGGTTCTGGGTGATCTTCATCGACGGACACATCGGACTACGGGCATCAAAGTTGAAGCATAAACCGTTGCCGTTACACTCCATCGCACCGCGCCACTGCTGGCGTACCGCAATGGGGATCTGCCGATCGAATGTACCGCGCTTCACCGCGTCCACTTTCATCATCGGCGCATCGAGACCTTCTGGCGGGCAAATCTTCCCTGGGTTGAGTCGGTTATGCGGGTCAAATGCCGCTTTCACTTTGCGCAGTTCTGCAAAAAGTTCCTCACCGAAAAACGCCGGGCTGTATTCAGCGCGAAAACCTTTGCCGTGCTCGCCCCACAACAAACCACCGTATTTCGCAGTCAGCGCCACCACGTCATCAGAGATTTGCTT  >  NC_000913/1766778‑1767251

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