Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,492,0820TG38.8% 79.6 / 16.8 49V490G (GTG→GGG) glpDsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (9/21);  new base G (19/0);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAC  >  REL606/3492009‑3492156
                                                                         |                                                                          
gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtgg                                                                          <  2:1357911/76‑1 (MQ=255)
gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggg                                                                          >  1:904334/1‑76 (MQ=255)
  ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg                                                                        >  1:896614/1‑74 (MQ=255)
   cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGAGGCTGGGGGAGt                                                                       >  1:1883286/1‑76 (MQ=255)
     gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat                                                                     <  1:648509/76‑1 (MQ=255)
      tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGgata                                                                    >  1:187765/1‑76 (MQ=255)
       ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATAc                                                                   <  1:1371347/76‑1 (MQ=255)
       ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATAc                                                                   >  2:234912/1‑76 (MQ=255)
        gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg                                                                  <  2:664593/76‑1 (MQ=255)
         cGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACgg                                                                 >  1:531143/1‑75 (MQ=255)
           tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcag                                                               <  1:303474/76‑1 (MQ=255)
           tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcag                                                               <  2:1407098/76‑1 (MQ=255)
            cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc                                                              <  1:1792382/76‑1 (MQ=255)
              cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag                                                            <  2:1173015/76‑1 (MQ=255)
               acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGa                                                           <  1:694101/76‑1 (MQ=255)
               acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGa                                                           <  2:1766777/76‑1 (MQ=255)
                caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa                                                          <  1:854978/76‑1 (MQ=255)
                  aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGc                                                        >  1:883965/1‑76 (MQ=255)
                   aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGAAGCt                                                       >  1:649190/1‑76 (MQ=255)
                     cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGt                                                     >  2:1833336/1‑76 (MQ=255)
                       aGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGAAGCTGTCg                                                   >  1:1418489/1‑76 (MQ=255)
                       aGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCg                                                   <  1:1534427/76‑1 (MQ=255)
                       aGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCg                                                   >  1:1361659/1‑76 (MQ=255)
                        ggCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCCGTCGc                                                  >  2:796937/1‑76 (MQ=255)
                          cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTg                                                <  2:941207/76‑1 (MQ=255)
                          cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTg                                                >  1:315055/1‑76 (MQ=255)
                            tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGc                                              <  1:103392/76‑1 (MQ=255)
                              tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgt                                            <  1:1089673/76‑1 (MQ=255)
                              tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCCGGcgg                                            >  1:1772726/1‑75 (MQ=255)
                              tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCgggggg                                            >  1:483230/1‑73 (MQ=255)
                                gccaaaTGCGGATCAACAATCTCGTGTGCGTCAGTGGCTGGTGGGGGATACGGCGCAGAAGGGGTTGCTTGggtcg                                          >  1:905189/4‑76 (MQ=255)
                                gCTAAATGCGGATCAACAATCTCGTGTGCGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcg                                          >  1:905209/1‑76 (MQ=255)
                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgt                                         <  1:275362/76‑1 (MQ=255)
                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCCGGcgtcgt                                         >  1:1852151/1‑76 (MQ=255)
                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgt                                         >  1:977440/1‑76 (MQ=255)
                                  tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTg                                        >  2:283597/1‑76 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAAAAGCCGTCGCTGGCGACGTgt                                       >  1:251688/1‑75 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGa                                       >  1:455358/1‑76 (MQ=255)
                                    aaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGAt                                      >  1:1667472/1‑76 (MQ=255)
                                     aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAtt                                     >  1:870204/1‑76 (MQ=255)
                                      tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt                                    <  2:1768900/76‑1 (MQ=255)
                                       gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa                                   >  1:1151235/1‑76 (MQ=255)
                                        cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTaa                                  <  2:885404/76‑1 (MQ=255)
                                          gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCCGCAGAAGCTGTCGCTGGCGGCGGGATTTAATg                                >  2:163973/1‑76 (MQ=255)
                                              aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAg                            >  1:37437/1‑76 (MQ=255)
                                              aacaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAg                            <  2:1629141/76‑1 (MQ=255)
                                              aacaaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAg                            >  1:1278471/1‑76 (MQ=255)
                                               acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGGGATTTAATGTAAgg                           >  2:124596/1‑76 (MQ=255)
                                                caaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGt                          >  1:935502/1‑76 (MQ=255)
                                                 aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTg                         <  1:1237707/76‑1 (MQ=255)
                                                    ctcGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTAAGgggggc                      >  2:410679/1‑74 (MQ=255)
                                                         ttgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGATTTAATGTAAGGTGGGCGGGtg                 >  2:823856/2‑75 (MQ=255)
                                                         gtgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTc                 <  1:76880/76‑1 (MQ=255)
                                                          tgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCa                <  1:980639/76‑1 (MQ=255)
                                                              tCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGGGATTTTATGTAAGGGGGCCCGGGCAGCtt            >  2:591423/1‑76 (MQ=255)
                                                                    gCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCAGATTTTAATc      >  1:951235/1‑76 (MQ=255)
                                                                     cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt     >  2:110147/1‑76 (MQ=255)
                                                                       ggtggAGTATAGGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGa   >  2:1656881/1‑76 (MQ=255)
                                                                        gtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAc  <  1:1691497/76‑1 (MQ=255)
                                                                        gtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCCGGTTTCCATCTTAc  >  1:329419/1‑76 (MQ=255)
                                                                         |                                                                          
GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAC  >  REL606/3492009‑3492156

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