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→GGG) flcfluoride 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

CACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGT  >  NC_000913/657745‑658038
                                                                                                                                                        |                                                                                                                                                
cacaCTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTc                                                                                                                                                 >  2:767911/1‑151 (MQ=255)
cacaCTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTc                                                                                                                                                 >  1:696557/1‑151 (MQ=255)
     tGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGt                                                                                                                                             >  1:719195/1‑150 (MQ=255)
        aTCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTacca                                                                                                                                         <  2:125302/151‑1 (MQ=255)
        aTCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTacc                                                                                                                                          >  2:407905/1‑150 (MQ=255)
         tCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTacca                                                                                                                                         >  2:286678/1‑150 (MQ=255)
               ttCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCa                                                                                                                                   >  2:501233/1‑150 (MQ=255)
                           ctaacccATGCGAACCCTTTCCTTTTGATAAATCCCCCAATCGGTTTGGT‑TGTCAGCGCCCCCAA‑GGAATCCCCTGATTGCAGGGGTTTAAATCCAATCCTAAACACCC‑ATCTCCCCACGTTCCCCGTCCCACCGCCAATAAAAACTGc                                                                                                                        <  2:599927/144‑1 (MQ=255)
                              aaCCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTaaa                                                                                                                    >  2:50743/1‑150 (MQ=255)
                               aCCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCCCGCTTCCCGTACCACCGCCAATAAAAACTGCTaaa                                                                                                                    >  2:231757/1‑149 (MQ=255)
                                  aTGCGAATCCTATTCCTATGATGAATGCCCAAATCAGGTTTGC‑TGTCAGCGCCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCCCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAaagaag                                                                                                               <  1:508596/151‑1 (MQ=255)
                                      gAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGt                                                                                                           >  1:523486/1‑151 (MQ=255)
                                         tcctattcctatGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAac                                                                                                        >  2:47125/1‑151 (MQ=255)
                                          cctattcctatGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAac                                                                                                        >  2:332132/1‑150 (MQ=255)
                                           ctactcctatGATAAATCCCCCAATCGGGTTTGC‑TGTCCGCCTCCCAAACGGATTCCCGGAA‑GGCGGTGGTTAAAATCCTATTTTAAACACAC‑ATCGAGCCGCTCTACCCGGCCCGCCGCCAATAAAAACTGCTAAAAGAAGTTGTaa                                                                                                         <  2:721950/148‑1 (MQ=255)
                                            tattcttatAATAAATGCCCCAACCAGGTTTGC‑TTTCACGGCCCCCAGCGGAACCGCCTGA‑TGCAGTGGGTAAAACCGAATAATAAACACCA‑ATAGCCCCACGCTTCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTaa                                                                                                         <  2:108006/147‑1 (MQ=255)
                                                   atgatgAATGCCCCAATCAGTTTTGC‑TGCCAGCCCCCCCACCGGAATCGCCTGA‑TGCAGTGGGTTAAATCCAATACTAAACACCC‑ATCTCGCCCCTCTCCCCGCACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCa                                                                                                  <  2:21937/147‑1 (MQ=255)
                                                             gccacaataGGTTTGCTTGCCAGCGCCCCGACGGGAGTCGCCTGA‑TGCATGGGGTTAACCCGAAAACTAACCACCC‑ATCACCCCGCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTa                                                                                       <  2:683479/141‑1 (MQ=255)
                                                             ccccAATCGGTTTTGT‑GCCCAGCGTCCCCAACGAACTCCCCGGT‑TGCGGTGGGTAAAATCGAATACTAAACACCC‑ATCTCGCCCCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTattt                                                                                    <  2:344188/151‑1 (MQ=255)
                                                             ccccAATCAGTTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGGGGTTTAAATCGCATACTTAACAGCC‑ACCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTat                                                                                      <  1:398647/149‑1 (MQ=255)
                                                             cccaagtCAGGTTTGC‑TGTCAGCCCCCCAACCGAAAGCGCTTGA‑GGAAGGGGGTAAAACCAAATACTAAACAGCA‑ATCCCGCCCCGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatt                                                                                     <  2:479525/144‑1 (MQ=255)
                                                                 aagcgggTTTGC‑TGTCGGCGCCCCCAACGAACTCGCTTGG‑TGCGGGGGTTAAAAACGCATACAAACCGCCC‑ATCTCCCCTCTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatttga                                                                                  <  2:75157/144‑1 (MQ=255)
                                                                        tttGC‑TGTCAGCGCCCCCAACGGAATCGCTTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ACCTCGCCACGCTCCCCTTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTAt                                                                           <  1:841338/149‑1 (MQ=255)
                                                                        tttGC‑TATCAGCGCCCCCAACGACACCGCCTAA‑GGCAGGGGTTAAAAACGCATAATTAAAACCA‑ATATCGCCTCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATg                                                                          <  2:234837/150‑1 (MQ=255)
                                                                              tGTCAGCTCCCCCAACGGAATCCCTTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCTTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg                                                                      <  1:836402/149‑1 (MQ=255)
                                                                              tGTCAGCGCCCCAAACGGAATCGCCTGA‑TGCAGTGGTTAAAATCGCATACTAAACAGCC‑ATCCCCCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg                                                                      <  1:323987/149‑1 (MQ=255)
                                                                                            acaacacccccTGA‑CGCGGGGGTTTAAACCAAATAC‑TAACCGCCAATCTCCCAACTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATaa                                                                  <  2:108633/130‑1 (MQ=255)
                                                                                             cgaaaTCGCCTGA‑TCCAGTGGTTAAAATCAAATAATTAACACCA‑ATCTCGCCACGCTCCCCGAACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGc                                                         <  2:217854/144‑1 (MQ=255)
                                                                                              gaaaTCCCTTGA‑TGCGGTGGTTTAAATCGAATACTAAACAGCC‑TTCCCGCCCCGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCg                                                       <  1:768106/146‑1 (MQ=255)
                                                                                                aaTCTCGTGA‑TGCAGTGGGTTAAATCACATACTTAACACCC‑ATCCCGCCACTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGa                                                      <  1:147249/147‑1 (MQ=255)
                                                                                                     ccTGA‑TGCGGTGGGTTAAATCGCATCCTTAACAGCC‑ATCTCCCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGaataat                                                 <  1:711459/147‑1 (MQ=255)
                                                                                                           tGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTgg                                        >  2:652249/1‑151 (MQ=255)
                                                                                                              attGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCg                                      >  2:615940/3‑150 (MQ=255)
                                                                                                              aGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCg                                      >  2:740837/1‑150 (MQ=255)
                                                                                                                  ggTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAAttt                                 >  2:612620/1‑151 (MQ=255)
                                                                                                                      aaaTCGCATTCTTAGCATCC‑ATCTCGCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTAc                              <  1:286766/150‑1 (MQ=255)
                                                                                                                      aaaTCGAATACTAAACAGCC‑ATCCCGCCACGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTAc                              <  2:696376/150‑1 (MQ=255)
                                                                                                                       aaTCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACg                             >  1:269572/1‑150 (MQ=255)
                                                                                                                        aTCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGAc                           >  1:311514/1‑151 (MQ=255)
                                                                                                                                  tAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagg                  >  1:180802/1‑150 (MQ=255)
                                                                                                                                  aaaCACCC‑ATCACCCCCCGCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagg                  <  1:673089/149‑1 (MQ=255)
                                                                                                                                   aaCACCC‑CTCCCCCCCCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagga                 <  1:47133/150‑1 (MQ=255)
                                                                                                                                        cc‑ATCTCGCCGCTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTg             <  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)
                                                                                                                                                        |                                                                                                                                                
CACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGC‑TGTCAGCGTCCCCAACGGAATCGCCTGA‑TGCAGTGGGTTAAATCGCATACTTAACAGCC‑ATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGT  >  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: ‑