Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,383	T→A	Q218L (CAG→CTG)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,383	0	T	A	98.3%	230.1 / ‑3.0	59	Q218L (CAG→CTG)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (18/40);  total (18/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

AATTTCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGA  >  REL606/3497310‑3497454                                                                          |                                                                        aaTTTCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGa                                                                       <  2:1607566/76‑1 (MQ=255)  aTTTCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGAt                                                                      <  1:1304631/76‑1 (MQ=255)  aTTTCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGAt                                                                      <  1:913986/76‑1 (MQ=255)   tttCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATc                                                                     <  2:210741/76‑1 (MQ=255)     tCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGt                                                                   <  2:1042861/76‑1 (MQ=255)      cGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTa                                                                  <  2:873374/76‑1 (MQ=255)       gCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTaa                                                                 <  2:285965/76‑1 (MQ=255)       gCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTaa                                                                 <  1:761823/76‑1 (MQ=255)          aCTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCg                                                              <  2:76227/76‑1 (MQ=255)          aCTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCg                                                              <  2:584464/76‑1 (MQ=255)          aCTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCg                                                              <  1:1384010/76‑1 (MQ=255)           cTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGa                                                             <  1:458867/76‑1 (MQ=255)            tGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGCAAGCGAc                                                            >  1:1106837/1‑76 (MQ=255)               cTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTcc                                                         >  2:1834618/1‑76 (MQ=255)                 tGCGCACTCCACAAGCGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGAGCGTAAGCGACTCCCa                                                       <  2:1406308/76‑1 (MQ=255)                 tGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCa                                                       <  2:1418899/76‑1 (MQ=255)                 tGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCa                                                       <  2:672073/76‑1 (MQ=255)                      aCTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAAtc                                                  >  1:934581/1‑76 (MQ=255)                       cTCCACAACCGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAAtct                                                 <  1:713371/76‑1 (MQ=255)                       cTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAAtct                                                 >  2:1766972/1‑76 (MQ=255)                       cTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAAtct                                                 <  1:350905/76‑1 (MQ=255)                       cTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAAtct                                                 <  1:116668/76‑1 (MQ=255)                          cacaAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCtt                                              >  1:815068/1‑76 (MQ=255)                          cacaAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCtt                                              >  1:1289292/1‑76 (MQ=255)                          cacaAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCtt                                              <  2:1392345/76‑1 (MQ=255)                          cacaAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCtt                                              <  2:1392322/76‑1 (MQ=255)                          cacaAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCtt                                              <  2:267722/76‑1 (MQ=255)                          cacaAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCtt                                              >  2:896755/1‑76 (MQ=255)                             aaaCGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCgg                                           >  2:473122/1‑76 (MQ=255)                                cGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGttt                                        >  2:27198/1‑76 (MQ=255)                                  cAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCa                                      <  1:1336186/76‑1 (MQ=255)                                   aaCGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCaa                                     <  1:1166878/76‑1 (MQ=255)                                   aaCGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCaa                                     <  1:1143399/76‑1 (MQ=255)                                   aaCGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCaa                                     <  1:568622/76‑1 (MQ=255)                                      gtgtTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGTTTTCAATcc                                  <  2:1356291/76‑1 (MQ=255)                                      gtgtTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATcc                                  <  2:292708/76‑1 (MQ=255)                                       tgtTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCTACTCCCAGAATCTCTTCGGTTTCAATCCa                                 <  2:1084688/76‑1 (MQ=255)                                       tgtTGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCa                                 <  2:250961/76‑1 (MQ=255)                                         ttGGTCGCGTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAgc                               <  1:1171689/76‑1 (MQ=255)                                              cgcgTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtt                          <  1:692930/76‑1 (MQ=255)                                              cgcgTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtt                          <  2:847020/76‑1 (MQ=255)                                              cgcgTCGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGtt                          >  1:1886059/1‑76 (MQ=255)                                                   cGGTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGttttttt                     >  1:1585297/1‑76 (MQ=255)                                                    ggTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTa                    <  1:1093843/76‑1 (MQ=255)                                                    ggTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTa                    <  1:1834647/76‑1 (MQ=255)                                                     gTATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTAc                   >  2:1577872/1‑76 (MQ=255)                                                      tATCGTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTAcc                  >  2:297254/1‑76 (MQ=255)                                                      tATCGTAACCAGGAATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTAcc                  >  2:1176646/1‑76 (MQ=255)                                                          gTAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTcc              >  1:1359769/1‑76 (MQ=255)                                                           tAACCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCct             <  1:685441/76‑1 (MQ=255)                                                            aaCCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg            <  2:1094033/76‑1 (MQ=255)                                                            aaCCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg            <  1:549952/76‑1 (MQ=255)                                                            aaCCAGGGATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCctg            >  1:1851203/1‑76 (MQ=255)                                                                 gggATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGaa       <  1:1295846/76‑1 (MQ=255)                                                                 gggATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGaa       <  1:984990/76‑1 (MQ=255)                                                                  ggATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAAt      <  1:1576749/76‑1 (MQ=255)                                                                  agATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAAt      <  1:260536/75‑1 (MQ=255)                                                                   gATTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATg     <  2:899280/76‑1 (MQ=255)                                                                    aTTATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGCATGc    >  1:245614/1‑76 (MQ=255)                                                                     ttATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCg   >  1:1876094/1‑76 (MQ=255)                                                                      tATCAGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGa  >  1:669758/1‑76 (MQ=255)                                                                          |                                                                        AATTTCGCTACTGGCTTGCGCACTCCACAAACGCAACGTGTTGGTCGCGTCGGTATCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGGTTTCAATCCAGCGCGTTTTTTTACCTTCCTGCTGAATGCGA  >  REL606/3497310‑3497454

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