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

GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCT  >  REL606/3492009‑3492153
                                                                         |                                                                       
gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtgg                                                                       <  2:253005/76‑1 (MQ=255)
  ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggggg                                                                     >  2:804195/1‑74 (MQ=255)
     gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGtaa                                                                  >  1:374196/1‑75 (MQ=255)
     gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat                                                                  <  1:10678/76‑1 (MQ=255)
       ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATAc                                                                >  2:759610/1‑76 (MQ=255)
         cGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgc                                                              <  2:1417041/76‑1 (MQ=255)
         cGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACgc                                                              >  1:228894/1‑76 (MQ=255)
            cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACgcagc                                                           >  2:1521579/1‑76 (MQ=255)
              cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag                                                         <  1:31644/76‑1 (MQ=255)
              cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag                                                         <  1:349977/76‑1 (MQ=255)
              cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag                                                         <  1:742622/76‑1 (MQ=255)
              cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag                                                         <  1:529153/76‑1 (MQ=255)
              cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag                                                         <  1:673348/76‑1 (MQ=255)
                caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa                                                       <  2:350189/76‑1 (MQ=255)
                caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa                                                       <  2:350183/76‑1 (MQ=255)
                  aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATAAGCAGCAGAAGc                                                     >  2:1173097/1‑76 (MQ=255)
                    aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTg                                                   >  1:148719/1‑76 (MQ=255)
                      aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTc                                                 >  2:1653056/1‑76 (MQ=255)
                       aGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCg                                                <  2:1469311/76‑1 (MQ=255)
                         gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCt                                              >  2:508716/1‑76 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTTAGTGGCTGGGGGAGTATACGCGGCAGAAACGTTCGCCGGcg                                          >  1:1324604/1‑76 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTgggg                                          >  1:1404346/1‑74 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGcg                                          >  1:739874/1‑76 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGctggcgctggcg                                          >  1:637581/1‑76 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGcg                                          >  1:766171/1‑76 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTgggg                                          >  1:976330/1‑74 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTgggg                                          >  1:432541/1‑74 (MQ=255)
                                  tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg                                     <  1:1312002/76‑1 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGa                                    >  2:26276/1‑76 (MQ=255)
                                    aaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAt                                   <  1:1534042/76‑1 (MQ=255)
                                    aaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGGGAt                                   >  2:1255507/1‑76 (MQ=255)
                                       gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa                                <  2:790571/76‑1 (MQ=255)
                                       gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa                                <  1:133413/76‑1 (MQ=255)
                                        cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTaa                               >  1:60895/1‑76 (MQ=255)
                                            tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa                           <  1:25208/76‑1 (MQ=255)
                                                    ctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTTATGTAAGGTGGcc                   >  2:1236970/1‑76 (MQ=255)
                                                     tcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCa                  <  2:945736/76‑1 (MQ=255)
                                                     tcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCa                  <  2:345920/76‑1 (MQ=255)
                                                       gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAgg                <  1:612021/76‑1 (MQ=255)
                                                       gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAgg                <  1:1011821/76‑1 (MQ=255)
                                                       gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAgg                <  2:1097027/76‑1 (MQ=255)
                                                           gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAg            <  2:1041452/76‑1 (MQ=255)
                                                                     cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt  >  2:1177221/1‑76 (MQ=255)
                                                                         |                                                                       
GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCT  >  REL606/3492009‑3492153

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