Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,509	A→G	A391A (GCA→GCG)	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,509	0	A	G	100.0%	98.2 / NA	34	A391A (GCA→GCG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/19);  total (15/19)

TGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGA  >  NC_000913/1280‑1736                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      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‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑t                                                                                                                          >  1:99441/1‑251 (MQ=255)                                                                                         ttCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑tt                                                                                                                         <  1:216942/250‑1 (MQ=255)                                                                                          tCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTCCTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGGGGTAGGCGTTGGTATGCGCACCTTGCGTGGGTTCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGt                                                                                                                       >  2:18652/1‑251 (MQ=255)                                                                                                  gCTGATTACGCAATCATCTTGCGCATACAGGATCAGTTTCTGCGTTCCGCAAAGCGACTGGGTGCGAGCTGAACGGGCAATGCAGGAATAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCGTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCAGGCCCGCGCCAATATCCACA‑TTGTCGCCATTg                                                                                                               <  2:173914/251‑1 (MQ=255)                                                                                                         aCGCAATCATCTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGACGCTGGCGGTGAAGGAACGG‑CCGGCCATTATCTCGGTGGTAGGAGATGGTATTAGCACCTTGCGGGGGTGCTCGGAGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCCACA‑TTATCTCCATTGCTCAggg                                                                                                        >  1:110444/1‑251 (MQ=255)                                                                                                                catcTTCCGAATACAGTATCAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGCAACGG‑CTGGCCCTTATCCCGGTGGTAGGTGCTGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCACTATCAACA‑TTTTCGCCCTTTCTCAGGGAtcttct                                                                                                 >  1:139935/1‑251 (MQ=255)                                                                                                                     tcgcaatacggtaTCAGCTTCTGCGTTCCGCAAGGCGACTGGGTGCGAGCTGAACGGGCAATTCAGGAATAGTTCTCCCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACCGAACGT‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCACATATCAACA‑TTGTCGCCATTGCTCAGTGATCTTCTg                                                                                                <  2:226251/235‑1 (MQ=255)                                                                                                                               gtatccgTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAAAGG‑CTGGCCATTATATAGGTGGAAGGGGATGGTATGCGAACCCTGCGTGGGGTTTAGGCGAAATTTTTTGGCGCGCTGGGACGCGCCAATATCA‑CACTTGTCGCCAATGCTCGGGGCTCTTCTGAACCCTCaaattat                                                                                  >  2:264113/7‑246 (MQ=255)                                                                                                                                   cAGTTTCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGCAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATATCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCg                                                                              >  1:109980/1‑251 (MQ=255)                                                                                                                                       ttCTGCGTTCCGCAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTgg                                                                           <  1:68385/250‑1 (MQ=255)                                                                                                                                        tCTGCGTTCCGCAAAGCGACGGTGTGCGAGCTGAGCGGGTAATGCAGGAAGAGTGCTACCTGGAACTGAAAGAAGGCTTACTGGGGCCGCTGGCGTTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGAGCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCATTAGCAGCA‑TTGTCGCGATTGCGCAGGGATCTTCTGAACGCTCAATCTTTGTCGTGGt                                                                          <  2:181184/250‑1 (MQ=255)                                                                                                                                                         gactttggtcgtgCTGAACGGGCATTGCAGGAAGAGTGATTGCTGGAAGTGAACGAAGGCTTTCTGGAGCCGATGGAGGTGACGGACCGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCTCCTTGCTTGGGAGCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATACCGATGATGCGAc                                                        <  2:9183/239‑1 (MQ=255)                                                                                                                                                            tgggggCGAGCTGAACCGGCAATGCGGTATGAGGTCTACCTGGCACTGAGAGAAGGCTTACGGGGGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGGCGGCATTGCTCAGGGATCTTCTGAACTCTCAATCCCTGTCGTGGTAAATAACGATGATGCGACCa                                                      <  2:128877/245‑1 (MQ=255)                                                                                                                                                            gtggtgcgagctaatctggtactGCAGGAGGAGTGCTACCTTGAACTGAAAGAAGGCTTACTTGGTCCGCTGGCGGTGACGGAACGG‑CTCGCCATTATCTCGGTGGTAGGTGATGGTATGCGCCCCTTGCGTGGCAGCTCGGCGAAAGTCTTTGCCGCGCTGGCCCTCGCCAATATCAACA‑TTGTCGCCATCGCTCAGGGATCTGCTGAACGCTCAATCTCTGTCGCGGTAAATAACGATGATGCGACCa                                                      <  2:107179/228‑1 (MQ=255)                                                                                                                                                                 gCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCGGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCAGTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCg                                                <  1:194185/251‑1 (MQ=255)                                                                                                                                                                  cGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGATCTGAAAGAAGGCTTACTGGAGCCGCTGGCGGGGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGt                                               <  2:89263/251‑1 (MQ=255)                                                                                                                                                                                              tCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCaa                    <  1:203279/250‑1 (MQ=255)                                                                                                                                                                                                 aCCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAAt                   <  1:187916/248‑1 (MQ=255)                                                                                                                                                                                                      gAACTGAAAGAAGGCGTACTGGAGCCGCTGGCGGTGACGTAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGGTACTCAGCAGATGCTGTTCAATACCGATCa           <  2:27863/251‑1 (MQ=255)                                                                                                                                                                                                       aaCTGAAAGAAGGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCGTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGAGGCTGTGCAATTCCGATCAg          <  2:52144/251‑1 (MQ=255)                                                                                                                                                                                                               gaagGCTTACTGGAGCCGCTGGCGGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGCATGCGCCCCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCGCTGGCCCGCGCCAATTTCAACA‑TTGTCGCCATTGCTCAGGGCTCTTCTGAACGCTCAATCTCTGTCGTGGTCAATAACGGTGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCgc  >  1:133814/1‑250 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      TGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGG‑CTGGCCATTATCTCGGTGGTAGGTGATGGTATGCGCACCTTGCGTGGGATCTCGGCGAAATTCTTTGCCGCACTGGCCCGCGCCAATATCAACA‑TTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGA  >  NC_000913/1280‑1736

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