Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,407	C→T	S357S (AGC→AGT)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,407	0	C	T	93.5%	88.4 / NA	31	S357S (AGC→AGT)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/10);  minor base G (0/2);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCT  >  NC_000913/1166‑1641                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            aGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTCATGAAGACGCATTACCGGTCAAGGGAATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTCTTTCCGTGGAGCTGGTTACGCCATCATCTTCCGCATACAGTATCCGTTTc                                                                                                                                                                                                                                    >  1:254309/1‑251 (MQ=255)   aTCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTGCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGTCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTAGCAGTTTCTg                                                                                                                                                                                                                                  <  2:135321/251‑1 (MQ=255)    tCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGACTAAAGTATCAGTTTCTGc                                                                                                                                                                                                                                 >  1:52144/1‑251 (MQ=255)          gCCTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTGTCTGCGTTccg                                                                                                                                                                                                                           <  2:20979/251‑2 (MQ=255)            cTGATTAAAAATACCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATTCAGTATCAGTTTCTGCGTTccgc                                                                                                                                                                                                                          >  1:175886/1‑248 (MQ=255)                        aCCGGAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACtgtgt                                                                                                                                                                                                             >  2:203279/1‑251 (MQ=255)                            gAAATCCTCAAGCTCCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGa                                                                                                                                                                                                         <  2:142823/251‑1 (MQ=255)                                          ccAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCaa                                                                                                                                                                                            <  1:125756/250‑1 (MQ=255)                                                        gggTGCCAGCCTGGATGAAGACGAATTACCGGTCAAGGGCAGTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGTATGAAAGGGATGTTTGGCATGGCGGCGCGCGTGTTTGCAGCTATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTGCCGAAGACTGTATCAGTTTCTGCGTGCCTCAAAGCGACTGTGTGCGAGCTGAGCGGGCAATGCAGGAAGAGTTc                                                                                                                                                                              <  2:250914/249‑1 (MQ=255)                                                         ggTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGCTGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCtc                                                                                                                                                                            >  1:27863/1‑250 (MQ=255)                                                           tGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGAGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGCGCTGAACGGGCAATGCAGGAAGAGTTCTAcc                                                                                                                                                                          >  2:103020/1‑251 (MQ=255)                                                               acccGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCCCGCGCCCGTATTTCCGTGGTGGTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCCAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTgg                                                                                                                                                                       >  2:187916/3‑250 (MQ=255)                                                                                                                   tGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATCGAGGCGCGCGTGGTTGCCGCGATGACACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTAACAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGACGGGACGG‑CTGGCCAtt                                                                                                                 >  1:181184/1‑251 (MQ=255)                                                                                                                          gTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGCGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGCGGACCGGGCTGGC‑ATTTTCTCgg                                                                                                          >  1:135329/1‑251 (MQ=255)                                                                                                                                cGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGGATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCggtggt                                                                                                      >  1:65751/1‑249 (MQ=255)                                                                                                                                                        gttGGTTGGCATGGCGGCGCGCGGGTTTGCAGCGATGTCACGCGGCCGTACTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACCGTATCCGCTTCTTCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCCGGAAGTGTTCTACCTGGAACTGAACGAGGGCCTACTGGAGCCGCAGGCGGGGCAGGAACAG‑CCTGCCATTAACTCGGGGGTAggggagcgtattgcgcacatccgtg                                                                            >  2:72732/3‑228 (MQ=255)                                                                                                                                                           ggttggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGGGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGGCGGAACGG‑CGGGCCCTTTACTCGGGGGGAGGGGGTGGGTTGCGCACCTTGCGTgggg                                                                         >  2:128914/5‑250 (MQ=255)                                                                                                                                                               ggCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCTCACCGTGCGTGGTAtctc                                                                     <  1:79513/251‑1 (MQ=255)                                                                                                                                                                           cgcgtgattGCAGCACTGTCACGCGCCCGTATTTCAGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTTCGTGGGGTCTCGGCGAAATTCtt                                                         >  2:35778/8‑251 (MQ=255)                                                                                                                                                                              gtgtttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGc                                                      <  1:6113/248‑1 (MQ=255)                                                                                                                                                                                gtttGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGcc                                                     <  2:85857/249‑1 (MQ=255)                                                                                                                                                                                             tCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CCGGCCATTATCTCGGGGGTAGGTGATGGTATGCGCACCTTGCGTGGGGTCTCGGCGGAATTCTTTGCCGCGGTGGCCCccgc                                       >  2:216942/1‑251 (MQ=255)                                                                                                                                                                                               aCGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGGGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGcc                                      >  1:183362/1‑250 (MQ=255)                                                                                                                                                                                                 gcgcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGGATCAGTTTCTGCGTTCCGCAAAGCGACGGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGTCTTTCTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTGTGCCGCGCTGGCCCGCGCCAat                                   <  2:143534/251‑1 (MQ=255)                                                                                                                                                                                                       gTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGGATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAAc                             <  2:65751/251‑1 (MQ=255)                                                                                                                                                                                                         aTTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACAt                           >  1:99441/1‑251 (MQ=255)                                                                                                                                                                                                           ttCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACAtt                          <  1:216942/250‑1 (MQ=255)                                                                                                                                                                                                            tCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGGGGTAGGCGTTGGTATGCGCACCTTGCGTGGGTTCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACATTGt                        >  2:18652/1‑251 (MQ=255)                                                                                                                                                                                                                    gCTGATTACGCAATCATCTTGCGCATACAGGATCAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGCAGGAATAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCGTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCAGGCCCGCGCCAATATCCACATTGTCGCCATTg                <  2:173914/251‑1 (MQ=255)                                                                                                                                                                                                                           aCGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGACGCTGGCGGTGAAGGAACGG‑CCGGCCATTATCTCGGTGGTAGGAGATGGTATTAGCACCTTGCGGGGGTGCTCGGAGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCCACATTATCTCCATTGCTCAggg         >  1:110444/1‑251 (MQ=255)                                                                                                                                                                                                                                  catcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGCAACGG‑CTGGCCCTTATCCCGGTGGTAGGTGCTGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCACTATCAACATTTTCGCCCTTTCTCAGGGAtcttct  >  1:139935/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            AGATCCCTTGCCTGATTAAAAATACCGGAAATCCTCAAGCACCAGGTACGCTCATTGGTGCCAGCCGTGATGAAGACGAATTACCGGTCAAGGGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACATTGTCGCCATTGCTCAGGGATCTTCT  >  NC_000913/1166‑1641

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