Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2070TC37.0% 30.0 / 14.0 27S215G (AGT→GGT) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (16/1);  new base C (0/10);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 1.30e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556068‑3556343
                                                                                                                                            |                                                                                                                                        
accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑Ct                                                                                                                                          >  1:157170/1‑140 (MQ=255)
   aacgCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTccc                                                                                                                                       >  1:506248/4‑140 (MQ=255)
        aCTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCa                                                                                                                                  >  1:722697/1‑140 (MQ=255)
                 cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGa                                                                                                                         >  1:459944/1‑140 (MQ=255)
                      ccccGACAAAAAAAAAGCGTTCGGTATTTTTTTCCCTTTCGACTTAAAGCGAGCGGTTTTTCCCCGGCCCCGGGTCGCGGGGCGGTTTATGAATTTTCTGTTCA‑GGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCa                                                                                                                    <  2:852769/138‑1 (MQ=255)
                        ccGACGACAAAAAAGCGTGCGGTAATTTTTTCCCTTTCGCCTTCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCGGTTTATAAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACg                                                                                                                  <  2:933958/140‑1 (MQ=255)
                          gacaaaaaaaagGCGTTCGGTGATTTTTTCATTTTCGACTTCAGGCGAGAGGGTATTCCCCGCCCCCTGGTCGCGGGGCAGTTTAGAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                <  2:506248/134‑1 (MQ=255)
                            cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                              >  1:1415518/1‑140 (MQ=255)
                            cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                              >  2:180961/1‑140 (MQ=255)
                                     ggCGTTGGTTTATATTTTCCCTTTCCCCTTAAAGCGAGACGGTATTACCCGGCCCCGGGCGGCGCGGCGGGTTATGAATTTTCTTTTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                     <  2:723383/139‑1 (MQ=255)
                                        tttCGGTTATTTTTTCTTTTCCGACTTAAAGCGAGACGTTATTCCCCGCCCCCTGGTCGCGCGGCGGTTTATGAATTTTCTTTTCAGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                  <  2:722697/139‑1 (MQ=255)
                                         ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                 >  1:1041261/1‑140 (MQ=255)
                                                      tccttttCGACTTCAAGCGAGCCGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATTTTCTTTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                    <  2:861095/136‑1 (MQ=255)
                                                              gACTTAAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                            <  2:157170/140‑1 (MQ=255)
                                                                      gCGAGAGGGTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                    <  2:459944/140‑1 (MQ=255)
                                                                             ggTATTCCCCGCCCCCTGGCCGCGCGGCAGGTTATGATTTTTCTGTTCA‑GGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                             <  2:1322143/140‑1 (MQ=255)
                                                                                         ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                 >  2:48414/1‑140 (MQ=255)
                                                                                         ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                 >  1:1114798/1‑140 (MQ=255)
                                                                                         ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGC‑CAATGTTCCcgcg                                                >  1:1258594/1‑140 (MQ=255)
                                                                                                   cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                       >  1:61125/1‑140 (MQ=255)
                                                                                                   cgcgGCAGGTTATAAATATTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCACCAAg                                       <  1:180961/140‑1 (MQ=255)
                                                                                                          ggTTATGAATTTTCTGTTAAGGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                <  2:61125/140‑1 (MQ=255)
                                                                                                                      tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGac                    >  1:489335/1‑139 (MQ=255)
                                                                                                                      tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                    >  1:225254/1‑140 (MQ=255)
                                                                                                                               gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg              >  2:1500289/1‑137 (MQ=255)
                                                                                                                               gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg              <  1:1500289/137‑1 (MQ=255)
                                                                                                                               gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc           >  2:167873/1‑140 (MQ=255)
                                                                                                                               gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc           >  1:624305/1‑140 (MQ=255)
                                                                                                                                   gaaaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgaga       <  2:671009/138‑1 (MQ=255)
                                                                                                                                        aa‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc  >  2:1072635/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                        
ACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556068‑3556343

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: ‑