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	98.3%	217.3 / ‑5.0	58	F133L (TTC→CTC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (22/35);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGAG  >  REL606/3497567‑3497711                                                                         |                                                                         tGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACNTCGCCAGAGAATCAAGGaggca                                                                       <  2:257859/76‑1 (MQ=255)     cATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAgcc                                                                   <  1:1392124/76‑1 (MQ=255)      aTACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAggcg                                                                  >  1:632636/1‑76 (MQ=255)        cccGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAgccgcc                                                                <  1:1587412/75‑1 (MQ=255)        aCCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAgccgcc                                                                >  1:1604824/1‑76 (MQ=255)          cGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAg                                                              >  1:1563528/1‑76 (MQ=255)           gTAGTCATAGCGGATGCCGTAACCGCCCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGa                                                             <  1:1527407/76‑1 (MQ=255)            tAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGAc                                                            <  2:629782/76‑1 (MQ=255)            tAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGAc                                                            >  1:1752112/1‑76 (MQ=255)              gTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGt                                                          <  1:202677/76‑1 (MQ=255)                cATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTcc                                                        >  2:566044/1‑76 (MQ=255)                cATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTcc                                                        >  1:279265/1‑76 (MQ=255)                  tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGCCGTCCCa                                                      >  2:791177/1‑76 (MQ=255)                  tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCa                                                      <  2:1713801/76‑1 (MQ=255)                  tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCa                                                      <  1:879360/76‑1 (MQ=255)                     cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGc                                                   <  2:844788/76‑1 (MQ=255)                     cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGc                                                   >  2:1031077/1‑76 (MQ=255)                     cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAACCAAGGCGGCAAGCCGCCAGACGTCCCAGGc                                                   >  2:1245702/1‑76 (MQ=255)                        aTGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCAc                                                >  1:1634811/1‑76 (MQ=255)                         tGCCGTAGCCGGGCCCCGTCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCAcc                                               <  1:264128/76‑1 (MQ=38)                          gCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCg                                              >  1:873606/1‑76 (MQ=255)                            cGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAGGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGtt                                            >  2:1088893/1‑76 (MQ=255)                             gTACCCGCCCCCCGGCACCCCTACCGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                           <  1:626837/76‑1 (MQ=37)                             gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                           >  1:764575/1‑76 (MQ=255)                             gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                           >  1:764595/1‑76 (MQ=255)                             gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                           <  2:358108/76‑1 (MQ=255)                             gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                           >  2:71981/1‑76 (MQ=255)                             gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                           <  1:268617/76‑1 (MQ=255)                             gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGCTa                                           >  1:28466/1‑76 (MQ=255)                              tAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTAc                                          >  2:1818977/1‑76 (MQ=255)                              tAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTAc                                          >  2:1691550/1‑76 (MQ=255)                                  cgcgCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAg                                      <  1:745295/76‑1 (MQ=255)                                     gcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGcc                                   <  1:1808513/76‑1 (MQ=255)                                     gcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGcc                                   >  1:811863/1‑76 (MQ=255)                                       cccGGCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCg                                 <  1:1109125/76‑1 (MQ=255)                                          ggCAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGAt                              <  2:1461452/76‑1 (MQ=255)                                            cAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCa                            <  2:1844301/76‑1 (MQ=255)                                            cAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCa                            <  2:610621/76‑1 (MQ=255)                                            cAACCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCa                            <  2:777884/76‑1 (MQ=255)                                             aaCCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCAt                           <  2:839016/76‑1 (MQ=255)                                             aaCCCTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCAt                           <  1:1203280/76‑1 (MQ=255)                                               cccTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCAttt                         <  1:285569/76‑1 (MQ=255)                                                ccTAACGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCAtttt                        >  2:1587674/1‑76 (MQ=255)                                                   aaCGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTttctt                     >  1:1274256/1‑76 (MQ=255)                                                    aCGTCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTttcttc                    <  2:684337/76‑1 (MQ=255)                                                      ctCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTtcat                  <  2:661416/75‑1 (MQ=255)                                                       tCGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTtcatc                 <  1:302416/76‑1 (MQ=255)                                                       ccGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCA‑TTTCTtcatca                <  2:638037/75‑1 (MQ=37)                                                        cGCCAGAGAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTtcatca                >  2:6529/1‑76 (MQ=255)                                                        cGCCAGAGAAACAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTtcatca                <  2:1223374/76‑1 (MQ=255)                                                            agagAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca            <  1:1290418/76‑1 (MQ=255)                                                            agagAATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca            <  1:1180194/76‑1 (MQ=255)                                                             gagaATCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAg           <  2:959912/76‑1 (MQ=255)                                                                 aTCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGctct       <  2:1856180/76‑1 (MQ=255)                                                                 aTCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGctct       >  2:1190486/1‑76 (MQ=255)                                                                  tCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCtt      <  2:402432/76‑1 (MQ=255)                                                                  tCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCtt      <  2:1621091/76‑1 (MQ=255)                                                                  tCAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCtt      <  2:1358478/76‑1 (MQ=255)                                                                  acAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCtt      <  2:378634/75‑1 (MQ=255)                                                                   cAAGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTc     <  2:1638218/76‑1 (MQ=255)                                                                    aaGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCg    <  2:1222417/76‑1 (MQ=255)                                                                    aaGGAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCg    <  1:1592018/76‑1 (MQ=255)                                                                     aggaggCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCga   >  2:1335263/1‑76 (MQ=255)                                                                      ggAGGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCgag  <  2:1518260/76‑1 (MQ=255)                                                                         |                                                                         TGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGAG  >  REL606/3497567‑3497711

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