Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,492,0820TG34.3% 63.0 / 14.6 35V490G (GTG→GGG) glpDsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (3/20);  new base G (12/0);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 5.45e-07
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                                                                           >  2:837956/1‑76 (MQ=255)
 cccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtgggggg                                                                          >  2:804461/1‑75 (MQ=255)
  ccTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAg                                                                         <  1:771130/76‑1 (MQ=255)
   cTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGt                                                                        <  1:495056/76‑1 (MQ=255)
      tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata                                                                     <  1:92315/76‑1 (MQ=255)
      tgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGtata                                                                     <  2:405032/76‑1 (MQ=255)
       ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGtataa                                                                    >  2:685491/1‑75 (MQ=255)
        gCGTCGCACAAAACAAGGGATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACg                                                                   >  2:724195/1‑76 (MQ=255)
        gCGTCGCACAAAACAAGGGATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTAGCTGGGGGAGTATACg                                                                   >  2:724174/1‑76 (MQ=255)
           tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcag                                                                <  1:1222364/76‑1 (MQ=255)
              cacaAAACAAGGCATGTGGCTGAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGcagcag                                                             <  1:1055476/76‑1 (MQ=255)
                caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa                                                           <  1:863845/76‑1 (MQ=255)
                caAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGaa                                                           <  1:954400/76‑1 (MQ=255)
                 aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAg                                                          <  1:1118631/76‑1 (MQ=255)
                 aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAg                                                          >  1:830092/1‑76 (MQ=255)
                   aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCt                                                        <  1:651107/76‑1 (MQ=255)
                      aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTc                                                     >  2:44547/1‑76 (MQ=255)
                         gcgtgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCt                                                  <  1:733012/73‑1 (MQ=255)
                         gCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCt                                                  <  1:1043986/76‑1 (MQ=255)
                          cATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTg                                                 >  1:528491/1‑76 (MQ=255)
                           aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTgg                                                <  2:534383/76‑1 (MQ=255)
                            tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGc                                               <  2:468273/76‑1 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcg                                              <  1:1040994/76‑1 (MQ=255)
                                gCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTAAAAGCAACAAAAACTGTCGCTGGcgtcg                                           >  1:84260/1‑76 (MQ=255)
                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGcgtcgg                                          >  2:974115/1‑75 (MQ=255)
                                  tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg                                         >  2:1050588/1‑76 (MQ=255)
                                  tAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTg                                         <  2:1101682/76‑1 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGa                                        >  1:259529/1‑76 (MQ=255)
                                      tGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAttt                                     <  1:438244/76‑1 (MQ=255)
                                       gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTa                                    <  1:342379/76‑1 (MQ=255)
                                        cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTaa                                   <  1:569020/76‑1 (MQ=255)
                                        cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTaa                                   <  1:1121247/76‑1 (MQ=255)
                                           aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt                                <  1:53565/76‑1 (MQ=255)
                                           aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt                                <  2:22209/76‑1 (MQ=255)
                                           aTCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGt                                <  1:1007030/76‑1 (MQ=255)
                                              aacaaTCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTTGCTGGGGGCGTGGTTTTATGTAAg                             >  2:464811/1‑76 (MQ=255)
                                              aacaaTCTCGTGTGAGTCAGTGGCTGGGGGAATAAACGCAACAAAAACTGTCGCTGGCGGCGTGAATTAAAGTAAg                             >  2:194662/1‑76 (MQ=255)
                                                 aaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTg                          <  2:338805/76‑1 (MQ=255)
                                                   tctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGc                        <  1:478785/76‑1 (MQ=255)
                                                    ctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGGTTTTATGTAAGgggggc                       >  2:704001/1‑74 (MQ=255)
                                                    ctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTTATGTAAGGGGGcc                       >  1:147931/1‑76 (MQ=255)
                                                       gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGGTTTAATGTAAGGTGGCCAgg                    >  2:764240/1‑76 (MQ=255)
                                                       gtgtGAGTCAGTGGCTGGGGGGGGATACGCCGCAGAAGCTGTCGCTGGCGTCGGGGTTTAATGTCAGGGGGCCAgg                    >  2:764249/1‑76 (MQ=255)
                                                        tgtgAGTCAGTGGCTGGTGGAGTATACGAAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGt                   <  2:144379/76‑1 (MQ=255)
                                                                     cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt      >  2:11590/1‑76 (MQ=255)
                                                                         tggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAcc  <  1:31404/76‑1 (MQ=255)
                                                                         |                                                                           
GCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGACC  >  REL606/3492009‑3492157

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