Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,393,0610TG46.2% 36.3 / 22.7 39V211G (GTG→GGG) aidBisovaleryl CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (3/18);  new base G (18/0);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 2.74e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC  >  REL606/4392991‑4393136
                                                                      |                                                                           
cAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggg                                                                        >  1:1270705/1‑75 (MQ=255)
 aGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCa                                                                       <  1:246059/76‑1 (MQ=255)
  gggCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGggggggggg                                                                      >  1:1403848/1‑73 (MQ=255)
    gCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATaa                                                                    <  2:1429371/76‑1 (MQ=255)
      ggTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAAt                                                                  <  2:676709/76‑1 (MQ=255)
       gTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATg                                                                 >  1:1843843/1‑76 (MQ=255)
           cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttt                                                             <  2:1781563/76‑1 (MQ=255)
           cGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAAAGGttt                                                             >  1:255640/1‑76 (MQ=255)
             aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGttttt                                                           <  2:1158691/76‑1 (MQ=255)
             aTGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGttttt                                                           >  2:616059/1‑76 (MQ=255)
               gTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTtct                                                         <  2:1412401/76‑1 (MQ=255)
                  aTGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTTTCgg                                                      >  2:491875/1‑76 (MQ=255)
                   tGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGt                                                     >  2:1108036/1‑76 (MQ=255)
                     aGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTc                                                   <  2:1609980/76‑1 (MQ=255)
                       cAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCg                                                 <  1:1452790/76‑1 (MQ=255)
                         acacCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGGTAAATGGGTTTTTTCGGGTTCGcc                                               >  1:733124/1‑75 (MQ=255)
                          caccacCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGcca                                              >  2:898878/1‑74 (MQ=255)
                           accaccCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCaaa                                             <  2:1403848/76‑1 (MQ=255)
                                   tGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGCAAAGCGATgcg                                     >  2:1214577/1‑76 (MQ=255)
                                      agagCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGGTTTTTTCGGTTCCGCACAGCGATGCGCAt                                  >  2:658971/1‑76 (MQ=255)
                                        agCGTCTGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTTTTTGTTCCGCCAAGCGATGCGCATCt                                >  2:903515/1‑76 (MQ=255)
                                             cTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGc                           <  1:1315174/76‑1 (MQ=255)
                                             cTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGc                           <  1:1221440/76‑1 (MQ=255)
                                              tGGAAGATGGCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCt                          >  1:1475379/1‑76 (MQ=255)
                                              tGGAAGATGGCTCTTATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCt                          >  2:1197232/1‑76 (MQ=255)
                                               ggAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTg                         <  1:1158223/76‑1 (MQ=255)
                                                 aaGATGGCTCTTATCGGCTGGTGGGGGATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGgc                       >  2:1221913/1‑76 (MQ=255)
                                                 aaGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTggn                       <  1:387133/76‑2 (MQ=255)
                                                 aaGATGGCTCTTATCGGCTGGTGGGGCATAAAAGGTTTTTTTTGGTTCCGCAAAACGATGCGCATCTGGTGCTGgc                       >  2:150287/1‑76 (MQ=255)
                                                 aaGATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGGGCTggg                       >  2:258859/1‑75 (MQ=255)
                                                   gATGGCTCTTATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGGTCCGCAAAGCGATGCGCATCTGGTGCTGgcgc                     >  2:426456/1‑76 (MQ=255)
                                                     tGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAg                   <  1:67450/76‑1 (MQ=255)
                                                       gCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGAc                 <  1:388108/76‑1 (MQ=255)
                                                       gCTCTTATCGGCTGGGGGGGGATAAATGGTTTTTTTCGGTTCCGCAAAGCCATGCGCATCTGGGTCTGGCGCAGAc                 >  2:1344283/1‑76 (MQ=255)
                                                           ttATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACcgcg             >  2:137471/1‑76 (MQ=255)
                                                           ttATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCGGAGcgcg             >  2:831707/1‑76 (MQ=255)
                                                            tATCGGCTGGGGGGGGATAAAAGGTTTTTCTCGGTTCCGCAAAACGATGCGCATCTGGGGCTGGCGCAGACCGCgg            >  2:1419651/1‑76 (MQ=255)
                                                            tATCGGCTGGGGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCgg            >  2:1460437/1‑76 (MQ=255)
                                                            tATCGGCTGGGGGGGCATAAATGGGTTTTCTCGGTTCCGCCAAGCGATGCGCATCTGGTGCTGGCGCCGACCGCgg            >  2:1839386/1‑76 (MQ=255)
                                                              tCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggt          <  2:1315055/76‑1 (MQ=255)
                                                              tCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggt          <  1:1174622/76‑1 (MQ=255)
                                                               cGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtg         <  2:551954/76‑1 (MQ=255)
                                                                ggCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGggtgg        >  1:1119292/1‑76 (MQ=255)
                                                                      tggGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTc  <  2:524692/76‑1 (MQ=255)
                                                                      |                                                                           
CAGGGCGGTTCCGATGTTATGAGCAACACCACCCGTGCAGAGCGTCTGGAAGATGGCTCTTATCGGCTGGTGGGGCATAAATGGTTTTTCTCGGTTCCGCAAAGCGATGCGCATCTGGTGCTGGCGCAGACCGCGGGTGGTCTGTC  >  REL606/4392991‑4393136

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