Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,703,693	C→T	A206T (GCG→ACG)	ydgJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,703,693	0	C	T	97.1%	104.7 / NA	34	A206T (GCG→ACG)	ydgJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/17);  total (16/18)
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.98e-01

TACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCA  >  NC_000913/1703596‑1703789                                                                                                  |                                                                                                 taCCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTc                                                                                                <  2:317434/100‑1 (MQ=255)    cGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCgg                                                                                             <  2:1295646/100‑1 (MQ=255)       gCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGcc                                                                                          <  1:387359/100‑1 (MQ=255)            cTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAAc                                                                                     >  1:584501/1‑100 (MQ=255)            cTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAAc                                                                                     >  2:661649/1‑100 (MQ=255)                       ccAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATc                                                                          >  2:338122/1‑100 (MQ=255)                                    gtgtAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTg                                                             >  2:1102199/1‑100 (MQ=255)                                       tAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAcc                                                          >  2:52684/1‑100 (MQ=255)                                       tAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAcc                                                          <  1:734358/100‑1 (MQ=255)                                       tAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAcc                                                          <  1:918240/100‑1 (MQ=255)                                          aaTGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGt                                                       >  1:833855/1‑100 (MQ=255)                                              aCTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc                                                   >  2:1078115/1‑100 (MQ=255)                                              aCTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAAc                                                   <  2:833849/100‑1 (MQ=255)                                                tCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCa                                                 >  2:618084/1‑100 (MQ=255)                                                 cgccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCaa                                                <  1:669438/100‑1 (MQ=255)                                                 cgccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCaa                                                >  2:175736/1‑100 (MQ=255)                                                 cgccgcTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGACATGCTGACCGGTAAACCaa                                                <  2:799455/100‑1 (MQ=255)                                                          ggATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTaa                                       >  2:436811/1‑100 (MQ=255)                                                          ggATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTa                                        >  2:284897/1‑99 (MQ=255)                                                               ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                  >  1:382879/1‑100 (MQ=255)                                                               ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                  >  2:939136/1‑100 (MQ=255)                                                               ggACAAGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGCCATGCTGACCGGTAAACCAAAAAGCGTAATGGcc                                  <  2:682840/100‑1 (MQ=255)                                                                 aCAAGATGGCGGGGAAATAAACGGTCGACTGCGCTCCGAACCGTAACTTTGCCAAATCTACCGTCATGCTAAAAAAAAAAAAAAAAAGCGTAATGGcc                                  <  2:1278041/98‑1 (MQ=255)                                                                    aGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATc                             <  1:338125/100‑1 (MQ=255)                                                                    aGATGGCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATc                             <  1:357915/100‑1 (MQ=255)                                                                        ggCGTGGAAATAATCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCaaga                         <  2:1196285/100‑1 (MQ=255)                                                                                aataatCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTAAAAAAAAAAAAAAAAAGCGTAATGGCCTGATCAAGAAGATgtgg                 <  2:674987/100‑1 (MQ=255)                                                                                 ataatCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAAAAAAAAAAAAAAAAGCGTAATGGCCTGATCAAGAAGATgtggt                <  2:116545/100‑1 (MQ=255)                                                                                   aatCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGc              <  2:817156/100‑1 (MQ=255)                                                                                    atCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAACAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCt             >  2:427934/1‑100 (MQ=255)                                                                                    atCGGTCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCt             <  2:185543/100‑1 (MQ=255)                                                                                         tCGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATc        >  2:1027033/1‑100 (MQ=255)                                                                                          cGACTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGAAAAAAAAAAAAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCg       <  2:1250201/100‑1 (MQ=255)                                                                                            aCTGCGTTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTa     >  2:1224894/1‑100 (MQ=255)                                                                                               gcgtTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCa  >  1:401669/1‑100 (MQ=255)                                                                                                  |                                                                                                 TACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTCCGGGCCGTAACTGTGCCAAATCTACCGTCATGCTGACCGGTAAACCAAAAAGCGTAATGGCCTGATCAAGAAGATGTGGTGCTAAATCGTACCA  >  NC_000913/1703596‑1703789

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