Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,8940TC40.4% 53.7 / 19.4 48G15G (GGA→GGGflcfluoride efflux channel, dual topology membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/6);  new base C (0/19);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 2.56e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGT  >  NC_000913/657745‑658038
                                                                                                                                                        |                                                                                                                                                
cacaCTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTc                                                                                                                                                 >  2:767911/1‑151 (MQ=255)
cacaCTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTc                                                                                                                                                 >  1:696557/1‑151 (MQ=255)
     tGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGt                                                                                                                                             >  1:719195/1‑150 (MQ=255)
        aTCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTacca                                                                                                                                         <  2:125302/151‑1 (MQ=255)
        aTCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTacc                                                                                                                                          >  2:407905/1‑150 (MQ=255)
         tCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTacca                                                                                                                                         >  2:286678/1‑150 (MQ=255)
               ttCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCa                                                                                                                                   >  2:501233/1‑150 (MQ=255)
                           ctaacccATGCGAACCCTTTCCTTTTGATAAATCCCCCAATCGGTTTGGTTGTCAGCGCCCCCAAGGAATCCCCTGATTGCAGGGGTTTAAATCCAATCCTAAACACCCATCTCCCCACGTTCCCCGTCCCACCGCCAATAAAAACTGc                                                                                                                        <  2:599927/144‑1 (MQ=255)
                              aaCCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTaaa                                                                                                                    >  2:50743/1‑150 (MQ=255)
                               aCCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCCCGCTTCCCGTACCACCGCCAATAAAAACTGCTaaa                                                                                                                    >  2:231757/1‑149 (MQ=255)
                                  aTGCGAATCCTATTCCTATGATGAATGCCCAAATCAGGTTTGCTGTCAGCGCCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCCCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAaagaag                                                                                                               <  1:508596/151‑1 (MQ=255)
                                      gAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGt                                                                                                           >  1:523486/1‑151 (MQ=255)
                                         tcctattcctatGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAac                                                                                                        >  2:47125/1‑151 (MQ=255)
                                          cctattcctatGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAac                                                                                                        >  2:332132/1‑150 (MQ=255)
                                           ctactcctatGATAAATCCCCCAATCGGGTTTGCTGTCCGCCTCCCAAACGGATTCCCGGAAGGCGGTGGTTAAAATCCTATTTTAAACACACATCGAGCCGCTCTACCCGGCCCGCCGCCAATAAAAACTGCTAAAAGAAGTTGTaa                                                                                                         <  2:721950/148‑1 (MQ=255)
                                            tattcttatAATAAATGCCCCAACCAGGTTTGCTTTCACGGCCCCCAGCGGAACCGCCTGATGCAGTGGGTAAAACCGAATAATAAACACCAATAGCCCCACGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTaa                                                                                                         <  2:108006/147‑1 (MQ=255)
                                                   atgatgAATGCCCCAATCAGTTTTGCTGCCAGCCCCCCCACCGGAATCGCCTGATGCAGTGGGTTAAATCCAATACTAAACACCCATCTCGCCCCTCTCCCCGCACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCa                                                                                                  <  2:21937/147‑1 (MQ=255)
                                                             gccacaataGGTTTGCTTGCCAGCGCCCCGACGGGAGTCGCCTGATGCATGGGGTTAACCCGAAAACTAACCACCCATCACCCCGCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTa                                                                                       <  2:683479/141‑1 (MQ=255)
                                                             ccccAATCGGTTTTGTGCCCAGCGTCCCCAACGAACTCCCCGGTTGCGGTGGGTAAAATCGAATACTAAACACCCATCTCGCCCCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTattt                                                                                    <  2:344188/151‑1 (MQ=255)
                                                             ccccAATCAGTTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGGGGTTTAAATCGCATACTTAACAGCCACCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTat                                                                                      <  1:398647/149‑1 (MQ=255)
                                                             cccaagtCAGGTTTGCTGTCAGCCCCCCAACCGAAAGCGCTTGAGGAAGGGGGTAAAACCAAATACTAAACAGCAATCCCGCCCCGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatt                                                                                     <  2:479525/144‑1 (MQ=255)
                                                                 aagcgggTTTGCTGTCGGCGCCCCCAACGAACTCGCTTGGTGCGGGGGTTAAAAACGCATACAAACCGCCCATCTCCCCTCTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatttga                                                                                  <  2:75157/144‑1 (MQ=255)
                                                                        tttGCTGTCAGCGCCCCCAACGGAATCGCTTGATGCAGTGGGTTAAATCGCATACTTAACAGCCACCTCGCCACGCTCCCCTTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTAt                                                                           <  1:841338/149‑1 (MQ=255)
                                                                        tttGCTATCAGCGCCCCCAACGACACCGCCTAAGGCAGGGGTTAAAAACGCATAATTAAAACCAATATCGCCTCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATg                                                                          <  2:234837/150‑1 (MQ=255)
                                                                              tGTCAGCTCCCCCAACGGAATCCCTTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCTTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg                                                                      <  1:836402/149‑1 (MQ=255)
                                                                              tGTCAGCGCCCCAAACGGAATCGCCTGATGCAGTGGTTAAAATCGCATACTAAACAGCCATCCCCCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg                                                                      <  1:323987/149‑1 (MQ=255)
                                                                                            acaacacccccTGACGCGGGGGTTTAAACCAAATACTAACCGCCAATCTCCCAACTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATaa                                                                  <  2:108633/130‑1 (MQ=255)
                                                                                             cgaaaTCGCCTGATCCAGTGGTTAAAATCAAATAATTAACACCAATCTCGCCACGCTCCCCGAACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGc                                                         <  2:217854/144‑1 (MQ=255)
                                                                                              gaaaTCCCTTGATGCGGTGGTTTAAATCGAATACTAAACAGCCTTCCCGCCCCGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCg                                                       <  1:768106/146‑1 (MQ=255)
                                                                                                aaTCTCGTGATGCAGTGGGTTAAATCACATACTTAACACCCATCCCGCCACTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGa                                                      <  1:147249/147‑1 (MQ=255)
                                                                                                     ccTGATGCGGTGGGTTAAATCGCATCCTTAACAGCCATCTCCCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGaataat                                                 <  1:711459/147‑1 (MQ=255)
                                                                                                           tGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTgg                                        >  2:652249/1‑151 (MQ=255)
                                                                                                              attGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCg                                      >  2:615940/3‑150 (MQ=255)
                                                                                                              aGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCg                                      >  2:740837/1‑150 (MQ=255)
                                                                                                                  ggTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAAttt                                 >  2:612620/1‑151 (MQ=255)
                                                                                                                      aaaTCGCATTCTTAGCATCCATCTCGCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTAc                              <  1:286766/150‑1 (MQ=255)
                                                                                                                      aaaTCGAATACTAAACAGCCATCCCGCCACGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTAc                              <  2:696376/150‑1 (MQ=255)
                                                                                                                       aaTCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACg                             >  1:269572/1‑150 (MQ=255)
                                                                                                                        aTCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGAc                           >  1:311514/1‑151 (MQ=255)
                                                                                                                                  tAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagg                  >  1:180802/1‑150 (MQ=255)
                                                                                                                                  aaaCACCCATCACCCCCCGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagg                  <  1:673089/149‑1 (MQ=255)
                                                                                                                                   aaCACCCCTCCCCCCCCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagga                 <  1:47133/150‑1 (MQ=255)
                                                                                                                                        ccATCTCGCCGCTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTg             <  2:592557/149‑1 (MQ=255)
                                                                                                                                               cGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCt     >  1:728598/1‑151 (MQ=255)
                                                                                                                                               cGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCt     >  1:479547/1‑151 (MQ=255)
                                                                                                                                                gCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTg    >  1:172710/1‑151 (MQ=255)
                                                                                                                                                 ccACGCTCCCCGTACCGCCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCt     <  1:408006/149‑1 (MQ=255)
                                                                                                                                                  cACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGt  >  2:148941/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
CACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGT  >  NC_000913/657745‑658038

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