Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077794,062,1140TG42.9% 84.8 / 29.1 63G142G (GGT→GGG) gntRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (9/27);  new base G (26/1);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.37e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGCCTCGAATCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGATGACCACTGCCATTA  >  NC_007779/4062017‑4062214
                                                                                                 |                                                                                                    
aaCGCCTCGAATCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTAtt                                                                                                   <  1:785304/101‑1 (MQ=255)
aaCGCCTCGAATCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGAtt                                                                                                   >  3:56242/1‑101 (MQ=255)
     cTCGAATCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCgt                                                                                              <  3:811605/101‑1 (MQ=255)
        gAATCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAGGTGGCGGGTATTCCCgtggt                                                                                           >  1:520900/1‑101 (MQ=255)
        gAATCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCgtggt                                                                                           <  3:8381/101‑1 (MQ=255)
         aaTCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCgtggtg                                                                                          <  4:4115/101‑1 (MQ=255)
         aaTCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATGCCCgtggtg                                                                                          <  4:85895/101‑1 (MQ=255)
          aTCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTccccgggggg                                                                                         >  3:1428608/1‑94 (MQ=255)
          aTCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGtggggg                                                                                         >  1:617291/1‑101 (MQ=255)
              aTGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTTCCGTGGGGGAACt                                                                                     >  1:112418/1‑101 (MQ=255)
              aTGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGGACt                                                                                     >  1:126657/1‑101 (MQ=255)
                gCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGGACTgg                                                                                   >  1:1190495/1‑100 (MQ=255)
                 ctctCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGAt                                                                                  <  3:1272000/101‑1 (MQ=255)
                 ctctCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGGACTGAt                                                                                  >  3:1404858/1‑101 (MQ=255)
                  tctcCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGGACTGATg                                                                                 >  1:1374921/1‑101 (MQ=255)
                     ccTGGCCTATAGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGGACTGATGGAc                                                                              >  4:90062/1‑101 (MQ=255)
                     ccTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGGACTGATGGAc                                                                              >  4:36359/1‑101 (MQ=255)
                      cTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACa                                                                             <  3:641606/101‑1 (MQ=255)
                      cTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACa                                                                             <  3:555365/101‑1 (MQ=255)
                            tatCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGt                                                                       <  3:1282008/101‑1 (MQ=255)
                             atCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTc                                                                      <  1:1458651/101‑1 (MQ=255)
                              tCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCg                                                                     <  3:26172/101‑1 (MQ=255)
                               cGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGc                                                                    <  2:159461/101‑1 (MQ=255)
                               cGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGc                                                                    <  1:1191797/101‑1 (MQ=255)
                                gACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTTCCGTGGTGGGACTGATGGGCAGCAAGTTGcc                                                                   >  3:1367695/1‑101 (MQ=255)
                                gACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGTGGAACTGATGGACAGCAAGTCGcc                                                                   >  1:1535967/1‑101 (MQ=255)
                                 aCGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCa                                                                  <  1:63854/101‑1 (MQ=255)
                                 aCGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCCTGGGGGGACTGATGGACAGCAAGTCGCCa                                                                  >  1:1090025/1‑101 (MQ=255)
                                     ccTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTTCCGTGGGGGAACTGGTGGACAGCAAGGCGCCATGcc                                                              >  1:1396619/1‑101 (MQ=255)
                                      cTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCt                                                             <  1:1009390/101‑1 (MQ=255)
                                       tGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCACGccct                                                            >  3:1226475/1‑99 (MQ=255)
                                       tGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGAACTGATGGACAGCAAGTCGCCATGCCtt                                                            >  3:781291/1‑101 (MQ=255)
                                        gATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGTGGAACTGATGGACAGCAAGGCGCCATGCCTTg                                                           >  1:937761/1‑101 (MQ=255)
                                         aTCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGa                                                          <  3:468518/101‑1 (MQ=255)
                                          tCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGat                                                         <  3:1440418/101‑1 (MQ=255)
                                            cTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGTGGGACTGATGGACAGCAAGTCGCCATGCCTTGatat                                                       >  1:7113/1‑101 (MQ=255)
                                             tCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATc                                                      <  3:832168/101‑1 (MQ=255)
                                                ccGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGGACTGATGGACAGCAAGGCGCCATGCCTTGATATCGcc                                                   >  1:858970/1‑101 (MQ=255)
                                                   aaCGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCGCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTc                                                <  3:1297875/101‑1 (MQ=255)
                                                      gTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCCCGGGGGGACTGATGGACAGCAAGAGGCCCCGGCCTTATATTGCCGTCGGt                                             >  3:280738/1‑101 (MQ=255)
                                                      gTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGt                                             >  4:56052/1‑101 (MQ=255)
                                                       tACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGGACTGATGGACAGCAAGTCGCCCTGCCCTGATATTGCCGTCgggt                                            >  3:292206/1‑99 (MQ=255)
                                                        aCCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGttt                                           <  1:1027980/101‑1 (MQ=255)
                                                           cacaCGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGGGGGGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGa                                        >  4:49180/1‑101 (MQ=255)
                                                                gccgcGCACCTTAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGAACTGATGGGCAGCAAGTCGCCATGCCTTGATATCGCCGTCGGGTTTGATAACt                                   >  1:83846/1‑101 (MQ=255)
                                                                   gcgcACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTg                                <  3:1100334/101‑1 (MQ=255)
                                                                    cgcACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGACGCCATGCCCTGATATCGCCGTCGGTTTTTATAAATTTGa                               >  3:516773/1‑101 (MQ=255)
                                                                      cACCTTACAGATGATTGCAGTGGCGGGGATTCCCGGGGGGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAag                             >  4:49851/1‑101 (MQ=255)
                                                                      cACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAag                             <  3:1548663/101‑1 (MQ=255)
                                                                         cTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGGGGAACCGATGGGAAGCAAAGCGCCATGCCTTGATATCGCCGTCGGTTTTGATAAATTTGAagcag                          >  1:200802/1‑101 (MQ=255)
                                                                           tAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAgcagca                        <  3:1425078/101‑1 (MQ=255)
                                                                           tAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAgcagca                        <  1:846340/101‑1 (MQ=255)
                                                                           tAAAGATGATTGAAGTGGCGGGGATTCCCGTGGGGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAgcagca                        >  3:1357963/1‑101 (MQ=255)
                                                                              aGATGATTGAAGTGGCGGGTATTCCCGTGGGGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGc                     >  1:1117697/1‑101 (MQ=255)
                                                                              aGATGATTGAAGTGGCGGGGATTCCCGTGGGGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGGTTTGATAACTTTGAAGCAGCACGc                     >  1:1189677/1‑101 (MQ=255)
                                                                                atgatTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCa                   <  1:1095530/101‑1 (MQ=255)
                                                                                atgatTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCa                   <  3:340619/101‑1 (MQ=255)
                                                                                   atTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGAt                <  1:796734/101‑1 (MQ=255)
                                                                                   atTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGAt                <  1:230061/101‑1 (MQ=255)
                                                                                    ttGAAGTGGCGGGTATTCCCGTGGGGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGATg               >  1:59190/1‑101 (MQ=255)
                                                                                    ttGAAGTGGCGGGGCTTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGATg               >  4:154555/1‑101 (MQ=255)
                                                                                      gAAGTGGCGGGTATTCCCGTGGGGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGATGAc             >  1:637762/1‑101 (MQ=255)
                                                                                             cGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGATGACCACTGcc      >  1:828427/1‑101 (MQ=255)
                                                                                                 tATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGATGACCACTGCCatta  >  3:177410/1‑101 (MQ=255)
                                                                                                 tATTCCCGTGGAGGAACTGTTGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGATGACCACTGCCatta  >  1:718786/1‑101 (MQ=255)
                                                                                                 |                                                                                                    
AACGCCTCGAATCCATGCTCTCCTGGAATATCGACGGCCTGATCCTCACCGAACGTACCCACACGCCGCGCACCTTAAAGATGATTGAAGTGGCGGGTATTCCCGTGGTGGAACTGATGGACAGCAAGTCGCCATGCCTTGATATCGCCGTCGGTTTTGATAACTTTGAAGCAGCACGCCAGATGACCACTGCCATTA  >  NC_007779/4062017‑4062214

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