Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091334,5080GC34.4% 56.2 / 16.9 32G70A (GGA→GCA) caiFcai operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (16/5);  new base C (0/11);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 6.77e-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
                                                                   |                                                                 
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:516567/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:124069/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:733993/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:473486/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:299765/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  1:281156/1‑68 (MQ=255)
atTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCgggg                                                                   >  2:832761/1‑68 (MQ=255)
  tCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTTAGATGATCCCTAATAAGCTGGAAGGGCGGGGAt                                                                 >  1:927915/1‑68 (MQ=255)
  tCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGAt                                                                 >  1:33988/1‑68 (MQ=255)
  tCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGAt                                                                 >  1:634411/1‑68 (MQ=255)
  tCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCAt                                                                 <  1:439699/68‑1 (MQ=255)
  tCTGTCGGAAGTCACAGAAAGAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGAt                                                                 >  2:155588/1‑68 (MQ=255)
                   aaaTAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTg                                                <  2:296600/68‑1 (MQ=255)
                    aaTATGCTGCTATGTTAAGATTATCCCAAATATGCTGTAAGGGCGGGCATGCCAGTGTCAGCGACTgg                                               <  2:281158/68‑1 (MQ=255)
                    aaTATGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTgg                                               <  1:117794/68‑1 (MQ=255)
                    aaTAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTgg                                               <  1:372003/68‑1 (MQ=255)
                    aaTAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTgg                                               <  2:507673/68‑1 (MQ=255)
                    aaTAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGCATGCCAGTGTCAGCGACTgg                                               <  1:124069/68‑1 (MQ=255)
                           tGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGTGGATGCCAGTGTCAGCGACTGGTTAAAGt                                        <  2:288564/68‑1 (MQ=255)
                             cGAAGTTATGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTgg                                      <  2:765679/68‑1 (MQ=255)
                              gAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGt                                     >  2:37556/1‑68 (MQ=255)
                                      gatgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTgg                                               >  2:519584/1‑50 (MQ=255)
                                      gatgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTgg                                               <  1:519581/50‑1 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  2:762833/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  2:329942/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  2:399254/1‑68 (MQ=255)
                                        tgatCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGAt                           >  2:794839/1‑68 (MQ=255)
                                                  taaGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAAttt                 >  1:698492/1‑68 (MQ=255)
                                                    aGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTAc               >  2:190933/1‑68 (MQ=255)
                                                     gcaggAAGGGCGGGGATGCCAGTGTCAGAGACTGGATAAAGTGGTCGATATCGATGAGCAAATTTAcg              >  2:867946/4‑68 (MQ=255)
                                                     gCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTAcg              >  2:678414/1‑68 (MQ=255)
                                                      cTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTAcgc             >  1:296000/1‑68 (MQ=255)
                                                             ggCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTg      >  1:211001/1‑68 (MQ=255)
                                                               cGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  1:190933/68‑1 (MQ=255)
                                                               cGGGCATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  2:539128/68‑1 (MQ=255)
                                                               cGGGCATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  1:516564/68‑1 (MQ=255)
                                                               cGGGCATGCCAGTGTCAGCGACTGGTTAAAGTGGGCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  2:786598/68‑1 (MQ=255)
                                                               cGGGCATGCCAGCGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  2:787890/68‑1 (MQ=255)
                                                               cGGGCACGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTgcg    <  2:265779/68‑1 (MQ=255)
                                                                 gtgATGCCAGTGTCAGCGACTTGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTGCGCa  <  1:867942/66‑1 (MQ=255)
                                                                   |                                                                 
ATTCTGTCGGAAGTCACAGAAATAAGCTGCGAAGTTAAGATGATCCCTAATAAGCTGGAAGGGCGGGGATGCCAGTGTCAGCGACTGGTTAAAGTGGTCGATATCGATGAGCAAATTTACGCGCGCCTGCGCA  >  NC_000913/34441‑34573

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