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.

ACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATG  >  NC_000913/4166613‑4166898
                                                                                                                                                                                                                                              |                                                  
acacGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑Aggggg                                      >  1:71252/1‑250 (MQ=11)
  acGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGAcc                                   >  1:337592/1‑251 (MQ=31)
    gTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGATTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCtt                                 <  1:533607/251‑1 (MQ=14)
      aaTTCATTACGCAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGAAACTGCCTGATGGGGGGGGAT‑AAATACTGGAAACGGGGGCTAATACCGCA‑TAACGTGGCAAAACCAAAG‑Aggggggccgtcg                               >  2:760757/1‑244 (MQ=11)
      aaTTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAAAGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑ACCTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGCCCTTCg                               >  1:302765/1‑251 (MQ=14)
       aTTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGGTGGAGGGGGAT‑AACTACTGGAAAAGGTAGCTAATACCGCA‑TAACGTCGCAAAACCCAAG‑ATGGGGACCTTCgg                              >  2:370540/1‑251 (MQ=9)
         tCATTACGCAGTTTAATTCTTTGAGCGGCAACCTTTTAAGTTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTCACAGGAAGAAGCTTGCTTCTTTGCTGACGAGAGGAGGGCGGGGGGGGGAAGTCTGGGAAACTGCCTGGGGGGAGGGGAT‑CACTGCAGGAAACGGTGGGTAATACCGCC‑TCACGTCTCAAGACAAAAG‑AGGGGGACCTTCgggg                            >  2:443289/1‑250 (MQ=2)
         tCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGAACAAAG‑AGGGGGACCTTCGGGc                            >  1:203522/1‑251 (MQ=14)
           aTTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGTCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGGTGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCct                          >  1:573927/1‑251 (MQ=14)
           aTTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCct                          <  1:790902/251‑1 (MQ=31)
            agaCGAAGTTGAATTCTTGTAGCGTCAAACTTTTAAATTGAATAGTTTGAGCATGGCTCAGAATGAACGATGGCGGCAGGCCTTACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTTTTTGCTGACGAGTTGCTGACGGGGGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAACACCGCA‑TAACGTCGCAAGACCAAAG‑ATGTGGACCTTCGGGCctc                         <  2:476564/249‑1 (MQ=11)
               cGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTg                      <  2:774572/251‑1 (MQ=31)
               cGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGATGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGAAGGGAGAGGAATGTCTGGGAAAATGCCTGATGGAGGGGGAG‑AACTAATGGACACGGTAGCTAATAACGCA‑TAACGTCGCAAGACCAAAGCAGGGGGACATTGGGGcccttg                       >  2:464406/1‑247 (MQ=11)
                gAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAAACCAAAG‑AGGGGGACCTTCGGGCCTCTTGc                     <  2:677787/251‑1 (MQ=14)
                gAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCGCAGATTGAACGCTGGCGGCAGGCCTTACACATGCAAGTTGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGTCGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGTGCCTCTTGc                     <  2:299598/251‑1 (MQ=9)
                 aaGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGCGTTTGATCATGGCTCAGATTGAACGCTGCCGACAGGACTAAAACATGCAAGGCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGGATGGCGGACGGGTGAGTAATGCTTGGGAAAACACCCGATGGTGGGGGATGAAC‑ACTGGAAACGGTCGCTAA‑ACCACACAAACGGAGCCAGACAAAAG‑AGAGAGACCTTCGGGCCCCTTGcc                    >  2:374283/1‑251 (MQ=2)
                  aGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGGGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAAAGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGcc                    >  1:806176/1‑250 (MQ=14)
                   gTTTAATTCGTTGAGCGGCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCAt                  <  2:834238/251‑1 (MQ=14)
                     ttAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATc                 <  1:842110/250‑1 (MQ=31)
                     ttAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTAAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGACGAAGCTTGCTTCTTTGCTGACGCGTGGCGGACGGGGGAGTAATATCTGGGAAACTGCCTGCTGGCGGGGGAT‑AACTACTGGAAACGGTAACTAATACCGAA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCAGGCCTCTTGCCATCg                >  1:112315/1‑251 (MQ=11)
                       aaTTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCgg               >  1:648660/1‑250 (MQ=31)
                        aTTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGATCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAA‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCGTGCCATCGGAt             >  1:871615/1‑251 (MQ=9)
                          tCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAAACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCGGAtgt           >  1:478530/1‑251 (MQ=14)
                            tttGAGCGTCAAACTTTTAAATTGAGGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTACCAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGGCGGGGGAGAAATGTCTGGGAAACTGCCCGATGGAGGGGGAT‑ACCAACTGGAAACGGGAGCTAATACCGTG‑TCACGACGCAAAACCCAAG‑AAGGGGGCATTCGGGCCTCTTGCCCTCGGAtgtgg         >  2:136750/1‑250 (MQ=11)
                            tttGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAAACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCGGGTGTGc         >  1:286419/1‑251 (MQ=14)
                             ttGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGGTTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGCAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGcc        >  1:343176/1‑251 (MQ=14)
                             ttGAGCGTCAAACTTTTAAAGTGAAGAGTTTGATCATGGCGCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACGGCCTTATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACGGCA‑GAACTTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCGGAtgtg          <  1:33286/249‑1 (MQ=11)
                               gAGCGTCAAACTTTTAAATTGAATAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTTACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACTGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCTGGCCTCTTGCCATCGTATGTGCCCa      <  2:857858/251‑1 (MQ=9)
                                  cGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGccc       >  1:227898/1‑247 (MQ=11)
                                  cGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGAAACAGGAAGCAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGGGGGGGAT‑AACTACTGGAAACGGTTGCTAATACCGCA‑TAACGTCGCAAGACCACAG‑AGGGGGACCTACTGGCCTCTTGCCATCGGATGTGCCCAg     >  1:844316/1‑249 (MQ=11)
                                   gTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACGGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAAACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATg  <  1:486412/251‑1 (MQ=11)
                                                                                                                                                                                                                                              |                                                  
ACACGTAATTCATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAT‑AACTACTGGAAACGGTAGCTAATACCGCA‑TAACGTCGCAAGACCAAAG‑AGGGGGACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATG  >  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: ‑