Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,984	T→G	K18Q (AAG→CAG)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,984	0	T	G	100.0%	181.8 / NA	44	K18Q (AAG→CAG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (19/25);  total (19/25)

AGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATT  >  REL606/3497911‑3498058                                                                          |                                                                           aGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGaa                                                                          <  1:1066251/76‑1 (MQ=255) aGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGaa                                                                          >  1:845840/1‑76 (MQ=255)     aTTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGc                                                                      >  2:1287241/1‑76 (MQ=255)       tCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCtt                                                                    >  2:859026/1‑76 (MQ=255)        cATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTc                                                                   >  2:13762/1‑76 (MQ=255)         aTGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCt                                                                  <  1:1441827/76‑1 (MQ=255)         aTGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCt                                                                  <  1:576575/76‑1 (MQ=255)          tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTa                                                                 <  2:1688794/76‑1 (MQ=255)          tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTa                                                                 <  2:198610/76‑1 (MQ=255)           gTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTAc                                                                >  1:204518/1‑76 (MQ=255)            tttATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACg                                                               >  1:368082/1‑76 (MQ=255)                        gacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGc                                                   <  2:1475805/76‑1 (MQ=255)                        gacCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGc                                                   >  2:1524555/1‑76 (MQ=255)                          ccGGGTCCTTTCCAATCATAAACACCAGCTTGTAAGCGATAGAGTGCAGAAGGGCTTCTACGCTAAGCGTGGGCGa                                                 <  2:572836/76‑1 (MQ=255)                             ggTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGa                                              <  1:1350013/76‑1 (MQ=255)                              gTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGaa                                             >  1:1608723/1‑76 (MQ=255)                                ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAata                                           <  1:948308/76‑1 (MQ=255)                                 cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGTGATGAatat                                          >  2:799734/1‑76 (MQ=255)                                 cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAatat                                          <  2:562362/76‑1 (MQ=255)                                 cTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAatat                                          <  1:1362586/76‑1 (MQ=255)                                    tCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTa                                       <  1:32922/76‑1 (MQ=255)                                    tCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTa                                       >  1:1229849/1‑76 (MQ=255)                                      cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaaa                                     <  1:946412/76‑1 (MQ=255)                                      cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaaa                                     <  2:1096380/76‑1 (MQ=255)                                         tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACgg                                  <  1:468231/76‑1 (MQ=255)                                          cGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGa                                 >  1:1427947/1‑76 (MQ=255)                                          cGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGa                                 <  2:1330483/76‑1 (MQ=255)                                          cGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGa                                 >  2:1598818/1‑76 (MQ=255)                                          cGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGa                                 >  1:291603/1‑76 (MQ=255)                                           gTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAg                                >  2:1627158/1‑76 (MQ=255)                                             aaaCATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCa                              >  1:1123927/1‑76 (MQ=255)                                             aaaCATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCa                              <  1:516574/76‑1 (MQ=255)                                             aaaCATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCa                              <  1:167208/76‑1 (MQ=255)                                                   cAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCata                        <  2:1171402/76‑1 (MQ=255)                                                    aGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCatat                       <  2:420426/76‑1 (MQ=255)                                                    aGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCatat                       <  1:722101/76‑1 (MQ=255)                                                        tGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGc                   >  2:805708/1‑76 (MQ=255)                                                          tAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGt                 <  2:219445/76‑1 (MQ=255)                                                             gCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTc              <  1:335802/76‑1 (MQ=255)                                                               gATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCt            >  2:1157516/1‑76 (MQ=255)                                                                  agagTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa         >  1:987900/1‑76 (MQ=255)                                                                  agagTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa         <  2:1704095/76‑1 (MQ=255)                                                                    agTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaaa       <  1:887419/76‑1 (MQ=255)                                                                      tGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACt     >  2:37316/1‑76 (MQ=255)                                                                         ttAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAtt  >  2:1455072/1‑76 (MQ=255)                                                                          |                                                                           AGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATT  >  REL606/3497911‑3498058

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