Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,043,4260TA21.0% 161.7 / 11.2 63I200F (ATT→TTT) yihOpredited transporter
Reads supporting (aligned to +/- strand):  ref base T (31/18);  new base A (0/13);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGC  >  REL606/4043351‑4043496
                                                                           |                                                                      
tACTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAt                                                                        >  1:761601/1‑76 (MQ=255)
tACTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAt                                                                        <  1:1475751/76‑1 (MQ=255)
tACTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAt                                                                        >  2:1618650/1‑76 (MQ=255)
 aCTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATa                                                                       <  2:1154079/76‑1 (MQ=255)
 aCTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATa                                                                       >  1:1875959/1‑76 (MQ=255)
  cTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaa                                                                      >  2:1644180/1‑76 (MQ=255)
   tGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                     <  2:795855/76‑1 (MQ=255)
   tGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                     <  1:1467175/76‑1 (MQ=255)
   tGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                     >  2:145563/1‑76 (MQ=255)
    ggCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAg                                                                    >  1:1398990/1‑76 (MQ=255)
    ggCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAg                                                                    >  1:523564/1‑76 (MQ=255)
        ttATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgccgc                                                                >  2:107150/1‑76 (MQ=255)
        ttATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgccgc                                                                <  1:1607114/76‑1 (MQ=255)
        ttATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgccgc                                                                >  1:318662/1‑76 (MQ=255)
         tATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgccgcc                                                               >  2:2214394/1‑76 (MQ=255)
          aTGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCa                                                              >  1:577434/1‑76 (MQ=255)
          aTGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCa                                                              <  1:667208/76‑1 (MQ=255)
          aTGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCa                                                              <  2:96089/76‑1 (MQ=255)
           tGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCaa                                                             <  2:413809/76‑1 (MQ=255)
              ccGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATgg                                                          >  2:630166/1‑76 (MQ=255)
                 gTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGgaga                                                       <  2:277126/76‑1 (MQ=255)
                 gTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGgaga                                                       <  2:1441602/76‑1 (MQ=255)
                   tGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAAc                                                     >  2:2029972/1‑76 (MQ=255)
                    ggTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACa                                                    >  1:558530/1‑76 (MQ=255)
                     gTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCAAAAAAATAAAGCCGCCAATGGAGAACAt                                                   <  2:1057491/76‑1 (MQ=255)
                     gTGCTGTGCCCACATAATGCTCTTTGACGCCTCTGAAGCACAGCATCAAAAAAATAAAGCCGCCAATGGAGAACAt                                                   <  1:1682975/76‑1 (MQ=255)
                      tGCGGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATa                                                  <  1:470827/76‑1 (MQ=255)
                         tgtgTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGt                                               <  2:318662/76‑1 (MQ=255)
                          gtgtCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTg                                              <  1:1971435/76‑1 (MQ=255)
                           tgtCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGc                                             <  1:782572/76‑1 (MQ=255)
                           tgtCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGc                                             <  1:1565041/76‑1 (MQ=255)
                            gtCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGcc                                            >  1:964455/1‑76 (MQ=255)
                               cacaTAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGca                                         >  1:2110888/1‑76 (MQ=255)
                                gcaTAATGCTCGTGGGCGCTTCTGTAGCACAGCTTAATAAAAAGAAATCCGCCAATGGAGAACATAAGTGCCGcac                                        <  1:948546/75‑1 (MQ=255)
                                   tAATGNTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAgg                                     >  2:22662/1‑76 (MQ=255)
                                   tAATGCTCTTTGACGCCTCTGTAGCACAGCATAATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAgg                                     <  1:903838/76‑1 (MQ=255)
                                    aaTGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGc                                    >  1:2188587/1‑76 (MQ=255)
                                    aaTGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGc                                    >  1:921738/1‑76 (MQ=255)
                                      tGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCAt                                  >  2:521863/1‑76 (MQ=255)
                                      tGCTCTTGGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCAt                                  <  2:1734609/76‑1 (MQ=255)
                                       gCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATa                                 <  2:1598311/76‑1 (MQ=255)
                                        ctttttGACGCCTCTGTAGCACAGCATCAAAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATaa                                <  1:741189/73‑1 (MQ=255)
                                         tctTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAAc                               >  1:1166577/1‑76 (MQ=255)
                                         tctTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAAc                               >  1:396921/1‑76 (MQ=255)
                                         tctTTGACGCCTCTGTAGCACAACATAAAAAAAATAAGGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAAc                               <  1:734730/76‑1 (MQ=255)
                                          ctTTGCCGCCTCTGTAGCACAGCATCATAAAAAGAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAAcc                              <  2:755161/76‑1 (MQ=255)
                                          ctTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAAcc                              >  2:1135655/1‑76 (MQ=255)
                                            ttGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGa                            >  1:1686340/1‑76 (MQ=255)
                                             tGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGAc                           >  2:993949/1‑76 (MQ=255)
                                             tGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGAc                           >  1:1165662/1‑76 (MQ=255)
                                             tGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGAc                           >  1:5517/1‑76 (MQ=255)
                                             tGACGCCTCTGTAGCACAGCATCAAAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGAc                           <  2:106545/76‑1 (MQ=255)
                                              gACGCCTCTGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACg                          <  1:707631/76‑1 (MQ=255)
                                               aCGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACgg                         <  1:164683/76‑1 (MQ=255)
                                                 gCCTCTGAAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGtt                       <  1:2154034/76‑1 (MQ=255)
                                                  ccTCTGTAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTg                      <  1:2006966/76‑1 (MQ=255)
                                                  ccTCTGTAGCACAGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTg                      <  1:1522403/76‑1 (MQ=255)
                                                    tctGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAg                    >  2:923007/1‑76 (MQ=255)
                                                     ctGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGt                   <  2:1478033/76‑1 (MQ=255)
                                                     ctGAAGCAGAGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGt                   <  2:1105955/76‑1 (MQ=255)
                                                        tAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAg                <  1:2178759/76‑1 (MQ=255)
                                                         aGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGa               >  1:1256146/1‑76 (MQ=255)
                                                             caGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGaaaaa           >  1:2036618/1‑76 (MQ=255)
                                                             caGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGaaaaa           >  1:1538166/1‑76 (MQ=255)
                                                              aGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGaaaaaa          <  1:1713883/76‑1 (MQ=255)
                                                              aGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGaaaaaa          <  1:548792/76‑1 (MQ=255)
                                                                catcatAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGc        <  2:1871940/76‑1 (MQ=255)
                                                                 atcatAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGATAAAAGCg       >  1:1366754/1‑76 (MQ=255)
                                                                 atcaaAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCg       <  1:534584/76‑1 (MQ=255)
                                                                      aaaaaaaaAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGc  <  2:2029292/76‑1 (MQ=255)
                                                                      aaaaaTAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGc  >  2:1976780/1‑76 (MQ=255)
                                                                           |                                                                      
TACTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGC  >  REL606/4043351‑4043496

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