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

ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC  >  REL606/3492007‑3492157
                                                                           |                                                                           
acgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggt                                                                             >  2:58889/1‑76 (MQ=255)
 cgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtg                                                                            <  1:1149182/76‑1 (MQ=255)
  gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtgg                                                                           <  1:1498204/76‑1 (MQ=255)
    ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                                         <  1:1398484/76‑1 (MQ=255)
       gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat                                                                      <  2:428199/76‑1 (MQ=255)
       gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtat                                                                      <  2:1413789/76‑1 (MQ=255)
         ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATAc                                                                    <  1:170195/76‑1 (MQ=255)
          gCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg                                                                   <  2:1568037/76‑1 (MQ=255)
            gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgca                                                                 <  2:474261/76‑1 (MQ=255)
              cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc                                                               <  2:1515928/76‑1 (MQ=255)
                cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag                                                             <  2:806759/76‑1 (MQ=255)
                cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACgcaggcg                                                             >  2:280222/1‑73 (MQ=255)
                cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGcagcag                                                             >  1:1249865/1‑76 (MQ=255)
                cacaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGcagcag                                                             >  1:661875/1‑76 (MQ=255)
                  caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGaa                                                           >  1:933314/1‑76 (MQ=255)
                   aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAg                                                          >  2:1191327/1‑76 (MQ=255)
                    aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGc                                                         <  1:481986/76‑1 (MQ=255)
                      aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTg                                                       <  2:1943714/76‑1 (MQ=255)
                        aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTc                                                     >  2:1250920/1‑76 (MQ=255)
                            cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTTGccg                                                 >  2:1385187/1‑74 (MQ=255)
                            cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGccg                                                 >  2:611026/1‑74 (MQ=255)
                            cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCCGAAGCTGTCGCTg                                                 >  2:1385178/1‑76 (MQ=255)
                            cATGTGGCTAAATGCGGATCAACAATCTCGTGGGAGTCAGCGGCTGGAGGCGTATTCGCAGCACAAACCGTCGcgg                                                 >  2:611006/1‑74 (MQ=255)
                              tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATAAGCAGCAGAAGCTGTCGCTggg                                               >  1:1931111/1‑75 (MQ=255)
                               gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCCGGcg                                              >  1:1762083/1‑76 (MQ=255)
                               gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGGAGCAGAAGGCGGCGCTgggg                                              >  1:328894/1‑74 (MQ=255)
                                tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTg                                                 <  2:1793613/72‑1 (MQ=255)
                                tgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTg                                                 >  1:1793613/1‑72 (MQ=255)
                                  gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcg                                           <  2:1137708/76‑1 (MQ=255)
                                  gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGgcggcg                                           >  1:232590/1‑76 (MQ=255)
                                   cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgt                                          <  2:1249865/76‑1 (MQ=255)
                                   cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGggtggt                                          >  2:1135067/1‑76 (MQ=255)
                                    tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGGGGCGTg                                         >  1:253415/1‑76 (MQ=255)
                                        tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGttt                                     >  2:687360/1‑76 (MQ=255)
                                          cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGGGGCGGGGTTTaa                                   >  1:1302394/1‑76 (MQ=255)
                                           ggATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAAt                                  <  1:516966/76‑1 (MQ=255)
                                            gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATg                                 <  1:627782/76‑1 (MQ=255)
                                            gATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATg                                 <  1:1436695/76‑1 (MQ=255)
                                             aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt                                <  1:1319069/76‑1 (MQ=255)
                                              tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGatttaattta                               >  2:1458514/1‑76 (MQ=255)
                                               caacaaTCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTaa                              >  2:523859/1‑76 (MQ=255)
                                                 acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAgg                            <  1:1609516/76‑1 (MQ=255)
                                                 acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAgg                            <  1:1578114/76‑1 (MQ=255)
                                                  caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGt                           <  1:228977/76‑1 (MQ=255)
                                                   aaTCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGAAGCTGTCGCTGGGGTCGGGATTTAATGTAAgggg                          >  2:1772142/1‑74 (MQ=255)
                                                   aaTCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTg                          >  2:1076153/1‑76 (MQ=255)
                                                     tctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGc                        <  1:1835164/76‑1 (MQ=255)
                                                       tcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCa                      <  2:676843/76‑1 (MQ=255)
                                                          tgtgAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGATTTAATGTAAGGGGGCCAggg                   >  2:1617187/1‑75 (MQ=255)
                                                                 cAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGAttt            <  2:522620/76‑1 (MQ=255)
                                                                    tGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCaa         >  2:1842256/1‑76 (MQ=255)
                                                                     ggCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAAt        >  1:1339003/1‑76 (MQ=255)
                                                                      gCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATc       >  1:982628/1‑76 (MQ=255)
                                                                       cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCAGATTTTAATCt      >  1:1726199/1‑76 (MQ=255)
                                                                        tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg     <  1:579896/76‑1 (MQ=255)
                                                                           tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc  >  2:349994/1‑76 (MQ=255)
                                                                           |                                                                           
ACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC  >  REL606/3492007‑3492157

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: ‑