Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,674,839	A→G	9.2%	K2K (AAA→AAG)	hom →	homoserine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,674,839	0	A	G	9.2%	142.7 / 5.0	54	K2K (AAA→AAG)	hom	homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (21/28);  new base G (2/3);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

AGGTCATTCCGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGAAGTGGCACAGATCGCCATGCGCTCGCAGAACCCG  >  NC_002947/1674692‑1674980                                                                                                                                                     |                                                                                                                                               atgTCATTCCGGGCTACCAGCCGGCACCACAGGTCCCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAAc                                                                                                                                               <  2:66444/147‑1 (MQ=255)   gTCATTCCGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCg                                                                                                                                             <  8:106620/149‑1 (MQ=255)   gTCATTCCGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCg                                                                                                                                             >  4:144772/1‑149 (MQ=255)      aTTCCGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGt                                                                                                                                           >  1:61912/1‑148 (MQ=255)       ttCCGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTc                                                                                                                                          >  1:184917/1‑148 (MQ=255)          cGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCGAGCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCaaa                                                                                                                                       <  1:66544/149‑1 (MQ=255)          cGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGt                                                                                                                                           <  1:185859/144‑1 (MQ=255)          cGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGt                                                                                                                                           >  2:185859/1‑144 (MQ=255)          cGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCaaa                                                                                                                                       <  4:139726/148‑1 (MQ=255)          cGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAg                                                                                                                                      >  3:180426/1‑149 (MQ=255)          cGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGGTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCaaa                                                                                                                                       <  6:12613/148‑1 (MQ=255)             gCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGT‑ATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAgg                                                                                                                                  >  3:252173/1‑149 (MQ=255)               tACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAgg                                                                                                                                  >  1:94391/1‑148 (MQ=255)                        gCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCtgt                                                                                                                           <  1:55319/146‑1 (MQ=255)                        gCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCtgt                                                                                                                           >  2:55319/1‑146 (MQ=255)                           ccacAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCATTGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTg                                                                                                                     <  2:294786/149‑1 (MQ=255)                           ccacAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTg                                                                                                                     >  7:336468/1‑149 (MQ=255)                           ccacAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTg                                                                                                                     <  8:335373/149‑1 (MQ=255)                            cacaGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTgg                                                                                                                    <  3:155651/149‑1 (MQ=255)                               aGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGt                                                                                                                                           <  2:211644/123‑1 (MQ=255)                               aGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGt                                                                                                                                           >  1:211644/1‑123 (MQ=255)                                          cTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGc                                                                                                      <  2:82281/149‑1 (MQ=255)                                             gCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGa                                                                                                   <  3:53229/149‑1 (MQ=255)                                               aCTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAAc                                                                                                 <  3:208939/149‑1 (MQ=255)                                                   gCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCtt                                                                                              <  4:180426/148‑1 (MQ=255)                                                   gCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCtt                                                                                              >  5:69731/1‑148 (MQ=255)                                                   gCCAAGT‑ATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCtt                                                                                              <  4:252173/147‑1 (MQ=255)                                                     cAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCaa                                                                                           >  2:3880/1‑149 (MQ=255)                                                     cAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCaa                                                                                           >  2:287715/1‑149 (MQ=255)                                                      aaGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAAt                                                                                          <  1:287715/149‑1 (MQ=255)                                                      aaGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAAt                                                                                          <  5:218764/149‑1 (MQ=255)                                                            tCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCCCAGTGAAGCCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACtt                                                                                    >  4:91238/1‑149 (MQ=255)                                                            tCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAGCCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACtt                                                                                    >  2:92223/1‑149 (MQ=255)                                                                           cgtcgtCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCgag                                                                     >  2:288681/1‑149 (MQ=255)                                                                               gtcATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGgaga                                                                 <  3:144772/149‑1 (MQ=255)                                                                                      ggTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCccg                                                          <  2:94391/149‑1 (MQ=255)                                                                                            ggCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGc                                                   >  5:299140/1‑149 (MQ=255)                                                                                              cAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGCTTGCCCG‑CCGTGCCg                                                 >  7:326604/1‑149 (MQ=255)                                                                                                 gTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGGCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGgc                                              <  8:336468/149‑1 (MQ=255)                                                                                                  tttCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGgcg                                             <  2:61912/149‑1 (MQ=255)                                                                                                    tCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGgcgc                                            <  6:69731/148‑1 (MQ=255)                                                                                                        ggTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTAt                                       >  8:261502/1‑149 (MQ=255)                                                                                                        ggTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTAt                                       >  1:161653/1‑149 (MQ=255)                                                                                                            ggCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGGTGGGGACCGTCGGTGGCGGAACCTTCAATGTACATCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGaa                                   <  2:76055/149‑1 (MQ=255)                                                                                                                   gtgtCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGAAGTGGcac                            <  8:129647/149‑1 (MQ=255)                                                                                                                          gCCGATCAATGGGGAGTTAACAGTGAAACCGGTCAAATTAGGCATCAGTCGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGAGCAACGCCGAGGAGATTGCCCT‑CCGTCCCGGGAGCGGTATTGAAGTGGAACAGATCg                      >  2:214664/1‑149 (MQ=255)                                                                                                                              ttCAATGGGGAGTTCACAGTGAAGCCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGAAGTGGCACAGATCGCCAt                  <  4:353454/149‑1 (MQ=255)                                                                                                                               tCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACg                                                                          <  7:60545/93‑1 (MQ=255)                                                                                                                               tCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACg                                                                          >  8:60545/1‑93 (MQ=255)                                                                                                                                     gggAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGAAGTGGCACAGATCGCCATGCGCTc            <  7:261502/148‑1 (MQ=255)                                                                                                                                      ggAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGAAGTGGCACAGATCGCCATGCGCTCGc          <  2:161653/149‑1 (MQ=255)                                                                                                                                             acaGTGAAGCCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCGCCCGTGCCGGGCGCGGTATTGAAGTGGCACAGATCGCCATGCGCTCGCAGAAcc    <  1:92223/149‑1 (MQ=255)                                                                                                                                             acaGTGAAGCCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGAAGTGGCACAGATCGCCATGCGCTCGCAGAAcc    <  3:91238/148‑1 (MQ=255)                                                                                                                                               aGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGAAGTGGCACAGATCGCCATGCGCTCGCAGAACCCg  >  2:160626/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                               AGGTCATTCCGGGCTACCAGCCGGCACCACAGGTCGCCAAGCTGGCACTTGCCAAGTAATCCAGTTCAGTACGCCGTCGTCATGGGGTGACGGCAAGTTTCACGGTCGGCAAACGTGTCG‑GCCGTTCAATGGGGAGTTCACAGTGAAACCGGTCAAAGTAGGCATCTGTGGGTTGGGGACCGTCGGTGGCGGAACCTTCAATGTACTTCAGCGCAACGCCGAGGAGATTGCCCG‑CCGTGCCGGGCGCGGTATTGAAGTGGCACAGATCGCCATGCGCTCGCAGAACCCG  >  NC_002947/1674692‑1674980

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