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

GGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556073‑3556345
                                                                                                                                      |                                                                                                                                          
ggCACTGCCGCGCTTTTCACCGACGAAAAAAAGGCGTTCGGTGATTTTTCCCCTTTCCCTTCCAAGCGAGACGGTATTCCCCGCCCCCGGCTGGCGCGGCAGTTTATGAATTTTCTGTGCATGGAGGGAAAACCCCCCCg                                                                                                                                       <  2:1206970/140‑1 (MQ=255)
                   ccGAAGAAAAAAAAGGTTCCGGTATTTTTTTCTCTTCCGTCTAAAGGGGAGCGGGTTTTCCCGGCCCCCTGGCCGCGGGGCGGTTTAAAAATATTGTTTTATTGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACg                                                                                                                    <  2:181052/140‑1 (MQ=255)
                     gacgaaAAAAAAGCGTTGGGTGATTTTTTCACTTTCGACTTCAAGCGAGAGGTTATTCCCCGGCCCCGGGCCGCGGGGCGGTTTATAAATTTTCTGTTCAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                  <  2:769520/140‑1 (MQ=255)
                      acgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                 >  1:181052/1‑140 (MQ=255)
                       cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                >  1:544815/1‑140 (MQ=255)
                         acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                              >  1:610391/1‑140 (MQ=255)
                                ggcgtggggTGTTTTTTTCACTTCCGCCTTCAGGGGAGAGGTTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                       <  2:77757/133‑1 (MQ=255)
                                ggCGTTGGGTTATTTTTTCTCTTTCGACTTAAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                       <  2:1269682/139‑1 (MQ=255)
                                  gcttCGGTTATTTTTTCTCTTCCGCCTTAAACCGAGAGGTTATTCCCCGCCGCCTGGTCGCGCGCCAGTTTAAGAATTTTTTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                                     <  2:863307/138‑1 (MQ=255)
                                       ggTGTTTTTTTCACTTTCGACTCCAAGCGAGGCGTTTTTCCCCGGCCCCTGGCGGCGCGGCAGTTTATAAATTTTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                <  2:108251/140‑1 (MQ=255)
                                       ggTGATTTTTTCCCTTTCGACTTAAGGCGAGACGGTTTTCCCCGGCCCCGGGCCGCGGGGCAGTTTATAAATATTCTTTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                <  2:1065045/140‑1 (MQ=255)
                                       ggTGATTTTTTCACTTTCGACTTAAAGCGAGCCGGTATTCCCCGGCCCTTGGCCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                <  2:108440/140‑1 (MQ=255)
                                          aagtttttCACTTTCGACTTCAGGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                             <  2:233000/136‑1 (MQ=255)
                                                 tCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                      >  2:738396/1‑140 (MQ=255)
                                                        cGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGactg                                                                               >  1:233000/1‑137 (MQ=255)
                                                          aCTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                             >  1:138796/1‑140 (MQ=255)
                                                                 gCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATACTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                      >  2:776969/1‑140 (MQ=255)
                                                                      agggTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                 <  1:1335556/138‑1 (MQ=255)
                                                                        ggTATTCCCCGGCCCCTGGCCGCGGGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGca                                                               <  2:610391/140‑1 (MQ=255)
                                                                                  ggCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                     >  2:824679/1‑140 (MQ=255)
                                                                                    ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:245617/1‑140 (MQ=255)
                                                                                    ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:1038721/1‑140 (MQ=255)
                                                                                    ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:1229851/1‑140 (MQ=255)
                                                                                    ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTAATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:884737/1‑140 (MQ=255)
                                                                                           tCGCGCGGCAGTTTATAAATATTCTGTCCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                            <  1:107104/140‑1 (MQ=255)
                                                                                             gcgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCaa                                          >  2:1367744/1‑140 (MQ=255)
                                                                                                    gggTTATAAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   <  2:138796/139‑1 (MQ=255)
                                                                                                                ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  1:1018828/1‑140 (MQ=255)
                                                                                                                ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       <  1:245617/140‑1 (MQ=255)
                                                                                                                ttCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  2:46951/1‑140 (MQ=255)
                                                                                                                 tCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  2:1295663/1‑140 (MQ=255)
                                                                                                                    tttCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                   <  2:354358/139‑1 (MQ=255)
                                                                                                                    gTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                   >  2:584311/1‑140 (MQ=255)
                                                                                                                     ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  1:14699/1‑140 (MQ=255)
                                                                                                                     ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  1:1423261/1‑140 (MQ=255)
                                                                                                                     ttCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  1:1213054/1‑140 (MQ=255)
                                                                                                                     ttCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  <  1:1229851/140‑1 (MQ=255)
                                                                                                                      tCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc                 >  2:637177/1‑140 (MQ=255)
                                                                                                                          gCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             >  2:1090579/1‑140 (MQ=255)
                                                                                                                            aGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCga           <  2:1213054/140‑1 (MQ=255)
                                                                                                                             gggAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag          >  1:354358/1‑140 (MQ=255)
                                                                                                                                    aCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg   <  1:1367744/140‑1 (MQ=255)
                                                                                                                                     cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  <  2:14699/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                          
GGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556073‑3556345

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