Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273595,0350TG71.9% ‑0.1 / 23.2 33V289G (GTG→GGG) cusACu(+)/Ag(+) efflux RND transporter permease subunit CusA
Reads supporting (aligned to +/- strand):  ref base T (0/9);  new base G (23/0);  total (24/9)
Fisher's exact test for biased strand distribution p-value = 3.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGC  >  NZ_CP009273/594948‑595123
                                                                                       |                                                                                        
gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                        <  2:922030/90‑1 (MQ=255)
        gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCgg                                                                                >  1:977175/1‑89 (MQ=255)
             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTGCCGGGATCCgg                                                                           >  1:615053/1‑90 (MQ=255)
             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTcccg                                                                           >  2:1122375/1‑88 (MQ=255)
             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCgg                                                                           >  2:52624/1‑90 (MQ=255)
             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCgg                                                                           >  2:466803/1‑90 (MQ=255)
             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGGGCCgg                                                                           >  1:314581/1‑90 (MQ=255)
             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTCCgg                                                                           >  2:1018105/1‑90 (MQ=255)
             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATcccg                                                                           >  2:1065208/1‑88 (MQ=255)
               tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCa                                                                         >  2:13983/1‑90 (MQ=255)
                ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCaa                                                                        >  2:305488/1‑90 (MQ=255)
                  aGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaa                                                                      <  2:166867/90‑1 (MQ=255)
                    aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc                                                                    <  2:237569/90‑1 (MQ=255)
                       ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCCGCAAAAAAGcc                                                                 >  2:1011688/1‑90 (MQ=255)
                       ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTCCGGTCCGGCAAAAACAcc                                                                 >  2:803162/1‑90 (MQ=255)
                       ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCAAAAACGcc                                                                 >  2:230880/1‑90 (MQ=255)
                        gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAACGccc                                                                >  1:402653/1‑90 (MQ=255)
                        gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGccc                                                                >  2:409820/1‑90 (MQ=255)
                                  gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGAt                                                      >  2:54475/1‑90 (MQ=255)
                                  gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGAt                                                      >  2:45268/1‑90 (MQ=255)
                                        gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                <  2:1166220/90‑1 (MQ=255)
                                                 cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                       <  2:977175/90‑1 (MQ=255)
                                                  gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCAGGGAAAAAGCCCCAGCAGTGGGCGCCGCCGTGGAggtc                                      >  1:809118/1‑88 (MQ=255)
                                                  gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGGAAAAAAGCCCGAGAAGGGATCGCCGCCGCGAAGGAc                                      >  2:987210/1‑90 (MQ=255)
                                                  gAACTAAACAGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATCGCCCCCGGGGAgggc                                      >  2:862589/1‑88 (MQ=255)
                                                       aaaCGGCGAAGGCGAAGTGGCGGGCGGGGGGGAGATCCCGCGATCCGGCCAAAAACCACGAGAAGTGATCGCCGCCGCGAAGGACAAACt                                 >  2:134833/1‑90 (MQ=255)
                                                                 ggcgaagTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAACGCCCCAGAAGGGATCGCCGCCGGGAAGGGAAAACTGGAAAAGCTg                       >  1:445691/1‑90 (MQ=255)
                                                                     aagTGGCGGGCGGGGGGGGGATCCTTGGATCCCGCGAAAACGCCCTTGGATTGGGCGCCGGCGGGGAGGGCGAAATGGGAAAGCTGaaaa                   >  1:431934/1‑90 (MQ=255)
                                                                           cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc             <  2:402653/90‑1 (MQ=255)
                                                                           cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc             <  2:255095/90‑1 (MQ=255)
                                                                           cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc             <  2:724076/90‑1 (MQ=255)
                                                                           cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc             <  1:999249/90‑1 (MQ=255)
                                                                            gggcgggGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCCAGAAGGGATCGCCGCCGGGAAGGGCAAAATGGAAAAGCTGAAAAGTAGTCt            >  2:1202431/1‑90 (MQ=255)
                                                                                  ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa      <  2:809118/90‑1 (MQ=255)
                                                                                   gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa     <  2:384188/90‑1 (MQ=255)
                                                                                   gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa     <  2:144069/90‑1 (MQ=255)
                                                                                      gtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc  <  1:276395/90‑1 (MQ=255)
                                                                                       |                                                                                        
GCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGC  >  NZ_CP009273/594948‑595123

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