Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091334,5080GC25.7% 85.2 / 13.9 36G70A (GGA→GCA) caiFcai operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (22/4);  new base C (0/9);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.01e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTGCGCA  >  NC_000913/34441‑34573
                                                                   |                                                                 
atTCTGTCGGAAGTCGCAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:36418/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAGGCTGCGAAGTTAAGATGATCACTAATAAGCTGGAAGGGCgggg                                                                   >  2:1141312/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:237459/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:940383/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:728510/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:668421/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGATAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:239202/1‑68 (MQ=255)
  tCTGTCGGAAGTCACAGATATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGAt                                                                 >  2:896564/1‑68 (MQ=255)
  tCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGAt                                                                 >  1:948193/1‑68 (MQ=255)
  tCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGAt                                                                 >  2:91218/1‑68 (MQ=255)
        ggATGTCACTGAAATAAGCTGCGAAGTTAATATTATCCCTAATAAGCTGGAAGGGCGGGCATGCCatt                                                           <  1:668423/68‑3 (MQ=255)
            gTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCa                                                       >  1:555107/1‑68 (MQ=255)
                   ctaTAAGCTGCGAAGTTATGATGATCCCTACTAAGCTGGAAGGGCGGGCACGCCAGTGCCAGCGACTg                                                <  2:483201/66‑1 (MQ=255)
                   aaaTATGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTg                                                <  2:529646/68‑1 (MQ=255)
                   aaaTAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTg                                                >  2:672344/1‑68 (MQ=255)
                    aaTAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTgg                                               >  2:500626/1‑68 (MQ=255)
                         gcagCGAAGTGAAGATGGTCCATGAAAAGCTGAAAGGGCGGGGATGCCTGAGTCAGCGACAGGTTaaa                                          >  2:1138929/4‑68 (MQ=255)
                         gCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTaaa                                          >  1:54129/1‑68 (MQ=255)
                         gCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCAAGTGTCAGCGACTGGTTaag                                          >  2:115131/1‑67 (MQ=255)
                           tGCGATGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGt                                        <  2:606561/68‑1 (MQ=255)
                                  ttAAGTTTATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGat                                 <  2:920754/68‑1 (MQ=255)
                                      gatgatCCCTAATAAGCTGGAAGGGAGGGGATGCAAGTGTCAGAGACTGGTTTAAGTGGTCGatatag                             >  1:1073844/1‑66 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  2:120667/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  2:744262/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  1:437226/1‑68 (MQ=255)
                                            cccTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGc                       <  1:1078205/68‑1 (MQ=255)
                                                 atatGCTGGAAGGGCGTTCATGCCAGTGTCAGCTACTGGTTAAAGTGGGCGATATCGATGAGCATAtt                  <  1:304417/68‑1 (MQ=255)
                                                 ataaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAAtt                  >  1:61661/1‑68 (MQ=255)
                                                  taaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTgg                                               <  1:909230/38‑1 (MQ=255)
                                                  taaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTgg                                               <  1:751615/38‑1 (MQ=255)
                                                  taaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTgg                                               >  2:909227/1‑38 (MQ=255)
                                                  taaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTgg                                               >  2:751613/1‑38 (MQ=255)
                                                  taaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATTTCGATGAGCAAAttt                 >  1:77088/1‑68 (MQ=255)
                                                   aaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAGGTGGTCGATATCGATGAGCAAAtttg                >  2:310258/1‑67 (MQ=255)
                                                   aaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTa                >  1:570838/1‑68 (MQ=255)
                                                     gCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTAcg              >  2:157394/1‑68 (MQ=255)
                                                     gCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTAcg              >  2:2891/1‑68 (MQ=255)
                                                         caaGGGCGGGCATGCCAGTGTCAGCGACTGGGTAAAGTGGCCGATATCGACGAGCAAATCTAcgcgcg          <  2:815746/67‑1 (MQ=255)
                                                          aaGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAGTTTCCgcgcgc         >  1:805156/1‑68 (MQ=255)
                                                              gCGGGCATGCCAGTGGCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgc     <  2:253515/68‑1 (MQ=255)
                                                               cGGGTATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  2:863527/68‑1 (MQ=255)
                                                               cGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATTTCGATGAGCAAATTTACGCGCGCCTgcg    <  2:77088/68‑1 (MQ=255)
                                                               cGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  2:364409/68‑1 (MQ=255)
                                                                 gggATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTGCGCa  <  2:61661/68‑1 (MQ=255)
                                                                   |                                                                 
ATTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTGCGCA  >  NC_000913/34441‑34573

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: