Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,515,9090AC54.9% 59.5 / 106.1 71V260G (GTG→GGG) ypfIpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (22/10);  new base C (1/38);  total (23/48)
Fisher's exact test for biased strand distribution p-value = 9.60e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-10
Rejected as consensus: Frequency below/above cutoff threshold.

GTTAACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATT  >  REL606/2515834‑2515984
                                                                            |                                                                               
gTTAACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGacca                                                                                 <  1:2952770/76‑1 (MQ=255)
  tAACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGaccacc                                                                               <  1:1189814/76‑1 (MQ=255)
  tAACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGaccacc                                                                               <  1:517561/76‑1 (MQ=255)
   aacaacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGAccacca                                                                              <  1:4894199/76‑1 (MQ=255)
    acaacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccag                                                                             <  2:3685617/76‑1 (MQ=255)
    acaacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccag                                                                             <  2:3915840/76‑1 (MQ=255)
    acaacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccag                                                                             <  2:6520452/76‑1 (MQ=255)
     caacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagc                                                                            <  1:2284223/76‑1 (MQ=255)
      aacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcc                                                                           <  2:4735215/76‑1 (MQ=255)
      aacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcc                                                                           <  1:573452/76‑1 (MQ=255)
       acGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcca                                                                          <  2:3467830/76‑1 (MQ=255)
        cGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagccag                                                                         <  1:2984876/76‑1 (MQ=255)
           tCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCCGGTATGGCTGCGGGTTCATCCAACCCCAGCCAGCCg                                                                     >  1:4730430/1‑76 (MQ=255)
            ggagaAAAGCCCGAAACCGGTAGATCCACCAGTGGCGTTGGTAGGCCTGCGGCTTCATCGCCCCCCAGCCAGTCgg                                                                    <  2:304984/75‑1 (MQ=255)
            cGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGGGGCTTTCTCGGCCCCCCACCAGTCgg                                                                    >  1:3045015/1‑76 (MQ=255)
            aggggAAGACGCGACACCAGTAGATGCAAAATTGGCGGTGGTAGGGCGGCGTCTTCATCGCCCCCCAGCCAGTCgg                                                                    <  2:6490858/73‑1 (MQ=255)
             ggagAAAACGGCAATACAAGTGGATGAAACAGGGGGGCAGGTGTGGCTGCGGCTTCATCGCCCCCCAGCCAGTCgg                                                                    <  2:3533895/73‑1 (MQ=255)
             gagaAAACCGCGATACCAGTTGATGCAACAAGGGCGCAGTTAGGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGc                                                                   <  2:56983/76‑1 (MQ=255)
             gagaAAAACGCGATCCCTGTTGATGCAACATTGGCGCAGTGAGGGCTGCGGTTTCATCCACCCCCAGCCAGTCGGc                                                                   <  2:4118392/76‑1 (MQ=255)
             gAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGGGGCTTTCTCTACCACCCGCCAAGCGGc                                                                   >  2:5101456/1‑76 (MQ=255)
               ggAAACGGCGATACCAGTTGAGGAAAAAAGGGCGCAGGTATGGCGGCGGCTTCATCGACCCCCAGCCAGTCGGCtt                                                                 <  2:650091/76‑1 (MQ=255)
                 aaaaCGCGATACCGTTGGATGCAAAATTGCGGAAGGTTTGGTGGCGGCTACATCCCCCCCCAGCCAGTCGGCTTGc                                                               <  2:6732491/76‑1 (MQ=255)
                        cgATCCGAGTTGATGCAACAATGGCGCAGGTATGCCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCg                                                       <  2:2839941/76‑1 (MQ=255)
                        caaTACCGGTTGATGCAACAAGGGCGCAGGTTGCGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCg                                                       <  2:4957911/74‑1 (MQ=255)
                           tcccAGTGGATGCAACAATGGCGCAGTTAGGTCGGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTg                                                    <  2:1629472/74‑1 (MQ=255)
                            aCCAGTTGATGCAACAAGGGCGCAGGTATGGCTGCGTCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTgg                                                   <  1:4348574/76‑1 (MQ=255)
                            aCCAGTTGAGGCAACAAGGGCGGAGGTAGGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTgg                                                   <  2:3181454/76‑1 (MQ=255)
                               atttGATGCAACAAGGGCGCAGTTATGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTa                                                <  2:4559594/74‑1 (MQ=255)
                               aGTTGATGCAACAAGGGCGCGGGTATGCCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTa                                                <  2:3947446/76‑1 (MQ=255)
                                gTTGATGCAACAATGGCGCAGGTAGGGCTGCGTTTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTaa                                               <  1:3764371/76‑1 (MQ=255)
                                gTTGAGCCAACAATGGCGCGGGTTGGGCTGCGTCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTaa                                               <  1:6948083/76‑1 (MQ=255)
                                 ttGATGCAACAATGGCGCAGGTATGGCTGCGGCTTTATCGACCACCAGCCAGTCGGCCTTCTCCCAGCTGGCtcaa                                              >  1:5201794/1‑73 (MQ=255)
                                 ttGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAAc                                              <  2:6899612/76‑1 (MQ=255)
                                 ttGAGGAAAAAAGGGCGGGGGTAGGGCGGCGGCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAAc                                              <  1:2223970/76‑1 (MQ=255)
                                   gATGCAACAATGGCGCGGGTTTGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACaa                                            <  1:1506785/76‑1 (MQ=255)
                                   gATGCAACAAGGGCGCGGGTGGGCCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACaa                                            <  1:3862784/76‑1 (MQ=255)
                                     ggaaacaaGGGCGCGGGTATGGCTGCGGCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAgg                                          <  2:4433833/73‑1 (MQ=255)
                                      gCAACAATGGCGCGGGTATGGCGGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGc                                         <  2:4252460/76‑1 (MQ=255)
                                      gCAACAAGGGCGGGGGTATGGCTGCGCCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGc                                         <  2:6072447/76‑1 (MQ=255)
                                      gCAACAAGGGCGCGGGTTCGCCGGCGGCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGc                                         <  2:5589307/76‑1 (MQ=255)
                                       aaacaaTGGCGCGGTTATGGCTGCGGTTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCa                                        <  1:432284/75‑1 (MQ=255)
                                        aacaaTGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCCTGCTCATCGCCGGCTAACAAAGCAt                                       >  1:6828835/1‑76 (MQ=255)
                                          aagtGGCGCAGGTGTGCGTGCGTCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATcc                                    <  1:4947464/73‑1 (MQ=255)
                                           aagggCGCAGGTGTGGCGGCGTCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCg                                   <  1:1051282/73‑1 (MQ=255)
                                           aaTGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTTCTCCTCCCCGGGCTAACAAGGCAAcc                                    >  1:1338578/1‑76 (MQ=255)
                                            gtGGCGCGGGTAGGGCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCgg                                  <  2:3447983/75‑1 (MQ=255)
                                            aTGGCGCAGTTTTGGCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCgg                                  <  2:691875/76‑1 (MQ=255)
                                               gcgcAGGTATGGCTGGGGCTTCATCCACCCCCCGCCCGTCGGCGTTCTCACTCCCGGGGTAAAAAGGCACCCgggg                                >  2:6532344/1‑74 (MQ=255)
                                                   aGGTATGGCTGCGGCTTTATCCCCCAACCGCCAGCCGGGCTTCCTCACCGCCGGTGAACAAGAAAGCCGCCGCaga                            >  2:5724194/1‑74 (MQ=255)
                                                    gggaTGCCTGCGGTTTCCTCACCCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATaa                          <  2:4921264/73‑1 (MQ=255)
                                                    ggTATGGCTGCGGCTTTCTCGACCACCCGCCAGGCGGCGTTCTCCTCCCCGGGCTAAAAAGGACTCCGCGGCaaca                           >  2:1988752/1‑74 (MQ=255)
                                                     gTATGGCTGGGCCTTCTCCGCCCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATaaa                         <  1:1674248/76‑1 (MQ=255)
                                                     gTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGGTAAAAAGGGATCCGGGGCCATaaa                         >  1:4308160/1‑76 (MQ=255)
                                                     gTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAAAAAGGCATCCGGCGCAATaaa                         >  1:2583743/1‑76 (MQ=255)
                                                      ttgggCTGCGGCTTCATCCCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAg                        <  1:4023879/73‑1 (MQ=255)
                                                      gagggcggcggcTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAg                        <  2:1221084/73‑1 (MQ=255)
                                                       aTGGCTGCGGCTTTATCGACCACCAACCAGTCGGCTTGCTCCTCCCTGGCGAACAAAGGATCCCGGGGCaaaaaaa                        >  2:6087286/1‑75 (MQ=255)
                                                       aTGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTTCTCATCGCTGGCTAAAAAAGCATCCGGCGCaaaaaaac                       >  1:2581980/1‑74 (MQ=255)
                                                       aTGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCCTCGCTGGCTAACAAAGCCTCCCGCGCAATaaaac                       >  2:3398240/1‑74 (MQ=255)
                                                       aTGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATaaaac                       >  1:6550743/1‑74 (MQ=255)
                                                       aTGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAAAAAGGCCTCCCGCGCAATAAAGc                       >  1:5667450/1‑76 (MQ=255)
                                                        tGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAAAAAAGCACCCGGCGCCACAAAACg                      >  2:3632306/1‑76 (MQ=255)
                                                        gggcggcGCCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCg                      <  1:3854254/75‑1 (MQ=255)
                                                        gggCTGCGGCTCCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCg                      <  1:2138969/75‑1 (MQ=255)
                                                         ggCTGGGTCTTCATCGCCCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGa                     <  2:854715/76‑1 (MQ=255)
                                                          gcggcgcCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGaa                    <  1:4096821/73‑1 (MQ=255)
                                                          gCTGCGTCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGaa                    <  2:252596/76‑1 (MQ=255)
                                                          gCTGCGGCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGaa                    <  1:5690144/76‑1 (MQ=255)
                                                           gggCGGCTTCATCCCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGaaa                   <  1:2968717/74‑1 (MQ=255)
                                                           ccgcgcCTTCATCCACCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGaaa                   <  1:4136011/74‑1 (MQ=255)
                                                           cTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCAACGCCGGCTAACAAGGCATCCGGCGCAATAAAACGaaa                   >  2:1409819/1‑76 (MQ=255)
                                                           cTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCCATAAAGCGaaa                   >  1:4119050/1‑76 (MQ=255)
                                                            tGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAAc                  >  1:3434732/1‑76 (MQ=255)
                                                             gCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACt                 >  1:6597948/1‑76 (MQ=255)
                                                              ggccTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACtt                <  1:4865551/73‑1 (MQ=255)
                                                               ggCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAAATtc               >  2:2573481/1‑76 (MQ=255)
                                                                gCTTCCTCGACCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTtct              <  1:2338893/76‑1 (MQ=255)
                                                                gCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTtct              >  1:2600815/1‑76 (MQ=255)
                                                                gCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTtcc              >  1:2914042/1‑75 (MQ=255)
                                                                gCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAAGCATCCGGCGCAATAAAGCGAAACtttt              >  1:282511/1‑74 (MQ=255)
                                                                  ttCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCg            >  1:3542424/1‑76 (MQ=255)
                                                                  ttCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCg            >  1:1690569/1‑76 (MQ=255)
                                                                  ttCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCg            >  2:562255/1‑76 (MQ=255)
                                                                  ttCATCGACCACCAGCCAGCCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGTGCAATAAAGCGAAACTTCTCg            >  2:3817125/1‑76 (MQ=255)
                                                                   tCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGc           >  2:1833625/1‑76 (MQ=255)
                                                                    cATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGcc          <  2:6216216/76‑1 (MQ=255)
                                                                     aTCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGccc         >  1:5001424/1‑76 (MQ=255)
                                                                       cGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGc       >  1:5708552/1‑76 (MQ=255)
                                                                        gACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCa      >  1:5790562/1‑76 (MQ=255)
                                                                         accaccAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGcca     >  1:6671342/1‑74 (MQ=255)
                                                                         accaccAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCaa     >  1:5895568/1‑76 (MQ=255)
                                                                           caccaGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAAt   >  1:3210616/1‑76 (MQ=255)
                                                                            accaGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCCAAtt  >  2:1889831/1‑76 (MQ=255)
                                                                            |                                                                               
GTTAACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAAATT  >  REL606/2515834‑2515984

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