Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,166,8490GA22.6% 27.6 / 10.3 31noncoding (191/1542 nt)rrsB16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (15/9);  new base A (5/2);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATG  >  NC_000913/4166613‑4166898
                                                                                                                                                                                                                                              |                                                  
acacGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAggggg                                      >  1:71252/1‑250 (MQ=11)
  acGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGAcc                                   >  1:337592/1‑251 (MQ=31)
    gTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGATTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCtt                                 <  1:533607/251‑1 (MQ=14)
      aaTTCATTACGCAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGAAACTGCCTGATGGGGGGGGATAAATACTGGAAACGGGGGCTAATACCGCATAACGTGGCAAAACCAAAGAggggggccgtcg                               >  2:760757/1‑244 (MQ=11)
      aaTTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAAAGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATACCTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGCCCTTCg                               >  1:302765/1‑251 (MQ=14)
       aTTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGGTGGAGGGGGATAACTACTGGAAAAGGTAGCTAATACCGCATAACGTCGCAAAACCCAAGATGGGGACCTTCgg                              >  2:370540/1‑251 (MQ=9)
         tCATTACGCAGTTTAATTCTTTGAGCGGCAACCTTTTAAGTTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTCACAGGAAGAAGCTTGCTTCTTTGCTGACGAGAGGAGGGCGGGGGGGGGAAGTCTGGGAAACTGCCTGGGGGGAGGGGATCACTGCAGGAAACGGTGGGTAATACCGCCTCACGTCTCAAGACAAAAGAGGGGGACCTTCgggg                            >  2:443289/1‑250 (MQ=2)
         tCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGAACAAAGAGGGGGACCTTCGGGc                            >  1:203522/1‑251 (MQ=14)
           aTTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGTCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGGTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCct                          >  1:573927/1‑251 (MQ=14)
           aTTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCct                          <  1:790902/251‑1 (MQ=31)
            agaCGAAGTTGAATTCTTGTAGCGTCAAACTTTTAAATTGAATAGTTTGAGCATGGCTCAGAATGAACGATGGCGGCAGGCCTTACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTTTTTGCTGACGAGTTGCTGACGGGGGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAACACCGCATAACGTCGCAAGACCAAAGATGTGGACCTTCGGGCctc                         <  2:476564/249‑1 (MQ=11)
               cGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTg                      <  2:774572/251‑1 (MQ=31)
               cGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGATGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGAAGGGAGAGGAATGTCTGGGAAAATGCCTGATGGAGGGGGAGAACTAATGGACACGGTAGCTAATAACGCATAACGTCGCAAGACCAAAGCAGGGGGACATTGGGGcccttg                       >  2:464406/1‑247 (MQ=11)
                gAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAAACCAAAGAGGGGGACCTTCGGGCCTCTTGc                     <  2:677787/251‑1 (MQ=14)
                gAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCGCAGATTGAACGCTGGCGGCAGGCCTTACACATGCAAGTTGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGTCGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGTGCCTCTTGc                     <  2:299598/251‑1 (MQ=9)
                 aaGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGCGTTTGATCATGGCTCAGATTGAACGCTGCCGACAGGACTAAAACATGCAAGGCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGGATGGCGGACGGGTGAGTAATGCTTGGGAAAACACCCGATGGTGGGGGATGAACACTGGAAACGGTCGCTAAACCACACAAACGGAGCCAGACAAAAGAGAGAGACCTTCGGGCCCCTTGcc                    >  2:374283/1‑251 (MQ=2)
                  aGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGGGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAAAGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGcc                    >  1:806176/1‑250 (MQ=14)
                   gTTTAATTCGTTGAGCGGCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCAt                  <  2:834238/251‑1 (MQ=14)
                     ttAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATc                 <  1:842110/250‑1 (MQ=31)
                     ttAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTAAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGACGAAGCTTGCTTCTTTGCTGACGCGTGGCGGACGGGGGAGTAATATCTGGGAAACTGCCTGCTGGCGGGGGATAACTACTGGAAACGGTAACTAATACCGAATAACGTCGCAAGACCAAAGAGGGGGACCTTCAGGCCTCTTGCCATCg                >  1:112315/1‑251 (MQ=11)
                       aaTTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCgg               >  1:648660/1‑250 (MQ=31)
                        aTTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGATCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAAAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCGTGCCATCGGAt             >  1:871615/1‑251 (MQ=9)
                          tCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAAACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGAtgt           >  1:478530/1‑251 (MQ=14)
                            tttGAGCGTCAAACTTTTAAATTGAGGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTACCAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGGCGGGGGAGAAATGTCTGGGAAACTGCCCGATGGAGGGGGATACCAACTGGAAACGGGAGCTAATACCGTGTCACGACGCAAAACCCAAGAAGGGGGCATTCGGGCCTCTTGCCCTCGGAtgtgg         >  2:136750/1‑250 (MQ=11)
                            tttGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAAACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGGTGTGc         >  1:286419/1‑251 (MQ=14)
                             ttGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGGTTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGCAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGcc        >  1:343176/1‑251 (MQ=14)
                             ttGAGCGTCAAACTTTTAAAGTGAAGAGTTTGATCATGGCGCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACGGCCTTATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACGGCAGAACTTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGAtgtg          <  1:33286/249‑1 (MQ=11)
                               gAGCGTCAAACTTTTAAATTGAATAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTTACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACTGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCTGGCCTCTTGCCATCGTATGTGCCCa      <  2:857858/251‑1 (MQ=9)
                                  cGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGccc       >  1:227898/1‑247 (MQ=11)
                                  cGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGAAACAGGAAGCAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGGGGGGGATAACTACTGGAAACGGTTGCTAATACCGCATAACGTCGCAAGACCACAGAGGGGGACCTACTGGCCTCTTGCCATCGGATGTGCCCAg     >  1:844316/1‑249 (MQ=11)
                                   gTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACGGGAAACGGTAGCTAATACCGCATAACGTCGCAAAACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATg  <  1:486412/251‑1 (MQ=11)
                                                                                                                                                                                                                                              |                                                  
ACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATG  >  NC_000913/4166613‑4166898

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