Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779598,7990TG42.3% 21.3 / 12.9 26V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (3/12);  new base G (11/0);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 4.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATA  >  NC_007779/598709‑598899
                                                                                          |                                                                                                    
aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGAGCCTg                                                                                            >  1:1048168/1‑101 (MQ=255)
       cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCg                                                                                     <  3:403452/101‑1 (MQ=255)
               cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGTCCTGCGATCCGGCCAAAAc                                                                             >  3:381203/1‑101 (MQ=255)
                           ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGTCCTGCGGCCCGGCAAAAACGCCCGAGAAGTg                                                                 >  1:1549892/1‑101 (MQ=255)
                               cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGGGCTGGGGTGCGGGAAAAACGGCCGAGAAGTGaatg                                                             >  3:885051/1‑98 (MQ=255)
                                 cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                           <  1:305600/101‑1 (MQ=255)
                                        gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGa                                                    <  1:46150/101‑1 (MQ=255)
                                            cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                <  3:765476/101‑1 (MQ=255)
                                                     gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGa                                       <  1:197873/101‑1 (MQ=255)
                                                     gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGa                                       <  1:605248/101‑1 (MQ=255)
                                                      ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGaa                                      >  1:1313120/1‑101 (MQ=255)
                                                          aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGACAAACTGGAAACGc                                  >  1:1147376/1‑101 (MQ=255)
                                                            tAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTg                                >  3:1172552/1‑101 (MQ=255)
                                                               aCGGCGAAGGGGAAGTGGCGTTCGGGGTGTTGAGCTTGCTATCTGTCAAAAGCACCCGAGAAGTAATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                             <  1:959241/101‑1 (MQ=255)
                                                               aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCCAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                             >  3:301060/1‑101 (MQ=255)
                                                                 ggcgaaggcgaaGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAag                           <  3:486286/101‑1 (MQ=255)
                                                                 ggcgaaggcgaaGTGGCGGGCGGGGGGGTGAGCCCTTGATCTGGCAAAAAAACCCGGGAAGTGATAGCCGCCGGGAATGACAAACTGGAAACgatttacat                           >  3:1567610/1‑93 (MQ=255)
                                                                  gcgaaggcgaagTGGCGGGCGGGGGGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAagt                          >  3:926134/1‑101 (MQ=255)
                                                                        gcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                    <  1:597735/101‑1 (MQ=255)
                                                                         cgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGc                   >  1:730602/1‑101 (MQ=255)
                                                                             gTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGa               >  1:19623/1‑101 (MQ=255)
                                                                                  gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCg          >  3:1102169/1‑101 (MQ=255)
                                                                                    gcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg        <  1:226367/101‑1 (MQ=255)
                                                                                     cgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg       <  3:1546615/101‑1 (MQ=255)
                                                                                      ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGga      >  3:16157/1‑101 (MQ=255)
                                                                                       gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgag     <  3:732991/101‑1 (MQ=255)
                                                                                        ggtggtGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga    >  3:325310/1‑101 (MQ=255)
                                                                                          tggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa  >  1:338389/1‑101 (MQ=255)
                                                                                          tggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa  >  3:1407614/1‑101 (MQ=255)
                                                                                          |                                                                                                    
ATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATA  >  NC_007779/598709‑598899

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