Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,000,353	C→T	G235E (GGG→GAG)	yigE ←	DUF2233 domain‑containing protein YigE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,000,353	0	C	T	95.8%	68.0 / NA	24	G235E (GGG→GAG)	yigE	DUF2233 domain‑containing protein YigE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/11);  minor base G (0/1);  total (12/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.96e-01

TGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGC  >  NC_000913/4000107‑4000594                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  tGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCACCGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGcc                                                                                                                                                                                                                                               >  1:37382/1‑251 (MQ=255)      ggtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCttt                                                                                                                                                                                                                                            <  2:76010/249‑1 (MQ=255)       gtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCat                                                                                                                                                                                                                                         <  2:199547/251‑1 (MQ=255)       gtggCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCat                                                                                                                                                                                                                                         >  1:86876/1‑251 (MQ=255)          gcggcATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCAtata                                                                                                                                                                                                                                      <  1:220159/251‑1 (MQ=255)                                                     ggTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACg                                                                                                                                                                                           >  1:67514/1‑251 (MQ=255)                                                           gcgcACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCa                                                                                                                                                                                       <  2:7614/249‑1 (MQ=255)                                                                       tGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGAGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGc                                                                                                                                                                          >  2:48534/1‑250 (MQ=255)                                                                                                         aaCGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCctgct                                                                                                                                       <  1:80198/251‑1 (MQ=255)                                                                                                           cGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCAACTGCGCCATATTCAGTATCTGACACAAAACTATCGTTTTAACCTTTACGCTCGCCGGAACTCATGGTGACAAAAGGATCGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAACTCGTGCCATCAAGATAAAGCAGATGCTCAACGTTCAGTTTCGCTTTGGCATAACAGTCTAAATCCTAAAAATTTGATGCCTGCTgg                                                                                                                                     >  2:130925/1‑251 (MQ=255)                                                                                                                            aGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGc                                                                                                                    <  1:21694/251‑1 (MQ=255)                                                                                                                            aGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCGCGACGGAAATCATGGTGACAAAAGGATAGCGATGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTAAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTGGCCTGCTGGCTCAACAAAAACACGGc                                                                                                                    <  2:40982/251‑1 (MQ=255)                                                                                                                             gggCGTGTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAAATATCGTTTTAACCTGTCCGCTCGTCGGAAATCATGGTGACAAAATGATAGCGTTGCCACGGGATTGCCGCGCCTTTCATATACATATGTGAACTCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGc                                                                                                                    <  1:130667/250‑1 (MQ=255)                                                                                                                                 gTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTcc                                                                                                               >  2:189596/1‑251 (MQ=255)                                                                                                                                               cGGCATCTGCGCCATCTTCAGTATCTGCCACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGAGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAACCACGGCGTTCCCATGTTTATTAATc                                                                                                 >  2:157235/1‑251 (MQ=255)                                                                                                                                                                 cAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAagcag                                                                                                                                                                                                     <  1:54707/133‑1 (MQ=255)                                                                                                                                                                 cAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAagcag                                                                                                                                                                                                     >  2:54707/1‑133 (MQ=255)                                                                                                                                                                         gACACAAAACTATCGTTTTAACCTTTCCGCTCGACGTAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCtt                                                                         <  1:74715/249‑1 (MQ=255)                                                                                                                                                                            acaAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCGACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGc                                                                    >  2:220273/1‑251 (MQ=255)                                                                                                                                                                               aaaCTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGAAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGAc                                                                 >  1:237272/1‑251 (MQ=255)                                                                                                                                                                                             cccTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAATTCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGGTTCGCTTTGGAATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGACTAc                                                   >  2:177473/2‑251 (MQ=255)                                                                                                                                                                                                 ttCCGCTCGACGGAAATCATGGTGACAAAAGGATGTCGGTGCCATGGTATTGCCTCGCCTGTCATATACATATTTGAAAGCGTGCCATCAGGATAAAGCAGCTGCTCAACGTTCATTTTCGCTTTGGCATAACAGGCAAAATCATAAAAAGTTGGGGCCGGCTGGCTCAACAAAAACACGGCGGTCCCATGTTTGTTAATTCCAACACCGTTACGAATTTTGCTTGATTCTAGGTTGGGATGAATAcgc                                                 <  2:244948/249‑1 (MQ=255)                                                                                                                                                                                                  tCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCTCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCgg                                               <  1:78402/250‑1 (MQ=255)                                                                                                                                                                                                                                              gggAATGCCTCGCCTTTCATATCCATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGCCGTTGGGATGAATACGCGGTTTAATTACACCGTTTTCCATCAACATTGGCCATGACTGAACCGc  >  2:201889/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  TGGCTGGTGGCGGCATACGCCCGGCTGGAGTCGGTGTAACGTTGCCGGATGCGGTGCTGCGCACCTTATTTGGCCTAAAAAATCATTCAGATTCAATAAATTGCAACGTCATGTAGGCCGGATAGGGCGTTTACGCCGCATCCGGCATCTGCGCCATCTTCAGTATCTGACACAAAACTATCGTTTTAACCTTTCCGCTCGACGGAAATCATGGTGACAAAAGGATAGCGTTGCCACGGGATTGCCCCGCCTTTCATATACATATGTGAAATCGTGCCATCAAGATAAAGCAGCTGCTCAACGTTCAGTTTCGCTTTGGCATAACAGGCAAAATCATAAAAATTTGTTGCCTGCTGGCTCAACAAAAACACGGCGTTCCCATGTTTATTAATCCCAACACCGTTACGAATTTTGCTTGAGGCGACGTTGGGATGAATACGCGGATTAATTACACCGTTTTCCATCAACATTGGCCCTGACTGCACCGC  >  NC_000913/4000107‑4000594

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