Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,498,251	G→A	L413F (CTT→TTT)	glgA ←	glycogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,498,251	0	G	A	100.0%	167.5 / NA	47	L413F (CTT→TTT)	glgA	glycogen synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/28);  total (19/28)

CGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTA  >  REL606/3498181‑3498323                                                                       |                                                                         cGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTtct                                                                     >  2:124694/1‑76 (MQ=255)  gAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACCNACNCACTGGCGACGCCANCTGCAAAGTtctc                                                                    >  2:1451939/1‑76 (MQ=255)  gAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTtctc                                                                    <  2:767494/76‑1 (MQ=255)   aaTAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCa                                                                   >  2:1416841/1‑76 (MQ=255)     tAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAag                                                                 >  2:137452/1‑76 (MQ=255)      aGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAaga                                                                <  1:443882/76‑1 (MQ=255)       gCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAagag                                                               <  1:714884/76‑1 (MQ=255)           ttAACAGCGTCCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACa                                                           <  2:555366/75‑1 (MQ=255)           gTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACa                                                           <  1:1164297/76‑1 (MQ=255)            tAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGCCAACGCCATCTGCAAGGTTCCCAAGAGACCAg                                                          >  2:957298/1‑76 (MQ=255)             naCAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTNCTCAAGAGAACAGt                                                         >  2:231493/2‑76 (MQ=255)             aaCAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGt                                                         >  1:1080816/1‑76 (MQ=255)             aaCAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGt                                                         >  2:1033469/1‑76 (MQ=255)              aCAGCGACCAGGCANTACTATCTTCAAATACAAANCCACTTGCNACGCCATCTGCAAAGTTCTAAAGAGAACAgcc                                                        <  2:1331835/76‑3 (MQ=255)               cAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCa                                                       <  1:1048475/76‑1 (MQ=255)               cAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCa                                                       <  1:559640/76‑1 (MQ=255)                aGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAg                                                      <  1:933212/76‑1 (MQ=255)                   gACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGaaa                                                   >  1:1147615/1‑76 (MQ=255)                        ggCATTGCTATCTTCAAAGACAAACCCACTGGCGACGCCACCTGCGAAGTTCTCAAGAGAACAGTCAGAAACCGTa                                              <  1:1483313/76‑1 (MQ=255)                         gCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTAt                                             <  1:550397/76‑1 (MQ=255)                           aTTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCa                                           <  1:666578/76‑1 (MQ=255)                           aTTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCa                                           <  1:472414/76‑1 (MQ=255)                            ttACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAg                                          >  2:663862/1‑76 (MQ=255)                                tATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAg                                      <  2:566478/76‑1 (MQ=255)                                 aTCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCTGTAAGc                                     <  1:96450/76‑1 (MQ=255)                                 aTCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGc                                     >  2:1037876/1‑76 (MQ=255)                                 aTCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGc                                     >  2:691680/1‑76 (MQ=255)                                  tCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGcc                                    <  1:985377/76‑1 (MQ=255)                                   cTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGccc                                   >  1:1171789/1‑76 (MQ=255)                                      cAAAGACAAACCCACTGGCGACGCCATCTGCAAAGTNCTCAACACAATAGACGGAAACCGTATCAGCAAGCNCAcc                                >  1:182814/1‑76 (MQ=255)                                       aaaGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCg                               >  2:808233/1‑76 (MQ=255)                                         aGACAAACCCACTGGCGACGCCATTTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGt                             >  2:792979/1‑76 (MQ=255)                                         aGACAAACCCACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGt                             <  2:12835/76‑1 (MQ=255)                                                 ccACTGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCa                     <  2:172305/76‑1 (MQ=255)                                                     tGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTaa                 <  1:685550/76‑1 (MQ=255)                                                     tGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTaa                 <  2:1087884/76‑1 (MQ=255)                                                     tGGCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTaa                 >  1:474553/1‑76 (MQ=255)                                                      ggCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAAc                <  1:1023150/76‑1 (MQ=255)                                                      ggCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAAc                >  1:1208902/1‑76 (MQ=255)                                                      ggCGACGCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGGATCAGCAAGCCCACCGGTGCGCCGCACTAAc                <  1:1133945/76‑1 (MQ=255)                                                          aCGCCATCTGCANAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCa            >  1:18127/1‑76 (MQ=255)                                                            gCCATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGc          >  2:180066/1‑76 (MQ=255)                                                             ccATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcg         >  1:1093053/1‑76 (MQ=255)                                                              cATCTGCAAAGTTCTCAGGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgt        <  1:177517/76‑1 (MQ=255)                                                              cATCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgt        <  1:1343846/76‑1 (MQ=255)                                                               aTCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgta       <  2:474553/76‑1 (MQ=255)                                                               aTCTGCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgta       <  2:654654/76‑1 (MQ=255)                                                                   gCAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgtaccgt   <  2:283812/76‑1 (MQ=255)                                                                    cAAAGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGcgtaccgta  <  2:172510/76‑1 (MQ=255)                                                                       |                                                                         CGAATAGCCCGTAACAGCGACCAGGCATTACTATCTTCAAAGACAAACCCACTGGCGACGCCATCTGCAAGGTTCTCAAGAGAACAGTCAGAAACCGTATCAGCAAGCCCACCGGTGCGCCGCACTAACGGCAGCGTACCGTA  >  REL606/3498181‑3498323

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