Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,638	A→C	F133C (TTC→TGC)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,638	0	A	C	100.0%	78.0 / NA	21	F133C (TTC→TGC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (11/10);  total (11/10)

TGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCG  >  REL606/3497567‑3497709                                                                        |                                                                        tGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGgcagca                                                                     >  2:765288/1‑76 (MQ=255) tGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGgcagca                                                                     <  2:488472/76‑1 (MQ=255)   aaCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaag                                                                   >  1:756524/1‑76 (MQ=255)       tACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAgccgc                                                               >  1:341917/1‑76 (MQ=255)          cgaaGTCATAGCGGATGCCCTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAg                                                            >  2:573498/4‑76 (MQ=255)           gTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAgc                                                           >  2:432670/1‑75 (MQ=255)                      ggATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGcc                                                >  2:468619/1‑76 (MQ=255)                          gCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCg                                            <  2:457096/76‑1 (MQ=255)                           ccGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCNGACGTCCCAGGCCACCGt                                           >  1:619428/1‑76 (MQ=255)                             gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                         >  2:587370/1‑76 (MQ=255)                               aaCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAcc                                       >  1:340588/1‑76 (MQ=255)                                    cgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGc                                  <  1:511129/76‑1 (MQ=255)                                           gCACCCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATc                           <  1:489300/76‑1 (MQ=255)                                            cAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGATGCCCGGATCa                          <  1:653624/76‑1 (MQ=255)                                            cAACCCAAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCa                          >  1:121647/1‑76 (MQ=255)                                                   aaCGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTttctt                   >  2:416906/1‑76 (MQ=255)                                                        cGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTtcatca              <  2:472329/76‑1 (MQ=255)                                                             gagaATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAg         <  2:756520/76‑1 (MQ=255)                                                               gaATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGct       <  1:587370/76‑1 (MQ=255)                                                                aaTCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGctc      >  1:61870/1‑76 (MQ=255)                                                                   cAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTc   <  1:147573/76‑1 (MQ=255)                                                                    aaGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCg  >  1:293337/1‑76 (MQ=255)                                                                    aaGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCg  <  1:444678/76‑1 (MQ=255)                                                                    aaGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCg  <  2:743290/76‑1 (MQ=255)                                                                        |                                                                        TGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCG  >  REL606/3497567‑3497709

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