Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,498,013	G→A	S8L (TCA→TTA)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,498,013	0	G	A	98.9%	354.2 / ‑3.5	88	S8L (TCA→TTA)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (39/48);  total (39/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACG  >  REL606/3497940‑3498085                                                                          |                                                                         gTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATaaa                                                                        <  2:464327/76‑1 (MQ=255)  tCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGCGGGCGATAAat                                                                       >  1:864933/1‑76 (MQ=255)   ccTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAata                                                                      <  2:1357461/76‑1 (MQ=255)    cTTTCCAATCGTAAACATCAGCTTGTGAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAatat                                                                     <  2:1130845/76‑1 (MQ=255)      ttCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGt                                                                   >  2:516977/1‑76 (MQ=255)       tCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTa                                                                  <  2:751013/76‑1 (MQ=255)        ccAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTaa                                                                 <  1:771455/76‑1 (MQ=255)        ccAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTGAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTaa                                                                 <  2:1373923/76‑1 (MQ=255)         cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTaaa                                                                <  1:943333/76‑1 (MQ=255)         cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTaaa                                                                <  2:663261/76‑1 (MQ=255)         cAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTaaa                                                                <  2:763285/76‑1 (MQ=255)          aaTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAAc                                                               >  2:150834/1‑76 (MQ=255)           aTCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACg                                                              <  2:1697101/76‑1 (MQ=255)            tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACgg                                                             <  2:1832565/76‑1 (MQ=255)            tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACgg                                                             <  1:1887984/76‑1 (MQ=255)            tCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACgg                                                             <  2:421528/76‑1 (MQ=255)             cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGa                                                            <  2:85378/76‑1 (MQ=255)             cGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGa                                                            >  2:1946543/1‑76 (MQ=255)              gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAg                                                           <  2:1987400/76‑1 (MQ=255)              gTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAg                                                           >  1:63923/1‑76 (MQ=255)               tAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGc                                                          >  2:1020609/1‑76 (MQ=255)                aaaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCa                                                         <  2:1907229/76‑1 (MQ=255)                 aaCATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCAt                                                        <  2:1680655/76‑1 (MQ=255)                   cATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTc                                                      <  1:580480/76‑1 (MQ=255)                      cAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCata                                                   <  2:1840474/76‑1 (MQ=255)                        gCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCAtata                                                 <  1:1040872/76‑1 (MQ=255)                         cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAg                                                >  1:1425166/1‑76 (MQ=255)                         cTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAg                                                <  2:48304/76‑1 (MQ=255)                          ttGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAgg                                               <  2:1433807/76‑1 (MQ=255)                           tGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGc                                              <  1:26531/76‑1 (MQ=255)                           tGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGc                                              <  1:1445902/76‑1 (MQ=255)                           tGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGc                                              >  1:958935/1‑76 (MQ=255)                            gTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCg                                             >  1:1165715/1‑76 (MQ=255)                            gTAAGCGATAAAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCg                                             <  1:382391/76‑1 (MQ=255)                              aaGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGtt                                           >  2:1714543/1‑76 (MQ=255)                               aGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGttt                                          <  1:732450/76‑1 (MQ=255)                               aGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGttt                                          <  2:1038024/76‑1 (MQ=255)                                gCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTc                                         <  2:1551820/76‑1 (MQ=255)                                gCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTc                                         <  2:571428/76‑1 (MQ=255)                                    tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                     >  1:175359/1‑76 (MQ=255)                                    tAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                     >  1:1549633/1‑76 (MQ=255)                                     agagTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                    <  2:507577/76‑1 (MQ=255)                                     agagTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                    <  1:1947276/76‑1 (MQ=255)                                     agagTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                    >  2:628844/1‑76 (MQ=255)                                     agagTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                    >  2:109635/1‑76 (MQ=255)                                     agagTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                    >  1:1515727/1‑76 (MQ=255)                                       agTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaaa                                  >  1:811045/1‑76 (MQ=255)                                       agTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCCGaaaa                                  >  2:404526/1‑76 (MQ=255)                                         tGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACt                                >  2:1413347/1‑76 (MQ=255)                                         tGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACt                                <  1:522197/76‑1 (MQ=255)                                            ttAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAtt                             >  1:32198/1‑76 (MQ=255)                                            ttAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAtt                             >  1:307762/1‑76 (MQ=255)                                            ttAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAtt                             >  1:362861/1‑76 (MQ=255)                                            ttAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAtt                             >  1:1967347/1‑76 (MQ=255)                                            ttAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAtt                             >  1:1361634/1‑76 (MQ=255)                                             tAAGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                            >  1:16985/1‑76 (MQ=255)                                              aaGGGCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                           <  1:2044376/76‑1 (MQ=255)                                                gggCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa                         <  1:358003/76‑1 (MQ=255)                                                 ggCTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAg                        >  1:1184753/1‑76 (MQ=255)                                                   cTTCTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCg                      <  1:823121/76‑1 (MQ=255)                                                      cTACGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATa                   >  2:606334/1‑76 (MQ=255)                                                        aCGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGt                 >  2:53700/1‑76 (MQ=255)                                                         cGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTa                >  2:2064762/1‑76 (MQ=255)                                                         cGCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTa                >  2:393301/1‑76 (MQ=255)                                                          gCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaa               <  1:779882/76‑1 (MQ=255)                                                          gCTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaa               <  2:853762/76‑1 (MQ=255)                                                           cTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa              <  2:835907/76‑1 (MQ=255)                                                           cTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa              >  2:569141/1‑76 (MQ=255)                                                           cTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa              >  2:1262770/1‑76 (MQ=255)                                                           cTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa              <  1:1457074/76‑1 (MQ=255)                                                           cTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa              <  2:1202109/76‑1 (MQ=255)                                                           cTAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTaaa              <  1:2020232/76‑1 (MQ=255)                                                            tAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg             >  2:1309668/1‑76 (MQ=255)                                                            tAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg             >  2:829728/1‑76 (MQ=255)                                                            tAAGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg             >  2:735422/1‑76 (MQ=255)                                                             aaGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGc            >  2:89139/1‑76 (MQ=255)                                                             aaGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGc            >  1:1797326/1‑76 (MQ=255)                                                             aaGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGc            >  2:575577/1‑76 (MQ=255)                                                              aGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCt           <  2:918360/76‑1 (MQ=255)                                                              aGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCt           >  2:783331/1‑76 (MQ=255)                                                              aGCGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCt           <  2:1726211/76‑1 (MQ=255)                                                                cGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGcnta         <  2:1452834/76‑4 (MQ=255)                                                                cGTGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa         >  1:1803975/1‑76 (MQ=255)                                                                  tGGGCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg       <  1:2033059/76‑1 (MQ=255)                                                                   gggCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg      >  1:487186/1‑76 (MQ=255)                                                                     gCGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTa    <  2:468552/76‑1 (MQ=255)                                                                      cGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTAc   <  2:74090/76‑1 (MQ=255)                                                                      cGATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTAc   <  1:1177756/76‑1 (MQ=255)                                                                       gATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACg  <  1:490879/76‑1 (MQ=255)                                                                       gATAAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACg  <  2:1222002/76‑1 (MQ=255)                                                                          |                                                                         GTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACG  >  REL606/3497940‑3498085

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