Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,693	C→T	E115K (GAA→AAA)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,693	0	C	T	96.3%	195.7 / ‑1.5	54	E115K (GAA→AAA)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (23/29);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGC  >  REL606/3497621‑3497763                                                                         |                                                                       tCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTttatc                                                                     <  1:820506/76‑1 (MQ=255)     cAGAGAATCAAGGAAGCAAGCCGCCAGTCGTCCCAGGCCGCCGTTACCGAGGCCCGGATCATTTTCTTTatcaatc                                                                 <  1:582248/76‑1 (MQ=255)     cAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTatcaatc                                                                 <  2:188424/76‑1 (MQ=255)     cAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTatcaatc                                                                 >  2:492693/1‑76 (MQ=255)      cgagaACCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTatcaatca                                                                <  2:706033/75‑1 (MQ=255)      agagAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTatcaatca                                                                <  1:1404089/76‑1 (MQ=255)           aTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGctct                                                           <  1:503604/76‑1 (MQ=255)           aTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGctct                                                           <  2:1160621/76‑1 (MQ=255)            tCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCtt                                                          <  1:1445448/76‑1 (MQ=255)             cAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTc                                                         <  2:1576128/76‑1 (MQ=255)              aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCg                                                        <  2:1121782/76‑1 (MQ=255)              aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCg                                                        <  1:120639/76‑1 (MQ=255)              aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCg                                                        <  2:692644/76‑1 (MQ=255)               aGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCgg                                                       >  2:567413/1‑75 (MQ=255)               aGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCga                                                       <  2:1358028/76‑1 (MQ=255)               aGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCga                                                       >  2:1029532/1‑76 (MQ=255)               aGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCga                                                       >  2:1105298/1‑76 (MQ=255)                ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCgag                                                      >  1:339829/1‑76 (MQ=255)                ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCgag                                                      <  2:22288/76‑1 (MQ=255)                ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCACTTTCTTTTTCAATTAGCTCTTccgt                                                      >  1:1286021/1‑73 (MQ=25)                     caagcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTa                                                 <  2:89170/76‑1 (MQ=255)                      aagcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTaa                                                >  1:1452933/1‑76 (MQ=255)                          cgccAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAAcccc                                            >  2:151743/1‑76 (MQ=255)                           gccAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCa                                           >  1:206824/1‑76 (MQ=255)                           gccAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCa                                           <  1:655690/76‑1 (MQ=255)                           gccAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCa                                           <  2:389209/76‑1 (MQ=255)                            ccAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCAt                                          <  2:255368/76‑1 (MQ=255)                              aGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCg                                        <  1:868373/76‑1 (MQ=255)                              aGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCg                                        >  2:1134334/1‑76 (MQ=255)                                  gTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTc                                    >  1:1031522/1‑76 (MQ=255)                                  gTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTc                                    >  1:1071626/1‑76 (MQ=255)                                    cccAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                  >  1:279485/1‑76 (MQ=255)                                        ggCCACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTgc                              <  1:1377951/76‑1 (MQ=255)                                           cACCGTTTCCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGcc                           >  2:1320863/1‑76 (MQ=255)                                           cACCGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGcc                           >  2:1219671/1‑76 (MQ=255)                                             ccGTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCt                         <  2:1286433/76‑1 (MQ=255)                                              ccttACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTg                        <  2:135595/74‑1 (MQ=255)                                              cGTGACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTg                        <  2:912159/76‑1 (MQ=255)                                               gTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGt                       <  1:1022651/76‑1 (MQ=255)                                               gTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGt                       <  1:911025/76‑1 (MQ=255)                                               gTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGt                       >  1:195981/1‑76 (MQ=255)                                               gTTACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGt                       <  1:811716/76‑1 (MQ=255)                                                 tACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTa                      >  2:465426/1‑75 (MQ=255)                                                 tACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTa                      <  1:465427/75‑1 (MQ=255)                                                 tACCGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTAc                     >  1:616079/1‑76 (MQ=255)                                                    cGAGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATc                  >  2:276964/1‑76 (MQ=255)                                                      aGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCtt                >  2:381566/1‑76 (MQ=255)                                                      aGGCCCGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCtt                <  2:691503/76‑1 (MQ=255)                                                          ccGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTa            <  2:1435187/76‑1 (MQ=255)                                                           cGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCGCATCGCTTCCAGTGCGCCCTGTACATCTTCGTaa           >  1:1317602/1‑76 (MQ=255)                                                           cGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTaa           >  1:1188052/1‑76 (MQ=255)                                                           cGGATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTaa           <  1:1115612/76‑1 (MQ=255)                                                             gATCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAAt         <  2:1452751/76‑1 (MQ=255)                                                              aTCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAAtt        >  1:556453/1‑76 (MQ=255)                                                              aTCATTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCACCGCTTCCAGCGCGCCCCGTACATCTTCGTAAAtt        >  1:950355/1‑76 (MQ=255)                                                                 aTTTTCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCt     <  2:428185/76‑1 (MQ=255)                                                                  ttttCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTa    >  1:1134358/1‑76 (MQ=255)                                                                   tttCTTTATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAg   >  2:1106302/1‑76 (MQ=255)                                                                    ttctttATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGc  >  2:137610/1‑76 (MQ=255)                                                                         |                                                                       TCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGC  >  REL606/3497621‑3497763

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