Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	939	G→A	V201V (GTG→GTA)	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	939	0	G	A	92.9%	70.5 / ‑2.0	30	V201V (GTG→GTA)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (24/2);  total (26/4)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

GCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTA  >  NC_000913/707‑1183                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      gCCGTGGCGAGAAAATGTCGATCGCCATTATGTCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCCGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGGGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGGACTTGGACGCAACGGTTcc                                                                                                                                                                                                                                     >  1:35869/1‑251 (MQ=255)            aaaaTGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTActccgc                                                                                                                                                                                                                          >  1:189956/1‑248 (MQ=255)                       cGCCATTATGGCCGGCGGATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGCCTActccccggcggtgccgg                                                                                                                                                                                                               >  2:63015/1‑237 (MQ=255)                         ccATTAAGGCCGGCGTATTAGACGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGTTATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGATTCACTGCCGGTAATGACAAAGGCGACCTGGTGTTACTTTGACGCAACGGTTCCGACTACTCCGCGGCGGTgccggct                                                                                                                                                                                                             >  1:302491/1‑251 (MQ=255)                                       gTATTAGAAGCGCGCGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTCCTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGATGCCTGTTTACgcgc                                                                                                                                                                                               >  1:21671/1‑251 (MQ=255)                                              aaGCGCGAGGTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGACATTACCTCGAATCTACCGTCGAAATTGCTGAGTCCCCCCGCCGTATTGCGCCAAGTCGCATTCCGGCTGCTCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGCCTACTACGCGGCGGAGCTGGCTGCCAGTTTATGCGACGActgt                                                                                                                                                                                        >  1:313276/1‑251 (MQ=255)                                                      ggTCACAACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGGTGGGTGCCCTTTTTCGCGCCCGTTGTTGCgaggt                                                                                                                                                                                >  2:175749/1‑249 (MQ=255)                                                          acaACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTACACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTAATCCGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGctagcttttg                                                                                                                                                                            >  2:142420/1‑242 (MQ=255)                                                          acaACGTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTgg                                                                                                                                                                            >  1:318425/1‑251 (MQ=255)                                                               gTTACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTggacggg                                                                                                                                                                       >  1:183809/1‑250 (MQ=255)                                                                ttACCGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCAACCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGCTGGCAGGTTTCACCGCCGGTCATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCAGCGGCGGGTCTGGCTGCCTGTTTTAGCGCCGGttgttggcggtattgtgaggac                                                                                                                                                                      >  2:466144/1‑235 (MQ=255)                                                                              ccGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATCTTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCTCCGATTGTTGCGAGATTTGGATGG‑ACGTTGTCGGGGTCtt                                                                                                                                                       >  1:46057/1‑250 (MQ=255)                                                                               cGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTTGTTACCCGCCGCCTTTTGCGTGATTTGGACGG‑GCGtttgcggagtgttt                                                                                                                                                      >  2:184590/1‑239 (MQ=255)                                                                                   cGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGGGGCGGCTGCCTGTTTACGCGCCGATTTTTGCGGGGTTGTGACGGGACGTTGAGGGGGtgtttacc                                                                                                                                                   >  2:332005/1‑244 (MQ=255)                                                                                       aaaCTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCACGTCGCATTCCGGCTGATCCCATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCCGTTTCCCCGCCGGTTGTTGCGGGGTTTGGACGG‑ACGTTGACGGGGGCTTTACCTGCgc                                                                                                                                              >  1:163257/1‑250 (MQ=255)                                                                                            gctgGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCGGGGTGCCTGTTTTCGCGCCGATTGGTGCGGGGATTGGGCGG‑ACGTTGACGGGGGCTATACCTGCGGCCcgc                                                                                                                                         >  2:126821/1‑251 (MQ=255)                                                                                                                           gTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTActccgcg                                                                                                                                                                                                                         >  1:415229/1‑137 (MQ=255)                                                                                                                           gTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTc                                                                                                          >  1:98955/1‑251 (MQ=255)                                                                                                                                             aCCCGCCGTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGGGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGGTGACGGGGGCTATACCTGGGACCCGCGTCAGGTGCCCGATGCGGGGTTGGTGAAGTCGATGGCCTACCAGGAAGc                                                                                        >  1:127743/1‑251 (MQ=255)                                                                                                                                                cgcgttgtggcgtaaagcctctttCCGGCGGATCACATGGTGCTGATGGCAGTTTTCCCCGTCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGAt                                                                                     <  1:85033/229‑1 (MQ=255)                                                                                                                                                   cgTATTGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATggg                                                                                  >  1:83974/1‑250 (MQ=255)                                                                                                                                                       ttGCGGCAAGTCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGAGGGTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGGGCtt                                                                              >  1:304090/1‑251 (MQ=255)                                                                                                                                                                               aTCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTAAGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCCGCGCTTAAGttctt                                                      >  1:239051/1‑251 (MQ=255)                                                                                                                                                                                                     tttCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGGCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTAcccccc                                >  1:457735/1‑250 (MQ=255)                                                                                                                                                                                                            gccgGTAATGAAAAAGGCGAACTGGAGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATAACCCCCATCGccc                          >  1:352514/1‑250 (MQ=255)                                                                                                                                                                                                               ggTAATGAAAAAGGCGAACTGGTGGTCCTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGtt                      <  2:83974/251‑1 (MQ=255)                                                                                                                                                                                                               ggTAATGAAAAAGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGCAGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGGTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCCAAAGTTCTTCACCCCCGCACCATTACCCCCCTCGCCCAGtt                      >  1:25293/1‑251 (MQ=255)                                                                                                                                                                                                                      aaaaaGGCGAACTGGTGGTACTTGGACGCACCGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCCGTTTACGCGCCGGTTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGGGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGGTGGGGCTTTCCTACTTCGGCGCTTAAGTTCTTCACCCCCGCACCCTTTCCCCCCTCGccccggtccacgtc               >  2:230236/1‑240 (MQ=255)                                                                                                                                                                                                                        aaaGGCGAACTGGTGGTACTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATccc             <  1:391195/251‑1 (MQ=255)                                                                                                                                                                                                                            gCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCCGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGc         <  1:142420/251‑1 (MQ=255)                                                                                                                                                                                                                                    gtggtaCTTGGACGCAACGGTTCCGACTACTCCGCGGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGCGATTTGGGCGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGGCCTACCAGGAAGCGGTGGAGCTTTCCTAATTCGGCGCTAAAGTTCTTCACCCCCGCACCCTTACCCCCATCGCCCAGTTCCCGATCCCTTGCCTGATTa  >  2:12345/1‑250 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      GCCGTGGCGAGAAAATGTCGATCGCCATTATGGCCGGCGTATTAGAAGCGCGCGGTCACAACGTTACTGTTATCGATCCGGTCGAAAAACTGCTGGCAGTGGGGCATTACCTCGAATCTACCGTCGATATTGCTGAGTCCACCCGCCGTATTGCGGCAAGCCGCATTCCGGCTGATCACATGGTGCTGATGGCAGGTTTCACCGCCGGTAATGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGG‑ACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTA  >  NC_000913/707‑1183

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