Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947648,2400TG62.2% ‑2.2 / 25.9 38intergenic (+6/+85)acoR/accCacetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base T (5/9);  new base G (0/23);  total (5/33)
Fisher's exact test for biased strand distribution p-value = 4.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGCTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGG‑CCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCG  >  NC_002947/648107‑648373
                                                                                                                                      |                                                                                                                                      
tGCTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTgg                                                                                                                                      >  1:1737745/1‑136 (MQ=255)
 gCTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTggg                                                                                                                                     >  1:1007152/1‑136 (MQ=255)
  cTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTgggg                                                                                                                                    >  2:202117/1‑136 (MQ=255)
    ggACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTgggggg                                                                                                                                  >  1:4037124/1‑134 (MQ=255)
    ggACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTgggggg                                                                                                                                  >  1:44323/1‑134 (MQ=255)
       cgcgCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGG‑Ccgcg                                                                                                                              <  1:904513/136‑1 (MQ=255)
       cgcgCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGG‑Ccgcg                                                                                                                              <  2:1097251/136‑1 (MQ=255)
         ggtttaTCCGCCATCGGTGTACGCGGGCTTTGGGGGCAGGGGCTTTGGGATTTTC‑GCGGGCGACGATTTATCGGCGTGTGAATGACCCCCGCATCGACATGCCGCGTTTTACGGGGTAGCGCCGTGGGG‑CCGCGAc                                                                                                                            <  2:3187191/131‑1 (MQ=255)
             tagccGCCATGGCTGGACGCGGGCTTTGGGGGCCGGGGCTTTGGGAATTCC‑GCGGGCGACGCTTTATCGGCGAGTGCATCACCACCGCATCCACATGCCGCGTATGACGGGGTAGCCCCGGGGGG‑ccgcgccgcgc                                                                                                                        <  2:11734/133‑2 (MQ=255)
                                             ccggggCTTTGGGCATTCC‑GCGGGCGACGCTTTAGCGGCGAGTGCATGCCCACCGCAGCGACAGGCCGCGGATGATGGGGTAGCGCCCGGGGG‑CCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                        <  2:2403203/133‑1 (MQ=255)
                                             ccggggCTTGGCGCTTTCG‑GCGGGCCCCGCTTTCGCGGGGCGGGCTTGCGCCCCGCACCGACTGGCCGCGGTTGAGGGGGTCGCCCCGGGGGG‑CCGCGCGGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                        <  2:3251229/133‑1 (MQ=255)
                                             caggggCTTTGGG‑ATTTCGGGGGGCTACGTTTTATGGGCCTGTCCTTGAGCACCGCATCGACATGCGGCGTATTATGGGTTAGCGCCGGGGGG‑CCGCGCCGCGCCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                        <  2:3321388/133‑1 (MQ=255)
                                             cACGGTGTTGGGGAATTTC‑GCGGCCGACGATTTATCGGCGTGGGCATGACCACCGAATCTAAAGGCCGCGGATGAAGGGTTCGCGCCCTGGGG‑CGGCTACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                        <  2:2698460/136‑1 (MQ=255)
                                             cACGGGGTTTGGCATTTTC‑GCGGGCGCCGCTTTATCGGCGTGTGCATCAGCACCGAATCGACATGCCGCGTATTAAGGGGTAGCGCCCGGGGG‑CCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                        <  2:1115060/136‑1 (MQ=255)
                                             cACGGGGTTGGGGCTTTTC‑GCGGGGGACGCTTTATGGGCGTTGGCTTGAGCACCGCTTCGACTTGCGGGGGTTGAAGGGTCACCGCCCTGGGG‑CCGGTCCGGGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                        <  2:2673421/136‑1 (MQ=255)
                                             cACGGGCTTTGGGGATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATTAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCGGGGG‑CCTCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                        <  1:3976423/136‑1 (MQ=255)
                                             cACGGGCTTTGGCCTTTGC‑GGGGGCGACGCTTTATCGGCGTGTCCATCCCCCCCGCATCGACTGGCCGCGGATGACGGGGTCGCCCCCGGGGG‑CCGCGCCGCGCCCCCAGATTCTGGCAGGCAATAAAAAACCCCg                                                                                        <  2:2205947/136‑1 (MQ=255)
                                                     ttGGCCTTTGC‑GCGGGCCCCGCTTTTGCGGCGAGGGCATGACCCCCGCATCCACATGCCGCGGATGACGGGGTAGCGCCCGGGGG‑CCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAggg                                                                                <  2:3638234/136‑1 (MQ=255)
                                                        ggCTTTGC‑GGGGGCTCCGCTTTAGGGGCGCGTCCATCCCCACCGCACCGACCGGCGGGGTTTGAGGGGGTAGCCCCGGGGGG‑CCGCGCCGGGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCg                                                                             <  2:4577263/136‑1 (MQ=255)
                                                                 gcgGCCGACTCTTTCTCGGCGTGTCCATTCCCCCCGAATCTACA‑GCCCCGGTTGGAGGGGTACCCCCCCGGGGCCCGCTCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTg                                                                     <  2:883639/136‑1 (MQ=255)
                                                                      cgccGCTTTATCGGCGCGTGCATCCCCCCCGCTCCGCCAGGCCGCGTTTGAGGGGGTCGCGCCGGGGGG‑CCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACa                                                                <  2:16062/133‑1 (MQ=255)
                                                                      cccggCTTTAGCGGCGAGTGCATGACCCCCCCATCGACAGGCCGCGGATGACGGGGTCGCGCCCGGGGG‑CCGCGACGCGCCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACa                                                                <  2:3607912/132‑1 (MQ=255)
                                                                      cGACGTTTTAGCGGCGTGTCCATGCGCACCGCATCTACTTCCCGCGTTGGACGGGGTAGCGCCGGGGGG‑CCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACa                                                                <  2:401369/136‑1 (MQ=255)
                                                                      cGACGCTTTAACGCCGAGGGCACGACCACCGCATCCACATGCCGCGCATGAAGGGCCCCCCCCCTGGGG‑CCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACa                                                                <  2:3114758/136‑1 (MQ=255)
                                                                        ccGCTTTAGCGGCGAGGCCAGCCCCACCGCACCCCCACGCCGCGTATGACGGGGCCGCCCCCGGGGG‑CCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATc                                                              <  2:3650823/135‑1 (MQ=255)
                                                                           cttttgcGGCGTGTGCATCACCACCCATTCCACTTGCCGGGCATTACGGGGTAGCCCCCGGGGG‑CCGCCACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATcaca                                                           <  2:2991896/130‑1 (MQ=255)
                                                                            ttttccGGCGTCTGCTTGCGCACCGCATCCACCTGCCGTGTTTTATGGGGTCCCCCGCGGGGG‑CCGCGCGGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCacac                                                          <  2:2118876/131‑1 (MQ=255)
                                                                                 gggCGTGTGCATCCCCCCCCCATCCACAGCCCGCGTTTTATGGGGTAGCGCCGGGGGG‑CCGCGAGGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTacac                                                     <  1:1065952/135‑1 (MQ=255)
                                                                                 gggCGCGTGCTTGAGCCCCGCAGCGACAGCCCGCGTATGCGGGGGGCGCGCCGGGGGG‑CCGGGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTacac                                                     <  2:976931/135‑1 (MQ=255)
                                                                                          caggAGCACCGCATCCACAGGCCGCTGATCTTGGGTGAGCGCCCGGGGG‑CCGCGTCGCGCCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTcc                                            <  2:4105058/133‑1 (MQ=255)
                                                                                           atccccACCGCATCGACTGGCGGGGTATGAGGGGGGTGCGCGGGGGGG‑CCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCt                                           <  2:4579963/131‑1 (MQ=255)
                                                                                                   cGCATCCACATGCCGCGTTTTATGGGTTCGCGCGGGGGGG‑CCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGcc                                   <  1:1176732/136‑1 (MQ=255)
                                                                                                     aaTCGACTTGCGGCGGTTTATGGGTTAGAGCCGGGGGG‑CCGCGACGGGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAg                                 <  2:1737745/135‑1 (MQ=255)
                                                                                                     aaTCGACATGCCGCGGATGACGGGGTAGAGCCCTGGGG‑CCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAg                                 <  2:1180215/135‑1 (MQ=255)
                                                                                                     aaTCCACATCCCGCGTATGACGGGGACGAGCCCGGGGC‑CCGCGACGCGCCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAg                                 <  2:2233712/135‑1 (MQ=255)
                                                                                                                   ggATGACGGGGTAGAGCCCGGGGG‑CCGCGCCGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGt                   <  2:2214173/136‑1 (MQ=255)
                                                                                                                     ttGATGGGGTAGACCCCTGGGG‑CCGCCACGGGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAg                 <  2:2705385/135‑1 (MQ=255)
                                                                                                                                    ccTGGGG‑CCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgccg  <  2:1821370/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
TGCTGGACGCGCTGATCCGCCATCGCTGGAAGCCGGCTGATGCGGCACGGGCTTTGGGCATTTC‑GCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGG‑CCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCG  >  NC_002947/648107‑648373

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