Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	790,832	Δ1 bp	coding (198/1047 nt)	galT ←	galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	790,832	0	G	.	96.8%	133.7 / ‑3.2	31	coding (198/1047 nt)	galT	galactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (6/24);  minor base A (1/0);  total (7/24)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACCCATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCC  >  NC_000913/790587‑791079                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        ccccATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTa                                                                                                                                                                                                                                                   <  1:93551/250‑1 (MQ=255)  cccATGGGTACGTTTTCCCCAGGTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTAt                                                                                                                                                                                                                                                  <  1:119157/251‑1 (MQ=255)           aCGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCGTGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGGCATTAGTGTAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGGCTGTc                                                                                                                                                                                                                                          <  2:143740/250‑1 (MQ=255)                                      ccTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCCCCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAg                                                                                                                                                                                                               <  1:198354/250‑1 (MQ=255)                                      ccTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAg                                                                                                                                                                                                               <  1:198338/250‑1 (MQ=255)                                            aGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTg                                                                                                                                                                                                        >  2:76382/1‑251 (MQ=255)                                                                      gcAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGttt                                                                                                                                                                               <  1:1209/250‑1 (MQ=255)                                                                      gcAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTtg                                                                                                                                                                              <  1:69555/251‑1 (MQ=255)                                                                              gAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCg                                                                                                                                                                       <  2:103076/250‑1 (MQ=255)                                                                                       cAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTgcgc                                                                                                                                                             >  1:157363/1‑251 (MQ=255)                                                                                             tttACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGGTTCCTGCGCCCCCt                                                                                                                                                        <  1:47704/250‑1 (MQ=255)                                                                                               tACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACGCCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGc                                                                                                                                                      <  1:94476/250‑1 (MQ=255)                                                                                                       tCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGGGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGGGCGCCGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGgcgc                                                                                                                                             >  2:205487/1‑251 (MQ=255)                                                                                                         cGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCtt                                                                                                                                           <  1:60543/251‑1 (MQ=255)                                                                                                                          cccGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCGCATTACCGGCGCAGAGGGATCAAGCTGGAGCGGGCGCTGGTAACACCGGTTGGGCTGGCGTTTCCTGCGGCGCCTGCGAGGGGGGGTCAGCGCGGGGCGATGaaa                                                                                                                          <  2:122576/251‑1 (MQ=255)                                                                                                                                      cgcgcgcgCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCAtt                                                                                                              <  1:93015/251‑1 (MQ=255)                                                                                                                                                      aCGCATCAGCGGATCGTGACTTTCTGGCTCATTTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAATTCCCGGTGTAATCTGG‑TTTTTATCGCCTGTCACCCGCACATTTCCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACTCCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGTCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGtt                                                                                              <  1:114950/251‑1 (MQ=255)                                                                                                                                                         cATCAGCGGATCGTGACTTTCTGGCGCATCTGGAGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTa                                                                                           <  2:192544/251‑1 (MQ=255)                                                                                                                                                              gCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAgcggc                                                                                      <  2:143651/251‑1 (MQ=255)                                                                                                                                                               cGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAgcggcg                                                                                     <  1:144986/251‑1 (MQ=255)                                                                                                                                                                                             gtCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGt                                                       <  1:180800/251‑1 (MQ=255)                                                                                                                                                                                                 gACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATg                                                   <  1:173088/251‑1 (MQ=255)                                                                                                                                                                                                  aCATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGGAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATgg                                                  >  2:2348/1‑251 (MQ=255)                                                                                                                                                                                                     tCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCg                                               >  1:55828/1‑251 (MQ=255)                                                                                                                                                                                                      cAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGt                                              >  2:160088/1‑251 (MQ=255)                                                                                                                                                                                                          gCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCt                                          <  1:171622/251‑1 (MQ=255)                                                                                                                                                                                                             gCAAAGGCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAGACCAGAATCTATTGCCCTGGGAGCGTTTTGTTGCGGCGATGTGGATGATCGACGGGATGAAAATGCGTCATCGACGTTCCTTaa                                       <  2:111808/251‑1 (MQ=255)                                                                                                                                                                                                                           gAAAACGTAAATCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCCCCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACTCCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCGGTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGAGTAAATTGCGTCATGGTCGTGCCTTAATCGTGATCTCCCtg                         <  2:29605/251‑1 (MQ=255)                                                                                                                                                                                                                                 gTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATg                   <  2:72268/251‑1 (MQ=255)                                                                                                                                                                                                                                        caggTGTAATCGAGATTTT‑ATCGCCTGTCACCCGCACACTACCTGCGCAGAAGACGCAATATGGATCGTGCGCAGGTAACCACTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCACCATCCATTGCCCGGTGAGCGGGCTGTGGCGGCGATGGGGACGATCAACGGGATTAACTTGAGTCCTGGTCGGTCCTTAATCAGGATATCCCTGTGGATGGCGTGAc            >  2:240561/3‑251 (MQ=255)                                                                                                                                                                                                                                             gtAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCAGTGCCCGGTGAGCGGGTTGTAGCGGCGATGGGGATGGTTAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGTGAGATCCCTATGGACGGCGTGACtgcc        <  2:148059/250‑1 (MQ=255)                                                                                                                                                                                                                                                    gggTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgccagtgcc  >  2:35585/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        ACCCATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGCC  >  NC_000913/790587‑791079

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: ‑