Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,0850TC35.0% 10.6 / 11.3 20S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (8/5);  new base C (1/6);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 7.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCC  >  NC_000913/70841‑71323
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
cGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGATCCATCCACCGATGGATAATCGGGAACGCGAGGCTTGACAGTACATCAGCGATCAACTGGCAGACAGCAATTTTGATATCGCCACAGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCg                                                                                                                                                                                                                                           >  1:88209/1‑251 (MQ=255)
   tAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGACGCGATTAACGAGTCGCTCCATCCACCGCTGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCCATTTTGGTCTCGCCCGCTTCTCACCGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTCGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATCAGCCAGGCGaa                                                                                                                                                                                                                                         >  1:102158/1‑250 (MQ=255)
                     tGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGACGCGAGGACCAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCAccac                                                                                                                                                                                                                       <  1:5100/250‑4 (MQ=255)
                                cgcATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCATTACATCAGCGATCACCTTGCAGACAGCAATGTTGATATCGCCATCGGCGCCCATGATGTTTGCTTTTCGCCGTCGCGGCTGTCACAGCTTGTCCGTCAGCAGTTAGGGATTAGCGTCTTAAGCCGGCGCGAGGACCAACGCATCAGCCAGGCGAAGCTGCtttt                                                                                                                                                                                                                                <  1:89248/230‑1 (MQ=255)
                                  cATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCACCGCCCGGATGCCTATCg                                                                                                                                                                                                         <  1:49976/251‑1 (MQ=255)
                                     ngAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTTTCCGTACATCAGCGATCACCTGGCCGACAGCAATTTTGCTATCGACAGCGTCGCACCGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGCTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGcc                                                                                                                                                                                                       >  1:20697/2‑250 (MQ=255)
                                                                                                               aCCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAgn                                                                                                                            <  1:33627/251‑2 (MQ=255)
                                                                                                                           gCAATTTTGATATCGCCAGCGTCGCACAGTATGTTTGCTTGTCGCCTTCGCGNCTGNCACANCTNTTCCGCCAGCATTTANGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGNNGNNNCNTCTTTTGa                                                                                                                                                                                                                              <  2:120527/141‑1 (MQ=255)
                                                                                                                           gCAATTTTGATATCGCCAGCGTCGCACAGTATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGACAGCAGTTAGGGATTAGAGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGa                                                                                                                                                                                                                              >  1:120527/1‑141 (MQ=255)
                                                                                                                                            aGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTga                                                                                                 >  1:76883/1‑248 (MQ=255)
                                                                                                                                                  gCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGAGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTCTTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt                                                                                        >  1:31671/1‑251 (MQ=255)
                                                                                                                                                                               tCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCGGGCGCGAGGACCAACGCATTAGTCGTGCGAATCTGCTTTTGAGCTGTACCCGGATGTCTGGCGCCACCGTCTGTCGCAAGGTTTTTTTTAACGATCGACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                                                                                                       <  1:116007/192‑1 (MQ=255)
                                                                                                                                                                                cTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTCGCCGTCAAGTTGTCATAAtt                                                          >  1:156698/1‑251 (MQ=255)
                                                                                                                                                                                    cacaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgc                                                                                                                                                                                                                                       >  1:41841/1‑75 (MQ=255)
                                                                                                                                                                                      caTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGTGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAAc                                                    <  1:159197/251‑1 (MQ=255)
                                                                                                                                                                                           tttCCGCCATCAGTTAGGGATTAGCGTCTTAAGCTGGCGCTAGGACCAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTGGCAATGTTGGTTTTTACGATTAACGCTATTTCTCTCGGGGATTGAAAAAATGCACCGGGGCCAGGCCTAGCGAGTATCGGTCCTTTGTTGAAGAAAAAGTGAATGATGGAGCCGTCAAGTTGTCATAATTGGTAACg                                                   <  1:159168/247‑1 (MQ=255)
                                                                                                                                                                                                             gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCCCCGTCGGTCGCAATGTTGTTTTTGCCGATCAACTCTATTTCTCGCGAGTATTTAAAAAAtn                                                                                                                                         <  1:51721/144‑2 (MQ=255)
                                                                                                                                                                                                                                cgcgAGGACCAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGGTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGTCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTGGACGGAGTAGCATAGGGTTTg          <  1:67335/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aCCAACGCATTAGTCAGGAGAAACTGCTTTTGAGCACTCCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCAAGCACGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACACTTGCCGGCTTGACGGAGTAGCATAGGGTTTGCAGCATc   >  1:59582/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATcn  <  1:46904/251‑2 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
CGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCC  >  NC_000913/70841‑71323

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