Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273595,0320TG46.7% 36.6 / 57.9 60V288G (GTG→GGG) cusACu(+)/Ag(+) efflux RND transporter permease subunit CusA
Reads supporting (aligned to +/- strand):  ref base T (1/31);  new base G (28/0);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 2.80e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAG  >  NZ_CP009273/594943‑595121
                                                                                         |                                                                                          
tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                            <  1:2556948/90‑1 (MQ=255)
tCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                            <  1:180200/90‑1 (MQ=255)
  tGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGgggggggggg                                                                                          >  1:494831/1‑90 (MQ=255)
   gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggg                                                                                         >  2:2135506/1‑89 (MQ=255)
   gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggg                                                                                         >  2:16805/1‑89 (MQ=255)
     gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggggg                                                                                       >  2:1469158/1‑87 (MQ=255)
     gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                       <  1:872233/90‑1 (MQ=255)
                  ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTTCTGCGGTCCgg                                                                         >  2:1574840/1‑90 (MQ=255)
                  ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGTCTGCGATcccg                                                                         >  1:732731/1‑88 (MQ=255)
                  ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCgg                                                                         >  2:742737/1‑90 (MQ=255)
                  ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTACCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTcccg                                                                         >  2:147625/1‑88 (MQ=255)
                        gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGAGCCGGCgaaaa                                                                   >  1:1535955/1‑90 (MQ=255)
                         aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc                                                                  <  2:1290976/90‑1 (MQ=255)
                         aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAAc                                                                  <  2:1382707/90‑1 (MQ=255)
                            ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGATGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGcc                                                               >  2:829053/1‑90 (MQ=255)
                            ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCcaaaaaacc                                                               >  1:2843556/1‑86 (MQ=255)
                            ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCaaaaaaacc                                                               >  2:2365925/1‑86 (MQ=255)
                            ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCCGCcaaaaaacc                                                               >  1:320560/1‑86 (MQ=255)
                            ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                            <  2:1279165/62‑1 (MQ=255)
                             gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCCCGCAAAAAAGccc                                                              >  1:548901/1‑90 (MQ=255)
                             gCCCGGAGAGGCGCCGGGGAATTCCCGACCAAACCGGCAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCg                                                             <  2:769818/90‑1 (MQ=255)
                                 ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgaga                                                          <  2:743809/90‑1 (MQ=255)
                                 ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCgaga                                                          >  1:996875/1‑90 (MQ=255)
                                  gagaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGGTCCGGCAAAAAAGCCCGGGaa                                                         >  1:2584068/1‑90 (MQ=255)
                                    gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTCCGGTCCCGCAAAAAAGCCCGCGAAgg                                                       >  1:1960067/1‑89 (MQ=255)
                                       gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCCAAAAAGCCCGAGAAGGGAt                                                    >  1:492239/1‑90 (MQ=255)
                                       gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGGACCCGCGAAAAAGCCCGAGAAggggt                                                    >  1:388171/1‑86 (MQ=255)
                                        cgcgGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                   <  1:2125155/90‑1 (MQ=255)
                                          ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc                                                 <  1:2263712/90‑1 (MQ=255)
                                            ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGGTCCGGCAAAAACGCCCCAGAAGGGATcgccg                                               >  2:2943864/1‑90 (MQ=255)
                                            ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCCGCGAAAAAGCCCGAGAAGGGGTcgccg                                               >  1:2177439/1‑90 (MQ=255)
                                            ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCCGGCAAAAAAGCCCGAGAAGGGGTcgccg                                               >  1:2928349/1‑90 (MQ=255)
                                             gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                              <  1:1704054/90‑1 (MQ=255)
                                             gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                              <  1:1574840/90‑1 (MQ=255)
                                             gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                              <  2:2643332/90‑1 (MQ=255)
                                             gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                              <  1:16805/90‑1 (MQ=255)
                                             gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                              <  2:705090/90‑1 (MQ=255)
                                             gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                              <  1:829053/90‑1 (MQ=255)
                                                 aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGGTCCGGCGAAAAAACCCCCGAAGGGGGCgccgcccgg                                          >  1:2643332/1‑87 (MQ=255)
                                                 aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTg                                          >  2:149215/1‑90 (MQ=255)
                                                  ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGGCAAAAAAGCCCGGGAAGGGAGCGCCGCCGTGa                                         >  2:1796133/1‑90 (MQ=255)
                                                  ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCCAAAAACCCCCGGAAGGGATCgccgcccgga                                         >  2:1704054/1‑86 (MQ=255)
                                                   tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                        <  2:1621844/90‑1 (MQ=255)
                                                    gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGGGAAg                                       >  2:1342568/1‑90 (MQ=255)
                                                       gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGAc                                    <  1:2041737/90‑1 (MQ=255)
                                                       gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGAc                                    >  2:46709/1‑90 (MQ=255)
                                                          cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                 <  2:492239/90‑1 (MQ=255)
                                                                 gcgaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                          >  2:1113279/1‑90 (MQ=255)
                                                                  cgaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAc                         <  2:641903/90‑1 (MQ=255)
                                                                    aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc                       <  2:2630518/90‑1 (MQ=255)
                                                                    aaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGc                       <  2:2843556/90‑1 (MQ=255)
                                                                      ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                     <  1:621511/90‑1 (MQ=255)
                                                                      ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                     <  1:3002323/90‑1 (MQ=255)
                                                                      ggcgaagTGGCGGGCGGGGGGGGGGTCCTGCGAACCGGCAAAAAAGCCCGAGAAGGGGTCGCCGCCGGGGAGGACAAACAGGAAAAGCTg                     >  1:2167176/1‑90 (MQ=255)
                                                                         gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                  <  2:294053/90‑1 (MQ=255)
                                                                           agTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                <  2:940108/90‑1 (MQ=255)
                                                                              ggcgggcgGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag             <  2:392965/90‑1 (MQ=255)
                                                                              ggcgggcgGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag             <  2:1879665/90‑1 (MQ=255)
                                                                              ggcgggcgGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagtag             <  1:1979003/90‑1 (MQ=255)
                                                                                cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc           <  2:2504691/90‑1 (MQ=255)
                                                                                cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc           <  2:2727975/90‑1 (MQ=255)
                                                                                cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc           <  2:1920671/90‑1 (MQ=255)
                                                                                cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTc           <  1:588393/90‑1 (MQ=255)
                                                                                cgggcgggGGGGGGATCCTGCGATCCGGCCAAAAAGCCCGAGGAGGGATCGCCGCCGTGGAGGACAAAATGGGAAAGCTGAAAAAAAGTc           >  1:1869412/1‑90 (MQ=255)
                                                                                   gcgggGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGGCAAAAAGGAAAAGCTGAAAAAAAGTCTGc        >  2:2572080/1‑90 (MQ=255)
                                                                                     ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCg      <  1:2406503/90‑1 (MQ=255)
                                                                                      gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg     >  1:1648714/1‑90 (MQ=255)
                                                                                       ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa    <  2:732731/90‑1 (MQ=255)
                                                                                       ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa    <  1:411911/90‑1 (MQ=255)
                                                                                       ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa    <  2:1501394/90‑1 (MQ=255)
                                                                                       ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa    <  2:1960067/90‑1 (MQ=255)
                                                                                       ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa    <  2:22736/90‑1 (MQ=255)
                                                                                         tggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg  <  2:134504/90‑1 (MQ=255)
                                                                                         tggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg  <  2:500399/90‑1 (MQ=255)
                                                                                         |                                                                                          
TCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAG  >  NZ_CP009273/594943‑595121

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