Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,496,783	C→T	G418D (GGT→GAT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,496,783	0	C	T	100.0%	168.5 / NA	43	G418D (GGT→GAT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/15);  total (28/15)

AGAGTGCAGTTCCGATACACCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGACCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGCAAGGTTT  >  REL606/3496710‑3496855                                                                          |                                                                         agagTGCAGTTCCGATACACCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCg                                                                        >  2:717485/1‑76 (MQ=255)   agTGCAGTTCCGATACACCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGtt                                                                      >  2:1394376/1‑76 (MQ=255)          ttCCGATACACCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATtca                                                               >  1:764040/1‑76 (MQ=255)          ttCCGATACACCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATtca                                                               >  1:1559980/1‑76 (MQ=255)           tCCGATACACCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATtcat                                                              <  2:1229434/76‑1 (MQ=255)             cGATACACCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATtcatca                                                            >  1:52434/1‑76 (MQ=255)                   acCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATc                                                      >  2:1563211/1‑76 (MQ=255)                   acCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATc                                                      >  2:990413/1‑76 (MQ=255)                    ccGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCg                                                     <  1:1356105/76‑1 (MQ=255)                    ccGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCg                                                     >  2:1204210/1‑76 (MQ=255)                    ccGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCg                                                     <  2:2913/76‑1 (MQ=255)                    ccGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCg                                                     >  1:9764/1‑76 (MQ=255)                     cGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCTATACGCACACGACGATCGTTGGATTCATCAATGATCGa                                                    <  2:564035/76‑1 (MQ=255)                     cGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGa                                                    <  1:321141/76‑1 (MQ=255)                     cGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGa                                                    >  1:682552/1‑76 (MQ=255)                     cGTTAACTTGGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGa                                                    <  2:767870/76‑1 (MQ=255)                      gTTAACTTTGTGGCCCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGAc                                                   >  2:143858/1‑76 (MQ=255)                            tttGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGt                                             >  1:197033/1‑76 (MQ=255)                            tttGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGt                                             >  2:327500/1‑76 (MQ=255)                            tttGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGt                                             >  1:159994/1‑76 (MQ=255)                              tgtgGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTcc                                           >  2:1077055/1‑76 (MQ=255)                                  gCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCcagc                                       >  2:176462/1‑76 (MQ=255)                                   cTCACAACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCcagca                                      <  2:1112756/76‑1 (MQ=255)                                     cacaACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGa                                    >  1:1092346/1‑76 (MQ=255)                                     cacaACAACCGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGa                                    >  2:97983/1‑76 (MQ=255)                                          cacccGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGt                               <  2:210070/73‑1 (MQ=255)                                             ccGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATc                            >  2:688071/1‑76 (MQ=255)                                              cGCCAGCCAGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCg                           >  1:84350/1‑76 (MQ=255)                                               gccagccaGGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGt                          <  1:802613/76‑1 (MQ=255)                                                    ccagGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGCTCGGTATCTTTCGGa                     >  2:132442/1‑76 (MQ=255)                                                    ccagGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGa                     >  2:132447/1‑76 (MQ=255)                                                     cagGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGAt                    >  2:1550472/1‑76 (MQ=255)                                                     cagGCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGAt                    >  2:1161330/1‑76 (MQ=255)                                                        gCCATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACt                 <  2:1477662/76‑1 (MQ=255)                                                          cATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGt               >  2:632591/1‑76 (MQ=255)                                                          cATACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGt               >  1:875239/1‑76 (MQ=255)                                                           aTACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGtt              <  1:1431966/76‑1 (MQ=255)                                                           aTACGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGtt              >  2:525057/1‑76 (MQ=255)                                                              cGCACACGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCt           <  1:92459/76‑1 (MQ=255)                                                                cacaCGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGc         <  1:273775/76‑1 (MQ=255)                                                                cacaCGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGc         >  2:1659818/1‑76 (MQ=255)                                                                cacaCGACGATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGc         <  2:1358283/76‑1 (MQ=255)                                                                   acgacgATCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGCAAg      <  2:1685746/76‑1 (MQ=255)                                                                     gacgatCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGCAAGGt    <  1:1310448/76‑1 (MQ=255)                                                                     gacgatCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGCAAGGt    >  2:1586317/1‑76 (MQ=255)                                                                       cgatCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGCAAGGttt  <  1:938925/76‑1 (MQ=255)                                                                       cgacCGTTGGATTCATCAAGGGTCGACGCCCGTCCCAGCAGATCGGTATCGTTCAAATACTGTTCCTGCATGGttt  <  2:1307691/76‑1 (MQ=255)                                                                          |                                                                         AGAGTGCAGTTCCGATACACCGTTAACTTTGTGGCTCACAACAACCGCCAGCCAGGCCATACGCACACGACGACCGTTGGATTCATCAATGATCGACGCCCGTCCCAGCAGATCGGTATCGTTCGGATACTGTTCCTGCAAGGTTT  >  REL606/3496710‑3496855

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