Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,223,5970CT29.5% 176.3 / 39.0 88R11K (AGA→AAA) iclRDNA‑binding transcriptional repressor IclR
Reads supporting (aligned to +/- strand):  ref base C (17/45);  new base T (26/0);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 2.88e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTG  >  NC_000913/4223549‑4223646
                                                |                                                 
aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTc                                                   <  1:3489312/49‑1 (MQ=255)
aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTc                                                   <  1:2572243/49‑1 (MQ=255)
aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTc                                                   <  1:5709994/49‑1 (MQ=255)
aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTc                                                   <  1:1668155/49‑1 (MQ=255)
 aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTg                                                 <  1:8097743/50‑1 (MQ=255)
 aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTg                                                 <  1:153304/50‑1 (MQ=255)
 aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTg                                                 >  1:5554294/1‑50 (MQ=38)
 aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTg                                                 >  1:5725824/1‑50 (MQ=38)
  aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGc                                                <  1:760473/50‑1 (MQ=255)
  aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTgg                                                >  1:3763480/1‑49 (MQ=38)
   gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGcc                                               <  1:6834958/50‑1 (MQ=255)
    aCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcg                                              >  1:4828880/1‑50 (MQ=38)
     cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCcgc                                             <  1:1146999/50‑1 (MQ=255)
     cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGTcgc                                             >  1:6895622/1‑50 (MQ=38)
     cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCcgc                                             >  1:7597003/1‑50 (MQ=38)
      tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcgcg                                            >  1:6658865/1‑50 (MQ=38)
       gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCCCGt                                           >  1:7608458/1‑50 (MQ=38)
        aaCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcgcgct                                          >  1:4673084/1‑48 (MQ=37)
         aCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGttt                                         >  1:5433212/1‑50 (MQ=38)
           cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcg                                       <  1:6775865/50‑1 (MQ=255)
           cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcg                                       <  1:2015457/50‑1 (MQ=255)
           cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGgtttg                                       >  1:6596963/1‑48 (MQ=25)
           cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGgtttg                                       >  1:8133371/1‑48 (MQ=37)
           cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTcg                                       >  1:3698930/1‑50 (MQ=39)
            tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTtgc                                      >  1:4767192/1‑50 (MQ=25)
              tCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcgc                                      <  1:4234470/48‑1 (MQ=255)
                cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCggg                                   <  1:248691/49‑1 (MQ=255)
                cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGa                                  <  1:3229893/50‑1 (MQ=255)
                cAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTGCGCGGTTCGCggg                                   >  1:1447648/1‑49 (MQ=25)
                 aGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCgggga                                 >  1:7585229/1‑48 (MQ=255)
                 aGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCgggg                                  >  1:5446809/1‑48 (MQ=255)
                  gTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTTGCGGGGAt                                >  1:2600125/1‑50 (MQ=38)
                  gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCGCGGTTCGCGGGGAt                                >  1:3385052/1‑50 (MQ=255)
                  gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAAt                                <  1:3238814/50‑1 (MQ=255)
                  gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAAt                                <  1:1528631/50‑1 (MQ=255)
                   tCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATg                               <  1:1369607/50‑1 (MQ=255)
                   tCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATg                               <  1:2998199/50‑1 (MQ=255)
                    cGCTGGTGCGGTGGCAACGGCGGGGTTTTTTTCCCGTTTTGCGGGGATgg                              >  1:1963089/1‑50 (MQ=16)
                    cGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTCGCGGGGATgg                              >  1:4303269/1‑50 (MQ=255)
                     gCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCGCGGTTCGCGGGGATtgg                             >  1:2659997/1‑50 (MQ=25)
                     gCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATggg                             <  1:7876877/50‑1 (MQ=255)
                     gCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATggg                             <  1:4403070/50‑1 (MQ=255)
                     gCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATgg                              <  1:4153162/49‑1 (MQ=255)
                     gCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGGAAggg                             >  1:517991/1‑50 (MQ=255)
                     gCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGCTTCGCGGGAATggg                             >  1:3792885/1‑50 (MQ=37)
                     gCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGATggg                             >  1:5953791/1‑50 (MQ=255)
                     gCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTCGCGGGGAAggg                             >  1:1637315/1‑50 (MQ=255)
                      cTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGt                            <  1:6477979/50‑1 (MQ=255)
                      cTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTCCGCGGGAATGGGt                            <  1:7673466/50‑1 (MQ=255)
                      cTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGt                            >  1:4886117/1‑50 (MQ=255)
                      cTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGATGGGt                            >  1:1798908/1‑50 (MQ=255)
                        ggTGCGGTGGCAACGGCGGGTTTTTTTCCGCGTTTCCCGGGGATAGGGGc                          >  1:4222346/1‑50 (MQ=19)
                         gTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg                           <  1:5890881/48‑1 (MQ=255)
                         gTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGc                          <  1:2814850/49‑1 (MQ=255)
                         gTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCg                         <  1:880248/50‑1 (MQ=255)
                         gTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCg                         <  1:111783/50‑1 (MQ=255)
                          tGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGa                        <  1:2870279/50‑1 (MQ=255)
                            cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGAcc                      <  1:259592/50‑1 (MQ=255)
                            cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGAcc                      <  1:3332734/50‑1 (MQ=255)
                             ggTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCa                     <  1:1432568/50‑1 (MQ=255)
                              gTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCAt                    <  1:7454228/50‑1 (MQ=255)
                              gTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCAt                    <  1:5973117/50‑1 (MQ=255)
                                 gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGa                  <  1:4404083/49‑1 (MQ=255)
                                 gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGa                  <  1:2887654/49‑1 (MQ=255)
                                 gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAc                 <  1:1593132/50‑1 (MQ=255)
                                 gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAc                 <  1:5440991/50‑1 (MQ=255)
                                    aCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt              <  1:4894489/50‑1 (MQ=255)
                                    aCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt              <  1:3568093/50‑1 (MQ=255)
                                     cggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtc             <  1:1710235/50‑1 (MQ=255)
                                      ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct            >  1:1689691/1‑50 (MQ=255)
                                       gcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtctc           <  1:136963/50‑1 (MQ=255)
                                       gcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtctc           <  1:837013/50‑1 (MQ=255)
                                        cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct          >  1:3857229/1‑50 (MQ=255)
                                        cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct          >  1:2298915/1‑50 (MQ=255)
                                        cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct          <  1:7989724/50‑1 (MQ=255)
                                        cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct          >  1:8238211/1‑50 (MQ=255)
                                         gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt         >  1:1660802/1‑50 (MQ=255)
                                         gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt         >  1:130813/1‑50 (MQ=255)
                                         gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt         >  1:844830/1‑50 (MQ=255)
                                          ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttt        <  1:765380/50‑1 (MQ=255)
                                          ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttt        >  1:3040590/1‑50 (MQ=255)
                                          ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttt        >  1:5055051/1‑50 (MQ=255)
                                          ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttt        <  1:8078027/50‑1 (MQ=255)
                                            ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttttt      <  1:6731052/50‑1 (MQ=255)
                                            ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttttt      <  1:1700070/50‑1 (MQ=255)
                                            ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttttt      >  1:5559726/1‑50 (MQ=255)
                                             tttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtttttt     >  1:3145805/1‑50 (MQ=255)
                                              ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc    >  1:7292687/1‑50 (MQ=255)
                                              ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc    >  1:2490530/1‑50 (MQ=255)
                                               tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc    <  1:2789163/49‑1 (MQ=255)
                                               tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc    <  1:482721/49‑1 (MQ=255)
                                               tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt   >  1:700868/1‑50 (MQ=255)
                                                cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt   <  1:1776619/49‑1 (MQ=255)
                                                cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt   <  1:5490392/49‑1 (MQ=255)
                                                cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt   <  1:2269911/49‑1 (MQ=255)
                                                cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTg  >  1:3470749/1‑50 (MQ=255)
                                                |                                                 
AAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTG  >  NC_000913/4223549‑4223646

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