Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091334,5080GC35.7% 63.6 / 32.7 42G70A (GGA→GCA) caiFcai operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (21/6);  new base C (0/15);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.10e-06
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
                                                                   |                                                                 
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:970818/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:1137458/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:137395/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:42831/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:986056/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:863391/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:8174/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCACTAATAAGCTGGAAGGGCgggg                                                                   >  1:639843/1‑68 (MQ=255)
  tCTGTCGGAAGTCACAGATATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGAt                                                                 >  1:257298/1‑68 (MQ=255)
   cTGGCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATg                                                                <  2:807277/68‑1 (MQ=255)
            gTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTcg                                                       >  2:875417/1‑67 (MQ=255)
            gTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCa                                                       >  2:155584/1‑68 (MQ=255)
            gTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCa                                                       >  2:850064/1‑68 (MQ=255)
               acaTAAATAAGCTGCGATGTTATGATGACCCCTAATATGCTGGATGGCCGGGCATGCCAGTGTCAGCg                                                    <  2:395300/68‑1 (MQ=255)
                 aGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGAc                                                  <  2:35285/68‑1 (MQ=255)
                   aaaTAAGCTGCGAAGTTTAGATGATCCCTAATAAGCTGGATGGGCGGGCATGCCAGTGTCAGCGACTg                                                <  2:68011/68‑1 (MQ=255)
                   aaaTAAGCTGCGAAGTTAAGATGATCCCTAATAATCTTGAAGGGCGGGCATGCCAGTGTCAGCGACTg                                                <  1:39242/68‑1 (MQ=255)
                   aaaTAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGTAAGGGCGGGCATGCCAGTGTCAGCGACTg                                                <  1:714590/68‑1 (MQ=255)
                   aaaTAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTg                                                <  1:1192272/68‑1 (MQ=255)
                       aaGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTa                                            >  1:183549/1‑68 (MQ=255)
                        aGCTGCGAAGTTAAGATGATCCCTAATAAGCTGTAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTaa                                           <  1:879237/68‑1 (MQ=255)
                         gCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTaaa                                          >  2:648569/1‑68 (MQ=255)
                         gCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTaaa                                          >  2:1162799/1‑68 (MQ=255)
                         gCTGAGAAGTTAAGATGATCCCTAATAAGCGGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTtaa                                          >  2:763073/1‑68 (MQ=255)
                           tGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGt                                        <  2:176360/68‑1 (MQ=255)
                             ctaaGTTAAGACGATCCCTAATAAGCTGGACGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGTgg                                      <  2:437567/66‑1 (MQ=255)
                             cGAAGTTAAGATGATCCCTAATATGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGTgg                                      <  2:877559/68‑1 (MQ=255)
                             cGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGCgg                                      <  2:183549/68‑1 (MQ=255)
                             cGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGTgg                                      <  1:657764/68‑1 (MQ=255)
                             cGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGTgg                                      <  1:171054/68‑1 (MQ=255)
                                  ttAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGGCAGCGACTGGTTAAAGTGGTCGat                                 <  2:1054669/68‑1 (MQ=255)
                                      gatgatCCCTATTAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCg                             >  1:567366/1‑68 (MQ=255)
                                      gatgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCg                             >  2:488978/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCGGCGACTGGTGAAAGTGGTCGATATCGAt                           >  2:1034055/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTTAAGTGGTCGATATCGAt                           >  1:1058608/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  2:606802/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  1:487041/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  2:725313/1‑68 (MQ=255)
                                          atCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGa                         <  1:178564/68‑1 (MQ=255)
                                             ccTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCa                      >  2:1045912/1‑68 (MQ=255)
                                                aataaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAAt                   <  2:921532/68‑1 (MQ=255)
                                                  tacgCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTTGATAAAGTGGTCGAAATCGATGGGCAAAttt                 >  1:271388/4‑68 (MQ=255)
                                                  taaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGat                                 <  2:1183783/52‑1 (MQ=255)
                                                  taaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGat                                 >  1:1183785/1‑52 (MQ=255)
                                                   aaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTa                >  1:1172284/1‑68 (MQ=255)
                                                     gCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTAcg              >  1:1137464/1‑68 (MQ=255)
                                                         gacgggcgggcAGCCCAGTGTCAGCGACTGTGTAAAGTGGTCGATATCGATGAGCAAACTTAcgcgcg          <  2:1137456/65‑1 (MQ=255)
                                                         gAAGGGCGGGCATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGACGAGCAAATTTAcgcgcg          <  2:676417/68‑1 (MQ=255)
                                                          aaGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACgcgcgc         >  2:713032/1‑68 (MQ=255)
                                                              gCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgc     <  1:1045913/68‑1 (MQ=255)
                                                               cGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  2:863391/68‑1 (MQ=255)
                                                                 gggATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTGCGCa  <  2:1137462/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: