Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,8940TC38.5% 55.7 / 13.1 39G15G (GGA→GGG) flcfluoride efflux channel, dual topology membrane protein
Reads supporting (aligned to +/- strand):  ref base T (19/5);  new base C (0/15);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 7.71e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGC  >  NC_000913/657745‑658041
                                                                                                                                                     |                                                                                                                                                   
cacaCTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTc                                                                                                                                                    >  1:446984/1‑151 (MQ=255)
 acacTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTcc                                                                                                                                                   >  1:810454/1‑151 (MQ=255)
          cAATGTTCGTCATCCTGCTGAACCATGCGCATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACGTAACAGCCATCTCGCCACGCTTCCCGTaccac                                                                                                                                           <  1:430471/150‑1 (MQ=255)
                   tCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATaa                                                                                                                                  >  2:185186/1‑150 (MQ=255)
                     aTCCTGCTGAACCATGCGAACCCTATCCCTATGATAAATGCCCCAACCAGGTTGCCTGTCAGCGCCCCCAACGAAATCGCCTGATGCAGTGGGTTAAATCGAATACTTAACAGCCATCTCCCCACGCTTCCCGTACCACCGCCAATaa                                                                                                                                  <  1:480073/148‑1 (MQ=255)
                           ctgACCCATGCGAATCCTATCCCTATGATGAATGCCCCAATCGGTTTTGCTGTCAGCGCCCCCAACGGAACCGCCTGATGCAGTGGGTTAAATCGCATACTAAACAGCCATCCCCCCACGCTTCCCGCCCCACCGCCAATAAAAACTGc                                                                                                                           <  1:95012/149‑1 (MQ=255)
                           ctgAACCATGCGAACCCTTTTCCTATGATAAATGCCCCAATCAGTTTTGCTGTCAGCGCCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCCCCCCACGCTCCCCGTCCCACCGCCAATAAAAACTg                                                                                                                            <  2:71115/148‑1 (MQ=255)
                              acccATGCGACTCCTATCCCTATAATAACTGCCCCAAACGGGTTTGCTGCCAGCCCCCCCAACGGAACCCCCTTGTGCAGTGGTTAAAATCGAATACTAAACAGCCATCTCGCCACGCTCCCCGTCCCACCGCCAATAAAAACTGCTa                                                                                                                         <  1:502594/146‑1 (MQ=255)
                              aaCCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTaaa                                                                                                                       >  1:472411/1‑150 (MQ=255)
                                  aTGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAaag                                                                                                                     >  1:581135/1‑148 (MQ=255)
                                    gCGAATCCTATTCTTATGATGAACCCCCCAATCAGTTTTGCTGCCAGCGACCCCAACGGAATCGCCTAAGGCAGTGAATAAAATCACAAAATAAACCGCCATCTCGCCACGTTCCCCGAACCACCGCCAATAAAAACTGCTAAaag                                                                                                                     <  1:183122/146‑1 (MQ=255)
                                     cGATTCCTATTCCTATGATGAATGCCCCAAGCGGGTTTGCTGTCAGCTTCCCCACCGGAACCGCCTGATGTATTGGTTTAAATCGCATATTTAACACCAATCCCGCCCCGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGt                                                                                                                 <  1:584232/149‑1 (MQ=255)
                                      gAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACAGCCAATAAAAACTGCTAAAAGAAGTTGt                                                                                                              >  2:402489/1‑151 (MQ=255)
                                       aaTCCTATTCCTATGATGAATGCCCCAATCAGTTTTGCTGTCAGCGTCCCCAACGAAATCGCCTGATGCAGTGGGTAAAATCGCATATTTAACAGCCATCTCGCCCCGTTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTg                                                                                                               <  2:1766/149‑1 (MQ=255)
                                                cttaagaagaACGCCCAAATCAGGTTGGCTGTAAGCTCCCCCACCGGAACCGCCTAAAGCAGTGGTTTAAAACGAATAATAACCACCAACCGCGCCCCTCTCCCCGACCCCCCGCCAATAAAAACTGCTAAAAGAAGTTGTAaca                                                                                                          <  1:772895/137‑1 (MQ=255)
                                                           tGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTat                                                                                         >  2:393683/1‑151 (MQ=255)
                                                             cccaAATCAGGTTTGCTGTCAGCGCCCCCACCGAAATCGCCTGATGCAGTGGGTAAAATCGAATATTAAACACCCATCTCGCCACTCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTattt                                                                                       <  1:106405/151‑1 (MQ=255)
                                                                 aaccAGTTTTGCTTCCAGCCCCCCCACCGTATCCCCCGGTTGCATTGGTTTAACTCCAATACTAAACACCCATCCCCCCCCCCTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatttga                                                                                     <  2:138583/146‑1 (MQ=255)
                                                                  aTCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatttgattt                                                                                  >  2:471958/1‑151 (MQ=255)
                                                                    cAGTTTTGCTGTCAGCGCCCCCAACGGACCCGCCTGATGCAGGGGTTTAAACCGAATAATAAACACCCATCTCCCCACCCTCCCCTTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTatttgatttg                                                                                 <  1:801034/150‑1 (MQ=255)
                                                                    cAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGt                                                                                >  1:234381/1‑151 (MQ=255)
                                                                     aGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGt                                                                                >  2:60536/1‑150 (MQ=255)
                                                                        tttGCTGCCAGCGCCCCAACCGAAACCGCCTGATGCAGTGGGTTAAACCGCATACTTAACACCAATCTCGCCACTCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGt                                                                                <  1:771591/147‑1 (MQ=255)
                                                                            cTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTg                                                                         >  1:642552/1‑150 (MQ=255)
                                                                                 aGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAg                                                                    >  1:701096/1‑150 (MQ=255)
                                                                                   cccccccAACGGAACCGCCTGAGGCAGTGGTTAAAATCAAATAATAAACAGCCATCCCGCCACGTTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAgtgt                                                                 <  1:612647/148‑1 (MQ=255)
                                                                                   cGTCCCCAACGAAATCGCCTGATGCAGTGGGTTAAATCGAATATTAAACAGCCATCTCGCCACGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATa                                                                      <  1:241663/146‑1 (MQ=255)
                                                                                    gcccccAACGGAATCGCCTGATGCAGTGGTTTAAACCGCCTATTTAACAGCCATCCCGCCTCTCTCCCCGTACCACCGCCAATAAAATCTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGa                                                                        <  2:290486/141‑1 (MQ=255)
                                                                                          accGGATTCGCCTGAGACAGTGGGTTAAACCGAAAATTAAACAGCCATCGCGCCCCCCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCg                                                          <  1:618549/149‑1 (MQ=255)
                                                                                             ggAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGaat                                                       >  2:161135/1‑151 (MQ=255)
                                                                                                          gCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGc                                          >  2:370581/1‑151 (MQ=255)
                                                                                                               gggTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAAt                                      <  2:262354/150‑1 (MQ=255)
                                                                                                                ggTAAAATCCCATACTTAACGTCCATCCCGCCACGTTCCCCGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAAtt                                     <  2:795588/150‑1 (MQ=255)
                                                                                                                          cATTCTTAACAGCCACCTCGCCCCGCTCCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGActct                           <  2:370374/150‑1 (MQ=255)
                                                                                                                               ttAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagg                     >  1:787005/1‑151 (MQ=255)
                                                                                                                                 aaCAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACagg                     >  1:118506/1‑149 (MQ=255)
                                                                                                                                            cGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCt        >  2:592779/1‑151 (MQ=255)
                                                                                                                                            cGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCt        >  2:14586/1‑151 (MQ=255)
                                                                                                                                                  gCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGc  >  2:70538/1‑151 (MQ=255)
                                                                                                                                                   tttCCGGTCCCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGtt    <  2:203763/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCAACGGAATCGCCTGATGCAGTGGGTTAAATCGCATACTTAACAGCCATCTCGCCACGCTTCCCGTACCACCGCCAATAAAAACTGCTAAAAGAAGTTGTAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGC  >  NC_000913/657745‑658041

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