Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,515,9090AC48.9% 25.5 / 34.3 45V260G (GTG→GGG) ypfIpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (14/9);  new base C (0/22);  total (14/31)
Fisher's exact test for biased strand distribution p-value = 6.81e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAA  >  REL606/2515837‑2515981
                                                                         |                                                                          
aacaacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGAccacca                                                                         <  2:3696786/76‑1 (MQ=255)
 acaacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccag                                                                        <  1:2610133/76‑1 (MQ=255)
   aacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcc                                                                      <  1:2433761/76‑1 (MQ=255)
   aacGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGCCGCAGGTATGGCTGCGGCTTCATCGACCAccagcc                                                                      <  2:1562581/76‑1 (MQ=255)
    acGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcca                                                                     <  1:2166013/76‑1 (MQ=255)
    acGTTCGAGGAAAACGCGATACCAGTTGAGGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagcca                                                                     <  2:3992043/76‑1 (MQ=255)
     cGTTCGAGTAAAACGCGATCCCAGTTGATGCAACAATGGGGCAGGTATGGCTGCGGCTCCATCGACCAccagccag                                                                    <  2:2704345/76‑1 (MQ=255)
     cGTTCGAGGAAACCGCGATCCCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCAccagccag                                                                    <  1:2073816/76‑1 (MQ=255)
           agaaaacCGCAACACCAGTTGATGCAAAAATGGCGCAGGTATGGCGGCGCCTTCATCGACCCCCAGCCAGTCGGCt                                                             <  2:2600680/73‑1 (MQ=255)
            ggAAACCCCGATCCCAGTTGATGCAACAAGGGCGCAGGTTTGGCGGCGGCTTCATCGCCCCCCAGCCAGTCGGCtt                                                            <  2:1358715/76‑1 (MQ=255)
            ggAAAACGCTATCCCAGTTGATGCAACAAGGGCGCAGGTAGGCCTGCGCCTTCACCGCCCCCCAGCCAGTCGGCtt                                                            <  2:2002979/76‑1 (MQ=255)
               aaaCGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTTCTCTACCACCCGCCAGTGGGCGTgtt                                                         >  1:2751897/1‑74 (MQ=255)
                  cgcgACACCATTTGATGCAACAAGGGCGCGGGTAGGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCAt                                                      <  1:4440026/76‑1 (MQ=255)
                      atcccATTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCt                                                  <  2:2858372/73‑1 (MQ=255)
                      aTACCAGTGGATGCAAAAAGGGCGCAGGTATGGCTGCGCCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCt                                                  <  2:3050233/76‑1 (MQ=255)
                       tACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTg                                                 <  1:3022381/76‑1 (MQ=255)
                       tACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCCTTATTGACCACCAACCAGGCGGCCTTGCACTTCGCt                                                  >  1:4943090/1‑76 (MQ=255)
                       tACCAGTTGATGCAACAAGGGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTg                                                 <  1:3704639/76‑1 (MQ=255)
                       acccAGTTGATCCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTg                                                 <  1:2347648/74‑1 (MQ=255)
                       acccAGTTAATGCAACAATGGCGCGGGTAGGGCTGCGTCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTg                                                 <  1:4905397/74‑1 (MQ=255)
                         ccAGTTGAGCCACAATGGGCGCAGGTTTGCCGGCGTCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGc                                               <  1:4177546/76‑1 (MQ=255)
                          cAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCt                                              <  2:4675110/76‑1 (MQ=255)
                              tgatgaAACAAGGGCGCAGGTATGGCTGCGGCTTCATCGACCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACa                                         <  2:718631/76‑1 (MQ=255)
                               gaggCAACAAGGGCGCAGGTATGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACaa                                        <  1:1508493/73‑1 (MQ=255)
                               gATGCAAAAAGGGCGCGGGTAGGGCTGCGGCTTCATCGCCCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACaa                                        <  1:3246139/76‑1 (MQ=255)
                               gATCCAACAGCGGCGCAGGTGTGGTGGCGCCTTCACCGACCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACaa                                        <  2:1819830/76‑1 (MQ=255)
                                 ggCAACAGGGGCGCGGTTGTGGCTGCGGCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAgg                                      <  1:1361733/75‑1 (MQ=255)
                                          gggCGCAGGTAGGGCGGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGgc                              <  2:2148673/75‑1 (MQ=255)
                                             cgcgGGTATGGCGGCGCCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCa                           <  1:3451779/76‑1 (MQ=255)
                                                gggTGTGGTTGCGCCTTCCCCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATa                        <  1:914762/75‑1 (MQ=255)
                                                gggTAGGGCGGCGGCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATa                        <  1:2024234/75‑1 (MQ=255)
                                                  ttATGGCGGCGCTTTCAACCCCCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATaaa                      <  1:270247/75‑1 (MQ=255)
                                                  gTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTTCTCATCGCTGGCTAACAAAGGATCCCGCGCAATaaa                      >  1:4436648/1‑76 (MQ=255)
                                                  gTATCGCTGCGGCTTCGTCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATaaa                      <  1:305283/76‑1 (MQ=255)
                                                   tttGGCTGCGCCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAg                     <  2:1637176/74‑1 (MQ=255)
                                                   ttgggCTGCGCCTCCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAg                     <  1:95727/73‑1 (MQ=255)
                                                   tATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTTCTCCTCGCTGGGTTAAAAAGCATCCCGCGCCATaaaa                     >  1:2017659/1‑75 (MQ=255)
                                                   tATGGCGGGGGCTTCCTCGCCCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAg                     <  2:1064041/76‑1 (MQ=255)
                                                     tGGCTGCGGCTTCATCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCg                   <  2:4899165/76‑1 (MQ=255)
                                                      ggCTGCGGCTTCATCCACCACCAGCCAGTCGGCTTGCTCATCGCCGGCGAACAAAGGCACCCGGGCAATAAAGCGa                  >  2:1329825/1‑76 (MQ=255)
                                                       gcggCGGCTTCACCCCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGaa                 <  1:4614099/73‑1 (MQ=255)
                                                       gCTGCGCCTTCAGCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGaa                 <  1:3774394/76‑1 (MQ=255)
                                                         tGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGGATCCGGCGCCATAAAGCGGAAc               >  1:1839335/1‑76 (MQ=255)
                                                         tGCGGCTTCATCCACCACCAGCCAGTCGGCTTGATCATCGCTGGCTAACAAGGGATCCGGCCGCATAAAACGaaaa               >  1:3514069/1‑75 (MQ=255)
                                                         tGCGGCTTCACCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAAc               <  1:2856372/76‑1 (MQ=255)
                                                             gCTTCGCCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTtct           <  2:3304998/76‑1 (MQ=255)
                                                             gCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGTTAACAAGGCATCCGGCGCAATAAAGCGGAACTtct           >  1:2925931/1‑76 (MQ=255)
                                                               ttCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAAATTCTCg         >  1:2875563/1‑76 (MQ=255)
                                                               ttCATCGACCACCAGCCAGTCAGCTTGGTCATCGATGGCTAACGACGGATCCCGAGCAAAAAAGTGAAACTtctct         >  2:4902746/1‑75 (MQ=255)
                                                                 cATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGcc       >  1:1018821/1‑76 (MQ=255)
                                                                  aTCGACCACCAGCCAGTCGGCTTGTCCATGGCTGGCTAACAAGGCATCCGGGGCAATAGAGCGAAACTTTTCGccc      >  2:794515/1‑76 (MQ=255)
                                                                   tCGCCCCCCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCg     <  2:1961667/76‑1 (MQ=255)
                                                                    cGACCACCCGCCAGTCGGCTTGCTCATCGCTGGCTAACACAGCATCCGGCTCAAATAAGAGAAACTTCTCGccctc    >  2:4975555/1‑74 (MQ=255)
                                                                    cGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGc    >  1:2706907/1‑76 (MQ=255)
                                                                    cGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGc    >  1:20901/1‑76 (MQ=255)
                                                                    cGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGc    >  2:524353/1‑76 (MQ=255)
                                                                      accaccAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCaa  >  1:3211739/1‑76 (MQ=255)
                                                                         |                                                                          
AACAACGTTCGAGGAAAACGCGATACCAGTTGATGCAACAATGGCGCAGGTATGGCTGCGGCTTCATCGACCACCAGCCAGTCGGCTTGCTCATCGCTGGCTAACAAGGCATCCGGCGCAATAAAGCGAAACTTCTCGCCCGCAA  >  REL606/2515837‑2515981

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