Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,329	G→A	A332T (GCG→ACG)	thrA →	bifunctional aspartokinase I/homeserine dehydrogenase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,329	0	G	A	100.0%	124.7 / NA	35	A332T (GCG→ACG)	thrA	bifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/12);  total (23/12)

TTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATA  >  REL606/1262‑1402                                                                    |                                                                          ttCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGgcgcgc                                                                   >  1:841672/1‑76 (MQ=255)    cAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCgtg                                                                >  2:1290498/1‑76 (MQ=255)          gAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCa                                                          >  1:280913/1‑76 (MQ=255)           aataatATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAg                                                         >  2:487514/1‑76 (MQ=255)             taataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCg                                                       >  1:1119757/1‑76 (MQ=255)              aataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGa                                                      >  1:990150/1‑76 (MQ=255)              aataTGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGCGTTTGCAGCgg                                                      >  1:990152/1‑75 (MQ=255)               atatGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTNTGCAGCGAt                                                     <  1:35698/76‑1 (MQ=255)                   ggCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCa                                                 >  2:1065280/1‑76 (MQ=255)                       aTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCAcgcg                                             >  1:831552/1‑76 (MQ=255)                       aTGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCAcgcg                                             <  2:477567/76‑1 (MQ=255)                          ctCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGccc                                          <  2:806857/75‑1 (MQ=255)                           tCAGCGTGTCCGGCCCGGGGATGAAAGGGAGGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCg                                         <  1:519648/76‑1 (MQ=255)                              gCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTAt                                      <  2:792002/76‑1 (MQ=255)                                gTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGAGTGTTTGCAGCGATGTCCCGCGCCCGTAttt                                    >  1:1133578/1‑76 (MQ=255)                                 tttCCGGCCCGGGGATGAAAGGGATGGATGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTc                                   >  1:1293751/1‑76 (MQ=255)                                  tgccggccggGGGATGAAAGGGAGGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTcc                                  <  2:719865/74‑1 (MQ=255)                                   tCCGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCg                                 >  2:849414/1‑76 (MQ=255)                                   tCCGCCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCTGCGATGTCACGCGACCGTGTTGNCg                                 >  2:1244582/1‑76 (MQ=255)                                    ccGGCCCGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCgt                                <  2:1069545/76‑1 (MQ=255)                                     cGGCCCGGGGATGAAAGGGATGGTTGGCATGACNGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCgtg                               >  2:1121923/1‑76 (MQ=255)                                       gCCGGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCgtggt                             <  1:886773/76‑1 (MQ=255)                                         ccGGGGATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGc                           >  2:722831/1‑76 (MQ=255)                                           ggggATGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTg                         >  1:1005925/1‑76 (MQ=255)                                               aTGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTa                     <  2:950340/76‑1 (MQ=255)                                                tGAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTAc                    >  2:655377/1‑76 (MQ=255)                                                 gAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTTCg                   >  1:660243/1‑76 (MQ=255)                                                 gAAAGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACg                   >  1:528884/1‑76 (MQ=255)                                                  aaaGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGc                  >  1:393665/1‑76 (MQ=255)                                                   aaGGGATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCa                 >  2:838340/1‑76 (MQ=255)                                                       gATGGTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAatca             >  1:954401/1‑76 (MQ=255)                                                          ggTTGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCt          <  2:860692/76‑1 (MQ=255)                                                             tGGCATGTCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTcc       <  2:251196/76‑1 (MQ=255)                                                             tGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTcc       >  2:699653/1‑76 (MQ=255)                                                             tGGCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGNATTTCCGTGGTGCTGATTACGCAATCATCTTcc       <  2:1114350/76‑1 (MQ=255)                                                              ggCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCg      <  2:393665/76‑1 (MQ=255)                                                               gCATGACGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGa     <  2:856995/76‑1 (MQ=255)                                                                  tGGCGGCGCTCGTGTTTGCAGCGATGTAACGCGCCCGTATTTCCGTGGTGCTGATTACGCAAGCATCTTCCGAATa  <  2:743739/76‑1 (MQ=255)                                                                    |                                                                          TTCCAATCTGAATAATATGGCAATGTTCAGCGTTTCCGGCCCGGGGATGAAAGGGATGGTTGGCATGGCGGCGCGCGTGTTTGCAGCGATGTCACGCGCCCGTATTTCCGTGGTGCTGATTACGCAATCATCTTCCGAATA  >  REL606/1262‑1402

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