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

CCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556069‑3556345
                                                                                                                                            |                                                                                                                                          
ccccGGCAGTGCCGCTTTTTTCCCCGACGAAAAAAAAGCGTTCGGTGATTTTTTCTCTTTCGACTTAAAGCGAGACGGTATTCCCCGCCCCCT‑GGTCGCGGGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAA‑CCTc                                                                                                                                           <  2:1009055/137‑1 (MQ=255)
    ggCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCg                                                                                                                                       >  1:755931/1‑140 (MQ=255)
                               aaaaaaGTGTTGGGTTATTTTTTCACTTTCGACTTCAAGCGAGAGGGTATTCCCCGCCCCCT‑GGTCGCGCGGCGGGTTATGAATATTTTTTTCATGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCaa                                                                                                            <  2:1354264/140‑1 (MQ=255)
                                    aGCGTTCGGTGATTTTTTCCCTTTCGACTTAAAGCGAGACGGTATTCCCCGGCCCCT‑GGCCGCGGGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                       <  2:982336/140‑1 (MQ=255)
                                        ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                   >  1:1505348/1‑140 (MQ=255)
                                        ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                   >  2:1306991/1‑140 (MQ=255)
                                           ggTGATTTTTTCACTTTCGACTTAAAGCGAGACGTTATTCCCCGCCCCCG‑GGTCGCGGGGCGGGTTATAAATTTTCTGTTCATGAAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                <  2:674993/140‑1 (MQ=255)
                                             tGATTTTTTCACTTTCGCCTTCAGGGGAGACGTTATTCCCCGGCCCCG‑GGTCGCGCGGCAGGTTATGAATTTTCTGTTCAGGAAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcaca                                                                                              <  2:729488/140‑1 (MQ=255)
                                                ttttttCCCTTTCGACTTAAAGCGAGACGGTATTCCCCGGCCCCG‑GGCCGCGCGGCAGTTTATGAATATTCTGTTCAGGAAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:1505348/138‑1 (MQ=255)
                                                ttttttCATTTTCGTCTACAAG‑GGGGCGGTTTTCCCCCGCCCCTGGGTCGCGCGGCAGGTTATAATTTTTCTGTTCATGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:1067174/138‑1 (MQ=255)
                                                         tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                  >  1:1751097/1‑140 (MQ=255)
                                                          ttcaatttCAAGCGAGACGTTTTTGCCCGGCCCCG‑GGCCGCGCGGGATTTTATGAATTTTCTGTTCAGGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                 <  2:681101/134‑1 (MQ=255)
                                                           tCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                >  2:115341/1‑140 (MQ=255)
                                                                  cAAGCGAGACGTTATTCCCCGCCCCCT‑GGCCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGGAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                         <  2:1499290/140‑1 (MQ=255)
                                                                       gagaCGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAA‑CCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  1:1478340/140‑1 (MQ=255)
                                                                         aaCGATATTCCCCGCCCCCT‑GGGCGCCCGGCAGTTTATGAATATTCTGTTCA‑GGAGGGAAAACCCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                  <  2:841542/139‑1 (MQ=255)
                                                                          cgggTATCCCCCGCCCCCT‑GGGCGCGCGGGAGGTTATAAATATTTTTTTCATGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACACGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                 <  2:419897/138‑1 (MQ=255)
                                                                          aCGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                 >  1:837939/1‑140 (MQ=255)
                                                                                      ggCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                     >  2:1033211/1‑140 (MQ=255)
                                                                                      ggCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                     >  1:1680885/1‑140 (MQ=255)
                                                                                      ggCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                     >  2:205509/1‑140 (MQ=255)
                                                                                        ccccT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:583034/1‑140 (MQ=255)
                                                                                        ccccT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  1:1278841/1‑140 (MQ=255)
                                                                                        ccccT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  1:1148138/1‑140 (MQ=255)
                                                                                        ccccT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACACGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccctg                                                   >  1:419897/1‑138 (MQ=255)
                                                                                                gcgcgggGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                            <  1:1730947/139‑1 (MQ=255)
                                                                                                        cattttATAAATTTTCTGTTCATGAAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  2:837939/136‑1 (MQ=255)
                                                                                                        agggTTATGAATATTCTGTCCATGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  1:1426734/138‑1 (MQ=255)
                                                                                                         aGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTACACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   >  2:823547/1‑140 (MQ=255)
                                                                                                          gtttATGAATTTTCTGTTCATGCAGGAAAAA‑CCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                  <  2:1278841/138‑1 (MQ=255)
                                                                                                                     ttCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  1:720751/1‑140 (MQ=255)
                                                                                                                      tCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  2:791786/1‑140 (MQ=255)
                                                                                                                          ttCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  1:1125011/1‑140 (MQ=255)
                                                                                                                          ttCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  2:559560/1‑140 (MQ=255)
                                                                                                                          ttCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  2:234851/1‑140 (MQ=255)
                                                                                                                          ttCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  1:1442887/1‑140 (MQ=255)
                                                                                                                           tCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAc                 >  2:1349195/1‑140 (MQ=255)
                                                                                                                             aTGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg               >  2:1484331/1‑140 (MQ=255)
                                                                                                                             aTGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg               >  1:735663/1‑140 (MQ=255)
                                                                                                                                    gAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc        >  1:782225/1‑140 (MQ=255)
                                                                                                                                        a‑aCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc    >  1:1078018/1‑140 (MQ=255)
                                                                                                                                           cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  <  1:583034/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCT‑GGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑ACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556069‑3556345

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