Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,580	C→A	E134* (GAA→TAA)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,580	0	C	A	96.9%	94.1 / NA	32	E134* (GAA→TAA)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (15/16);  total (15/17)
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

ACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGA  >  NC_000913/789344‑789819                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                aCTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTTCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGgaag                                                                                                                                                                                                                                   <  1:116647/251‑1 (MQ=255)      tgatgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGATGCGGAAGAACtt                                                                                                                                                                                                                              >  2:195419/1‑251 (MQ=255)          gatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAAcc                                                                                                                                                                                                                          >  1:23535/1‑251 (MQ=255)          gaggacgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGAGACTTACAGTGAAGCGGAAGAACTTAAcc                                                                                                                                                                                                                          <  2:150410/248‑1 (MQ=255)                        cacCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCTGCAGATGATAAAGCTGCTGCACTACGGTTCCGACCGCGACTTACCGTGAAGCGGAAGAACTTAACCCGGCCCCCTGCGGc                                                                                                                                                                                                            >  1:97329/1‑251 (MQ=255)                                            gCTTTGGTCCCCAGGGAGCGGCAATCGATCAGCAATGCAGGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTCCAGCGTACAAACTGGTTGACTGCTTCCTGACCGTTAAGCGCGATGTGTGCGCCGGCCAGCGGCACATGATAAAGCTGCTGCAAGACGTGTCCGACCGCGACCTACAGTGAAGCGGAAGCACGGAACCCGGCACCCTGCGGcacaatcccaatgagcacca                                                                                                                                                                                         <  2:110034/250‑17 (MQ=255)                                                tGGTCCCCAGTGAGCGTCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTc                                                                                                                                                                                     <  1:115006/250‑1 (MQ=255)                                                 ggTCCCCAGTGAGCGGCAATCTATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTGAAGCGCGATGTGTGCGCCTGCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCTACCGCGACTTACAGTGAAGCGGAAGAACGGAACCCGGCACCCTGCGGCACATTTCCGCTGATCTCCATGTcc                                                                                                                                                                                    <  1:99887/250‑1 (MQ=255)                                                      ccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcgccg                                                                                                                                                                              >  2:26009/1‑251 (MQ=255)                                                      ccAGGGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcgccg                                                                                                                                                                              <  2:149673/251‑1 (MQ=255)                                                               ggCAAGCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGTCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACATTGAAGCTGAAGAACGTATCCCGGCACCCTGCGGCACATTGCCGCTGATCGCCATGTCCACGTCGCCGAAGCtg                                                                                                                                                                     <  2:112225/251‑1 (MQ=255)                                                                                                cGAGCGCGGAAATTAGCTGATCCAGGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgcca                                                                                                                                    <  2:6398/251‑1 (MQ=255)                                                                                                  aGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTACCCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCAccacgccacg                                                                                                                                  <  2:200903/251‑1 (MQ=255)                                                                                                                  taTCCATGAGCCCTCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAATCGCGATTTGTGCGCCGTCCAGCGTCAGATTATAAAGCTGCTGCATTACGGTGCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGTCCCCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCGGTTGTTACTCAGTTGCAGATTTTTCACCACGCCACGAACGTAGTTAGccc                                                                                                                    <  1:228744/248‑1 (MQ=255)                                                                                                                       aTGATCCCGCAGTTACAGCCTACAAACTGGGTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATa                                                                                                             <  1:4058/251‑1 (MQ=255)                                                                                                                          aTCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGtt                                                                                                          >  2:74796/1‑251 (MQ=255)                                                                                                                               gCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCAt                                                                                                     <  1:26009/251‑1 (MQ=255)                                                                                                                                            aCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCAGTCCAGTGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACAATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTATCCACGCCACGAACGTAGTTAGCCCATTGATAGGTTTCATGTGCGACAATg                                                                                          <  2:97329/249‑1 (MQ=255)                                                                                                                                                 cTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGCACGAAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCAt                                                                                   >  1:26245/1‑251 (MQ=255)                                                                                                                                                 cTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCAt                                                                                   >  1:2991/1‑251 (MQ=255)                                                                                                                                                           gCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACGTAACCCGGCTCCCTGCGGCACATTGCCGCTGATCACCATGGCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGa                                                                            <  2:191693/248‑1 (MQ=255)                                                                                                                                                               ccTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCCCCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcg                                                                     >  1:40823/1‑251 (MQ=255)                                                                                                                                                                  gACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgtcg                                                                  >  2:225150/1‑251 (MQ=255)                                                                                                                                                                             gcgATTTGTGCGCCGTCCAGCGGCCGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGCGGGAAAACTCGTCGAGCTGATTTTc                                                       >  2:63664/1‑251 (MQ=255)                                                                                                                                                                                                 cGGCAGATGATAAAGCTGCTGCAATAAGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATCATCGGCTGCCATCAcgc                                   >  2:252322/1‑251 (MQ=255)                                                                                                                                                                                                         gatAAAGCTGCTGCAACACCGGTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTa                           <  2:41822/251‑1 (MQ=255)                                                                                                                                                                                                                       aaTACGGTTCCGAACGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGCAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtgg             >  2:211177/1‑251 (MQ=255)                                                                                                                                                                                                                        aTACGGTTCCGACCGCGACTTACAGTGAAGCGGCAGAACTTAAACCGGAACCCTGCGGCACATTGCCGCTGCTCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtggt            >  1:222363/1‑251 (MQ=255)                                                                                                                                                                                                                            ggTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGcac        >  2:16203/1‑251 (MQ=255)                                                                                                                                                                                                                             gTTCAGACCGCGACTTACAGAGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCAGCCGACGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCCCGAACGAAGTTAGCACATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATAGGGTGCCATCACGCGAACTTTACGTTCATAGCGGGGTGcac        >  2:235492/1‑250 (MQ=255)                                                                                                                                                                                                                                ccGACCGCGCCTTCCAGTGAAGCTGAAGAACTTACCCCGGCACCCTGCGGCACATTGCCGCTGAGCACCATGTCCACGCCGCCGAATCTGTTGTTACGCAGTTGCAGATTTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGACCTTTACGGTCATCGCGTGGTGCACAACt    <  2:148849/251‑1 (MQ=255)                                                                                                                                                                                                                                  gACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGa  <  1:225150/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                ACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGA  >  NC_000913/789344‑789819

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