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→TTT) emrKtripartite 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.

GAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTGCAAT  >  NC_000913/2483024‑2483508
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
gAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACt                                                                                                                                                                                                                                            >  2:311788/1‑251 (MQ=255)
                                   ggTGCGATTGGGGCATCAGTTTTATCCAGTGAATCTAAAATGTCCCCTGGTCGAACGTGGTTCGTTTCTTTATTATTAACGACAGTGACACTCCCTGAGATTTGTGAAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATTTCTTCTAATTCCATTGACAAAT‑AGAAATAGCCACCTGAAACCGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa                                                                                                                                                                                                         <  2:516614/250‑1 (MQ=255)
                                   aGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt                                                                                                                                                                                                            >  2:70290/1‑248 (MQ=255)
                                      gggATAGTGGCATAAGTTTAATCCAGTGAAACTAAAATGTCACGTGGTCGAACGTAGTTCGTATTTTTATAATAAACAAAAGTGACACTACCGGAGACTTGTACAGAAATTGAACCTGCATCCCCCGTGAAATAGGCGTA‑ATCTGTACT‑AATCATGTCTTATAATTCAAGGGACAAAT‑AGGAAAAGGCACCTAAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATg                                                                                                                                                                                                      <  1:311788/249‑1 (MQ=255)
                                          gaggggcatatgTTTTATCCAGTGAAACTAAAATGTCCCTTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTTCCTGAGACTTGTGCAGAAATTGGATCTGCTTTCCCCGGAACATAGGCGTC‑ATCTGTACT‑AAACATGTCTTCTATTTCCATTGCCAAAT‑AGGCATAGGCACTTGAAAAGGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttt                                                                                                                                                                                                  <  1:230678/240‑1 (MQ=255)
                                                 aTCAGTTTTACCCAGTAAAACTAAAATGTCACCTTGTCGAACGTGGTTCGTATTTTTATGATTAACGACAGTGACACTCCCTGAGATTTGTGCAGAAATGGGATCTGCATCCCCGGTGACATAGGCGTC‑ATCTTTACT‑AATCATGTCTTCTATTCCCATTGACAAAT‑AGGCATAGGCACCGAAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTg                                                                                                                                                                                           <  1:518098/251‑1 (MQ=255)
                                                   cAGTTTTATCCAGTGAAACTAAAATGTCACTTTGTCAAACGTAGTTCTTATCTTTATGATTAACGACAGTGACCCTCCCTGATACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATGGGCTTC‑TTCTGTACT‑AATCATGTCTTCTAATTCCATTGCCCAAT‑AGGCATAGGCCCCTGAAAACGCAATAAATAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGa                                                                                                                                                                                          <  1:533854/250‑1 (MQ=255)
                                                     gTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑TTCTGTACA‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTAAAAACGCAATAAAAAAAATTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt                                                                                                                                                                                       <  1:329094/251‑1 (MQ=255)
                                                         tATCCAGAGAAACTAAAAGGTCACTTGGTAGAACGTATTTCGTATCTTTATGATTAAGAACAGTGTCACTCCCGGAGACTTTAGCAAAAATTGTATCGGCATGCCCCGTGAAATAGGCTAC‑ATCGTTACT‑AACCAGGTCTTAAAATTCAATGCACCAAT‑AGGCAAAGGCTCCTAAAACCAAAAAAAAGAATACTACCGCCAATAAAGAAAAGTATTTTCTTCGGTTAGAATGTTTTTTATTTGAATTaag                                                                                                                                                                                    <  2:397821/250‑2 (MQ=255)
                                                                            gtccctttGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACCCTACCTGAGATTTGTGCAAAAATTGGATCTCAATTCCCGGTGACAAGGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCCCCTGAAAACGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctctc                                                                                                                                                                <  1:263983/245‑1 (MQ=255)
                                                                                  tttgcgaaggtggTTCTTATCTTTATGATTACCGACGTTGACCCTCCCGGGGATTTGTGCAGAATTTGGATCTCCATTCCCCGTAACAAAGGCTTC‑ATCTG‑ACTTAATAATGTCTTCTAATCCAATTGACCAAA‑AGGAATAGGCACTTGAAACCGCAATAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc                                                                                                                                                           <  2:405488/238‑1 (MQ=255)
                                                                                                 cGTATCTTTATGATTACCGACAGTGACACTCCCTGAGACTTGTGCAAAAATTGGTTCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACA‑AATCATTTCTTCTAATTCCAGTGACCAAT‑AGACAAAGACACCTAAAACCGCAATAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAAttt                                                                                                                                            <  2:218526/250‑1 (MQ=255)
                                                                                                        ttATGATTAACGACAGTGACACTACCTGAGACTTGTGCCGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg                                                                                                                                    >  2:424955/1‑251 (MQ=255)
                                                                                                         tagaataaaCGACAGTGACACTACATGAGGCTTGTCAAAAAATTGGATCAGCATCCCCGGGGCCATAGG‑GTCCATCTTTACA‑AATCATGTTTTCTAATTCCATTGACAAAT‑AGGCATAGGCCCCTGAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATCTAtgaattg                                                                                                                                    <  1:541911/243‑1 (MQ=255)
                                                                                                                aaCGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaa                                                                                                                             >  1:377867/1‑250 (MQ=255)
                                                                                                                         gacatgaCCTGAGACTTGTGTAGAATTTGTATCTGCATCCCCCGTGACAGAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATACAATTAACCAAT‑AGGAATAGAAACATAAAAACGCAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcattc                                                                                                                    <  2:28517/244‑1 (MQ=255)
                                                                                                                                   gagaCTTGTGCAGAAATTGGATCTGCATTCCCCGTAACATAGGCGTC‑TTCTGTACA‑AATCATGTCTTCAAATTCCATTGACCAAA‑AGGCATAGGCACCTGAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCat                                                                                                           <  1:71553/249‑1 (MQ=255)
                                                                                                                                         tgggCAGAATTTGGATCTGTATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATGCCATTGCCCAAT‑AGGAATAGACACCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTGCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCAGTCATTCTTCATATCaaa                                                                                                     <  1:52432/246‑1 (MQ=255)
                                                                                                                                              aGAAATTGGATCTGCTTCCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTTTTCTAATTCCATTGACAAAT‑AGGCATAGGCACCTAAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTg                                                                                              <  1:34572/251‑1 (MQ=255)
                                                                                                                                                      gATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGtt                                                                                      >  2:500322/1‑251 (MQ=255)
                                                                                                                                                                   cGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATccc                                                                          >  2:394712/1‑250 (MQ=255)
                                                                                                                                                                               gTC‑ATCTTTACA‑AATCATGTTTTCGATTCCCATTG‑CCAATTGGGCAGAGGCCCCTAAAAAGGCAAAAAAAAAAACTGCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGTATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTagagcgggggggtgggag                                                             <  1:115070/251‑18 (MQ=255)
                                                                                                                                                                                       gTACT‑AATCATTTCTTCAAATACAATTGACAAAT‑AGGAAAAGCCACATAAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAt                                                      <  2:138672/251‑1 (MQ=255)
                                                                                                                                                                                       gTACT‑AATAATGTCTTCTAATTCCATTGACAAAT‑AGGCATAGGAACCTGAAAACGAAAAAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAt                                                      <  1:64416/251‑1 (MQ=255)
                                                                                                                                                                                                   gTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGt                                           >  2:243282/1‑251 (MQ=255)
                                                                                                                                                                                                       tctAATTCCATTGACCAAT‑AGGCAGAGGCACCTGAAAACCAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTg                                          <  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
GAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAACCGTGGCTTTTGCAAT  >  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: ‑