Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,333	T→G	S235R (AGT→CGT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,333	0	T	G	98.1%	197.8 / ‑3.4	53	S235R (AGT→CGT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (15/37);  total (15/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

TTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAT  >  REL606/3497260‑3497405                                                                          |                                                                         ttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTTGCACGcc                                                                        >  1:342669/1‑76 (MQ=255) ttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTTGCACGcc                                                                        >  2:906178/1‑76 (MQ=255) ttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCAcccc                                                                        >  2:906157/1‑76 (MQ=255) ttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGcc                                                                        <  1:366490/76‑1 (MQ=255)  tATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCca                                                                       <  1:1275392/76‑1 (MQ=255)  tATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCca                                                                       <  1:1393239/76‑1 (MQ=255)      ttCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACaaa                                                                   <  2:2203164/76‑1 (MQ=255)       tCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAAc                                                                  <  2:220970/76‑1 (MQ=255)        ccACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACGCCACcaaag                                                                 >  1:1643189/1‑74 (MQ=255)        ccACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACg                                                                 <  2:867329/76‑1 (MQ=255)         cACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCCCGCCACAAACGc                                                                >  1:1927561/1‑76 (MQ=255)          aCTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACCCCCCCAACGcc                                                               >  1:1126383/1‑75 (MQ=255)          aCTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGTGCACGCCCCAAACGCa                                                               >  1:1499441/1‑76 (MQ=255)           cTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCCCAAAAGCaa                                                              >  1:2076472/1‑76 (MQ=255)            tGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCGCCCCACAAACGCCAc                                                             >  2:2019586/1‑76 (MQ=255)            tGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAAc                                                             <  1:1861936/76‑1 (MQ=255)              ccGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACgt                                                           <  1:118248/76‑1 (MQ=255)               cgcgAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTTGCCCGCCGCAAACGCAAcccg                                                          >  1:248396/1‑73 (MQ=255)               cgcgAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACgtg                                                          >  2:1080712/1‑76 (MQ=255)                gcgAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTTGCCCGCCACCAAACCAAAgtgt                                                         >  2:90992/1‑76 (MQ=255)                 cgAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACGCCACAAACGCCACGCGtt                                                        >  2:1457377/1‑76 (MQ=255)                  gAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCCCGCCCCAAACCCAACGTGTTg                                                       >  1:1313165/1‑76 (MQ=255)                  gAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACGCCCCAAACGCAACGTGTTg                                                       >  1:1607990/1‑76 (MQ=255)                   aaGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACGCCCCCAAAGCAACGTGTTgg                                                      >  2:235665/1‑76 (MQ=255)                   aaGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACGCCCCAAACGCAACGCGTTgg                                                      >  2:1156347/1‑76 (MQ=255)                    agtagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGt                                                     >  1:1445451/1‑76 (MQ=255)                    agtagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGt                                                     <  2:1525227/76‑1 (MQ=255)                    agtagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACCCAACGGGTTGGt                                                     >  1:309832/1‑76 (MQ=255)                     gtagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTc                                                    >  2:2049101/1‑76 (MQ=255)                      tagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCCCAAACGCAACGGGTTGGGcg                                                   >  1:1622849/1‑76 (MQ=255)                      tagtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGGcg                                                   >  1:1900324/1‑76 (MQ=255)                       agtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTcgc                                                  >  1:2063081/1‑76 (MQ=255)                       agtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTcgc                                                  <  1:523215/76‑1 (MQ=255)                       agtCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTcgc                                                  <  1:523233/76‑1 (MQ=255)                         tCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGTGCACGCCACAAACGCAACGTGTTGGTCGCGt                                                <  2:630002/76‑1 (MQ=255)                         tCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGGGTTGGGcgcgg                                                >  1:1123596/1‑75 (MQ=255)                         tCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGCGTTGGTcgcgg                                                >  2:920289/1‑75 (MQ=255)                          cACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTc                                               <  2:798458/76‑1 (MQ=255)                           aCCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCg                                              <  2:1054680/76‑1 (MQ=255)                           aCCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCg                                              >  2:1080521/1‑76 (MQ=255)                            cccTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCgg                                             <  2:1575753/76‑1 (MQ=255)                            cccTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCgg                                             <  2:1447365/76‑1 (MQ=255)                             ccTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACGCCCCAAACGCAACGTGTTGGTCGCGGCGGt                                            >  1:234533/1‑76 (MQ=255)                             ccTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGt                                            <  1:1705038/76‑1 (MQ=255)                             ccTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGt                                            <  2:261735/76‑1 (MQ=255)                              cTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCggga                                           >  1:308696/1‑74 (MQ=255)                               tGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTAt                                          <  2:1305997/76‑1 (MQ=255)                                ggTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACTCAACGTGTTGGTCGCGTCGGTATc                                         <  2:656253/76‑1 (MQ=255)                                ggTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATc                                         <  2:1715762/76‑1 (MQ=255)                                ggTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATc                                         <  2:1389454/76‑1 (MQ=255)                                   tGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCCCAAAAACAAAGCGTTGGGCGCGGCGGGATCGTa                                      >  1:895326/1‑76 (MQ=255)                                   tGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTa                                      >  2:1883676/1‑76 (MQ=255)                                    gAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTaa                                     <  2:245526/76‑1 (MQ=255)                                     aaTTTACCGAGGTTAATTTCGCTACTGGCTTTCGCACGCCCCAAACGCCACGGGTTGGTCGCGGCGGGGTCGTAAc                                    >  1:710218/1‑76 (MQ=255)                                      aTTTACCGAGGTTAATTTCGCTACTGGCTTGCGCCCGCCACAAACGCAACGGGTTGGGCGCGGCGGGATCGTAAcc                                   >  1:1012457/1‑76 (MQ=255)                                      aTTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAAcc                                   <  2:1633678/76‑1 (MQ=255)                                       tttACCGAGGTTAATTTCGCTACTGGCTTGCGCCCGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCa                                  >  2:2210955/1‑76 (MQ=255)                                        ttACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAAAGCAACGTGTTGGTCGCGTCGGGATCGTAACCAg                                 >  1:1664477/1‑76 (MQ=255)                                         tACCGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAgg                                <  1:397096/76‑1 (MQ=255)                                           ccGAGGTTAATTTCGCTACTGGCTTGCGCACGCCCCAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGa                              >  1:1971611/1‑76 (MQ=255)                                            cGAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGAt                             >  2:2181999/1‑76 (MQ=255)                                             gAGGTTAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGAtt                            <  2:1078740/76‑1 (MQ=255)                                              aGGTTAATTTCGCTACTGGCTTTCGCACGCCCCAAACGCAACGTGTTGGTCGCGGCGGTATCGTAACCAGGGGTTa                           >  2:2048200/1‑76 (MQ=255)                                                 ttAATTTCGCTACTGGCTTGCGCACGCCACAAACGCCACGGGTTGGTCGCGGCGGGATCGTAACCAGGGATTATCt                        >  1:6755/1‑76 (MQ=255)                                                 ttAATTTCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGCCGCGTCGGTATCGTAACCAGGGATTATCt                        <  1:1287951/76‑1 (MQ=255)                                                  tAATTTCGCTACTGGCTTGCGCACGCCCCAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATTTg                       >  2:2065023/1‑76 (MQ=255)                                                     tttCGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATc                    <  1:153229/76‑1 (MQ=255)                                                        cGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTa                 <  1:49295/76‑1 (MQ=255)                                                        cGCTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTa                 <  2:755258/76‑1 (MQ=255)                                                          cTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAg               <  1:1911623/76‑1 (MQ=255)                                                          cTACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGGCGGGATCGTAACCAGGGAGTATCTGATCGGAAg               >  2:1775832/1‑76 (MQ=255)                                                           tACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGc              <  2:1040817/76‑1 (MQ=255)                                                           tACTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGc              <  2:255908/76‑1 (MQ=255)                                                            aCTGGCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCg             <  2:1605702/76‑1 (MQ=255)                                                               ggCTTGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACt          >  2:1529983/1‑76 (MQ=255)                                                                   tGCGCACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCa      >  2:1988824/1‑76 (MQ=255)                                                                    gcgcACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAg     <  1:1733518/76‑1 (MQ=255)                                                                     cgcACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGa    <  2:782179/76‑1 (MQ=255)                                                                      gcACGCCACAACCGCAACGTTTTGTTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGaa   <  1:1036224/76‑1 (MQ=255)                                                                      gcACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGaa   <  1:933609/76‑1 (MQ=255)                                                                      gcACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGaa   >  2:246662/1‑76 (MQ=255)                                                                      gcACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGaa   <  2:1457298/76‑1 (MQ=255)                                                                       cACGCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAt  >  2:1198989/1‑76 (MQ=255)                                                                          |                                                                         TTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAT  >  REL606/3497260‑3497405

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