Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	788,967	C→T	6.8%	G338D (GGC→GAC)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	788,967	0	C	T	6.8%	178.6 / 2.7	59	G338D (GGC→GAC)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base T (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

AGCACTGTCCTGCTCCTTGTGATGGTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATA  >  NC_000913/788833‑789113                                                                                                                                       |                                                                                                                                                   aGCACTGTCCTGCTCCTTGTGATGGTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATccgccg                                                                                                                                      >  2:309414/1‑149 (MQ=255) aGCACTGTCCTGCTCCTTGTGATGGTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATccgcc                                                                                                                                       <  7:408008/148‑1 (MQ=255)  gCACTGTCCTGCTCCTTGTGATGGTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATccgccgc                                                                                                                                     >  5:375300/1‑149 (MQ=255)   cACTGTCCTGCTCCTTGTGATGGTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATccgccgcc                                                                                                                                    >  5:9237/1‑149 (MQ=255)      tGTCCTGCTCCGTGTGATGGTTTACAAACGTAAAAAGTCTATTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGt                                                                                                                                 >  6:66474/1‑149 (MQ=255)         tattcTCCTTGTTATGGTTTTAAAAGTTAAAAAGTGTATTTATTTCGTTTTTTTGCTTAATATTGTTCAGCGAGAGCTTGCTGTACGGCTGGCTCCAGCTCTTGCGGGATCAGCTCGACGATTCAGCTGCCAAATCCGCCGCCGGTCa                                                                                                                               <  2:334122/144‑1 (MQ=255)          cTGCTCCTTGTGATGGTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCAt                                                                                                                              <  1:127148/148‑1 (MQ=255)                 ttGTGATGGTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGc                                                                                                                      >  5:151461/1‑149 (MQ=255)                        ggTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCttt                                                                                                               >  8:64663/1‑149 (MQ=255)                          tttACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGt                                                                                                             <  1:309414/149‑1 (MQ=255)                            tACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTc                                                                                                            <  6:9237/148‑1 (MQ=255)                             aCAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGc                                                                                                          >  6:173422/1‑149 (MQ=255)                              cAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGc                                                                                                          >  7:357087/1‑148 (MQ=255)                               aaaCGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAACTCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGcc                                                                                                         >  6:233716/1‑148 (MQ=255)                               aaaCGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGcc                                                                                                         >  3:188525/1‑148 (MQ=255)                               aaaCGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTACGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCCTGCGTACGCCACCTTTGTCGcc                                                                                                         >  5:22932/1‑148 (MQ=255)                                        aaGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGTCGCCAAATCCGCCGCCgg                                                                                                                                  <  8:394629/114‑1 (MQ=255)                                        aaGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGTCGCCAAATCCGCCGCCgg                                                                                                                                  >  7:394629/1‑114 (MQ=255)                                             tctTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATcaca                                                                                                  <  1:247921/141‑1 (MQ=255)                                             tctTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATcaca                                                                                                  >  2:247921/1‑141 (MQ=255)                                             tctTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGAc                                                                                          >  5:114171/1‑149 (MQ=255)                                             tctTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGAc                                                                                          >  4:388678/1‑149 (MQ=255)                                               tttAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAAt                                                                                                      >  5:217327/1‑135 (MQ=255)                                               tttAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAAt                                                                                                      <  6:217327/135‑1 (MQ=255)                                               tttAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGa                                                                                        >  6:114393/1‑149 (MQ=255)                                                   aTACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTc                                                                                    <  6:375300/149‑1 (MQ=255)                                                   aTACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTc                                                                                    <  7:64663/149‑1 (MQ=255)                                                      ccTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTAc                                                                                 <  7:7613/149‑1 (MQ=255)                                                      ccTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTAc                                                                                 <  6:345135/149‑1 (MQ=255)                                                       cTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACAGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCAGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATAACAGCTTTGACGATTTCTAc                                                                                 >  5:29443/1‑148 (MQ=255)                                                          tttttGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCaga                                                                             <  5:173422/149‑1 (MQ=255)                                                               gCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAgt                                                                           <  4:208650/146‑1 (MQ=255)                                                               gCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAgt                                                                           >  3:208650/1‑146 (MQ=255)                                                                 ttCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGTCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCaa                                                                      >  3:145128/1‑149 (MQ=255)                                                                   cATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAAtt                                                                    >  7:1775/1‑149 (MQ=255)                                                                   cATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAAt                                                                     <  3:388678/148‑1 (MQ=255)                                                                              gCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGTCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGt                                                         <  4:145128/149‑1 (MQ=255)                                                                                gACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGa                                                       <  8:357087/149‑1 (MQ=255)                                                                                 aCAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGAt                                                      >  5:169102/1‑149 (MQ=255)                                                                                      ttGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCg                                                  >  7:65640/1‑148 (MQ=255)                                                                                        gCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGaa                                                <  6:23410/148‑1 (MQ=255)                                                                                        gCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGa                                                 <  4:188525/147‑1 (MQ=255)                                                                                        gCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCAGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGa                                                 <  8:415327/147‑1 (MQ=255)                                                                                            tACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAgt                                                                           <  5:373786/117‑1 (MQ=255)                                                                                            tACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAgt                                                                           >  6:373786/1‑117 (MQ=255)                                                                                              cGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAatcatc                                          <  6:114171/148‑1 (MQ=255)                                                                                              cGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATcg                                         <  8:1775/149‑1 (MQ=255)                                                                                               ggcaggcaCCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATcgc                                        <  8:65640/149‑1 (MQ=255)                                                                                               ggcaggcaCCAGCTCTTACGGGATCAGCGCGACCATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATcgc                                        >  6:138500/1‑149 (MQ=255)                                                                                                       ccAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATagag                                <  5:323042/149‑1 (MQ=255)                                                                                                            tctTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCAt                            >  8:37498/1‑148 (MQ=255)                                                                                                              ttCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCAt                            >  6:400351/1‑146 (MQ=255)                                                                                                                  gggATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTc                     <  5:114393/149‑1 (MQ=255)                                                                                                                         gcgcGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCAt               >  6:379327/1‑148 (MQ=255)                                                                                                                          cgcgACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCgg                                                               <  8:256370/99‑1 (MQ=255)                                                                                                                           gcgACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCaa            <  5:233716/149‑1 (MQ=255)                                                                                                                            cgacgaTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAAc           >  7:58659/1‑149 (MQ=255)                                                                                                                                gaTACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGCGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGc       <  6:22932/149‑1 (MQ=255)                                                                                                                                     aGCCGCCAAATCCGCCGCCGGGCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATa  <  6:203927/149‑1 (MQ=255)                                                                                                                                       |                                                                                                                                                   AGCACTGTCCTGCTCCTTGTGATGGTTTACAAACGTAAAAAGTCTCTTTAATACCTGTTTTTGCTTCATATTGTTCAGCGACAGCTTGCTGTACGGCAGGCACCAGCTCTTCCGGGATCAGCGCGACGATACAGCCGCCAAATCCGCCGCCGGTCATGCGTACGCCACCTTTGTCGCCAATCACAGCTTTGACGATTTCTACCAGAGTGTCAATTTGCGGCACGGTGATTTCGAAATCATCGCGCATAGAGGCATGAGACTCCGCCATCAACTCGCCCATA  >  NC_000913/788833‑789113

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