Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,417,4950TA20.3% 176.5 / 10.8 59L24F (TTA→TTTemrKEmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component
Reads supporting (aligned to +/- strand):  ref base T (39/8);  new base A (0/12);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 1.13e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTAT  >  REL606/2417421‑2417569
                                                                          |                                                                          
tAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaataaata                                                                           >  1:1516092/1‑76 (MQ=255)
tAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaataaata                                                                           >  1:183857/1‑76 (MQ=255)
tAGGCGCCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCaataaata                                                                           >  1:183853/1‑76 (MQ=255)
 aGGCGTCATCTGTCCTAATCATGTCTTCTAATTCCATTGACCAATAGGCACAGGCACCTGAAAACGCAataaataa                                                                          >  2:868220/1‑76 (MQ=255)
 aGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAataaataa                                                                          >  2:868239/1‑76 (MQ=255)
 aGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAataaataa                                                                          >  2:1229417/1‑76 (MQ=255)
    cGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAAc                                                                       >  1:109483/1‑76 (MQ=255)
    cGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAAc                                                                       >  1:1616074/1‑76 (MQ=255)
     gTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACt                                                                      <  2:517300/76‑1 (MQ=255)
      tCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTa                                                                     <  1:979376/76‑1 (MQ=255)
      tCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTa                                                                     <  1:162798/76‑1 (MQ=255)
       tcccTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACTGAAAACGGCATAAATAAAaccacc                                                                   >  2:890987/4‑76 (MQ=255)
        aTCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTAcc                                                                   >  1:1474391/1‑76 (MQ=255)
         tCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCg                                                                  >  1:1128062/1‑76 (MQ=255)
         tCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCg                                                                  >  1:1128041/1‑76 (MQ=255)
         tCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCg                                                                  <  1:59053/76‑1 (MQ=255)
          cTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGc                                                                 <  1:1335837/76‑1 (MQ=255)
           tGTACTAATCATGTCTTCTAATTCCCTTGACCACTAGGCATAGGCACCTCACAACGCAATAAATACAACTAACGcc                                                                >  1:47637/1‑76 (MQ=255)
           tGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACt                                                                      <  1:1087449/70‑1 (MQ=255)
           tGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACt                                                                      >  2:1087449/1‑70 (MQ=255)
           tGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGcc                                                                >  1:47601/1‑76 (MQ=255)
           tGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCACAGGCACCTGAAAACGCAATAAATAAAACTACCGcc                                                                >  1:47606/1‑76 (MQ=255)
            gTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCa                                                               >  2:1195618/1‑76 (MQ=255)
                tAATCATGTTTTCAAATTCCATTGACAAATAGGCATAGGCCCCTAAAAACGAAAAAAAAAAAACTACCGCCAATaa                                                           <  2:1831450/76‑1 (MQ=255)
                 aaTCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATaaa                                                          >  2:297130/1‑76 (MQ=255)
                  aTCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAg                                                         >  1:201347/1‑76 (MQ=255)
                  aTCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAg                                                         >  1:201356/1‑76 (MQ=255)
                  aTCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAg                                                         >  1:361977/1‑76 (MQ=255)
                  aTCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAg                                                         >  1:480835/1‑76 (MQ=255)
                   tCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGa                                                        >  1:782469/1‑76 (MQ=255)
                       gTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAg                                                    >  2:439311/1‑76 (MQ=255)
                       gTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAg                                                    >  2:439321/1‑76 (MQ=255)
                         nttctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTa                                                  >  1:545308/2‑76 (MQ=255)
                         cttctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTa                                                  >  2:1773757/1‑76 (MQ=255)
                         cttctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTa                                                  >  1:1617004/1‑76 (MQ=255)
                         cttctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTa                                                  >  1:1616990/1‑76 (MQ=255)
                          ttctAATTCCATTGACAAATAGGCATAGGCACCAAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTAt                                                 <  1:1286633/76‑1 (MQ=255)
                          ttctAATCCCATTGACAAATAGGCATAGGCACCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTAt                                                 <  1:1090374/76‑1 (MQ=255)
                           tctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTAtt                                                >  2:360172/1‑76 (MQ=255)
                           tctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTAtt                                                <  2:1649719/76‑1 (MQ=255)
                           tctAATTCCATTGACCAATAGGCATAGGCACCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTAtt                                                <  2:1210118/76‑1 (MQ=255)
                             tAATTCCATTGACAAATAGGCATAGGCACCTAAAACCGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTAtttt                                              <  1:1472548/76‑1 (MQ=255)
                              aaTTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTttc                                             >  1:936406/1‑76 (MQ=255)
                                ttCCATTGACCAATAGGCATAGGCACCTGAAAACGAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTttctt                                           <  2:1459760/76‑1 (MQ=255)
                                 tCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTttcttc                                          >  2:543415/1‑76 (MQ=255)
                                   cATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTg                                        >  2:1784631/1‑76 (MQ=255)
                                   cATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTg                                        >  1:1806256/1‑76 (MQ=255)
                                       gACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAg                                    >  2:1538308/1‑76 (MQ=255)
                                       gACCAAAAGGAATAGGCACCTAAAAACGAAAAAAAAAAACCTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAg                                    <  2:461270/76‑1 (MQ=255)
                                       gACAAATAGGCAGAGGCCCCAAAAAACAAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAg                                    <  1:1292709/76‑1 (MQ=255)
                                         ccAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGaa                                  >  1:1256875/1‑76 (MQ=255)
                                         ccAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGaa                                  >  1:121445/1‑76 (MQ=255)
                                         ccAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTCTCTTCTGTTAGaa                                  >  1:798348/1‑76 (MQ=255)
                                         caaaTAGGCATAGGCACCTGAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGaa                                  <  2:571991/74‑1 (MQ=255)
                                          cAATAGGCATAGGCACCTGAAAACGCAATAAATAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAAt                                 <  2:1514509/76‑1 (MQ=255)
                                          cAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAAt                                 >  1:834740/1‑76 (MQ=255)
                                            ataggcataggcACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGt                               >  2:250358/1‑76 (MQ=255)
                                            ataggcatagccACCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGt                               <  2:543668/76‑1 (MQ=255)
                                             taggcataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtt                              >  2:705779/1‑76 (MQ=255)
                                              aggcataggcaCCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGttt                             <  1:315505/76‑1 (MQ=255)
                                                gcataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGttttt                           >  1:1313611/1‑76 (MQ=255)
                                                 cataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttttt                          >  2:1719045/1‑76 (MQ=255)
                                                 cataggcaCCTGAAAACGCAATAAATAAAACTACCGCAAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttttt                          >  2:1719015/1‑76 (MQ=255)
                                                  ataggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTa                         >  1:545312/1‑76 (MQ=255)
                                                   taggcaCGGAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAt                        <  1:630803/76‑1 (MQ=255)
                                                   taggcaCCTAAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAt                        <  1:1730780/76‑1 (MQ=255)
                                                    aggcaCCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAtt                       >  1:897019/1‑76 (MQ=255)
                                                    aggcaCCTAAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTAtt                       <  2:1880228/76‑1 (MQ=255)
                                                        cccTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGaa                   <  2:1762852/75‑1 (MQ=255)
                                                         ccTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAt                  >  1:787148/1‑76 (MQ=255)
                                                          ctaaaaaCGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                 <  1:1691325/73‑1 (MQ=255)
                                                          cTGAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                 <  1:1691309/76‑1 (MQ=255)
                                                           tGAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTa                <  2:109179/76‑1 (MQ=255)
                                                           gaaaaaCGAAAAAAATAAAACAACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTa                <  2:860385/74‑1 (MQ=255)
                                                            aaaaaCGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTaa               <  2:993756/75‑1 (MQ=255)
                                                            aaaaaCAAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTaa               <  2:794333/75‑1 (MQ=255)
                                                                aCGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTg           <  1:1891863/76‑1 (MQ=255)
                                                                aCGAAAAAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTg           <  1:1163528/76‑1 (MQ=255)
                                                                    aaaaaaaaaaaCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTcc       <  2:960726/73‑1 (MQ=255)
                                                                     ataaaaaaATCAGCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCa      <  2:451665/76‑1 (MQ=255)
                                                                     aaaaaaaaaaCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCa      <  1:297130/74‑1 (MQ=255)
                                                                     aaaaaaaaaaCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCa      <  1:572036/74‑1 (MQ=255)
                                                                      taaataaaACTACCGCCAATAAAGAAAAGTATTTTATTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCAc     >  2:133253/1‑76 (MQ=255)
                                                                       aaaaaaaaCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACt    <  1:635951/76‑1 (MQ=255)
                                                                       aaaaaaaaCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACt    <  1:505316/76‑1 (MQ=255)
                                                                        aataaaACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACta   >  1:878961/1‑76 (MQ=255)
                                                                        aataaaACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACta   >  1:26828/1‑76 (MQ=255)
                                                                         ataaaACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACtat  >  1:1713799/1‑76 (MQ=255)
                                                                          |                                                                          
TAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTAT  >  REL606/2417421‑2417569

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