Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,639	A→G	F133L (TTC→CTC)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,639	0	A	G	100.0%	164.2 / NA	44	F133L (TTC→CTC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/26);  total (18/26)

TTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGA  >  REL606/3497566‑3497710                                                                          |                                                                        ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGTGCACCGGCAACNCTAACGTCGCCAGAGAATCAAGGaggc                                                                       <  1:1498721/76‑1 (MQ=255)  tGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGaggca                                                                      <  1:549746/76‑1 (MQ=255)     aCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGaggcaagc                                                                   <  1:1243830/76‑1 (MQ=255)      naTACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAgcc                                                                  >  1:425794/2‑76 (MQ=255)       aTACCGTAGTCATAGCGGATGCCGTAANCGCGCCCCAGCCACCCTAACGTCGCCAGAGAATCAAGGAGGCAAgccg                                                                 >  1:1502402/1‑76 (MQ=255)         aCCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAgccgcc                                                               <  2:769182/76‑1 (MQ=255)          ccGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAAAAAATCAAGGAGGCAAgccgccg                                                              >  1:1119548/1‑75 (MQ=255)             tAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGAc                                                           >  2:881431/1‑76 (MQ=255)              aGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACg                                                          >  1:29263/1‑76 (MQ=255)              aGTCATAGCGGATGCCGTAACCGCGCCCCGGCAAACCTAACNTCGCCAGAGAATCAAGGANGCNAGCCGCCAGACg                                                          >  1:1396065/1‑76 (MQ=255)                  aTAGCGGATGCCGTAACCGCGCCCCGGCAACCCTCACGTCGNCTGAGAATCAAGGAGGCANGCNGCCAGACGTccc                                                      >  1:1392316/1‑76 (MQ=255)                  aTAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCGGCCAGACGTccc                                                      >  1:489185/1‑76 (MQ=255)                   tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCa                                                     <  1:1554110/76‑1 (MQ=255)                      cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGc                                                  >  1:956625/1‑76 (MQ=255)                       ggATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGcc                                                 >  2:1294220/1‑76 (MQ=255)                       ggATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGcc                                                 >  2:1149552/1‑76 (MQ=255)                        gATGCNGCAACCGCGCCCCGGCAACCCTACCGTCGCCAGAGAATCAAGGAGGCAAGCCGCCNGACGTCCCAGGCCa                                                >  2:1714087/1‑76 (MQ=255)                           gCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCg                                             <  1:609896/76‑1 (MQ=255)                            ccGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGt                                            >  1:1626441/1‑76 (MQ=255)                             cGTACCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGtt                                           <  1:997076/76‑1 (MQ=255)                             cGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGtt                                           <  2:883826/76‑1 (MQ=255)                              gTACCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                          <  2:1658116/76‑1 (MQ=255)                              gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                          <  2:818773/76‑1 (MQ=255)                              gTAACCGCGCCCCGGCAACCCTAACGTCACCAGAGAATAAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGCTa                                          >  2:1092559/1‑76 (MQ=38)                               tAACCGCGCCCCGGCACCCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTAc                                         <  1:1306244/76‑1 (MQ=255)                                aaCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTAcc                                        >  2:1073040/1‑76 (MQ=255)                                 aCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCg                                       <  2:861338/76‑1 (MQ=255)                                   cgcgCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAg                                     <  1:306052/76‑1 (MQ=255)                                    gcgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCCGGCCACCGTTACCGAgg                                    <  1:1703162/76‑1 (MQ=255)                                    gcgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                    <  1:1104689/76‑1 (MQ=255)                                    gcgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                    <  1:1351492/76‑1 (MQ=255)                                         ccGGCACCCCTCACGCCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCgg                               <  2:640616/76‑1 (MQ=37)                                           ggCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGAt                             <  2:1599019/76‑1 (MQ=255)                                           ggCAACCCTAACGGCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGNCACCGTTCACNAGGGCCGGAt                             <  2:1688627/76‑1 (MQ=25)                                              aaCCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCAt                          <  2:1114776/76‑1 (MQ=255)                                                 ccTAACGTCGCCAGNGAATCAAGGAGGCAAGCCGCCAGACGTTCCAGGCCACCGTTACCGCGGGCCGGATCAtttt                       <  1:1469301/76‑1 (MQ=255)                                                    aaCGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTttctt                    >  2:763937/1‑76 (MQ=255)                                                          gCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaa              <  2:1238231/76‑1 (MQ=255)                                                            cAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatc            <  2:1471742/76‑1 (MQ=255)                                                            cAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatc            <  2:888355/76‑1 (MQ=255)                                                            cAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatc            <  2:493358/76‑1 (MQ=255)                                                             tgagaATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca           >  2:1102337/2‑76 (MQ=255)                                                             tgagaATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca           >  2:192327/2‑76 (MQ=255)                                                             agagAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca           <  2:956623/76‑1 (MQ=255)                                                              gagaATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAg          <  1:1216424/76‑1 (MQ=255)                                                               agaATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGc         <  1:1092561/76‑1 (MQ=255)                                                                gaATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGct        <  2:1291604/76‑1 (MQ=255)                                                                 aaTCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGctc       >  2:721050/1‑76 (MQ=255)                                                                  aTCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGctct      >  2:1495808/1‑76 (MQ=255)                                                                    cAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTc    >  2:1263804/1‑76 (MQ=255)                                                                     aaGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCATTCAGCTCTCCg   <  1:1574514/76‑1 (MQ=37)                                                                     aaGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCg   >  2:625656/1‑76 (MQ=255)                                                                      aggaggCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCga  >  2:504707/1‑76 (MQ=255)                                                                          |                                                                        TTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGA  >  REL606/3497566‑3497710

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