Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,208,2300TA27.6% 32.5 / 16.8 29noncoding (84/1542 nt)rrsE16S ribosomal RNA of rrnE operon
Reads supporting (aligned to +/- strand):  ref base T (8/13);  new base A (0/8);  total (8/21)
Fisher's exact test for biased strand distribution p-value = 6.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGT  >  NC_000913/4208083‑4208264
                                                                                                                                                   |                                  
aCCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgctt                                   <  2:1165710/149‑1 (MQ=32)
 ccAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgcttg                                  >  1:474323/1‑149 (MQ=32)
  cAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgcttgc                                 >  2:1159930/1‑149 (MQ=32)
    aGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAgcttgctt                               <  1:1534189/149‑1 (MQ=32)
     gTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATCAAAGTCGAACGGTAACAGGAAGAAGCTTGcttc                              >  1:1067460/1‑149 (MQ=17)
       ctcAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGcttctt                            <  2:1468073/149‑1 (MQ=32)
        tcAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGCTCATGGCGCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCttt                           <  2:323934/149‑1 (MQ=17)
        tcAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttc                           <  1:381907/149‑2 (MQ=17)
        tcAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTATGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttc                           <  1:1715782/149‑2 (MQ=12)
        tcAAGAGAGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGTATCAGATTGAACTCTGGCGGCAGGCCTAACACATGCAAGTCGAAAGGTAACAGGAAGAAGCTTGCTTCttt                           >  1:1125232/1‑149 (MQ=11)
         cAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTg                          >  2:1490378/1‑149 (MQ=32)
         cAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTg                          <  2:248247/149‑1 (MQ=32)
         cAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcc                          <  1:1544944/149‑3 (MQ=12)
          aaGTGTGAACACGTAATTCATTACGAAGTTTAATTATTTGAGCATCAAAGTTTTAAATTGAAGAGTTTTATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcct                         <  2:574416/149‑4 (MQ=11)
          aaGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcct                         <  1:395179/149‑4 (MQ=12)
          aaGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcct                         <  1:649907/149‑4 (MQ=12)
          aaGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcct                         <  1:946710/149‑4 (MQ=12)
           agagTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAACTTGAAGAGCTTGATCATGGCCCAGATTGAACGCGGGCGGCAGGCCTAACACATGCAAGTCGAACGGGACCCGGAAGAAGCTTGCTTCTTTGCt                        <  2:557626/149‑1 (MQ=11)
           agagTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCAAGcttcttcctg                        <  1:1265318/149‑5 (MQ=12)
                    acGTAATTCATTTCGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGACGAGGTTGATCATGGCTCAGATTGCACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTgg               <  2:555787/149‑1 (MQ=255)
                    acGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTgg               >  2:500914/1‑149 (MQ=32)
                    acGTAATTCATTACGAAGTTTAATTCTTTGAGAATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGAGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAATCTAGCTTCTTTGCTTACGAGTgg               >  2:1700774/1‑149 (MQ=11)
                     cGTAATTCATTACGAAGATTAATTCTTTGAGAATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGc              <  2:1297758/149‑1 (MQ=17)
                        atttCATTACGAAGTTTATTTCTTTGAGCATCTAACTTTTAAATTGTAGAGTTTGATCATGGCGCAGATTGACCGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGa           <  1:303570/147‑1 (MQ=11)
                        aaTTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGa           <  1:330348/149‑1 (MQ=32)
                           tCATTACGAAGTTTACTTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTATGCTGACGAGTGGCGGACgg        <  2:1737861/149‑1 (MQ=17)
                           tCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACgg        <  2:38966/149‑1 (MQ=32)
                           tCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACgg        <  2:94038/149‑1 (MQ=32)
                            cATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACggg       >  2:477655/1‑149 (MQ=32)
                               tACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGa    >  2:501940/1‑149 (MQ=32)
                               tACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGGGGCGGACGGGTGa    <  2:316655/149‑1 (MQ=17)
                                 cGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGt  <  1:1365794/149‑1 (MQ=31)
                                 cGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGt  >  2:80785/1‑149 (MQ=31)
                                                                                                                                                   |                                  
ACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGT  >  NC_000913/4208083‑4208264

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