Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,032,423	G→A	V219V (GTC→GTT)	yedJ ←	putative HD superfamily phosphohydrolase YedJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,032,423	0	G	A	97.0%	106.4 / ‑3.6	33	V219V (GTC→GTT)	yedJ	putative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (20/12);  total (21/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 = 4.00e-01

CCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCA  >  NC_000913/2032196‑2032657                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           ccTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAAGTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACAGAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTGGTGAGCAACACCTTCATTCTCCCCCGCCAg                                                                                                                                                                                                                     <  2:147667/251‑1 (MQ=255)    acacCATTAAGCTGCGCCACCAACGTTTTAACATCAAAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAAGCGG‑TTACTGATATATTTTCCTTCAGAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCTGTTATGGCCACGCTCTCAGCCAGCGGATGAAAACGCCGCTACCCCCTTGTGACTAACACTTTCATCCTCCCCCGCCAGTTc                                                                                                                                                                                                                  <  2:345704/251‑1 (MQ=255)                   gcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGc                                                                                                                                                                                                   >  2:91128/1‑251 (MQ=255)                   gcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGc                                                                                                                                                                                                   >  1:329458/1‑251 (MQ=255)                       ccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTCTCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCAt                                                                                                                                                                                                >  2:104488/1‑250 (MQ=255)                        caACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATaa                                                                                                                                                                                              <  1:321602/251‑1 (MQ=255)                                         aaaTTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACTCCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATgc                                                                                                                                                                             <  1:111022/251‑1 (MQ=255)                                              tGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGt                                                                                                                                                                         <  1:289397/250‑1 (MQ=255)                                               gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACCCCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTtg                                                                                                                                                                       <  1:309364/251‑1 (MQ=255)                                                       gtctcCAGCATTTCCTGTCCGGCATGCCCGGTGCTGTATGAATCTG‑TTACTGATATATTTTCCGGAATAAATTTCACCGGCCATCGGATAAGGCCGACAGGTTAACATAATTTCGGCCATTGAGGGCCTCGCCCTCGCCCAAGGGATGAAAACTCGTCTGTCCCCTTGTGAGCAACCCCGTCAAGCCCCCCCGCCAGTTCTTCATGGAGCTTCTCCATAAACTCGACTAAAAAATGCGCGTTGTGCTgcg                                                                                                                                                                <  1:333993/249‑1 (MQ=255)                                                                 ttGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTg                                                                                                                                                     >  1:317770/1‑251 (MQ=255)                                                                   gCCTGTTCGTCATTCCCGGTGCTGTA‑GGATCGGGTTACAGTTATATTTGCCTTAATAGATTTCCCCGGCCACCTGAGATGGCCGACGGGTTAACATAATTTCAGCCAGTTAGTGCCACGCGCCCAGCCAGCGGATGAAAACGTATCCATCACGTTGGGCGCAACACCTTCATTCCCCCCCGCCAGTTCGGCGCTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGc                                                                                                                                                    <  2:304701/250‑1 (MQ=255)                                                                                 cccGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATAGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGc                                                                                                                                      >  2:187915/1‑250 (MQ=255)                                                                                    ggTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCCTgg                                                                                                                                  >  2:289912/1‑251 (MQ=255)                                                                                              gAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACCCCTGCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAg                                                                                                                        <  1:172585/251‑1 (MQ=255)                                                                                                        acggATATATTTTCCTGCATAGATTCCACCGGCCATCTGAGAGGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCGCAGCCAGCGGATGAAAACGCATCTATCCCCTGGTGATCAACACCTCCATTCTCCCCCGCCAGTTCGGCACTGAGCGTCGCCATAAACTCGACTAAAAAATGCGCGTTGtgc                                                                                                                                                                    <  2:147496/195‑1 (MQ=255)                                                                                                        aCTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCCTTTCGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTCAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCa                                                                                                                >  1:220174/1‑250 (MQ=255)                                                                                                        aCTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGtgc                                                                                                                                                                    >  1:147496/1‑198 (MQ=255)                                                                                                                atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTc                                                                                                        >  1:302251/1‑250 (MQ=255)                                                                                                                atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCt                                                                                                       >  1:48820/1‑251 (MQ=255)                                                                                                                          aTAGATGTCACCGGCCATCTGGGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGATCAACCCCTCCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATc                                                                                             <  2:38563/251‑1 (MQ=255)                                                                                                                                   aCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAACCATTTAGGGCAACGCGCTAAGCCAGCGGATGACAACGCATCTATCACCTTGTGATCAAAACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAACCATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGACTGGAAATGATCCAGCGCGt                                                                                     >  1:322736/1‑250 (MQ=255)                                                                                                                                     cGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGc                                                                                  >  1:6953/1‑251 (MQ=255)                                                                                                                                               agaTGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTatcatc                                                                        >  2:101170/1‑251 (MQ=255)                                                                                                                                                      cGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGttt                                                                                                                               >  2:212695/1‑189 (MQ=255)                                                                                                                                                      cGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGttt                                                                                                                               <  1:212695/189‑1 (MQ=255)                                                                                                                                                                  aTAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTAGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCg                                                     >  1:262330/1‑251 (MQ=255)                                                                                                                                                                   tAATTACAGACATTTAGGGCCACGCGATCATCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCACGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGGAGCGTTTATCATCCAGCGGGCGATGCTGTGCg                                                     >  2:63272/1‑250 (MQ=255)                                                                                                                                                                         cAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGAt                                              >  1:334173/1‑251 (MQ=255)                                                                                                                                                                                     gCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGCTCGGCACTGACCTTCGCCCTACACTCCACTAAAAAATCCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGCAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCa                                   >  1:46361/1‑250 (MQ=255)                                                                                                                                                                                                           gAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGcccc            >  2:279586/1‑251 (MQ=255)                                                                                                                                                                                                               aCGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGa         >  1:14795/1‑250 (MQ=255)                                                                                                                                                                                                                      tATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGCGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCa  >  1:310690/1‑250 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           CCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGG‑TTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCA  >  NC_000913/2032196‑2032657

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