Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,7910AG56.2% 3.3 / 16.1 17S182P (TCC→CCC) betIcholine‑inducible betIBA‑betT divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (1/6);  new base G (9/0);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-03

CCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑TCAGGGAATTAGC‑GCGGGTTTTA‑TCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCG‑CGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGC‑GGCAATTCGCGACGAAACTCGCTCACCAGATTCGAC  >  NC_000913/328646‑328933
                                                                                                                                                   |                                                                                                                                                  
ccGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑TCAGGGAATTAg                                                                                                                                               <  1:592774/151‑1 (MQ=255)
   ggTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑TCAGGGAATTAgc‑g                                                                                                                                            <  1:14950/150‑1 (MQ=255)
   ggTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑GCTGGGGAGTAggtg                                                                                                                                            >  1:135298/1‑148 (MQ=255)
   ggTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGGGATAAAGTGGCGGG‑‑ACGGGGGAtaattgg                                                                                                                                            >  2:267258/1‑145 (MQ=255)
             ctcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑TCAGGGAATTAGC‑GCGGGTTTTA‑Tc                                                                                                                                <  1:51369/151‑1 (MQ=255)
               cGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCTGGGGGTATTAGC‑GGGGGttgtcta                                                                                                                                 >  1:537789/1‑146 (MQ=255)
                                        ggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑TCAGGGGAGTTAC‑GGGGGGTTTA‑TCTAAGGGTTTGTGCCTCTGCGCcgcgc                                                                                                      >  1:509311/1‑150 (MQ=255)
                                        ggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑GCTGGGGATGAAA‑GGGGGGGGTT‑TCTAGCGGCTTGTCGCTCCGAGCcgcgcg                                                                                                     >  1:506455/1‑151 (MQ=255)
                                                        atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑GCGGGGGAATAAA‑GGGGGGTTTT‑TTCAACGGGTTGGCTGGCAGGGGCGCGCGCAGCCaaaaaacaaca                                                                                     >  2:105354/1‑144 (MQ=255)
                                                             aaGCGGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑TCAGGGAATTAGC‑GCGGGTTTTA‑TCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAg                                                                                <  1:267998/151‑1 (MQ=255)
                                                                                     cATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑TCAGGGGAATAGCGGGGGGTTTTT‑TCCAAGGGGTT‑TCGCCCAGGGGCGCGCGCAGGCCGAAACCCTCAAACATCGCCGGCCGGGcccggcccccct                                                         >  2:158548/1‑140 (MQ=255)
                                                                                        cACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑GCCGGGGACTAGC‑GCGGGGTTTT‑TTCACCGGGGTGGCGCCCCGCGCCCCCCGCCGCCCTAAACCATCAAACAACC‑cggcccggcccgcgcccccctccc                                                     >  1:603913/1‑136 (MQ=255)
                                                                                            tctcCGCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGCGGG‑‑TCGGGGAATTAGC‑GCGGGTTTTA‑TCCAGCGGTTGGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCG‑CGGCCAGGCCGTAGCCCGCTTCCTGTGcc                                                <  1:521615/151‑1 (MQ=255)
                                                                                                    gTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑GCCGGGGAGTAGA‑GGGGGGTTTT‑TCCAGCGGGTTGGCGCCCCGAGCCGCGCGCCGCCCCAATCCCTAAATCACCC‑CGGCCCGGCCGGAGCCCCCCCCCCCTCCCTGTTcgcg                                        >  2:410470/1‑151 (MQ=255)
                                                                                                               ggTAGATGCTGAGTGATAAAGTGGCGGG‑‑GCAGGGAAATAAC‑GGGGGGGTTATTCCACCGG‑TTGCGGCTCACAGCCGCGCGCAGCCCTAATCCATAAAACAACG‑CGGGCCCGCCGGCGACCGCCCCCTCCGCCTCTTTCCGC‑GGGAATTct                             >  2:432212/1‑149 (MQ=255)
                                                                                                                gTAGATGCTGAGTGATAAAGTGGCGGG‑‑TGAGGGGGTTAGA‑GCGGGGTTTT‑TTCACCGGGGTGTCGCCCAGGGCCGCGGGCCGCCATAATCCATCTATCATCG‑CCGCCCGGCCGCAGACCCCCTCCCGTTCCCCTTCGCCCGGGCAAttgc                             >  2:496242/1‑148 (MQ=255)
                                                                                                                                          g‑‑TCAGGGAATTGGC‑GCGGGTTTTA‑CCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCG‑CGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGC‑GGCAATTCGCGACGAAACTCGCTCACCAGATTCGAc  <  1:542669/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGGGTTAATGGTCTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGG‑‑TCAGGGAATTAGC‑GCGGGTTTTA‑TCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCG‑CGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGC‑GGCAATTCGCGACGAAACTCGCTCACCAGATTCGAC  >  NC_000913/328646‑328933

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