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

ACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  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)
                        gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACg                                                                                                                             <  2:550793/151‑1 (MQ=255)
                            tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGccc                                                                                                                          <  2:122283/150‑1 (MQ=255)
                             aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGCCTGC‑GGTCCGGCaaaaa                                                                                                                               >  1:258710/1‑144 (MQ=255)
                                      tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GATCCGGCAAAAAAGCCCGAGAAGTg                                                                                                                  >  2:544895/1‑148 (MQ=255)
                                       gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATcg                                                                                                              <  2:51389/151‑1 (MQ=255)
                                         aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGGGATcgcc                                                                                                            >  2:651681/1‑151 (MQ=255)
                                           aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATcgccgc                                                                                                          <  2:258766/151‑1 (MQ=255)
                                             ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                         <  2:264129/150‑1 (MQ=255)
                                               cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGC‑GATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTg                                                                                                      >  2:701848/1‑151 (MQ=255)
                                                    ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                                  <  2:84714/150‑1 (MQ=255)
                                                     cGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAACGCCCGCGAAGGGATCGCCGCCGTGGAgg                                                                                                  >  1:242935/1‑149 (MQ=255)
                                                          aTCTGCGCGATGCTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAAc                                                                                            <  2:496239/150‑1 (MQ=255)
                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GATCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTGGAGGACAAACTgg                                                                                         >  2:194586/1‑149 (MQ=255)
                                                               cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaa                                                                                       <  2:679870/150‑1 (MQ=255)
                                                                gcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAACGCCCGAGGAGGGATCGCCGCCGCGAAGGACAAACTGGaaa                                                                                      >  1:336238/1‑150 (MQ=255)
                                                                  gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGaaa                                                                                      >  2:41838/1‑148 (MQ=255)
                                                                        gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaaa                                                                             <  2:577282/151‑1 (MQ=255)
                                                                          gacggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GATTCGGCAAAAAAACGCGGGGAGTGGGCGGCGCCCGGAAGGACAAACTGGGAAAACAGaaaa                                                                             >  1:72714/4‑149 (MQ=255)
                                                                          gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGACCTAAACGGCGAAGGGGAAGTGGAGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                                                                           <  2:35309/151‑1 (MQ=255)
                                                                              gTCCAGATTGGCCCGGAGATGCGCCGAGGTATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGGGATCGCCGGGGTGAAGGAAAAAATGGaaaa                                                                                     >  2:709214/1‑136 (MQ=255)
                                                                                ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTTCTGC‑GGTCCCGCCAAAACGGCCGAGAAGTGATCGCCGCCGTGAAGGACAAACATGCAACGATGAAAAGGAGTCTg                                                                     >  1:526955/1‑151 (MQ=255)
                                                                                            ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                         <  2:458058/151‑1 (MQ=255)
                                                                                               gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                       <  1:550824/150‑1 (MQ=255)
                                                                                                aTGCGCCGGGGCATTGCCGAACTAAGCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GATCCCGCAAAAAAGACCGAGAAGGGATCGCCGCCGGGAAGGACAAACTGGAAACGCTGAAAAGAAGTCTGCCGGAAGGCGTgg                                                        >  1:110621/1‑148 (MQ=255)
                                                                                                        gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTaca                                             <  2:418039/151‑1 (MQ=255)
                                                                                                            aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGC‑GACCCGGCAAAAACGCCCGGGGAGGGATCGCCGCCGCGAAGGACAAACTGGAAAAGCTGAAAAGGAGTCTGCCGGAAGGCGGGGGGGTGGTTacaa                                            >  1:280034/1‑148 (MQ=255)
                                                                                                             ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCCCCGGACCCGGC‑CAAACGCCCGAGGAGCGCCAGCCGCCGTGCAGGGCCAACCAGGAACGCTGAAAAGTAGGCTGCCGGACGGAGGGGAGATAGGTccaaca                                          >  1:680729/1‑149 (MQ=255)
                                                                                                              tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GGACCGGCAAAAACGCCCGGGAAGGGAGCGCCGCCGTGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGCCGGGGGGGGGGGAGATAGTTataaca                                          >  1:48883/1‑148 (MQ=255)
                                                                                                              tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGC‑GATCCGGCCAAAACGCCCGAGAAGGGATAGCCGCCCTGAAGGGCAAAATGGAAACCCTGAAAAGAAGTCTGCCGGAAGGGGGGGAGGTAGTTACAACATa                                        >  1:572971/1‑150 (MQ=255)
                                                                                                                 cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGC‑GATCCGGCCAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACCAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGa                                     >  2:443443/1‑150 (MQ=255)
                                                                                                                 cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGGAGGGATCGGCGCCGTGAAGGGCAAACTGGAAAACCTGGAAAGTAGGCTGCCGGGAGGGGGGGGGATAGTTACAACATACGa                                     >  1:213645/1‑150 (MQ=255)
                                                                                                                                       gTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgcg              <  2:334623/151‑1 (MQ=255)
                                                                                                                                         ggcgggcgGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCgcgc             <  1:275752/150‑1 (MQ=255)
                                                                                                                                           cgggcgggGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCAt          >  2:513969/1‑151 (MQ=255)
                                                                                                                                            gggcgggGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc         <  2:335346/151‑1 (MQ=255)
                                                                                                                                                 gggTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACaa    <  1:564611/151‑1 (MQ=255)
                                                                                                                                                 gggTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa     >  1:571596/1‑150 (MQ=255)
                                                                                                                                                  ggtggtGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAc   >  2:535749/1‑151 (MQ=255)
                                                                                                                                                   gtggtgATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  <  1:598587/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGC‑GATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  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: ‑