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

GCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAA  >  NZ_CP009273/594946‑595120
                                                                                         |                                                                                     
gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt                                                                                       <  2:1740075/90‑1 (MQ=255)
  gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGa                                                                                     >  1:251390/1‑90 (MQ=255)
               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGATCCgg                                                                        >  2:667185/1‑90 (MQ=255)
               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATTCTGCGGACCgg                                                                        >  1:50315/1‑90 (MQ=255)
               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCgg                                                                        >  1:1716651/1‑90 (MQ=255)
               ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCgg                                                                        >  1:556517/1‑90 (MQ=255)
                      aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc                                                                 <  1:1649618/90‑1 (MQ=255)
                         ggCCCGGAGATGCGCCGGGGGATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGcc                                                              <  1:1165810/90‑1 (MQ=255)
                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGGTCCCGCAAAAAAGcc                                                              >  2:67560/1‑90 (MQ=255)
                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAACGcc                                                              >  1:479225/1‑90 (MQ=255)
                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGcc                                                              >  1:620502/1‑90 (MQ=255)
                          gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGccc                                                             >  1:912612/1‑90 (MQ=255)
                           cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCg                                                            <  1:1769650/90‑1 (MQ=255)
                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGGTCCGGCAAAAACGCCCGCGAAGGGag                                                   >  2:742695/1‑89 (MQ=255)
                                    gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCCGCAAAAACGCCCGAGAAGGGaa                                                   >  1:1433722/1‑89 (MQ=255)
                                       cgggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc                                                <  2:1671808/88‑1 (MQ=255)
                                       ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc                                                <  2:251390/90‑1 (MQ=255)
                                         ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGACCCCGGAAAAACACCCCAGAAGGGATcgccg                                              >  2:355831/1‑90 (MQ=255)
                                                       cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                <  2:998322/90‑1 (MQ=255)
                                                           aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTgg                            >  1:724176/1‑90 (MQ=255)
                                                              gcgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                         <  1:1619842/90‑1 (MQ=255)
                                                                      gaagTGGCGGGGGGGGCGGTGGTCCTGGGATCCGGGAAAACCGCCCGAGAAGGGATCGCCGCCGTGAAGGGCAAACGGGAAACGCTGaaa                 <  2:1440104/90‑1 (MQ=255)
                                                                      gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                 <  1:937937/90‑1 (MQ=255)
                                                                                    ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa   <  1:1444291/90‑1 (MQ=255)
                                                                                     gtggtgACCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGaa  <  2:408218/90‑1 (MQ=255)
                                                                                         |                                                                                     
GCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAA  >  NZ_CP009273/594946‑595120

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