Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,945	T→G	K31Q (AAG→CAG)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,945	0	T	G	98.2%	222.2 / ‑4.0	54	K31Q (AAG→CAG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (25/28);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

CCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATA  >  REL606/3497873‑3498017                                                                         |                                                                         ccACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTcc                                                                       <  1:1348098/76‑1 (MQ=255)    cGAGACGATCGCGCACAGGAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAAt                                                                    <  1:597919/76‑1 (MQ=255)         cGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTaa                                                               >  2:1363088/1‑76 (MQ=255)          gATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTaaa                                                              >  1:1672798/1‑76 (MQ=255)           aTCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAAc                                                             <  1:840640/76‑1 (MQ=255)            tCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACa                                                            <  1:1301469/76‑1 (MQ=255)            tCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACa                                                            <  2:1782985/76‑1 (MQ=255)              gcgcACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATc                                                          <  1:1004753/76‑1 (MQ=255)              gcgcACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATc                                                          <  1:1758808/76‑1 (MQ=255)                 cacaGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGc                                                       >  2:2020272/1‑76 (MQ=255)                 cacaGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGc                                                       <  2:1634247/76‑1 (MQ=255)                    aGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTg                                                    <  2:797274/76‑1 (MQ=255)                      cAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTa                                                  <  1:1680370/76‑1 (MQ=255)                       aaaTAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTaa                                                 <  2:775482/76‑1 (MQ=255)                       aaaTAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTaa                                                 >  1:347865/1‑76 (MQ=255)                          tAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAgag                                              <  1:1417918/76‑3 (MQ=255)                          tAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCg                                              <  1:2051978/76‑1 (MQ=255)                             cGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATa                                           <  2:1341910/76‑1 (MQ=255)                              gTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATag                                          <  2:636774/76‑1 (MQ=255)                               ttGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATaga                                         >  2:574183/1‑76 (MQ=255)                                tGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATaggg                                        >  2:1069835/1‑74 (MQ=255)                                  cGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTg                                      <  1:200016/76‑1 (MQ=255)                                     tCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCtt                                   <  1:2028194/76‑1 (MQ=255)                                     tCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCtt                                   >  1:1346711/1‑76 (MQ=255)                                      cAGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTa                                  <  2:874448/76‑1 (MQ=255)                                       aGCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTaa                                 <  1:1422652/76‑1 (MQ=255)                                        gCCATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAg                                <  1:2105918/76‑1 (MQ=255)                                          cATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAggg                              >  1:566717/1‑76 (MQ=255)                                          cATTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAggg                              <  2:581299/76‑1 (MQ=255)                                           aTTCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGc                             >  2:382007/1‑76 (MQ=255)                                            ttCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGGGAAAGAGTGCTTAAGGGCt                            >  1:1690887/1‑76 (MQ=255)                                             tCATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCtt                           >  2:635725/1‑76 (MQ=255)                                              cATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTc                          >  1:609317/1‑76 (MQ=255)                                              cATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTc                          >  2:153528/1‑76 (MQ=255)                                              cATGTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTc                          >  1:1496088/1‑76 (MQ=255)                                                 gTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGCGTGCTTAAGGGCTTCTAc                       >  2:455507/1‑76 (MQ=255)                                                 gTTTATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTAc                       >  1:1231853/1‑76 (MQ=255)                                                  tttATTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACg                      >  1:1971566/1‑76 (MQ=255)                                                     aTTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTa                   >  1:166024/1‑76 (MQ=255)                                                     aTTGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTa                   <  2:1731195/76‑1 (MQ=255)                                                      ttGGCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTaa                  <  2:1109105/76‑1 (MQ=255)                                                        ggCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGc                <  2:2091102/76‑1 (MQ=255)                                                        ggCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGc                <  2:164439/76‑1 (MQ=255)                                                         gCGACGACCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCg               <  2:615671/76‑1 (MQ=255)                                                              gacCGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGc          >  2:1611683/1‑76 (MQ=255)                                                                 cGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGAt       <  2:225228/76‑1 (MQ=255)                                                                 cGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGAt       >  2:1376411/1‑76 (MQ=255)                                                                 cGGGTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGAt       <  1:154310/76‑1 (MQ=255)                                                                   ggTCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGa     >  2:1194652/1‑76 (MQ=255)                                                                     tCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAat   >  2:1497626/1‑76 (MQ=255)                                                                     tCCTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAat   >  2:1332536/1‑76 (MQ=255)                                                                      ccTTTCACAGCGTAAACATCAGCTTGTAAGCGATAGATTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAata  <  2:985100/76‑1 (MQ=255)                                                                      ccTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAata  >  1:1647764/1‑76 (MQ=255)                                                                      ccTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAata  >  1:1399151/1‑76 (MQ=255)                                                                      ccTGTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAata  <  2:985109/76‑1 (MQ=255)                                                                         |                                                                         CCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGTAAGCGATAGAGTGCTTAAGGGCTTCTACGCTAAGCGTGGGCGATGAATA  >  REL606/3497873‑3498017

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