Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,450	G→T	L705L (CTG→CTT)	thrA →	bifunctional aspartokinase I/homeserine dehydrogenase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,450	0	G	T	100.0%	261.4 / NA	74	L705L (CTG→CTT)	thrA	bifunctional aspartokinase I/homeserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (44/30);  total (44/30)

AACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAA  >  REL606/2402‑2497                                                 |                                                aaCGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTcc                                               >  1:6836175/1‑51 (MQ=255) aaCGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTcc                                               >  1:5168081/1‑51 (MQ=255) aaCGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTcc                                               <  1:2598640/51‑1 (MQ=255)  aCGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTccc                                              <  1:860814/51‑1 (MQ=255)  aCGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTccc                                              <  1:2446278/51‑1 (MQ=255)    gggACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGc                                            >  1:6725419/1‑51 (MQ=255)    gggACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGc                                            >  1:6928036/1‑51 (MQ=255)    gggACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGc                                            >  1:3868053/1‑51 (MQ=255)       aCGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCaga                                         <  1:2534549/51‑1 (MQ=255)       aCGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCaga                                         <  1:3237737/51‑1 (MQ=255)        cGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCagag                                        >  1:2473504/1‑51 (MQ=255)        cGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCagag                                        <  1:2640828/51‑1 (MQ=255)         gTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGt                                       >  1:6009445/1‑51 (MQ=255)         gTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGt                                       >  1:754369/1‑51 (MQ=255)         gTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGt                                       >  1:1312939/1‑51 (MQ=255)         gTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGt                                       >  1:3827445/1‑51 (MQ=255)          tGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGtt                                      >  1:4054544/1‑51 (MQ=255)          tGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGtt                                      <  1:4961986/51‑1 (MQ=255)          tGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGtt                                      <  1:3545906/51‑1 (MQ=255)          tGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGtt                                      >  1:2203441/1‑51 (MQ=255)          tGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGtt                                      <  1:2225791/51‑1 (MQ=255)           gAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGttt                                     <  1:975145/51‑1 (MQ=255)           gAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGttt                                     <  1:4987434/51‑1 (MQ=255)           gAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGttt                                     >  1:749416/1‑51 (MQ=255)           gAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGttt                                     <  1:5573831/51‑1 (MQ=255)             aCTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTaa                                   <  1:7195997/51‑1 (MQ=255)             aCTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTaa                                   >  1:4586177/1‑51 (MQ=255)              cTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAAc                                  >  1:2284298/1‑51 (MQ=255)              cTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAAc                                  >  1:657190/1‑51 (MQ=255)               tGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACg                                 >  1:4504891/1‑51 (MQ=255)               tGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACg                                 >  1:2408663/1‑51 (MQ=255)               tGGAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACg                                 >  1:6820977/1‑51 (MQ=255)                ggAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGc                                >  1:1614204/1‑51 (MQ=255)                 gAGCTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCt                               >  1:5416920/1‑51 (MQ=255)                    cTGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAg                            >  1:5579819/1‑51 (MQ=255)                     tGGCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAgg                           >  1:4559646/1‑51 (MQ=255)                      ggCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAggg                          >  1:2814277/1‑51 (MQ=255)                      ggCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAggg                          <  1:5493084/51‑1 (MQ=255)                       gCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGt                         <  1:3657594/51‑1 (MQ=255)                       gCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGt                         <  1:1458823/51‑1 (MQ=255)                       gCGGATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGt                         <  1:7095378/51‑1 (MQ=255)                          gATATTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCt                               <  1:5347316/42‑1 (MQ=255)                           atatTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATg                     >  1:4412443/1‑51 (MQ=255)                           atatTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATg                     >  1:5781102/1‑51 (MQ=255)                           atatTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATg                     <  1:5879522/51‑1 (MQ=255)                           atatTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATg                     >  1:663413/1‑51 (MQ=255)                           atatTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATg                     >  1:7246906/1‑51 (MQ=255)                           atatTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATg                     >  1:1838233/1‑51 (MQ=255)                            tatTGAAATTGAACCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGt                    <  1:3792883/51‑1 (MQ=255)                               tgaaattgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGc                 >  1:5118647/1‑51 (MQ=255)                               tgaaattgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGc                 <  1:1935968/51‑1 (MQ=255)                               tgaaattgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGc                 >  1:6180478/1‑51 (MQ=255)                               tgaaattgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGc                 <  1:2685666/51‑1 (MQ=255)                               tgaaattgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGc                 >  1:6440797/1‑51 (MQ=255)                                gaaattgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGcc                <  1:2916990/51‑1 (MQ=255)                                 aaattgaaCCTGTGCTTCCCGCAGAGTTTAACGCTTAGGGTGATGTTGCCg               >  1:2059101/1‑51 (MQ=255)                                 aaattgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCg               >  1:5960734/1‑51 (MQ=255)                                   attgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCt             >  1:1179612/1‑51 (MQ=255)                                    ttgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCtt            >  1:624030/1‑51 (MQ=255)                                    ttgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCtt            <  1:2214484/51‑1 (MQ=255)                                    ttgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCtt            >  1:5914917/1‑51 (MQ=255)                                     tgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCttt           <  1:3982730/51‑1 (MQ=255)                                     tgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCttt           >  1:4670319/1‑51 (MQ=255)                                     tgaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCttt           <  1:6365423/51‑1 (MQ=255)                                      gaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCtttt          <  1:6396651/51‑1 (MQ=255)                                      gaaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCtttt          >  1:5867204/1‑51 (MQ=255)                                       aaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTa         <  1:5931168/51‑1 (MQ=255)                                       aaCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTa         <  1:708607/51‑1 (MQ=255)                                        aCCTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTAt        >  1:2502188/1‑51 (MQ=255)                                          cTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATgg      >  1:5581788/1‑51 (MQ=255)                                          cTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATgg      >  1:3498679/1‑51 (MQ=255)                                          cTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATgg      <  1:2031879/51‑1 (MQ=255)                                          cTGTGCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATgg      >  1:772045/1‑51 (MQ=255)                                           tgtgCTTCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGc     >  1:5892876/1‑51 (MQ=255)                                             tgcttcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGa   >  1:411193/1‑51 (MQ=255)                                              gcttcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGaa  <  1:2601280/51‑1 (MQ=255)                                              gcttcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGaa  >  1:3750458/1‑51 (MQ=255)                                                 |                                                AACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAA  >  REL606/2402‑2497

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