Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,174,6480TG34.0% 75.5 / 26.7 50V409G (GTG→GGG) dldD‑lactate dehydrogenase, FAD‑binding, NADH independent
Reads supporting (aligned to +/- strand):  ref base T (15/18);  new base G (17/0);  total (32/18)
Fisher's exact test for biased strand distribution p-value = 9.03e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCGATAAATACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGTTTTTACACCG  >  REL606/2174574‑2174723
                                                                          |                                                                           
gcgcGATAAATACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCtggtg                                                                            <  2:321507/76‑1 (MQ=255)
   cgATAAATACGAGCATCATCTGCTTTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGat                                                                         <  1:388455/76‑1 (MQ=255)
   cgATAAATACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGat                                                                         <  2:1579603/76‑1 (MQ=255)
       aaaTACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGgttttt                                                                     >  2:1145919/1‑76 (MQ=255)
       aaaTACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGattatt                                                                     >  2:1956182/1‑76 (MQ=255)
         aTACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTc                                                                   <  2:2053591/76‑1 (MQ=255)
           aCGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTTaa                                                                 >  1:1368410/1‑76 (MQ=255)
           aCGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTCaa                                                                 >  1:1307669/1‑76 (MQ=255)
            cGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTTTTTTaaa                                                                >  1:1832454/1‑76 (MQ=255)
              aGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGGTTATTTCAAACa                                                              >  2:1369134/1‑76 (MQ=255)
              aGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTCAAACa                                                              >  1:1766431/1‑76 (MQ=255)
                catcatCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAgg                                                            <  1:221403/76‑1 (MQ=255)
                 atcatcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTTAAAAAGGc                                                           >  1:1748292/1‑76 (MQ=255)
                 atcatcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTTAAAAAGGc                                                           >  2:179668/1‑76 (MQ=255)
                  tcatcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTTAAAAAGGcc                                                          >  2:1533046/1‑76 (MQ=255)
                  tcatcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTCAAACAggac                                                          >  2:1008877/1‑74 (MQ=255)
                    atcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTCAAAAAGGCCGa                                                        >  2:702728/1‑76 (MQ=255)
                    atcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAACAGGCCGa                                                        >  2:1197966/1‑76 (MQ=255)
                     tcTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTTAAAAAAGCCGaa                                                       >  2:1370958/1‑76 (MQ=255)
                      ctgctgTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAg                                                      <  1:1992851/76‑1 (MQ=255)
                      ctgctgTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAg                                                      <  2:833670/76‑1 (MQ=255)
                      ctgctgTTAAAAATGGCGGGCGATGGCGGGGGCGAAGCCAAATCGTGGCTGGGGGGTTTTTTCAAACAGGCCGAAg                                                      >  1:1939790/1‑76 (MQ=255)
                        gctgTTAAAAATGGCGGGCGATGGCGTGGGTGAAGCCAAATCGTGGCTGGGGGGTTATTTTAAACAGGCCGAAggg                                                    >  2:1433763/1‑75 (MQ=255)
                          tgTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTTAAACAGGGCGGAggggg                                                  >  1:1218510/1‑73 (MQ=255)
                           gTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGGTTATTTCAAACAGGCCGAAGGCggt                                                 >  1:1796037/1‑74 (MQ=255)
                               aaaaTGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGAGTATTTTAAACAGGCCGAAAGCGATTTCt                                             >  2:1843995/1‑76 (MQ=255)
                                  aTGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTg                                          >  1:2048649/1‑76 (MQ=255)
                                   tGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGt                                         <  1:224369/76‑1 (MQ=255)
                                   tGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGt                                         >  1:1179765/1‑76 (MQ=255)
                                   tGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGt                                         >  2:15725/1‑76 (MQ=255)
                                      cgggcgATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTg                                      <  1:1246359/76‑1 (MQ=255)
                                       gggcgATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGt                                     <  2:1406763/76‑1 (MQ=255)
                                        ggcgATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTa                                    <  1:991763/76‑1 (MQ=255)
                                         gcgATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTAc                                   <  1:1883219/76‑1 (MQ=255)
                                          cgATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACg                                  >  2:789044/1‑76 (MQ=255)
                                           gATGGCGTGGGCGCAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGc                                 >  2:1372242/1‑76 (MQ=255)
                                           gATGGCGTGGGCGAAGCCAAATCGTGGCTGGGGGATTATTTCAAACAGGCCGAAGGCGATTTTTTTGTCTGTACGc                                 >  1:1579258/1‑76 (MQ=255)
                                               gCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCgga                             <  2:194275/76‑1 (MQ=255)
                                                 gTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCggagg                           >  2:418392/1‑76 (MQ=255)
                                                  tGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGGTTTCTTTGTCTGTACGCCggagga                          >  2:900570/1‑76 (MQ=255)
                                                  tGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCggagga                          >  2:1403480/1‑76 (MQ=255)
                                                   gggCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGaa                         >  1:650071/1‑76 (MQ=255)
                                                   gggCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGaa                         <  2:543430/76‑1 (MQ=255)
                                                   gggCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGaa                         >  1:648867/1‑76 (MQ=255)
                                                    ggCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAg                        <  2:416259/76‑1 (MQ=255)
                                                      ngAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGgc                      >  1:1525181/2‑76 (MQ=255)
                                                      cGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGGAGgc                      >  1:738265/1‑76 (MQ=255)
                                                       gAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGgca                     <  2:331561/76‑1 (MQ=255)
                                                            cAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCaaa                >  2:1759695/1‑76 (MQ=255)
                                                            cAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCaaa                <  1:1145919/76‑1 (MQ=255)
                                                                 cGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCACAGCGtt           >  2:811178/1‑76 (MQ=255)
                                                                  gTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGttt          <  2:1332088/76‑1 (MQ=255)
                                                                   tGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGtttt         <  2:622566/76‑1 (MQ=255)
                                                                    ggCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGttttt        >  2:23458/1‑76 (MQ=255)
                                                                      cTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGTTTTTac      >  1:1585877/1‑76 (MQ=255)
                                                                          tggATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGTTTTTACACCg  >  1:766106/1‑76 (MQ=255)
                                                                          |                                                                           
GCGCGATAAATACGAGCATCATCTGCTGTTAAAAATGGCGGGCGATGGCGTGGGCGAAGCCAAATCGTGGCTGGTGGATTATTTCAAACAGGCCGAAGGCGATTTCTTTGTCTGTACGCCGGAGGAAGGCAGCAAAGCGTTTTTACACCG  >  REL606/2174574‑2174723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑