Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG50.0% 44.4 / 17.9 40V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (2/18);  new base G (20/0);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 3.35e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599428‑599725
                                                                                                                                                    |                                                                                                                                                      
aCGCACGAAGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGCGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTTGCTgggggggggg                                                                                                                                                      >  1:566239/1‑151 (MQ=255)
 cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgg                                                                                                                                                      >  2:359976/1‑150 (MQ=255)
      gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATcc                                                                                                                                                >  2:215845/1‑151 (MQ=255)
       acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATcc                                                                                                                                                >  1:40953/1‑150 (MQ=255)
        cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTg                                                                                                                                              <  1:5191/151‑1 (MQ=255)
               aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggggg                                                                                                                                                   >  1:330265/1‑136 (MQ=255)
                        gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACg                                                                                                                             <  2:550793/151‑1 (MQ=255)
                            tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGccc                                                                                                                          <  2:122283/150‑1 (MQ=255)
                             aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGCGGTCCGGCaaaaa                                                                                                                               >  1:258710/1‑144 (MQ=255)
                                      tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGTg                                                                                                                  >  2:544895/1‑148 (MQ=255)
                                       gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcg                                                                                                              <  2:51389/151‑1 (MQ=255)
                                         aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATcgcc                                                                                                            >  2:651681/1‑151 (MQ=255)
                                           aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                                                                          <  2:258766/151‑1 (MQ=255)
                                             ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                         <  2:264129/150‑1 (MQ=255)
                                               cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTg                                                                                                      >  2:701848/1‑151 (MQ=255)
                                                    ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                                  <  2:84714/150‑1 (MQ=255)
                                                     cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGCGAAGGGATCGCCGCCGTGGAgg                                                                                                  >  1:242935/1‑149 (MQ=255)
                                                          aTCTGCGCGATGCTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                            <  2:496239/150‑1 (MQ=255)
                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTGGAGGACAAACTgg                                                                                         >  2:194586/1‑149 (MQ=255)
                                                               cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                       <  2:679870/150‑1 (MQ=255)
                                                                gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGCGAAGGACAAACTGGaaa                                                                                      >  1:336238/1‑150 (MQ=255)
                                                                  gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGaaa                                                                                      >  2:41838/1‑148 (MQ=255)
                                                                        gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                             <  2:577282/151‑1 (MQ=255)
                                                                          gacggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATTCGGCAAAAAAACGCGGGGAGTGGGCGGCGCCCGGAAGGACAAACTGGGAAAACAGaaaa                                                                             >  1:72714/4‑149 (MQ=255)
                                                                          gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGACCTAAACGGCGAAGGGGAAGTGGAGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                           <  2:35309/151‑1 (MQ=255)
                                                                              gTCCAGATTGGCCCGGAGATGCGCCGAGGTATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGGGGTGAAGGAAAAAATGGaaaa                                                                                     >  2:709214/1‑136 (MQ=255)
                                                                                ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTTCTGCGGTCCCGCCAAAACGGCCGAGAAGTGATCGCCGCCGTGAAGGACAAACATGCAACGATGAAAAGGAGTCTg                                                                     >  1:526955/1‑151 (MQ=255)
                                                                                            ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                         <  2:458058/151‑1 (MQ=255)
                                                                                               gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                       <  1:550824/150‑1 (MQ=255)
                                                                                                aTGCGCCGGGGCATTGCCGAACTAAGCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCCCGCAAAAAAGACCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGAAGTCTGCCGGAAGGCGTgg                                                        >  1:110621/1‑148 (MQ=255)
                                                                                                        gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTaca                                             <  2:418039/151‑1 (MQ=255)
                                                                                                            aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGACCCGGCAAAAACGCCCGGGGAGGGATCGCCGCCGCGAAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGCGGGGGGGTGGTTacaa                                            >  1:280034/1‑148 (MQ=255)
                                                                                                             ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCCCGGACCCGGCCAAACGCCCGAGGAGCGCCAGCCGCCGTGCAGGGCCAACCAGGAACGCTGAAAAGTAGGCTGCCGGACGGAGGGGAGATAGGTccaaca                                          >  1:680729/1‑149 (MQ=255)
                                                                                                              tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGACCGGCAAAAACGCCCGGGAAGGGAGCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGGGGGGGGGGAGATAGTTataaca                                          >  1:48883/1‑148 (MQ=255)
                                                                                                              tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCCAAAACGCCCGAGAAGGGATAGCCGCCCTGAAGGGCAAAATGGAAACCCTGAAAAGAAGTCTGCCGGAAGGGGGGGAGGTAGTTACAACATa                                        >  1:572971/1‑150 (MQ=255)
                                                                                                                 cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCCAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACCAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGa                                     >  2:443443/1‑150 (MQ=255)
                                                                                                                 cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGGAGGGATCGGCGCCGTGAAGGGCAAACTGGAAAACCTGGAAAGTAGGCTGCCGGGAGGGGGGGGGATAGTTACAACATACGa                                     >  1:213645/1‑150 (MQ=255)
                                                                                                                                       gTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgcg              <  2:334623/151‑1 (MQ=255)
                                                                                                                                         ggcgggcgGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc             <  1:275752/150‑1 (MQ=255)
                                                                                                                                           cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt          >  2:513969/1‑151 (MQ=255)
                                                                                                                                            gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc         <  2:335346/151‑1 (MQ=255)
                                                                                                                                                 gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa    <  1:564611/151‑1 (MQ=255)
                                                                                                                                                 gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa     >  1:571596/1‑150 (MQ=255)
                                                                                                                                                  ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAc   >  2:535749/1‑151 (MQ=255)
                                                                                                                                                   gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  <  1:598587/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599428‑599725

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: