Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,403	A→G	Y356C (TAC→TGC)	thrA →	Bifunctional aspartokinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,403	0	A	G	100.0%	94.2 / NA	25	Y356C (TAC→TGC)	thrA	Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/10);  total (15/10)

GGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTG  >  NC_000913/1258‑1546                                                                                                                                                  |                                                                                                                                                ggCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATgcagca                                                                                                                                            >  2:490003/1‑151 (MQ=255)               gAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGtt                                                                                                                              >  2:516556/1‑151 (MQ=255)                    aCATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCcaca                                                                                                                         >  2:152515/1‑151 (MQ=255)                         gCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCg                                                                                                                    >  1:353150/1‑151 (MQ=255)                             tGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACtg                                                                                                                <  2:10186/151‑1 (MQ=255)                                    cGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCg                                                                                                          <  2:285517/150‑1 (MQ=255)                                       ttCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAg                                                                                                        >  2:415207/1‑149 (MQ=255)                                                  ggATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAAt                                                                                           >  1:520975/1‑151 (MQ=255)                                                           ggATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAag                                                                                   <  2:445316/150‑1 (MQ=255)                                                                gTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAag                                                                                   >  2:95806/1‑145 (MQ=255)                                                                gTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAag                                                                                   <  1:95804/145‑1 (MQ=255)                                                                     cATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCt                                                                         <  1:36614/150‑1 (MQ=255)                                                                          cggcgCGCGTCTTTGCAGCGATGTCTCGCGCCCGTATGTCCGTGTTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAAAGGGCAAGGCAGGAAGAGTTCTACCTGGAAc                                                                    <  1:432841/150‑1 (MQ=255)                                                                                   tCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAgg                                                          >  1:11498/1‑151 (MQ=255)                                                                                    cTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAgg                                                          >  2:162165/1‑150 (MQ=255)                                                                                            cGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTgg                                                  <  2:171054/150‑1 (MQ=255)                                                                                            cGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTgg                                                  <  1:281385/150‑1 (MQ=255)                                                                                              aTGTCGCGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAg                                                >  2:155801/1‑150 (MQ=255)                                                                                              aTGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAg                                                >  1:327118/1‑150 (MQ=255)                                                                                              aTGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAg                                                >  1:283682/1‑150 (MQ=255)                                                                                                 tCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCg                                             >  1:492656/1‑150 (MQ=255)                                                                                                      cgcCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCa                                       >  1:116537/1‑151 (MQ=255)                                                                                                                          tGATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCAt                   <  2:115116/151‑1 (MQ=255)                                                                                                                           gATTACGCAATCATCTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCAtt                  >  1:137871/1‑151 (MQ=255)                                                                                                                                    atcatcTTCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCggtg         >  2:124316/1‑151 (MQ=255)                                                                                                                                           tCCGAATGCAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTg  <  1:388511/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                GGCATTTCCAATCTGAATAACATGGCAATGTTCAGCGTTTCTGGTCCGGGGATGAAAGGGATGGTCGGCATGGCGGCGCGCGTCTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATACAGCATCAGTTTCTGCGTTCCACAAAGCGACTGTGTGCGAGCTGAACGGGCAATGCAGGAAGAGTTCTACCTGGAACTGAAAGAAGGCTTACTGGAGCCGCTGGCAGTGACGGAACGGCTGGCCATTATCTCGGTGGTAGGTG  >  NC_000913/1258‑1546

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