Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,292,4750TG52.0% 64.6 / 70.6 50T141P (ACC→CCC) mdtOmultidrug efflux transporter permease subunit MdtO
Reads supporting (aligned to +/- strand):  ref base T (4/20);  new base G (26/0);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 2.25e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-09
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATC  >  NZ_CP009273/4292388‑4292563
                                                                                       |                                                                                        
tcaGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCTAGCATGGCGGGGAAGGttt                                                                                        <  1:866607/90‑1 (MQ=255)
           ggATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGTCCCTaaaaa                                                                             >  1:238186/1‑88 (MQ=255)
                gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGCAAAAAGCGAc                                                                        >  1:528138/1‑90 (MQ=255)
                gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCCAAAAAGCGAc                                                                        >  1:473789/1‑90 (MQ=255)
                gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGAc                                                                        >  2:383501/1‑90 (MQ=255)
                  ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAAAAACGACaa                                                                      >  1:878890/1‑90 (MQ=255)
                  ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCGAAAAAACGACaa                                                                      >  2:287707/1‑90 (MQ=255)
                  ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACaa                                                                      >  2:351622/1‑90 (MQ=255)
                  ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAGCGACaa                                                                      >  2:156117/1‑90 (MQ=255)
                       acaacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAGCGACAAAgggg                                                                 >  2:143536/1‑88 (MQ=255)
                         aacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGAc                                                               <  2:138088/90‑1 (MQ=255)
                            gATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGTTTGCCCGTAAATAACGACAATggcggcggc                                                            >  2:678198/1‑90 (MQ=255)
                                  ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCGAAAAAACGACAAAGGCGGCGGCGgaaaa                                                      >  1:315189/1‑90 (MQ=255)
                                  ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTGCCGTAAAAAAGGGCAATGGCGGCGGGGgaaaa                                                      >  2:520507/1‑90 (MQ=255)
                                  ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAAAGGCGGCGGCGgaaaa                                                      >  2:468490/1‑90 (MQ=255)
                                    acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAATAGCGACAAAGGCGGCGGCGgaaaaaa                                                    >  2:407498/1‑88 (MQ=255)
                                    acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCCTAAAAAGCGACAAAGGCGACGGCGaaaaaaa                                                    >  2:809553/1‑88 (MQ=255)
                                     caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAACGACAAAGGCGACGGGGGAAAAAAc                                                   >  1:210327/1‑90 (MQ=255)
                                     caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCCAAAAAACGACAAAAGGGGCGGCGgaaaaaaa                                                   >  1:902913/1‑87 (MQ=255)
                                      aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc                                                  <  1:678198/90‑1 (MQ=255)
                                      aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc                                                  <  2:156721/90‑1 (MQ=255)
                                            aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGcccc                                            <  2:939718/90‑1 (MQ=255)
                                            aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCCGTAATAGGAACAAATGGGAGGGGGGAAAAGACCAAcccc                                            >  1:751938/1‑90 (MQ=255)
                                            aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAAACCAGcccc                                            >  2:103933/1‑90 (MQ=255)
                                                 cACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCGTAAAAAGCGACAAAGGGGGCGGCGGAAAAAACCAGCCCCCAGCg                                       >  1:67853/1‑90 (MQ=255)
                                                 cACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAGCGACAATGGCGACGGCGGAAAAGACCAGCCCCAAGCg                                       >  2:125912/1‑90 (MQ=255)
                                                   cGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGAt                                     <  1:577394/90‑1 (MQ=255)
                                                         cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc                               <  2:7231/90‑1 (MQ=255)
                                                         cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc                               <  2:504349/90‑1 (MQ=255)
                                                         cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc                               <  2:205022/90‑1 (MQ=255)
                                                         cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc                               <  1:383501/90‑1 (MQ=255)
                                                         cTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgc                               <  1:905992/90‑1 (MQ=255)
                                                          ttCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgcg                              <  2:645352/90‑1 (MQ=255)
                                                                  aGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCaaaa                      <  2:473789/90‑1 (MQ=255)
                                                                   gTCGAGCATGGCGGGGGAGGGTTTCCCGTAAAAAACGACAAAGGCGGCGGGGGAAAAAACCAGCCCCCAGCGATGGGGGGGCCTCaaaaa                     >  1:205022/1‑90 (MQ=255)
                                                                   gTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAATGGGGACGGCGAAAAAAACCAGCCCCCAACGATGGGGGCGCATCaaaaa                     >  2:864646/1‑90 (MQ=255)
                                                                   gTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAAAGGGGACGGCGAAAAAAACCAGCCCCCAGCGGAGGGGGGGCATCaaaaa                     >  1:346248/1‑90 (MQ=255)
                                                                   gTCGAGCATGGCGGGGAAGGTTTGCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAAAACAGCCCCAAACGATGGGGGGGCATCaaaaa                     >  2:94972/1‑90 (MQ=255)
                                                                   gTCGAGCATGGCGGGGAAGGGTTTCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCCAGCGATGGGGGCGCATCaaaaa                     >  2:682631/1‑90 (MQ=255)
                                                                   gTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCCAGCGAGGGGGGCGCATCaaaaa                     >  2:531961/1‑90 (MQ=255)
                                                                   gTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCAAACGATGGGGGGGGATCaaaaa                     >  2:829893/1‑90 (MQ=255)
                                                                      gAGCATGGCGGGGAAGGGTTGCCCGTAAATAACGACAATGGGGACGGCGAAAAAAAACAGCCCCCAACCATGGGGGCGGAGCaaaaaaaa                  >  2:540617/1‑87 (MQ=255)
                                                                        gCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGc                >  2:395282/1‑90 (MQ=255)
                                                                           tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc             <  1:540617/90‑1 (MQ=255)
                                                                           tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc             <  2:528138/90‑1 (MQ=255)
                                                                           tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc             <  1:459600/90‑1 (MQ=255)
                                                                                      gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc  <  2:232180/90‑1 (MQ=255)
                                                                                      gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc  <  2:552483/90‑1 (MQ=255)
                                                                                      gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc  <  1:463070/90‑1 (MQ=255)
                                                                                      gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc  <  1:287707/90‑1 (MQ=255)
                                                                                      gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc  <  1:130167/90‑1 (MQ=255)
                                                                                       |                                                                                        
TCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATC  >  NZ_CP009273/4292388‑4292563

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: