Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,726,6520AG44.4% ‑4.9 / 44.9 36noncoding (2533/2904 nt)rrlG23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (2/18);  new base G (3/13);  total (5/31)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTC  >  NC_000913/2726522‑2726775
                                                                                                                                  |                                                                                                                           
ctctactactacTAGCAGCCCCCCTCATTTCTAAATCTCACACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGcccc                                                                                                                   >  2:1868191/7‑141 (MQ=1)
ctctactACTAGTAGAAGCCCCCCTCATTTCTCAATCGCACACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGcccc                                                                                                                   >  2:1868165/7‑141 (MQ=1)
                       ctccggtctcCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:18487/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:1308430/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:1110429/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:110255/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:1308465/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:1446719/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:1814005/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:223557/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:223567/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:223580/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:582778/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCAACATCGa                                                                                            <  2:1409605/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCAACATCGa                                                                                            <  2:1409287/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:373409/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:373410/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:7909/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:1225364/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:963533/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:950811/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:920696/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:680607/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:613972/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  2:1128339/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:1677585/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:16411/135‑1 (MQ=1)
                       ctcaaatctcAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACACTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGa                                                                                            <  1:613993/135‑1 (MQ=1)
                         cacttCTCCAACGCCCACGTCAGTTAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCCGGATGTGATGAGCCGACATCGAgg                                                                                          <  1:903821/138‑1 (MQ=1)
                         caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg                                                                                          <  1:603097/138‑1 (MQ=0)
                         caattCTCCAACGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGACTTCAGCCCCAGGATGTGATGAGCCGACATCGAgg                                                                                          <  1:30775/138‑1 (MQ=0)
                                caaacgtccACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCaaa                                                                                   >  1:276405/8‑141 (MQ=0)
                                caaacgtccACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCaaa                                                                                   >  1:276360/8‑141 (MQ=0)
                                caaacgtccACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCGGCTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCaaa                                                                                   >  1:1093614/8‑141 (MQ=0)
                                                                                                  ttcaaTGGCGACCAGCCATACCCTCGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACAACGAGGAGCCAAACACCGCCGACGATATCAACCCTCGGGCGGTATCAGCCTGATAACCCCGGAGTACCGTGCATCCGtc                 <  2:1944328/138‑2 (MQ=1)
                                                                                                                 ccttACCGTTGGTTCCGACTTCAGCCCCAGGATGTGTTTAGGCGACATCGAGTTGCCAAACTCCGCCGTCGATTTGAACTCTTGGGCGGTATCAGCCTGTTATCACCGGAGTTCCTTTTTTCGGTTGAGCGTTGGCGCTTc  <  2:86416/138‑1 (MQ=1)
                                                                                                                                  |                                                                                                                           
CTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCCTTGGGACCTACTTCAGCCCCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTC  >  NC_000913/2726522‑2726775

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