Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,858,058	A→G	M259T (ATG→ACG)	ydjJ ←	putative zinc‑binding dehydrogenase YdjJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,858,058	0	A	G	96.3%	93.4 / ‑5.8	27	M259T (ATG→ACG)	ydjJ	putative zinc‑binding dehydrogenase YdjJ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (13/13);  total (14/13)
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

CGCGTTTGTTGTTAACTGACTCTTCAAATGCCTGTTGTACATCCCGATAATCGTAAATATGCGTCACCATCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCATTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCTGTTCCGCCATTGCCAGACGTTTTTCCAGCACATCAACGACGGCAATTTCCGTTGCTCCCAGGCATTTGCACGCTTGCAACGTCATCAAACCA  >  NC_000913/1857825‑1858294                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             tgcgTTTGTTGTTAACTGTGTCTTCAAGTGCCTGTTGTACATCCCGATAATCGTACATATGCGTCACCATCGATTTCACATCGAATCGCCCTGAGGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCATTACCAGCTAAGGTGcc                                                                                                                                                                                                                             >  2:370350/2‑251 (MQ=255)      ttgttgTTAACTGACTCTTCAAATGCCTGTTGTACATCCCGATAATCGTAAATATGCGTCACCATCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGttt                                                                                                                                                                                                                        >  1:250020/1‑251 (MQ=255)                         cAAATGCCTGTTGTACATCCCGATAATCGTAAATATGCGTCACCATCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAAccc                                                                                                                                                                                                     >  1:231918/1‑251 (MQ=255)                          aaaTTCCTGTTGTACATCCCGATAATCGTAAATATGCGTCAGCATCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCTGAACCCg                                                                                                                                                                                                    <  2:282859/251‑1 (MQ=255)                                                 aaTCGTAAATATGCGTCACCATCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTgc                                                                                                                                                                              <  1:116237/250‑1 (MQ=255)                                                 aaTCGTAAATATGCGTCACCATCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTgc                                                                                                                                                                              >  2:116237/1‑250 (MQ=255)                                                           atGCGTCACCATCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGGTTCTAAAACAATATCTGCGCCTAGGTCtt                                                                                                                                                                   <  2:172686/251‑1 (MQ=255)                                                                     aTCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTg                                                                                                                                                         >  2:117840/1‑251 (MQ=255)                                                                           ttCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGAc                                                                                                                                                    >  1:82267/1‑250 (MQ=255)                                                                              acaTCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAAGATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCg                                                                                                                                                <  2:250020/251‑1 (MQ=255)                                                                                  cGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGc                                                                                                                                             >  1:112699/1‑250 (MQ=255)                                                                                  cGAATCGCCCTGAAGAAATAGCTTCAAGCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCa                                                                                                                                            <  1:219155/251‑1 (MQ=255)                                                                                                  aaTAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTg                                                                                                                            <  1:254428/251‑1 (MQ=255)                                                                                                     aGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTgcgc                                                                                                                         >  2:222282/1‑251 (MQ=255)                                                                                                                tGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAAcc                                                                                                                                                                                                      >  1:312175/1‑163 (MQ=255)                                                                                                                tGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAAcc                                                                                                                                                                                                      <  2:312175/163‑1 (MQ=255)                                                                                                                             aTTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCaa                                                                                                 <  2:112699/251‑1 (MQ=255)                                                                                                                                ggCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCt                                                                                              <  2:231918/251‑1 (MQ=255)                                                                                                                                                    ggATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGTTACAGTACCAACAATCATAATTTTACCGCCGCTCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCTGTTCCGCCATTGCCAGACg                                                                           <  2:141204/250‑1 (MQ=255)                                                                                                                                                                                          gCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGAGTTACCAGATAAGGTGCCTGTTTCACGGTGACCGCAGAACCCGCTGtt                                                                                                                                                                                               >  2:205670/1‑96 (MQ=255)                                                                                                                                                                                                    cGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCTGTTCCGCCATTGCCAGACGTTTTTCCAGCACATCAACGACGGCAATTTCCGTTGCTCCCAGGCATTTg                          >  1:46834/1‑251 (MQ=255)                                                                                                                                                                                                         aCAGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCTGATCCGCCATTGCCAGACGTTTTTCCAGCACATCAACGACGGCAATTTCCGTTGCTCCCAGGCATTTGCACGc                     >  1:358162/1‑251 (MQ=255)                                                                                                                                                                                                           aGTACCAACAATCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCTGTTCCGCCATTGCCAGACGTTTTTCCAGCACATCAACGACGGCAATTTCCGTTGCTCCCAGGCATTTGCACGCtt                   <  1:117840/251‑1 (MQ=255)                                                                                                                                                                                                            gTACCAACAATCATAAGTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGGCAGCGTGCAATAGTGTCGTCTTTTGCGCCGTTAATAACCACTGTCGCCCCAAGCTGTTCCGCCATTGCCAGACGTTTTTCCAGCCCAGCAACGACGCCAATTTCCGTGGCTCCCAGGCATTTGCACGCTTg                  <  1:370350/251‑1 (MQ=255)                                                                                                                                                                                                                    aaTCATAATTTTACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCTGTTCCGCCATTGCCAGACGTTTTTCCAGCACATCAACGACGGCAATTTCCGTTGCTCCCAGGCATTTGCACGCTTGCAACGtca          <  2:99537/251‑1 (MQ=255)                                                                                                                                                                                                                            ttttACCGCCGCGCGTTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGCACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCTGTTCCGCCATTGCCAGACGTTTTTCCAGCACATCCACGACTGCAATTTCCGTTGCTCCCAGGCATTTGCACGCTTGCAACGTCATCAAAcc   >  2:349760/1‑250 (MQ=255)                                                                                                                                                                                                                            ttttACCGCCGCGCGTTACCAGATAAGGTGACTGTTTAACGGAGACCGCAGAACCCGATGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCCAGCTGATCCGCCATTGCCAGACGTTTTTCCAGCACATCAACGACGGCAATTTCCGTTGCTCCCAGGCATTTGCACGCTTGCAACGTCATCAAAccc  >  1:194643/1‑250 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             CGCGTTTGTTGTTAACTGACTCTTCAAATGCCTGTTGTACATCCCGATAATCGTAAATATGCGTCACCATCGATTTCACATCGAATCGCCCTGAAGAAATAGCTTCAATCGTGACCGGATAACGATTGGCATAGCGGAATACCGTCTGGATAGTGACTTCGCGATTGATTTTGAGGAAATTGATTGCCGAATCGCCGGGTACAGTACCAACAATCATAATTTTACCGCCGCGCATTACCAGATAAGGTGCCTGTTTAACGGTGACCGCAGAACCCGCTGTTTCGAAAACAATATCTGCGCCCATGTCTTCGGTAAATTGCTGACAGCGTGCAATAGTGTCTTCTTTTGCGCCGTTAATAACCACTGTCGCACCAAGCTGTTCCGCCATTGCCAGACGTTTTTCCAGCACATCAACGACGGCAATTTCCGTTGCTCCCAGGCATTTGCACGCTTGCAACGTCATCAAACCA  >  NC_000913/1857825‑1858294

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