Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,277,0620GC35.3% 80.8 / 15.0 51P3P (CCG→CCCfabR_1HTH‑type transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (11/22);  new base C (18/0);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 2.35e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGCCCATGTGGGGCTTACCGCCGGCT  >  ACB143/4276928‑4277192
                                                                                                                                      |                                                                                                                                    
cccTACCACCCAGGCTTCAACGACACCAAATCCCGCCCTGCGCCACTCCGCAAGGAAATGGCCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGt                                                                                                                                    <  8:140679/136‑1 (MQ=255)
   ttCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTTGTAACCGACTGGGGTGCAAGCGCCATTCCCCCCGGCaa                                                                                                                                 >  3:351754/1‑137 (MQ=255)
       aCGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCTTCTGGGGAGCAAGCGCCAATCCCCCCGTCAATGcc                                                                                                                             >  8:146666/1‑137 (MQ=255)
       aCGAAGGCTGCAACGACAGAAAATCCCGCCCTGCTCCACTGCTCAAGGAACTGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGcc                                                                                                                             <  3:130372/137‑1 (MQ=255)
        cgcaGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCGTAACCAACTGAGGCGCCAGACCCCTGCCCCCCGGCAATGccc                                                                                                                            >  5:332785/4‑136 (MQ=255)
            ggCTTCAACGACAGCAAATACCCCCCTGCTCCACTCCTCAAGCAACTGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAAt                                                                                                                        <  6:208590/137‑1 (MQ=255)
             gCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATg                                                                                                                       <  1:8887/137‑1 (MQ=255)
              cTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCGACTGAGGAGCAAGCGCCCAGCCCCCGGGCAAAGCCGCAATgc                                                                                                                      >  4:307570/1‑137 (MQ=255)
                tCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTATTAACCAACATAGGAGCAAGCCCCATGCCCCCCGTCAAATCCGCCATgcgc                                                                                                                    >  2:213061/1‑137 (MQ=255)
                  aaCGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCtg                                                                                                                  <  8:200647/137‑1 (MQ=255)
                    cGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTAACCAAATGTGGACCAAGCGACAAGCCCCCCGGCAATTCCGCAATTCGCtgtt                                                                                                               >  8:333587/1‑136 (MQ=255)
                    cGACAGCAAATACCGCCCTGCGCCACTGCTCACGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCtgtg                                                                                                               <  4:51476/137‑1 (MQ=255)
                          cAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTGACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACt                                                                                                         >  6:163418/1‑137 (MQ=255)
                                  gCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGGTGTTGACCTACTGAGGAGCCACCCCCAAGCCCCCGGGGGATGCCGCCAAGCGCGGTGGCCACATTgtggag                                                                                                 >  1:309617/1‑137 (MQ=255)
                                   cccTGCTCGAGTCCTCAAGTATTGGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGc                                                                                                <  5:251486/137‑1 (MQ=255)
                                     cTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCTACTGAGGCGCAAGCGCCCTGCCCCCCGGCAATGCCGCCATGCGCGGGGGCCACGTTGTGGAGCgc                                                                                              >  4:45852/1‑136 (MQ=255)
                                        cTCCAGTGCTCAAGAAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGt                                                                                           <  1:1240/137‑1 (MQ=255)
                                           cACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGACCTACTGAGGAGCAAGCGCCCTGCCCCCCGGCAATGCCGCAATGCGCTGTGGCAACATTGAGGAGCGACGTGaa                                                                                       >  4:153482/1‑136 (MQ=255)
                                            aCTGCTCCAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGAcc                                                                                      <  8:261906/137‑1 (MQ=255)
                                               ccTCAAGGAAATGGGCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCggg                                                                                   <  6:45414/136‑1 (MQ=255)
                                                cTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTGTTAACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCCGGc                                                                                  >  2:103397/1‑137 (MQ=255)
                                                cTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTACACGACTGAGGAGCAAGCCCCAAGCCCCCCGTAAATGCCGCAAAGCGCGGGTCCAACATTGAGGAGCGACGTGACCgggg                                                                                  >  7:95686/1‑136 (MQ=255)
                                                   aaGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATTCCGCAATGCGCGGTGCCAAATTTGAGGAGCGACGTGCCCGGGCCAt                                                                               >  8:117385/1‑137 (MQ=255)
                                                   aaGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTTTTCACCTAAAGAGGCGCAAGCGCCATGCCCCCCGGCAATTCCGCAATGCGGGGTGGCAACATTTGGGCGCGACGTGACCGGGCCaa                                                                               >  6:118166/1‑136 (MQ=255)
                                                   aaGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTCTTGACCGACTGAGGAGGAAGCGCCCTGCCCCCCGTCAAAGCCGCAATGCGCTGTGGCAACTTTTAGGAGCGCCCTGACCGGGccca                                                                               >  8:34223/1‑135 (MQ=255)
                                                     ggAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCCCCATGCCCCCGGGCAATGCCGCAAAGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATgg                                                                             >  5:160784/1‑137 (MQ=255)
                                                     ggAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCCACAGAGGAGCAACCCCCAATCCCCCCGTGAAATCCCCAAATCGCTGTGCCCAATTTGAGGAGCGACGTGACCGGGCCAttg                                                                             >  4:33383/1‑135 (MQ=255)
                                                     ggAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGGTGTTCACCGACGGAGGCGCGACCCCCCAGCCCCCCGGCACTGCCGCCACTCGCGGGGGCACCGTTTAGGGGCGACGTGGCCCGGGCCagg                                                                              >  2:191150/1‑135 (MQ=255)
                                                     ggAAATGGTCGCGGCAGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTCACCAAATTAGGAGCAAGCCCCATGCCCCCCGGAAAATCCGCAAAGCGCAGTTCCCAACTTTATGAGCGACGGGACCGGGGCATgg                                                                             >  8:357634/1‑137 (MQ=255)
                                                             tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTTTTGACCGAATGAGGGGCAAGCGCCCTGCCCCCCGGGAATGGCCCAATGCGCGGGGGCAAAGTTGAGGAGCGACCTGGCCCGGCCATGGCAATGgtt                                                                     >  6:40418/1‑137 (MQ=255)
                                                             tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTGCCCCTCGGTGGGGGAAGCGCCCTGCCCCCCGGCCAAGCCGCAATGCGGGGTGGCCACATTGAGGGGCGGCCTGACCGGGCCATGGGAATGttt                                                                     >  2:257240/1‑137 (MQ=255)
                                                             tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGTTCACCTACTGAGGAGCAAGCGCCAAGCCCCCGGGCAATGCCGCCATGCGCTGTGGCAAAGTTGAGGCGCGACGTGGCCGGGCCCTGGGATTGgtt                                                                     >  7:290164/1‑137 (MQ=255)
                                                             tCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGGTTGGTGACCTACTGAGGATCAAGCGCCAATACCCCCGTCCATGCCGCCCAGAGCTGTGGCAACTATTAGGAGCGACGTGGACGGGCCATGGCCtttttt                                                                     >  4:359625/1‑137 (MQ=255)
                                                                cggccggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCGCCGACTGAGGAGCAACCGCCATGCCCCCCGGCAATGCCCCAATGCGCGGTGCCAACTTTGTGGAGCGCCGGTACCGGGCCCTGGGATTGTTTGcc                                                                  >  8:129339/1‑137 (MQ=255)
                                                                  gccggcTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGa                                                                <  2:282119/137‑1 (MQ=255)
                                                                       gtggctggGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAgg                                                           <  7:129339/136‑1 (MQ=255)
                                                                       ctggctggGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAgg                                                           <  7:18306/137‑1 (MQ=255)
                                                                             ggggCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGATTGCCGAAAAGGGCTTTg                                                     <  7:54956/137‑1 (MQ=255)
                                                                                 cgcAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGCCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCa                                                 <  1:197281/137‑1 (MQ=255)
                                                                                            tCGCGGGTTGTTGACCTACTGAGGGGCAAGCGCCAAGCCCCCCGTCAATTCCGCAATGCGCGGTGCCCACGTTGAGGAGCGCCGTGACCGGGCCATGGGATTGTTTGGCGAAAAGGGCTTTGGCCAGGTCAGCATgc                                      >  4:131381/1‑137 (MQ=255)
                                                                                                      ttcACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTggc                            <  6:168205/137‑1 (MQ=255)
                                                                                                      ttcACCTACTGAGGAGCAAGCGCCAGGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGCCGGGACCGGGCCATGGCATTGTTTGCCGAAACGGGCTTTGGCCAGGTCATCATGCGCGAGCTggc                            <  2:221573/137‑1 (MQ=255)
                                                                                                          ccTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGcc                        <  3:153482/137‑1 (MQ=255)
                                                                                                               tGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGCCCAtgt                   >  1:62241/1‑137 (MQ=255)
                                                                                                                  ggagCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGCCCATGTggg                <  7:333587/137‑1 (MQ=255)
                                                                                                                   gagCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGa                                                                >  8:119731/1‑88 (MQ=255)
                                                                                                                   gagCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGa                                                                <  7:119731/88‑1 (MQ=255)
                                                                                                                      cAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGCCCATGTGGGGCtt            >  5:44932/1‑137 (MQ=255)
                                                                                                                           gcCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGCCCATGTGGGGCTTAccgc       >  1:289630/1‑137 (MQ=255)
                                                                                                                             cATGACCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCAGTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGCCCATGTGGTGCTTAccgccg     <  1:90285/137‑1 (MQ=255)
                                                                                                                                gCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGCCCATGTGGGGCTTACCGCCGGCt  <  3:362836/137‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
GCCTTCCACGAAGGCTTCAACGACAGCAAATACCGCCCTGCTCCACTGCTCAAGGAAATGGTCGCGGCCGGCTGGCTGGGGCGCAAGAGCGGTCGCGGTTTCTTCACCTACTGAGGAGCAAGCGCCATGCCCCCGGTCAATGCCGCAATGCGCTGTGCCAACTTTGAGGAGCGACGTGACCGGGCCATGGCATTGTTTGCCGAAAAGGGCTTTGGCCAGGTCAGCATGCGCGAGCTGGCGGCCCATGTGGGGCTTACCGCCGGCT  >  ACB143/4276928‑4277192

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