Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273595,0350TG57.9% 3.3 / 17.0 19V289G (GTG→GGG) cusACu(+)/Ag(+) efflux RND transporter permease subunit CusA
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base G (11/0);  total (15/4)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGG  >  NZ_CP009273/594946‑595122
                                                                                         |                                                                                        
gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt                                                                                          <  1:625628/90‑1 (MQ=255)
  gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                        <  2:391247/90‑1 (MQ=255)
               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATcccg                                                                           >  2:386988/1‑88 (MQ=255)
                               gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGGTCCGGCCAAAAACCCCCggaa                                                           >  1:170702/1‑87 (MQ=255)
                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAggggt                                                     >  2:670965/1‑86 (MQ=255)
                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCCGCGGTCCGGCAAAAAAGCCCGAGAAGGGAt                                                     >  2:604186/1‑90 (MQ=255)
                                         ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGGTcgccg                                                >  1:198378/1‑90 (MQ=255)
                                         ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGGCCCGGCAAAAAAGCCCGGGAAGGGATggccc                                                >  1:346953/1‑89 (MQ=255)
                                         ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGCTCCGGCAAAAAAACCCGAGAAGTGATcgccc                                                >  2:353687/1‑89 (MQ=255)
                                              aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCCCGATCCCGCCAAAAAACCCGAGAAGGGAGCGccgccgcg                                           >  1:375074/1‑88 (MQ=255)
                                               ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGa                                          >  2:323644/1‑90 (MQ=255)
                                                    gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc                                     <  1:336034/90‑1 (MQ=255)
                                                    gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCCGGGATCCGGCAAAAAAACCCCAGGAGTGAACGCCCCCCGGAAgggc                                     >  2:628947/1‑88 (MQ=255)
                                                     aaCTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGGGGAgggaa                                    >  1:322938/1‑87 (MQ=255)
                                                               cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                          <  1:353687/90‑1 (MQ=255)
                                                                        agTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGAAGGAGAAAATGGGAACGCTGaaaaa                 >  2:651372/1‑89 (MQ=255)
                                                                             cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                          <  1:783884/60‑1 (MQ=255)
                                                                             cgggcgggGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                          >  2:783884/1‑60 (MQ=255)
                                                                                   gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg      >  1:786954/1‑90 (MQ=255)
                                                                                   gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg      >  1:425281/1‑90 (MQ=255)
                                                                                    ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa     <  2:175017/90‑1 (MQ=255)
                                                                                    ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa     <  2:330823/90‑1 (MQ=255)
                                                                                    ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa     <  2:697122/90‑1 (MQ=255)
                                                                                      tggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg   <  2:623082/90‑1 (MQ=255)
                                                                                       ggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg  >  2:485729/1‑90 (MQ=255)
                                                                                         |                                                                                        
GCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGG  >  NZ_CP009273/594946‑595122

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