Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,009,1490GA25.6% 65.6 / 20.6 39intergenic (+46/+36)norM/PP_5263multidrug resistance protein NorM/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base A (7/3);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGG  >  NC_002947/6009003‑6009284
                                                                                                                                                   |                                                                                                                                         
tGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTGGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGt                                                                                                                                         >  6:184598/1‑149 (MQ=255)
       tgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCtgtg                                                                                                                                  <  5:184598/149‑1 (MQ=255)
               tGGCGGTGACCTTTGGCTTTGAGTGTCGGATGAACCGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCggg                                                                                                                          <  3:60472/149‑1 (MQ=255)
                ggCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGc                                                                                                                         <  8:60911/149‑1 (MQ=255)
                    cTGACGTTTGGCTTTGAGTGGCGGATGATACGGTTGTTGGGTATGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACT‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGc                                                                                                                     <  6:83165/149‑1 (MQ=255)
                          tttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCGGGCGt                                                                                                                       <  4:276315/141‑1 (MQ=25)
                          tttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCGGGCGt                                                                                                                       >  3:276315/1‑141 (MQ=25)
                           ttGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCGCCCCCA‑GGATTGGGGGGGGCCGGGGGGGGCGGGGGTGG‑‑Ccgcg                                                                                                              >  6:203015/1‑147 (MQ=255)
                                tttGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGc                                                                                                       >  7:72502/1‑149 (MQ=255)
                                            aTGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCgg                                                                                                                           <  4:341811/119‑1 (MQ=21)
                                            aTGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCgg                                                                                                                           >  3:341811/1‑119 (MQ=21)
                                                aaCGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGCGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGTCACCCTc                                                                                       >  5:135515/1‑149 (MQ=255)
                                                    gttgttGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAAc                                                                                   <  8:195787/149‑1 (MQ=255)
                                                             tacggCGGGGGGTCGTACAGAAGCTGCGGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCAGGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCa                                                                          <  3:92043/146‑1 (MQ=38)
                                                                ggCGGGGGTTCGTACAGATGGTGCTGTTTCGGCCTGACTTCCGGGCGTCTTCGCGGGCACGCCCGCTCCAACA‑GTATTGGGCGGTCCCTGTGGGAGCGGGCGGGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGgc                                                                       <  6:257660/149‑1 (MQ=38)
                                                                   ggggggGCGTACTGATGCTGCTGTTTCGGCCTGAGTTCCGGGCCTCTTCGCGGGAACGCCCGCTCCCCCA‑GGATTGCGCGGTGCCTGTGGGAGGGGGGGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAAGTGGCGGCCAGGCGcc                                                                    <  2:12427/149‑1 (MQ=255)
                                                                   gggggTTCTTACAGAAGCTGCTGTTTCGCCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGcc                                                                    <  4:79298/149‑1 (MQ=255)
                                                                            tACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAgg                                                                                               <  2:215354/113‑1 (MQ=255)
                                                                            tACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAgg                                                                                               >  1:215354/1‑113 (MQ=255)
                                                                                       ctgTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCGACAACAGGCGGCCAGGCGCCAGCAGCAAATACGCCATCa                                                 >  4:259748/1‑148 (MQ=255)
                                                                                         gTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGGTTGCGCGGTCCCGGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCGCCAACTGGCGGCCAGGCGCCGGCAGCAAATACGCCCTCAAt                                               >  6:173269/1‑148 (MQ=38)
                                                                                          tttCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTgg                                                                                                                                 <  6:2982/67‑1 (MQ=255)
                                                                                          tttCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTgg                                                                                                                                 >  5:2982/1‑67 (MQ=255)
                                                                                                      ttCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGCGCGGGCGTGCTGGGGCGAAGGGGCCCTAAAGGAACCCCTCCAACTGGCGGCCCGGCGCCGGCCGCAAATACGCCATCAATTCCCCTATCGGt                                   >  1:290606/1‑149 (MQ=37)
                                                                                                      ttCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACCCCCTCCAACTGGCGGCCCGGCGCCTGCGGCAAATACGCCATCAATTCCCCGATCGGTAg                                 >  5:140233/1‑149 (MQ=255)
                                                                                                             ccTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTct                           >  6:95050/1‑148 (MQ=21)
                                                                                                                    gcgGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGAc                     >  8:201500/1‑147 (MQ=255)
                                                                                                                    gcgGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCa                   <  3:259748/149‑1 (MQ=255)
                                                                                                                            gCCCGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCt           >  3:238391/1‑149 (MQ=25)
                                                                                                                            gCCCGCTCCCACA‑GTATGCCGCGGTCCCTGTGGGAGCGGGCGGGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCt           <  7:287186/149‑1 (MQ=255)
                                                                                                                              ccGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCaa                                                <  3:184151/110‑1 (MQ=18)
                                                                                                                              ccGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCaa                                                >  4:184151/1‑110 (MQ=18)
                                                                                                                              ccGCTCCCACAGGGACCGCGCAAT‑CCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTg         >  6:234534/1‑149 (MQ=37)
                                                                                                                                gCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCa       <  2:290606/149‑1 (MQ=255)
                                                                                                                                  tCCCACAGGGAATGCGCAAT‑CCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCAcc     >  2:292254/1‑149 (MQ=38)
                                                                                                                                   cccACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCt    >  3:139177/1‑149 (MQ=255)
                                                                                                                                    ccACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCAcc     <  5:173269/147‑1 (MQ=255)
                                                                                                                                    ccACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATAAATTCCCCTATCGGTAGCGGTCTGTTGACCAGGAAGCCTTGCACCTg   >  6:19351/1‑149 (MQ=255)
                                                                                                                                     caca‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTgg  <  7:201500/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
TGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACA‑GGATTGCGCGGTCCCTGTGGGAGCGGGCGTGC‑‑CCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGG  >  NC_002947/6009003‑6009284

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑