Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2070TC41.7% 33.7 / 17.5 36S215G (AGT→GGT) ‡rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (19/2);  new base C (0/15);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 2.44e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGA  >  NC_000913/3556074‑3556338
                                                                                                                                     |                                                                                                                                   
gCACTGCCGCCTTTTTCCCCGACGAAAAAAAAGCGTTGGGTGATTTTTTCACTTTCGACTTCAGGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTCATGGAGGAAAAACCCCCCCGc                                                                                                                               <  2:441946/140‑1 (MQ=255)
gCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGc                                                                                                                               >  1:1388858/1‑140 (MQ=255)
                ccccgccgACAAAAAGGCGTTCGGTGATTTTTTCCTTTTCGACTCCAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATATTCTTTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCa                                                                                                               <  2:692494/138‑1 (MQ=255)
                  ccGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATTTTCTGTTCATGCAGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACg                                                                                                             <  1:644878/140‑1 (MQ=255)
                  ccGACGAAAAAAAAGCGTTGGGTGATTTTTTCCCTTTCGACTTAAAGCGAGACGGTATTCCCCGCCCCCTGGTCGCGGGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACg                                                                                                             <  2:1203449/140‑1 (MQ=255)
                            aaagGCGTTGGTTGATTTTTCCCTTTTCGCCTCCAGGCGAGCCGTTTTTCCCCGCCCCCTGGTCGGGGGGCGGGTTATAAATTTTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                   <  2:794020/140‑1 (MQ=255)
                                 cGTTGGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                              <  2:962544/140‑1 (MQ=255)
                                   ttCGGTGATTTTTTCTTTTTCGCTTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGGGGCGGGTTATAAATTTTTTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                            <  2:529753/140‑1 (MQ=255)
                                   ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                            >  1:529753/1‑140 (MQ=255)
                                        tGATTTTTTCACTTCCGACTTCAAGCGAGACGGTATTCCCCGCCCCCTGGCCGCGCGGCGGTTTATAAATTTTTTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca                                                                                       <  2:439603/140‑1 (MQ=255)
                                        tGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca                                                                                       >  2:1418706/1‑140 (MQ=255)
                                          atatTTCCCCTTTCGACTTAAAGCGAGACGGTATTCCCCGCCCCCTGGGCGGGCGGCGGGTTATAAATTTTCTGTACATGCAGAAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                     <  2:186494/140‑1 (MQ=255)
                                              tttCTCTTTCGATTTCAAGGGAGACGTTTTTCCCCGGCCCCTGGCCGGGGGGCAGGTTATGAATATTTTGTTAATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                 <  2:1211/140‑1 (MQ=255)
                                                tcccTTTCGATTTCAGGCGAGACGGTATTCCCCGCCCCCTGGTCGCGGGGCAGTTTATAAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                               <  2:1256497/137‑1 (MQ=255)
                                                      tCGCCTTCAAGCGAGACGGTATCCCCCGGCCCCTGGTCGGGGGGCAGTTTATGAATATTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                         <  2:73349/140‑1 (MQ=255)
                                                             aaaGGGAGACGTTTTTCCCCGCCCCCTGGCCGCGCGGCAGTTTATAAATTTTCTTTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                  <  2:1388858/139‑1 (MQ=255)
                                                                gCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                               >  2:80729/1‑140 (MQ=255)
                                                                gCGAGACGGTATTACCCGGCCCCTGGGCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                               >  2:1463074/1‑140 (MQ=255)
                                                                     cggttATTCCCCGCCCCCGGGCCGCGGGGCAGTTTATAAATTTTCTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                          <  2:41005/136‑1 (MQ=255)
                                                                     aCGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                          <  1:1235649/140‑1 (MQ=255)
                                                                        gTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                       >  1:681226/1‑140 (MQ=255)
                                                                                  gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                             >  1:1313033/1‑140 (MQ=255)
                                                                                   ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                            >  2:1031879/1‑140 (MQ=255)
                                                                                   ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                            >  1:1306997/1‑140 (MQ=255)
                                                                                   ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                            >  2:468392/1‑140 (MQ=255)
                                                                                          tCGCGGGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                     <  2:681226/140‑1 (MQ=255)
                                                                                          tCGCGCGGCAGTTTATAAATATTCTGTAAAAGCGGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                     <  2:97133/140‑1 (MQ=255)
                                                                                          ccgcgcgGCAGTTTATAATTTTTCTGTTCAGGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                     <  2:713345/139‑1 (MQ=255)
                                                                                             cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                  >  1:97133/1‑140 (MQ=255)
                                                                                                   agtttATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                            <  2:855231/137‑1 (MQ=255)
                                                                                                   aGGTTATGAATATTCTTTCCATGCAGAAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                            <  2:1313033/140‑1 (MQ=255)
                                                                                                    gtttATAAATATTCTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                           <  1:361811/138‑1 (MQ=255)
                                                                                                        aTGAATTTTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt                       <  1:1314207/140‑1 (MQ=255)
                                                                                                              atTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                 >  1:593080/1‑140 (MQ=255)
                                                                                                                 cTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg              >  2:915427/1‑140 (MQ=255)
                                                                                                                  tGTTCATGCAGGGAAAAACTCCCTGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc             >  2:1460633/1‑140 (MQ=255)
                                                                                                                   tttCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc            <  2:83234/139‑1 (MQ=255)
                                                                                                                   tttCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc            <  1:1322594/139‑1 (MQ=255)
                                                                                                                   gTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc            >  1:415310/1‑140 (MQ=255)
                                                                                                                    ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa           >  2:1475271/1‑140 (MQ=255)
                                                                                                                    ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa           >  1:492348/1‑140 (MQ=255)
                                                                                                                     tCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc          >  2:887798/1‑140 (MQ=255)
                                                                                                                            gggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag   >  1:1089906/1‑140 (MQ=255)
                                                                                                                             aaaaaacCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgaga  <  2:976543/138‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
GCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGA  >  NC_000913/3556074‑3556338

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