Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,691,8440TG41.9% 57.9 / 15.9 43N226T (AAC→ACC) yiaYpredicted Fe‑containing alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (6/19);  new base G (18/0);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 2.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCT  >  REL606/3691771‑3691918
                                                                         |                                                                          
ggCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGttt                                                                          <  2:1559905/76‑1 (MQ=255)
ggCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGttt                                                                          <  2:261833/76‑1 (MQ=255)
ggCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGttt                                                                          <  2:619653/76‑1 (MQ=255)
  cGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTc                                                                        <  1:164032/76‑1 (MQ=255)
  cGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTc                                                                        <  2:393284/76‑1 (MQ=255)
   gAGGAACTGGGCATAAGCCATTGCTTCATGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCg                                                                       >  2:782067/1‑76 (MQ=255)
   gAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGgttttg                                                                       >  1:1073304/1‑74 (MQ=255)
      gAACTGGGCATAAGCCATTGCTTCACGCGCTTTCCCATTACTGCCCTCTTCAACGGCTAACGGCAGGTTTTCGGCa                                                                    <  2:1291076/76‑1 (MQ=255)
         cTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATc                                                                 >  2:13870/1‑76 (MQ=255)
          tGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCa                                                                >  1:1390350/1‑76 (MQ=255)
           gggCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCAt                                                               <  2:71816/76‑1 (MQ=255)
           gggCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCAt                                                               <  2:715422/76‑1 (MQ=255)
           gggCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCAt                                                               >  2:1784908/1‑76 (MQ=255)
             gCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATgg                                                             >  2:702125/1‑76 (MQ=255)
              cATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGt                                                            >  1:4199/1‑76 (MQ=255)
                 aaGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGGAAACCTGGTCAc                                                         >  2:1379283/1‑76 (MQ=255)
                  aGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACg                                                        <  1:345014/76‑1 (MQ=255)
                      aTTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCtt                                                    <  2:35390/76‑1 (MQ=255)
                      aTTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCtt                                                    >  1:704204/1‑76 (MQ=255)
                       gtGCTTCGCGCGCTGTCGCATTACTGCCATCTCCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCttt                                                   <  2:1888228/75‑1 (MQ=255)
                             cACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTTAGTGc                                             >  1:308280/1‑76 (MQ=255)
                                  gcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAg                                        >  2:810439/1‑76 (MQ=255)
                                  gcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAg                                        >  2:1817845/1‑76 (MQ=255)
                                    tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCg                                      >  1:430250/1‑76 (MQ=255)
                                      tCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTc                                    <  1:512569/76‑1 (MQ=255)
                                       cGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCa                                   >  2:463984/1‑76 (MQ=255)
                                       cGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCCGTGCACAAGCGTcc                                   >  1:659080/1‑75 (MQ=255)
                                       cGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCa                                   >  1:1712090/1‑76 (MQ=255)
                                       cGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGGCACGGCTTTCAGTGCACAAGCggcc                                   >  1:113818/1‑73 (MQ=255)
                                         cATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTTGGCAATCATGGGCACGGCTTTTAGTGCACAAGCGTCAgg                                 >  1:1115040/1‑75 (MQ=255)
                                           ttACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCCTGGGCACGGCTTTCAGGGCACAAGCGTCAGGGa                               >  2:623747/1‑76 (MQ=255)
                                           ttACTGCCAACTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCCGTGa                               >  2:623725/1‑76 (MQ=255)
                                               tGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCgg                           <  1:594841/76‑1 (MQ=255)
                                                gCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGc                          >  1:682738/1‑76 (MQ=255)
                                                     cTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTggc                     >  1:1370662/1‑76 (MQ=255)
                                                     cTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTggc                     >  1:1722990/1‑76 (MQ=255)
                                                            ggCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATa              <  2:1632856/76‑1 (MQ=255)
                                                             gCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAg             <  2:316634/76‑1 (MQ=255)
                                                             gCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAg             <  2:1787159/76‑1 (MQ=255)
                                                              cTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGa            <  1:401535/76‑1 (MQ=255)
                                                              cTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGa            >  1:477523/1‑76 (MQ=255)
                                                                aaCGGCAGGATTTTGGCCATCATGGTCACGGCTTTCTGTGCCCAAGCGGCCGTGATCGGGGTGGCGGCAGTATaaa          >  2:831577/1‑76 (MQ=255)
                                                                  cGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACa        <  2:72638/76‑1 (MQ=255)
                                                                     cAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATAc     <  1:340250/76‑1 (MQ=255)
                                                                      aGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACg    <  1:49411/76‑1 (MQ=255)
                                                                      aGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACg    <  1:200999/76‑1 (MQ=255)
                                                                      aGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACg    <  1:1142541/76‑1 (MQ=255)
                                                                        gTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCt  >  1:675730/1‑76 (MQ=255)
                                                                         |                                                                          
GGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCT  >  REL606/3691771‑3691918

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