Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,2680TA39.3% 24.7 / 10.1 29L24F (TTA→TTTemrKtripartite efflux pump membrane fusion protein EmrK
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base A (0/11);  total (9/20)
Fisher's exact test for biased strand distribution p-value = 3.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTGCAAT  >  NC_000913/2483024‑2483508
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
gAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACt                                                                                                                                                                                                                                            >  2:311788/1‑251 (MQ=255)
                                   ggTGCGATTGGGGCATCAGTTTTATCCAGTGAATCTAAAATGTCCCCTGGTCGAACGTGGTTCGTTTCTTTATTATTAACGACAGTGACACTCCCTGAGATTTGTGAAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATTTCTTCTAATTCCATTGACAAATAGAAATAGCCACCTGAAACCGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                                                                                                                                         <  2:516614/250‑1 (MQ=255)
                                   aGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                                                                                                                            >  2:70290/1‑248 (MQ=255)
                                      gggATAGTGGCATAAGTTTAATCCAGTGAAACTAAAATGTCACGTGGTCGAACGTAGTTCGTATTTTTATAATAAACAAAAGTGACACTACCGGAGACTTGTACAGAAATTGAACCTGCATCCCCCGTGAAATAGGCGTAATCTGTACTAATCATGTCTTATAATTCAAGGGACAAATAGGAAAAGGCACCTAAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATg                                                                                                                                                                                                      <  1:311788/249‑1 (MQ=255)
                                          gaggggcatatgTTTTATCCAGTGAAACTAAAATGTCCCTTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTTCCTGAGACTTGTGCAGAAATTGGATCTGCTTTCCCCGGAACATAGGCGTCATCTGTACTAAACATGTCTTCTATTTCCATTGCCAAATAGGCATAGGCACTTGAAAAGGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttt                                                                                                                                                                                                  <  1:230678/240‑1 (MQ=255)
                                                 aTCAGTTTTACCCAGTAAAACTAAAATGTCACCTTGTCGAACGTGGTTCGTATTTTTATGATTAACGACAGTGACACTCCCTGAGATTTGTGCAGAAATGGGATCTGCATCCCCGGTGACATAGGCGTCATCTTTACTAATCATGTCTTCTATTCCCATTGACAAATAGGCATAGGCACCGAAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTg                                                                                                                                                                                           <  1:518098/251‑1 (MQ=255)
                                                   cAGTTTTATCCAGTGAAACTAAAATGTCACTTTGTCAAACGTAGTTCTTATCTTTATGATTAACGACAGTGACCCTCCCTGATACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATGGGCTTCTTCTGTACTAATCATGTCTTCTAATTCCATTGCCCAATAGGCATAGGCCCCTGAAAACGCAATAAATAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGa                                                                                                                                                                                          <  1:533854/250‑1 (MQ=255)
                                                     gTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCTTCTGTACAAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTAAAAACGCAATAAAAAAAATTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                                                                                                                                       <  1:329094/251‑1 (MQ=255)
                                                         tATCCAGAGAAACTAAAAGGTCACTTGGTAGAACGTATTTCGTATCTTTATGATTAAGAACAGTGTCACTCCCGGAGACTTTAGCAAAAATTGTATCGGCATGCCCCGTGAAATAGGCTACATCGTTACTAACCAGGTCTTAAAATTCAATGCACCAATAGGCAAAGGCTCCTAAAACCAAAAAAAAGAATACTACCGCCAATAAAGAAAAGTATTTTCTTCGGTTAGAATGTTTTTTATTTGAATTaag                                                                                                                                                                                    <  2:397821/250‑2 (MQ=255)
                                                                            gtccctttGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACCCTACCTGAGATTTGTGCAAAAATTGGATCTCAATTCCCGGTGACAAGGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCCCCTGAAAACGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctctc                                                                                                                                                                <  1:263983/245‑1 (MQ=255)
                                                                                  tttgcgaaggtggTTCTTATCTTTATGATTACCGACGTTGACCCTCCCGGGGATTTGTGCAGAATTTGGATCTCCATTCCCCGTAACAAAGGCTTCATCTGACTTAATAATGTCTTCTAATCCAATTGACCAAAAGGAATAGGCACTTGAAACCGCAATAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc                                                                                                                                                           <  2:405488/238‑1 (MQ=255)
                                                                                                 cGTATCTTTATGATTACCGACAGTGACACTCCCTGAGACTTGTGCAAAAATTGGTTCTGCATTCCCCGTGACATAGGCGTCATCTGTACAAATCATTTCTTCTAATTCCAGTGACCAATAGACAAAGACACCTAAAACCGCAATAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAAttt                                                                                                                                            <  2:218526/250‑1 (MQ=255)
                                                                                                        ttATGATTAACGACAGTGACACTACCTGAGACTTGTGCCGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg                                                                                                                                    >  2:424955/1‑251 (MQ=255)
                                                                                                         tagaataaaCGACAGTGACACTACATGAGGCTTGTCAAAAAATTGGATCAGCATCCCCGGGGCCATAGGGTCCATCTTTACAAATCATGTTTTCTAATTCCATTGACAAATAGGCATAGGCCCCTGAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATCTAtgaattg                                                                                                                                    <  1:541911/243‑1 (MQ=255)
                                                                                                                aaCGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaa                                                                                                                             >  1:377867/1‑250 (MQ=255)
                                                                                                                         gacatgaCCTGAGACTTGTGTAGAATTTGTATCTGCATCCCCCGTGACAGAGGCGTCATCTGTACTAATCATGTCTTCTAATACAATTAACCAATAGGAATAGAAACATAAAAACGCAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcattc                                                                                                                    <  2:28517/244‑1 (MQ=255)
                                                                                                                                   gagaCTTGTGCAGAAATTGGATCTGCATTCCCCGTAACATAGGCGTCTTCTGTACAAATCATGTCTTCAAATTCCATTGACCAAAAGGCATAGGCACCTGAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCat                                                                                                           <  1:71553/249‑1 (MQ=255)
                                                                                                                                         tgggCAGAATTTGGATCTGTATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATGCCATTGCCCAATAGGAATAGACACCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTGCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCAGTCATTCTTCATATCaaa                                                                                                     <  1:52432/246‑1 (MQ=255)
                                                                                                                                              aGAAATTGGATCTGCTTCCCCCGTGACATAGGCGTCATCTGTACTAATCATGTTTTCTAATTCCATTGACAAATAGGCATAGGCACCTAAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTg                                                                                              <  1:34572/251‑1 (MQ=255)
                                                                                                                                                      gATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGtt                                                                                      >  2:500322/1‑251 (MQ=255)
                                                                                                                                                                   cGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATccc                                                                          >  2:394712/1‑250 (MQ=255)
                                                                                                                                                                               gTCATCTTTACAAATCATGTTTTCGATTCCCATTGCCAATTGGGCAGAGGCCCCTAAAAAGGCAAAAAAAAAAACTGCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGTATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTagagcgggggggtgggag                                                             <  1:115070/251‑18 (MQ=255)
                                                                                                                                                                                       gTACTAATCATTTCTTCAAATACAATTGACAAATAGGAAAAGCCACATAAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAt                                                      <  2:138672/251‑1 (MQ=255)
                                                                                                                                                                                       gTACTAATAATGTCTTCTAATTCCATTGACAAATAGGCATAGGAACCTGAAAACGAAAAAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAt                                                      <  1:64416/251‑1 (MQ=255)
                                                                                                                                                                                                   gTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGt                                           >  2:243282/1‑251 (MQ=255)
                                                                                                                                                                                                       tctAATTCCATTGACCAATAGGCAGAGGCACCTGAAAACCAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTg                                          <  2:347524/248‑1 (MQ=255)
                                                                                                                                                                                                                           aggcaaagccaCCAGAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTa                    <  1:243282/251‑1 (MQ=255)
                                                                                                                                                                                                                              cataggcaCCTGTAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAAc                 >  2:228800/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaCGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAgg                                                            >  1:352070/1‑194 (MQ=255)
                                                                                                                                                                                                                                             acgaaaaaaataaaATTACCTCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTGCAAt  <  2:199168/244‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
GAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTGCAAT  >  NC_000913/2483024‑2483508

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