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

GGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCTT  >  REL606/3691771‑3691919
                                                                         |                                                                           
ggCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGttt                                                                           <  2:210474/76‑1 (MQ=255)
ggCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGttt                                                                           <  2:387833/76‑1 (MQ=255)
  cGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTc                                                                         <  2:110214/76‑1 (MQ=255)
    aGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCgg                                                                       >  2:245544/1‑76 (MQ=255)
      gAACTGGGCATATGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCa                                                                     <  2:73105/76‑1 (MQ=255)
      gAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTTGGCa                                                                     >  2:562188/1‑76 (MQ=255)
          tGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCa                                                                 >  2:515361/1‑76 (MQ=255)
          tGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCa                                                                 <  2:93386/76‑1 (MQ=255)
          tGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCa                                                                 >  2:307697/1‑76 (MQ=255)
          tGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGGTTCGGCAAACa                                                                 >  1:483276/1‑76 (MQ=255)
           gggCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCAt                                                                >  2:36832/1‑76 (MQ=255)
             gCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATgg                                                              >  2:1290279/1‑76 (MQ=255)
              cATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGt                                                             <  1:278244/76‑1 (MQ=255)
               aTAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTc                                                            <  1:638456/76‑1 (MQ=255)
               aTAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTc                                                            >  1:863703/1‑76 (MQ=255)
                tAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCa                                                           <  2:745482/76‑1 (MQ=255)
                 aaGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCAc                                                          <  1:1393254/76‑1 (MQ=255)
                 aaGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAAACATGGTccc                                                          >  2:850317/1‑74 (MQ=255)
                   gCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACgg                                                        <  2:653275/76‑1 (MQ=255)
                    ccATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGAAGGGTTTTGGTCCTTATGGTGACggg                                                       >  2:160607/1‑75 (MQ=255)
                       ttGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCttt                                                    <  1:806636/76‑1 (MQ=255)
                       ttGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCttt                                                    >  1:559261/1‑76 (MQ=255)
                        tGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTc                                                   <  1:703051/76‑1 (MQ=255)
                        tGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTc                                                   <  2:140887/76‑1 (MQ=255)
                             cACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTTGGGGATCATGGTTACGGGTTTTAGTgg                                              >  1:263515/1‑75 (MQ=255)
                               cgcgCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGcac                                            <  2:1299756/76‑1 (MQ=255)
                               cgcgCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGcac                                            >  1:133204/1‑76 (MQ=255)
                                gcgcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGcaca                                           <  1:293305/76‑1 (MQ=255)
                                 cgcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGAGCACaa                                          >  2:1407319/1‑76 (MQ=255)
                                  gcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAg                                         <  1:787058/76‑1 (MQ=255)
                                    tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCg                                       >  2:885461/1‑76 (MQ=255)
                                    tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCCGCAATTATTGTCACGGCTTTCAGTGCACAAGCg                                       >  2:940427/1‑76 (MQ=255)
                                        gCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAg                                   <  1:245544/76‑1 (MQ=255)
                                           ttACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGa                                >  2:1026006/1‑76 (MQ=255)
                                            tACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGAt                               >  1:910124/1‑76 (MQ=255)
                                             aCTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATc                              >  1:66553/1‑76 (MQ=255)
                                                gCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGc                           <  2:1249019/76‑1 (MQ=255)
                                                 ccATCTTCAACGGCTAACGGCAGGTTTTCGGGAATCATGGTCACGGCTTTCAGAGCACAAGCGTCAGTGATCGGCg                          >  1:315605/1‑76 (MQ=255)
                                                  cATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGt                         <  2:1409246/76‑1 (MQ=255)
                                                   aTCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTg                        <  2:1391271/76‑1 (MQ=255)
                                                   aTCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTg                        <  2:47866/76‑1 (MQ=255)
                                                    tCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCAAGGGGTTCAGTGGACAAGCGGCAGTGATCGGCGCgg                       >  2:1004695/1‑76 (MQ=255)
                                                     cTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTggc                      <  1:311622/76‑1 (MQ=255)
                                                       tCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTggcgg                    >  1:1509364/1‑76 (MQ=255)
                                                         aaCGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTTAGTGCACAAGCGGCCGGGGTCGGCGggggggga                  >  2:207481/1‑74 (MQ=255)
                                                          aCGGCTAACGGCAGGTTTTCGGCAATCATGGGCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGGGGCGGCaa                 >  2:1311643/1‑76 (MQ=255)
                                                            ggCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATa               <  2:1203620/76‑1 (MQ=255)
                                                              cTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGa             <  1:89748/76‑1 (MQ=255)
                                                              cTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCACTAGa             >  1:389707/1‑76 (MQ=255)
                                                                  cGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACa         <  1:525927/76‑1 (MQ=255)
                                                                  cGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACa         <  1:265830/76‑1 (MQ=255)
                                                                   ggCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACAt        <  2:1393561/76‑1 (MQ=255)
                                                                    gCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAACAGAAACATa       <  1:176769/76‑1 (MQ=255)
                                                                      aGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACg     >  1:852462/1‑76 (MQ=255)
                                                                         ttttCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCtt  >  2:425117/1‑76 (MQ=255)
                                                                         ttttCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCtt  >  2:1225387/1‑76 (MQ=255)
                                                                         |                                                                           
GGCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCTT  >  REL606/3691771‑3691919

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: ‑