Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,881	C→T	R52H (CGT→CAT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,881	0	C	T	98.1%	199.7 / ‑5.5	52	R52H (CGT→CAT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/27);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGGTAGTAAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAAT  >  REL606/3497809‑3497951                                                                         |                                                                       gACAGGTAGTAAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAgatgat                                                                     <  2:335534/76‑1 (MQ=255)   cAGGTAGTAAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATcg                                                                   <  2:1115537/76‑1 (MQ=255)       tagtaACCCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAGCCAGCGCTCCACGAGATGATCGCGca                                                               <  1:357302/76‑1 (MQ=255)       tagtaAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGca                                                               <  1:875270/76‑1 (MQ=255)         gtaAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGcaca                                                             <  1:1541655/76‑1 (MQ=255)          taAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAg                                                            >  2:791766/1‑76 (MQ=255)          taAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAg                                                            <  1:571943/76‑1 (MQ=255)          taAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAg                                                            >  1:411465/1‑76 (MQ=255)          taAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAg                                                            >  1:259319/1‑76 (MQ=255)           aaaCCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGc                                                           >  2:1154084/1‑76 (MQ=255)            aaCCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCa                                                          >  1:1654498/1‑76 (MQ=255)               cTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAAt                                                       >  1:811466/1‑76 (MQ=255)                 gACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATaa                                                     <  1:205518/76‑1 (MQ=255)                 gACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATaa                                                     <  1:990832/76‑1 (MQ=255)                 gACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATaa                                                     >  2:308979/1‑76 (MQ=255)                 gACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATaa                                                     <  2:1470646/76‑1 (MQ=255)                   cGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACg                                                   >  2:1144499/1‑76 (MQ=255)                     aGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGtt                                                 <  1:443104/76‑1 (MQ=255)                     aGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGtt                                                 >  1:1472418/1‑76 (MQ=255)                      gTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTg                                                <  1:243126/76‑1 (MQ=255)                       tttCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGc                                               >  2:1658674/1‑76 (MQ=255)                       tttCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGc                                               <  1:595857/76‑1 (MQ=255)                       tttCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGc                                               <  1:185870/76‑1 (MQ=255)                        ttCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCg                                              >  1:1470994/1‑76 (MQ=255)                            tGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCa                                          <  1:440844/76‑1 (MQ=255)                            tGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCa                                          <  2:1161569/76‑1 (MQ=255)                             gCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAg                                         <  2:40524/76‑1 (MQ=255)                             gCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAg                                         <  2:1422008/76‑1 (MQ=255)                                  aaCTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCAtt                                    >  2:262702/1‑76 (MQ=255)                                   aCTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTc                                   >  1:1421183/1‑76 (MQ=255)                                    cTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCa                                  <  1:574878/76‑1 (MQ=255)                                      gggCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATg                                >  2:1479995/1‑76 (MQ=255)                                      gggCACGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATg                                >  2:700534/1‑76 (MQ=255)                                          aCGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTa                            <  2:1088018/76‑1 (MQ=255)                                          aCGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTa                            <  1:882901/76‑1 (MQ=255)                                           cGGTTTGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTAt                           <  2:955007/76‑1 (MQ=255)                                                tGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGcg                      <  1:504447/76‑1 (MQ=255)                                                tGAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGcg                      >  1:1266966/1‑76 (MQ=255)                                                 gAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGcga                     <  1:56070/76‑1 (MQ=255)                                                 gAACGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGcga                     <  2:1253270/76‑1 (MQ=255)                                                   aCGTAACCAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGcgacg                   <  1:406430/76‑1 (MQ=255)                                                          cAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGAGCGGGt            >  2:404340/1‑76 (MQ=255)                                                          cAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGt            >  2:1266177/1‑76 (MQ=255)                                                          cAGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGt            >  2:1147432/1‑76 (MQ=255)                                                           aGCGCTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTc           >  2:1540848/1‑76 (MQ=255)                                                              gcTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCtt        <  1:1678691/76‑1 (MQ=255)                                                               cTCCACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCttt       <  2:1054610/76‑1 (MQ=255)                                                               cTCCACGAGACGATCCGGCACAGCAAATAACGTTGCGTTCAGCCATTCAGGTTTAGTGGCGAAGACCGGGTCGttt       >  2:1590011/1‑76 (MQ=255)                                                                 ccACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTcc     >  1:360933/1‑76 (MQ=255)                                                                 ccACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTcc     >  1:1638450/1‑76 (MQ=255)                                                                 ccACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTcc     >  1:1600543/1‑76 (MQ=255)                                                                 ccACGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTcc     <  1:1415123/76‑1 (MQ=255)                                                                    cGAGATGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAAt  >  1:581174/1‑76 (MQ=255)                                                                         |                                                                       GACAGGTAGTAAACCTGACGAGTTTCTTGCGACAACTGGGCACGGTTTGAACGTAACCAGCGCTCCACGAGACGATCGCGCACAGCAAATAACGTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAAT  >  REL606/3497809‑3497951

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