Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,492,0820TG35.9% 88.1 / 17.3 53V490G (GTG→GGG) glpDsn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (6/28);  new base G (18/1);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 3.31e-08
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:536146/76‑1 (MQ=255)
gCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggggg                                                                          >  1:1285178/1‑76 (MQ=255)
    tgtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTTTGAGTCAGTGGCTGGTGGAGta                                                                      <  1:1047853/76‑1 (MQ=255)
    tgtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGta                                                                      >  1:1176666/1‑76 (MQ=255)
       ggCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATAc                                                                   <  1:375973/76‑1 (MQ=255)
        gggTCGGACAAAACGAGGCATGGGGCGAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg                                                                  <  1:524204/74‑1 (MQ=255)
        gcgcCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACg                                                                  <  1:524179/76‑1 (MQ=255)
         cGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgc                                                                 <  2:463139/76‑1 (MQ=255)
          gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACgca                                                                >  2:912638/1‑76 (MQ=255)
          gTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACgca                                                                >  2:236641/1‑76 (MQ=255)
            cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagc                                                              <  1:1026508/76‑1 (MQ=255)
            cGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACgcagc                                                              >  1:1144604/1‑76 (MQ=255)
             gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACgcaggg                                                             >  2:350327/1‑74 (MQ=255)
               acaAAACAAGGCATGTGGCTAGATTCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGa                                                           <  1:1238258/76‑1 (MQ=255)
                  aaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAgg                                                        >  2:994447/1‑75 (MQ=255)
                   aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGAAGCt                                                       >  2:815873/1‑76 (MQ=255)
                   aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCt                                                       <  1:594206/76‑1 (MQ=255)
                    aCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTg                                                      <  2:412190/76‑1 (MQ=255)
                     cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGt                                                     >  1:366072/1‑76 (MQ=255)
                     cAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGt                                                     >  2:862887/1‑76 (MQ=255)
                      aaGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTc                                                    <  1:564023/76‑1 (MQ=255)
                           aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTgg                                               <  1:298738/76‑1 (MQ=255)
                           aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTgg                                               <  2:473772/76‑1 (MQ=255)
                           aTGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTgg                                               >  2:1299861/1‑76 (MQ=255)
                            tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGc                                              <  1:922025/76‑1 (MQ=255)
                            tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGGCGGTGGc                                              >  2:1043759/1‑76 (MQ=255)
                             gtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGcg                                             <  2:25718/76‑1 (MQ=255)
                                gCTAAAGGCGGATCAACAATGTNNTGTGAGTCAGNNCCTGGTGGAGTATATGCAGCAGAAGCTGTCGCTGGcgtcg                                          <  1:1185453/76‑1 (MQ=255)
                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATAAGCAGGAGAAGCTTTCGCTGgcggcgg                                         >  1:256233/1‑75 (MQ=255)
                                 cTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCGGgcggcgg                                         >  1:750015/1‑75 (MQ=255)
                                   aaaTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTgg                                       >  2:429637/1‑75 (MQ=255)
                                     aTGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAtt                                     <  1:512490/76‑1 (MQ=255)
                                       gCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGGGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGGtttt                                   >  2:999118/1‑75 (MQ=255)
                                        cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGATTTaa                                  >  2:26089/1‑76 (MQ=255)
                                        cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGTATACGCAGCAGAAGCTGTCGCTGGCGGCGGGGtttta                                  >  2:218980/1‑74 (MQ=255)
                                        cGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGAtttta                                  >  2:1175874/1‑74 (MQ=255)
                                            tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa                              <  1:312796/76‑1 (MQ=255)
                                            tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa                              <  1:312817/76‑1 (MQ=255)
                                            tCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTa                              <  1:393368/76‑1 (MQ=255)
                                            tCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGGGATTTTATGTa                              >  1:817636/1‑76 (MQ=255)
                                              aacaaTCTCGTGTGAGTCAGTGGCTGGTGGGGGATACGCCGCAGAAGCTGTCGCTGGCGGCGTGGTTTTATGTAAg                            >  2:693864/1‑76 (MQ=255)
                                              aacaaTCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTAATGTAAg                            >  2:725962/1‑76 (MQ=255)
                                               acaaTCTCGTGTGAGTCAGTGGCTGGTGNAGTATACGCGTCCGAAGCAGTCGCTGGCGTCGTGATTTAATGTAAgg                           <  1:113601/76‑1 (MQ=255)
                                               acaaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAgg                           <  1:791670/76‑1 (MQ=255)
                                               aaaaTCTCGTGTGGGTCGGGGGCTGGGGGAGTATACGCAGCGGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAgg                           <  1:791669/74‑1 (MQ=255)
                                                caaTCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGt                          <  2:765120/76‑1 (MQ=255)
                                                 aaTCTCGTGTGACTCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGGCGCTGGCGGCGGGGTTTTATTTAAGGTg                         >  2:459978/1‑76 (MQ=255)
                                                    ctcGTGTGAGTCAGTGGGTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGcc                      <  1:141027/76‑1 (MQ=255)
                                                    ctcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGGGGcc                      >  1:1017586/1‑76 (MQ=255)
                                                     tcGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCa                     <  2:778882/76‑1 (MQ=255)
                                                       gtgtGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAAGCTGTCGCTGGCGGCGTGGTTTTATGTAAAGGGGGCAgg                   >  2:1054940/1‑76 (MQ=255)
                                                          tgAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCa                <  1:557371/76‑1 (MQ=255)
                                                           gAGTCAGTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAg               <  2:1152523/76‑1 (MQ=255)
                                                              tCAGTGGCTGGGGGAGGATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGGGGGCAGGGCCGAtt            >  2:472308/1‑76 (MQ=255)
                                                                aGTGGCTGGGGGGGGATACGCAGCAGAAACTGTCGCTGGGGGGGTAGTTTTATTTAAGGGGGNCAGGGGAGATTTc          >  2:1197676/1‑76 (MQ=255)
                                                                 gTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCa         <  1:182117/76‑1 (MQ=255)
                                                                 gTGGCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCa         <  1:906357/76‑1 (MQ=255)
                                                                  cgacTGGTGGAGTNTACGCAGCAGGAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCaa        <  1:174142/73‑1 (MQ=255)
                                                                   ggCTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAAt       <  2:739732/76‑1 (MQ=255)
                                                                     cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCt     >  2:988044/1‑76 (MQ=255)
                                                                     cTGGTGGAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATTTAATGTAAGGTGGCCAGGGCAGATTTTAATCt     >  1:289912/1‑76 (MQ=255)
                                                                      tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATNTAATGTAAGGTGGCCAGGTCAGATTTCAATCTg    <  2:1113548/76‑1 (MQ=255)
                                                                      tggtggAGTATACGCAGCAGAAGCTGTCGCTGGCGTCGTGATCTAATGTAAGGTGGCCAGGTCAGAGTTCAATCTg    <  2:1129079/76‑1 (MQ=255)
                                                                        gtggAGTATACGCAGCAGAAGCTGTCGNTGGCGTCGTGATTTAATGTAAGGTGGCCAGGTCAGATTTCAATCTGAc  >  1:1309895/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 < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: