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

ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC  >  REL606/3492007‑3492157
                                                                           |                                                                           
acgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggt                                                                             <  1:772806/76‑1 (MQ=255)
acgCACTGTGTCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggt                                                                             <  1:674444/76‑1 (MQ=255)
  gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTTAGTGGCtggtgg                                                                           <  2:409350/76‑1 (MQ=255)
     cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGCGTGAGTCAGTGGCtgggggggg                                                                        >  2:39701/1‑73 (MQ=255)
         ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAc                                                                    >  2:237163/1‑76 (MQ=255)
             tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAGTAATCTCGTGTGAGTCAGTGNCTGGTGGAGTATACgcag                                                                <  1:1154838/76‑1 (MQ=255)
              cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACgcagc                                                               >  1:480760/1‑76 (MQ=255)
                   aaaaCAAGGCATGTGGCTNAATGCGGATCAACAATCTCGCGTGAGTCAGTGGCTGGTGGGGGATACGCAGCAGAAg                                                          >  1:1154305/1‑76 (MQ=255)
                   aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCCGCAGAAg                                                          >  2:67847/1‑76 (MQ=255)
                     aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCt                                                        <  1:415063/76‑1 (MQ=255)
                      aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTg                                                       >  1:1077793/1‑76 (MQ=255)
                       cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGt                                                      <  1:748361/76‑1 (MQ=255)
                        aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGGTGTc                                                     >  1:559796/1‑76 (MQ=255)
                        aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGAGGCTGGGGGGGGATACGCAGCAGAAGCTGTc                                                     >  2:200051/1‑76 (MQ=255)
                         aGGCATGTGGCTAAATGTGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCg                                                    >  1:433457/1‑76 (MQ=255)
                             aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTgg                                                >  2:1075414/1‑76 (MQ=255)
                                tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgt                                             <  1:502795/76‑1 (MQ=255)
                                 ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtc                                            <  2:762489/76‑1 (MQ=255)
                                 ggCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGGCGATGGcggg                                            >  2:812158/1‑74 (MQ=255)
                                   cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGcggcgt                                          >  1:725719/1‑76 (MQ=255)
                                    tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGGGTCGTg                                         <  2:933889/76‑1 (MQ=255)
                                    tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGgcggcggg                                         >  1:991439/1‑74 (MQ=255)
                                    tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCCGGCGTCGTg                                         >  2:1180315/1‑76 (MQ=255)
                                     aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGgcggcgggg                                        >  1:616388/1‑73 (MQ=255)
                                      aaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTggt                                       >  2:985883/1‑74 (MQ=255)
                                        tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt                                     <  2:673830/76‑1 (MQ=255)
                                         gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa                                    <  1:834931/76‑1 (MQ=255)
                                          cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCCGGCGGCGTGGtttta                                   >  1:909147/1‑74 (MQ=255)
                                           ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCATAAGCTGTCGCTGGCGTCGTGATTTAAt                                  <  1:1187821/76‑1 (MQ=255)
                                            gATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCCGTCGCTGGCGGCGTGATTTAATg                                 >  1:729352/1‑76 (MQ=255)
                                               caacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTaa                              <  1:638207/76‑1 (MQ=255)
                                                 acaaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGCGGCTGGGCGCGTGGTTTTATGTAAgg                            >  1:57376/1‑76 (MQ=255)
                                                   aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTg                          <  1:39701/76‑1 (MQ=255)
                                                    aTCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGGGGCGGGGTTTTATGTAAGgggg                         >  1:901097/1‑76 (MQ=255)
                                                      ctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGcc                       >  2:412080/1‑76 (MQ=255)
                                                       tcGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTAAGGTGGCCa                      >  1:946862/1‑76 (MQ=255)
                                                           gtgAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGGAAGGGGGCCAgggc                  >  2:716820/1‑74 (MQ=255)
                                                                tCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAtt             <  1:917535/76‑1 (MQ=255)
                                                                 cAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGGGGCGGGATTTTATGTAAGGGGGCCCGGGCCGAttt            >  2:436393/1‑76 (MQ=255)
                                                                    tGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGNGATTTAATGTAAGGTGGCCAGGTCAGATTTCaa         <  2:1016300/76‑1 (MQ=255)
                                                                      gCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATc       <  1:916106/76‑1 (MQ=255)
                                                                       cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt      >  2:11901/1‑76 (MQ=255)
                                                                       cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt      <  1:622938/76‑1 (MQ=255)
                                                                        tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg     >  1:1001731/1‑76 (MQ=255)
                                                                        tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg     <  1:760945/76‑1 (MQ=255)
                                                                           tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc  <  1:1029966/76‑1 (MQ=255)
                                                                           |                                                                           
ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC  >  REL606/3492007‑3492157

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