Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,492,0820TG42.5% 57.0 / 23.3 48V490G (GTG→GGG) glpDsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (4/23);  new base G (20/0);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.27e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC  >  REL606/3492009‑3492157
                                                                         |                                                                           
gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggg                                                                           >  1:12408/1‑76 (MQ=255)
  ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg                                                                         >  1:81971/1‑74 (MQ=255)
     gtgGAGACGAACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat                                                                      <  1:874105/76‑1 (MQ=255)
      tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata                                                                     <  1:1554184/76‑1 (MQ=255)
        gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg                                                                   <  2:1354261/76‑1 (MQ=255)
        gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg                                                                   <  1:906727/76‑1 (MQ=255)
           tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACgcag                                                                >  1:44012/1‑76 (MQ=255)
           tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACgcag                                                                >  1:1085778/1‑76 (MQ=255)
            cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATAagcagc                                                               >  1:1210809/1‑76 (MQ=255)
            cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACgcagc                                                               >  1:313198/1‑76 (MQ=255)
            cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACgcagc                                                               >  1:1277454/1‑76 (MQ=255)
               acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGa                                                            <  2:901826/76‑1 (MQ=255)
               acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGa                                                            >  1:444456/1‑76 (MQ=255)
               acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGa                                                            >  1:1744893/1‑76 (MQ=255)
                caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa                                                           <  2:1777272/76‑1 (MQ=255)
                 aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAg                                                          >  1:1186652/1‑76 (MQ=255)
                   aaCAAGGCATGTGGCTAAATGCAGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCt                                                        <  2:129942/76‑1 (MQ=255)
                     cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAAGCAGCAGAAGccgg                                                      >  1:218076/1‑73 (MQ=255)
                      aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTc                                                     >  2:1421520/1‑76 (MQ=255)
                       aGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCg                                                    <  1:542462/76‑1 (MQ=255)
                        ggCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGc                                                   <  2:1442980/76‑1 (MQ=255)
                         gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCt                                                  >  1:971411/1‑76 (MQ=255)
                          cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTg                                                 <  1:974372/76‑1 (MQ=255)
                          cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTg                                                 <  2:1690027/76‑1 (MQ=255)
                            tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGAATACGCAGCCGAAACTGTCGCGGGc                                               >  2:1745396/1‑76 (MQ=255)
                                gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcg                                           <  2:1724580/76‑1 (MQ=255)
                                gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcg                                           <  2:654664/76‑1 (MQ=255)
                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgt                                          <  1:444016/76‑1 (MQ=255)
                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgt                                          <  1:1769936/76‑1 (MQ=255)
                                  tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg                                         <  2:747664/76‑1 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGa                                        <  2:266564/76‑1 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGa                                        <  2:1382964/76‑1 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGa                                        >  1:1382386/1‑76 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGCGGCGTgg                                        >  1:1139969/1‑75 (MQ=255)
                                     aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAtt                                      >  1:110438/1‑76 (MQ=255)
                                       gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGGGTCGGGGtttt                                    >  2:445618/1‑75 (MQ=255)
                                            tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTa                               >  2:902458/1‑76 (MQ=255)
                                            tCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGCGATTTAATGTa                               >  2:1604014/1‑76 (MQ=255)
                                             caacaaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTTATGTaa                              >  2:1086801/1‑76 (MQ=255)
                                                caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGATTTAATGTAAggg                           >  2:1166979/1‑75 (MQ=255)
                                                   tctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTAAGGGGGc                        >  1:205401/1‑76 (MQ=255)
                                                   tctcGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGGTGTCGCTGGCGCCGTGTTTTAATTTAAGGTGGc                        >  2:1130563/1‑76 (MQ=255)
                                                   tctcGTGTGAGTCAGTGGCTGGCGGAGGATACGCAGCAGAAGCTGTCGCTGGGGTTGTGGTTTCATGTAAGGTGGc                        >  2:1763213/1‑76 (MQ=255)
                                                     tcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCa                      <  2:1256477/76‑1 (MQ=255)
                                                      cGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAg                     <  1:1772916/76‑1 (MQ=255)
                                                      cGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGGGATTTAATGTAAGGGGGCCAg                     >  1:740258/1‑76 (MQ=255)
                                                       gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAgg                    <  2:1548456/76‑1 (MQ=255)
                                                         gtgAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGATTTAATGTAAGGGGGCCAgggg                  >  1:633975/1‑74 (MQ=255)
                                                          tgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCa                 <  1:446882/76‑1 (MQ=255)
                                                               cAGTGGCTGGTGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGATTTAATGTAAGGGGGCCCGGTCAGAttt            >  2:1716791/1‑76 (MQ=255)
                                                                   ggCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAAt        >  2:718933/1‑76 (MQ=255)
                                                                     cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt      <  1:1196832/76‑1 (MQ=255)
                                                                      tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg     <  1:505174/76‑1 (MQ=255)
                                                                      tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg     >  2:1253893/1‑76 (MQ=255)
                                                                        gtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAc   >  1:1244666/1‑76 (MQ=255)
                                                                         tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc  >  2:670162/1‑76 (MQ=255)
                                                                         tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc  <  2:46685/76‑1 (MQ=255)
                                                                         tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc  >  2:1565024/1‑76 (MQ=255)
                                                                         |                                                                           
GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC  >  REL606/3492009‑3492157

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