Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2050AC27.6% 46.8 / 10.2 29S215R (AGT→AGG) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base A (21/0);  new base C (0/8);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 2.33e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAC  >  NC_000913/3556070‑3556339
                                                                                                                                       |                                                                                                                                       
cTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcc                                                                                                                                     >  1:7398/1‑140 (MQ=255)
   ggCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCg                                                                                                                                 >  1:1006262/1‑140 (MQ=255)
          aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAgtgt                                                                                                                          >  2:793499/1‑140 (MQ=255)
          aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTAATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAgtgt                                                                                                                          >  1:1007104/1‑140 (MQ=255)
               cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGa                                                                                                                     >  1:688556/1‑140 (MQ=255)
                    ccccGACGACAAAAAAGCGTTGGGTGATTTTTCCCCTTTCGACTTCAGGCGAGACGTTATCCCCCGCCCCCTGGCCGCGGGGCAGGTTAAAAATTTTTTTTTCAGGGGGGGAAAACCCCC‑CCGCCAGTGTAGCGATTTCa                                                                                                                <  2:173698/138‑1 (MQ=255)
                    cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCa                                                                                                                >  1:312144/1‑140 (MQ=255)
                      ccGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGCCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAACCCCC‑CCGCCAGTGTAGCGATTTCACg                                                                                                              <  2:1121131/140‑1 (MQ=255)
                        acagaaaaaaagGCTTTGGGTGATTTTTTCTCTTTCGACTAAAAGCGAGCCGTTATTCCCCGCCCCCGGGCCGCGGGGGGGTTTATGATTATTCTTTCCAGGAAGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCt                                                                                                            <  2:962801/134‑1 (MQ=255)
                          cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCa                                                                                                          >  1:841363/1‑140 (MQ=255)
                          cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCa                                                                                                          >  1:893491/1‑140 (MQ=255)
                                aaaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                    >  1:748496/1‑140 (MQ=255)
                                       ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgctgt                                                                                             >  1:962801/1‑137 (MQ=255)
                                        tgggTGTTTTTTTCATTTTCGCTTTAAAGCGAGCCGTTATTCCCCGGCCCTTGGTCGCGCGGGGGGTTATAAATTTTCTGTTCATGAAGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcga                                                                                            <  2:748496/138‑2 (MQ=255)
                                          ggTGATTTTTTCACTTTCGCCTTAAGGCGAGCGGTTTTTCCCCGCCCCCGGGTCGCGGGGCAGTTTATAAATATTTTTTTCATGCGGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                          <  2:7398/140‑1 (MQ=255)
                                               ttttttcttttcCGCCTTCAAGCGAGACGTTATCCCCCGGCCCCTGGCCGGGCGGCGGTTTAAAAATGTTCTTTAGGTGGAGGAAAAAACGCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                     <  2:1183427/138‑1 (MQ=255)
                                               ttttttCACTTCCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGGGGCAGGTTATGAATTTTCTGTTCAGGCAGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                     <  2:525163/138‑1 (MQ=255)
                                                    tcccTTTCGACTTCAAGCGAGACGTTATTCCCCGCCCCCTGGTCGCGGGGCGGGTTATAAATTTTTTGTTCAGGCAGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                <  2:819398/137‑1 (MQ=255)
                                                    cctcTTCCGCCTTCAGGGGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGAAGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                <  2:688556/137‑1 (MQ=255)
                                                       ttttCGACTTAAAGCGAGACGGTATTCCCCGCCCCCTGGCCGCGGGGCAGTTTAGAAATATTCTTTTCAGGCAGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                             <  2:812293/139‑1 (MQ=255)
                                                        tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                            >  2:235088/1‑140 (MQ=255)
                                                         ttCGCCTTCAAGGGAGACGGTTTTCCCCGCCCCCGGGCCGGGGGGCGGTTTATAAATTTTCTGTTCATGGGGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                           <  2:893491/140‑1 (MQ=255)
                                                                   ggCGAGGCGGTATTCCCCGCCCCCTGGCCGCGCGGCGGTTTATGAATTTTCTGTTCATGCAGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                 <  2:839433/139‑1 (MQ=255)
                                                                                  cccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                  >  1:22447/1‑140 (MQ=255)
                                                                                      gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                              >  1:798166/1‑140 (MQ=255)
                                                                                                      gggtttATAAATTTTCTGTTCAGGCAGGAAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                              <  2:119275/136‑1 (MQ=255)
                                                                                                      cggtttATAAATTTTTTGTTCAGGCAGGGAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                              <  1:646335/136‑1 (MQ=255)
                                                                                                      cggtttATAAATATTCTTTTCATGAAGGGAAAACCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                              <  1:793499/136‑1 (MQ=255)
                                                                                                             taaaTATTCTGTTCATGCAGGGAAAACCCCC‑CCGCCAGTGTAGCGATTTAACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                       <  1:831076/138‑1 (MQ=255)
                                                                                                                   ttCTGTTCATGCAGGGAAAACCTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                              <  1:649839/127‑1 (MQ=255)
                                                                                                                   ttCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                              >  2:649839/1‑127 (MQ=255)
                                                                                                                   ttCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAcct                 >  1:119275/1‑138 (MQ=255)
                                                                                                                    tCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                >  1:333037/1‑140 (MQ=255)
                                                                                                                     cTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg               >  2:1019102/1‑140 (MQ=255)
                                                                                                                        ttCATGCGGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa            >  1:1301443/1‑140 (MQ=255)
                                                                                                                         tCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc           >  1:205771/1‑140 (MQ=255)
                                                                                                                             tcAGGAAAAAACCCCGCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg        <  2:964631/139‑1 (MQ=255)
                                                                                                                               aGGGAAAATCCCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCGCGGGCga     <  1:644380/140‑1 (MQ=255)
                                                                                                                                gggAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag    >  2:408571/1‑140 (MQ=255)
                                                                                                                                  gAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAATGATGCCACGGGCGAGAc  >  2:802789/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCC‑CCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAC  >  NC_000913/3556070‑3556339

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