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

TACTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAAC  >  REL606/4043351‑4043499
                                                                           |                                                                         
tACTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAt                                                                           >  2:1042086/1‑76 (MQ=255)
 aCTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATa                                                                          >  1:264951/1‑76 (MQ=255)
   tGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                        <  2:536214/76‑1 (MQ=255)
   tGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                        >  2:1265976/1‑76 (MQ=255)
   tGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATaaa                                                                        <  1:1566735/76‑1 (MQ=255)
     gCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgc                                                                      >  1:1253847/1‑76 (MQ=255)
     gCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgc                                                                      <  1:951438/76‑1 (MQ=255)
     gCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgc                                                                      <  1:461842/76‑1 (MQ=255)
     gCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgc                                                                      >  1:170369/1‑76 (MQ=255)
         tATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAgccgcc                                                                  <  2:1148610/76‑1 (MQ=255)
          aTGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCa                                                                 >  2:979184/1‑76 (MQ=255)
          aTGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCa                                                                 <  2:31957/76‑1 (MQ=255)
            gTCGGGTTGGTGTTGTGCCCAAATAATGCTCTTTGCCCCTTTTGAAGCACACCATAAAAAAAAAAAGGCCGCCAAt                                                               <  2:1428575/76‑1 (MQ=255)
            gTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAAt                                                               >  2:1352689/1‑76 (MQ=255)
            gTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAAt                                                               >  1:1459609/1‑76 (MQ=255)
             tCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATg                                                              >  2:223242/1‑76 (MQ=255)
              cgggTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATgg                                                             <  2:671498/74‑1 (MQ=255)
              cgggTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATAATAAAAAAAAAGCCGCCAATgg                                                             <  2:976921/74‑1 (MQ=255)
                ggTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGgag                                                           >  1:4407/1‑76 (MQ=255)
                 gTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGgaga                                                          >  2:1724240/1‑76 (MQ=255)
                  ttGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGaa                                                         <  2:1482212/76‑1 (MQ=255)
                   tGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAAc                                                        >  2:1333471/1‑76 (MQ=255)
                    ggTGCTGTGTCCACATAATGCTCTTTGCCGCCTCTGTAGCACAGCATAAAAAAAAAAAAGCCGCCAATGGAGAACa                                                       <  1:498264/76‑1 (MQ=255)
                     gTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACAt                                                      >  1:1652168/1‑76 (MQ=255)
                       gCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCACTGGAGAACATaa                                                    >  2:767599/1‑76 (MQ=255)
                       gCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATaa                                                    <  2:450894/76‑1 (MQ=255)
                        cTGTGCCCAAATATTCTTCTTTGCCCCCTCTGTAGCACAGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAg                                                   <  1:1225625/76‑1 (MQ=255)
                         tgtgTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGt                                                  >  1:1862985/1‑76 (MQ=255)
                         tgtgTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGt                                                  >  1:1182600/1‑76 (MQ=255)
                         tgtgTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGt                                                  >  2:1120511/1‑76 (MQ=255)
                          gtgtCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTg                                                 >  2:294936/1‑76 (MQ=255)
                           tgtCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGc                                                >  1:1069567/1‑76 (MQ=255)
                            gtCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGcc                                               >  2:1038024/1‑76 (MQ=255)
                                acaTAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGcac                                           >  1:739809/1‑76 (MQ=255)
                                 caTAATGCTCTTTGCCGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGcaca                                          <  2:458988/76‑1 (MQ=255)
                                 caTAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGcaca                                          <  2:274167/76‑1 (MQ=255)
                                  aTAATGCTCTTTGCCGCCTCTGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAg                                         <  1:1451806/76‑1 (MQ=255)
                                  aTAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAg                                         >  1:332824/1‑76 (MQ=255)
                                   tAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAgg                                        >  2:1098039/1‑76 (MQ=255)
                                   tAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAgg                                        >  1:957192/1‑76 (MQ=255)
                                    aaTGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGc                                       >  2:1937620/1‑76 (MQ=255)
                                    aaTGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGc                                       >  2:1431886/1‑76 (MQ=255)
                                    aaTCCTCTTTGCCGCCTCTGTAGCACAGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGc                                       <  1:691611/76‑1 (MQ=255)
                                      tGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCAt                                     <  2:1918713/76‑1 (MQ=255)
                                      tGCTCGTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCAt                                     <  2:1918737/76‑1 (MQ=255)
                                       gCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATa                                    >  2:1913175/1‑76 (MQ=255)
                                       gCTCTTTGACGCCTCTGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATa                                    <  1:284858/76‑1 (MQ=255)
                                         tctTTGACGCCTCTGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAAc                                  <  1:1067844/76‑1 (MQ=255)
                                          ctTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAAcc                                 <  2:589664/76‑1 (MQ=255)
                                           tttGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCg                                >  1:1302175/1‑76 (MQ=255)
                                           tttGACGCCTCTGAAGCACAGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCg                                <  1:645897/76‑1 (MQ=255)
                                            ttGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGa                               >  2:1103006/1‑76 (MQ=255)
                                            ttGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGa                               >  2:1602010/1‑76 (MQ=255)
                                             tGACGCCTCTGTAGCACAGCATAATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGAc                              <  2:1565165/76‑1 (MQ=255)
                                              gACGCCTCTGTAGCACAGCATAAAAAAAAGAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACg                             <  2:585011/76‑1 (MQ=255)
                                               aCGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACgg                            >  1:24930/1‑76 (MQ=255)
                                               aCGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACgg                            >  2:600392/1‑76 (MQ=255)
                                                  ccTCTGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTg                         <  2:1551565/76‑1 (MQ=255)
                                                   ctctGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGa                        <  2:79047/76‑1 (MQ=255)
                                                   ctctGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGa                        <  2:343311/76‑1 (MQ=255)
                                                   ctctGTAGCACAGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGa                        <  1:1662128/76‑1 (MQ=255)
                                                    tctGTAGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAg                       <  2:1417759/76‑1 (MQ=255)
                                                     ctGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGt                      >  1:1589184/1‑76 (MQ=255)
                                                       gaaGCACAGCATCAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCa                    <  2:376157/74‑1 (MQ=255)
                                                       gTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCa                    >  1:1693006/1‑76 (MQ=255)
                                                        tAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAg                   >  2:1966291/1‑76 (MQ=255)
                                                         aGCACAGCATAAAAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGa                  <  1:104235/76‑1 (MQ=255)
                                                             caGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGaaaaa              >  2:907789/1‑76 (MQ=255)
                                                              aGCATAATAAAAAAAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGaaaaaa             <  2:1540678/76‑1 (MQ=255)
                                                                  tcatAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGa         >  2:890249/1‑76 (MQ=255)
                                                                    atAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACt       <  2:679104/76‑1 (MQ=255)
                                                                     tAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTg      >  1:795889/1‑76 (MQ=255)
                                                                     aaaaaaaaaAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTg      <  1:355461/75‑1 (MQ=255)
                                                                     aaaaaaTAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTg      <  1:403781/75‑1 (MQ=255)
                                                                      aaaaaaaaGCCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGc     <  2:1733018/76‑1 (MQ=255)
                                                                      aaaaaaaaAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTgn     <  1:361623/76‑2 (MQ=255)
                                                                      aaaaaaaaAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGc     <  1:527505/76‑1 (MQ=255)
                                                                      aaaaaTAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGc     <  2:1851577/76‑1 (MQ=255)
                                                                      aaaaaTAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGc     >  1:187919/1‑76 (MQ=255)
                                                                        aaaTAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCaa   >  2:1479385/1‑76 (MQ=255)
                                                                         aaTAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAAc  >  2:914032/1‑76 (MQ=255)
                                                                           |                                                                         
TACTGGCTTTATGTCCGGTTGGTGCTGTGTCCACATAATGCTCTTTGACGCCTCTGTAGCACAGCATCATAAAAATAAAGCCGCCAATGGAGAACATAAGTGCCGCACAGGCATAACCGACGGTTGAGTCAGAAAAAAGCGACTGCAAC  >  REL606/4043351‑4043499

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