Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,537,549	C→A	R70S (CGC→AGC)	yfeK →	DUF5329 domain‑containing protein YfeK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,537,549	0	C	A	97.1%	108.3 / NA	34	R70S (CGC→AGC)	yfeK	DUF5329 domain‑containing protein YfeK
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (24/9);  total (25/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 = 9.84e-01

CAAAACAAGGATAAGGTAATTCAATGAAGAAAATCATTTGTCTGGTCATTACACTACTAATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAAT  >  NC_000913/2537319‑2537782                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          caAAACAAGGATAAGGTAATTCAATGAAGAAAATCATTTGTCTGGTCATTACACTACTAATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGc                                                                                                                                                                                                                         >  2:137606/1‑249 (MQ=255)       aaGGATAAGGTAATTCAATGAAGAAAATCATTTGTCTGGTCATTACACTACTAATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTAtt                                                                                                                                                                                                                 >  2:195130/1‑251 (MQ=255)               ggTAATTCAATGAAGAAAATCATTTGTCTGGTCATTACACTACTAATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGt                                                                                                                                                                                                         >  2:74730/1‑251 (MQ=255)                  aaTTCAATGAAGAAAATCATTTGTCTGGTCATTACACTACTAATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGc                                                                                                                                                                                                      >  1:164123/1‑251 (MQ=255)                              aaaaTCCTTTGTCTGGTCATTACACTACTACTGACACTCCCCGTTTACGCGAAGCTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCCGTTTATTGATAAGGTTGCTTCGTCGTCATc                                                                                                                                                                                          >  2:363467/1‑251 (MQ=255)                                           tGGTCATTACACTACTAATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGa                                                                                                                                                                               >  1:321202/1‑249 (MQ=255)                                            ggTCATTACACTACTAATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGaa                                                                                                                                                                              <  1:236972/249‑1 (MQ=255)                                                aTTACACTACTAATGACACTCCATGTTTAAGCGAAGTTAAATGCCAATGCAGAAGACCGCCTCAACGCACTGATGGAGGACTAAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCCGTTTATTGATAAGGTTGCTTCGTCGTAATCGATTACTGGGCAGccat                                                                                                                                                                         >  2:390292/1‑248 (MQ=255)                                                  tACACTACTAATGACACTCCCCGTTTACGCGCCGTTAACTGACCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGCAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGtat                                                                                                                                                                       >  2:57049/1‑250 (MQ=255)                                                                 acTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCg                                                                                                                                                       <  1:235599/251‑1 (MQ=255)                                                                 acTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCg                                                                                                                                                       >  2:313644/1‑251 (MQ=255)                                                                       cGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTAGATTGTGAAGATCCCCGTTAaga                                                                                                                                                 <  1:170368/251‑1 (MQ=255)                                                                                          gCCCATGAAGAAGCCCGCATCAACGACATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGcc                                                                                                                              >  2:186668/1‑251 (MQ=255)                                                                                          gCCCATGAAGAAGCCCGCATAAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAAAATACCTGATATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCCTTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAACGATGAGAACGCACAGcc                                                                                                                              >  1:273915/1‑251 (MQ=255)                                                                                                    aagCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACAt                                                                                                                    >  2:215118/1‑251 (MQ=255)                                                                                                         cGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGtt                                                                                                               >  1:124116/1‑251 (MQ=255)                                                                                                            ctCAACGCCATGCTGGAGGGCTTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAAt                                                                                                            >  2:172211/2‑251 (MQ=255)                                                                                                                      tGCTGGAGGGATTAGCACAGAACAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGCTTACTGGGCAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATTCGTTAATTGCGCAGACg                                                                                                  >  2:422177/1‑251 (MQ=255)                                                                                                                          ggAGGGATTAGCACAGAAAAAGGCTTTGATATTTGTGCGCAACGGTGATGAGCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGAt                                                                                               >  1:154133/1‑250 (MQ=255)                                                                                                                                                         tttGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTa                                                                     >  1:289871/1‑245 (MQ=255)                                                                                                                                                         tttGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTa                                                                     <  2:289871/245‑1 (MQ=255)                                                                                                                                                          ttGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGc                                                               <  1:74730/250‑1 (MQ=255)                                                                                                                                                            gtgcgaacgggTGATGCACATACCTGCTATGAAGCGCTTTCTCATCTGCGTCTGCAGATCGGCAATACCCGAAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCATCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGCGCGATGCGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGACATTAATCCGCttt                                                            >  2:339610/10‑251 (MQ=255)                                                                                                                                                            gtgCGCAACGGTGATGAACATCCCTGCTATGCAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCACTACCCGTAACCGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGACGCCGTATATTGTGAAGATCCACGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAACCGGTGACTGCGGAAGGAAATTAATCCGCttt                                                            >  2:471159/1‑251 (MQ=255)                                                                                                                                                                     ggTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGc                                                    >  2:473666/1‑250 (MQ=255)                                                                                                                                                                        gatgaACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCAtt                                                 <  1:538848/250‑1 (MQ=255)                                                                                                                                                                             aCATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGTATTTACagg                                           <  1:456964/251‑1 (MQ=255)                                                                                                                                                                                        aTGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACTCTGCCGAGCAGTTTATTGATAAGGTTGCGTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCGTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATg                                 <  1:97596/250‑1 (MQ=255)                                                                                                                                                                                             gCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGaaaaa                            >  2:168039/1‑250 (MQ=255)                                                                                                                                                                                              cGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACATCATGGACACTGCCGAGCAGTTTATGGATAAGGTTGCGTCGGCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAAGTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAAc                          <  1:471159/251‑1 (MQ=255)                                                                                                                                                                                                 tttCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCt                       >  1:176885/1‑251 (MQ=255)                                                                                                                                                                                                             gacTGAAGCTCGGCAATACCCGTAACAGCATTGACACAGCAGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGCTTACTGGGCAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCGGACGGATAACACGGGGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGCGGTAACATGAAAAAACGCTTTATTTATCACg           >  2:398191/3‑251 (MQ=255)                                                                                                                                                                                                                  aaGCTCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGTATACCACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGGAACATGAAAAAACGCTTTATTTATCACGGTGa       >  1:546900/1‑250 (MQ=255)                                                                                                                                                                                                                      tCGGCAATACCCGTAACAGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAAt  >  1:132242/1‑251 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          CAAAACAAGGATAAGGTAATTCAATGAAGAAAATCATTTGTCTGGTCATTACACTACTAATGACACTCCCCGTTTACGCGAAGTTAACTGCCCATGAAGAAGCCCGCATCAACGCCATGCTGGAGGGATTAGCACAGAAAAAGGATTTGATATTTGTGCGCAACGGTGATGAACATACCTGCTATGAAGCGGTTTCTCATCTGCGTCTGAAGCTCGGCAATACCCGTAACCGCATTGACACTGCCGAGCAGTTTATTGATAAGGTTGCTTCGTCGTCATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAAT  >  NC_000913/2537319‑2537782

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