Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2070TC43.8% 62.5 / 13.9 48S215G (AGT→GGT) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (22/5);  new base C (0/21);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 3.22e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA  >  NC_000913/3556059‑3556355
                                                                                                                                                       |                                                                                                                                                    
gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcc                                                                                                                                                    <  1:725141/151‑1 (MQ=255)
  cGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTcccc                                                                                                                                                  <  2:21358/151‑1 (MQ=255)
        caccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgt                                                                                                                                            >  2:909987/1‑151 (MQ=255)
         cccgcGGCACTGACGCGTTTTCCACCGACGAAAAAAAAGCGTTG‑GGTGA‑TATTTTCTTTTTCGACTTAAAGGGAGACGTTATTCCCCGCCCCCGGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAgtg                                                                                                                                           <  2:291117/147‑1 (MQ=255)
                      cgcgTTTTCCCCCGACGAAAAAAAGGCGTTC‑GGTGA‑TTTTTCCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATTTTTTGTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTc                                                                                                                                 <  1:541084/148‑1 (MQ=255)
                      ccctcTTTCCCCCGACAAAAAAAAGGCGTCC‑GGTAA‑TTTTTTCCCTTCCGACTCCAAGCGAGAGGTTATCCCCCGCCCCCTGGCCCCGCGGCGGTTTATAAATTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTc                                                                                                                                 <  2:931736/144‑1 (MQ=255)
                       gctcTTTCCACGGACGAAAAAAAAGCG‑TC‑GGTTATTTTTTTCCCTTTCGACTAAAAGCGAGCCGTTATTCCCCGCCCCCTGGCCGCGGGGAAGTTTAGGAATATTCTTTTCAAGCAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTc                                                                                                                                 <  2:552723/144‑1 (MQ=255)
                       gcgcTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGc                                                                                                                             >  1:877924/1‑151 (MQ=255)
                         gcTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGAt                                                                                                                                    >  1:874032/1‑142 (MQ=255)
                              tCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCa                                                                                                                       >  1:700867/1‑150 (MQ=255)
                                  cgacgaCAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                  >  2:173524/1‑151 (MQ=255)
                                  cgacgaCAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                  >  1:570266/1‑151 (MQ=255)
                                     cgacAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                                               >  1:367739/1‑151 (MQ=255)
                                       acAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                             >  1:446185/1‑151 (MQ=255)
                                       aaaaaaaaGCGTTC‑GGTGA‑TTTTTCCCCTTTCGACTTAAAGCGAGACGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                              <  1:105634/148‑1 (MQ=255)
                                         aaaaaaGGGGTCGGGTGA‑TTTTTTCCTTTTCGACTAAAAGCGGGACTTTTTCCCCCGCCCCTTGGCCGCGCGGCGGTTTATGAATTTTCTGTTCAGGAAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                   <  2:553473/144‑1 (MQ=255)
                                          aaaaaGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                           >  1:916123/1‑150 (MQ=255)
                                                   tC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                  >  2:593195/1‑150 (MQ=255)
                                                       ggtattttttccccTTACGCCTTCAAGCGAGACGTTTTTCCCCGCCCCCGGGTCGCGGGGGAGTTTATGAATTTTTTTTTCATGGAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                  <  2:726285/135‑1 (MQ=255)
                                                            tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                          >  1:798944/1‑151 (MQ=255)
                                                                ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                                      >  2:1039024/1‑151 (MQ=255)
                                                                     ttccGACTTCAGGGGAGACGGTTTTCCCCGGCCCCAGGCCGCGGGGCAGTTAAGGAATTTTCTGTTCAGGCGGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                    <  2:567272/145‑1 (MQ=255)
                                                                      tcggcctaaaggggagaGGGTATCCCCCGGCCCCGGGGCCGGCGGAAGTTTATGATTTTTTTGTTCAGGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                 <  1:155110/137‑1 (MQ=255)
                                                                            ttaaaGCGAGAGGTTATACCCCGCCCCCGGGCCGCGCGGAGGGTTATAAATTTTCTGTTGATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                           <  1:910028/147‑1 (MQ=255)
                                                                               aaGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                        >  2:932443/1‑150 (MQ=255)
                                                                                   gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                   >  2:316510/1‑151 (MQ=255)
                                                                                     gaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCGCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTc                                                                 <  1:158158/151‑1 (MQ=255)
                                                                                      acgttatccccgggCCCCGGGGCGCGCGCCTGTTTATAAATTTTTTTTTCTTGCAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                   <  1:587039/136‑1 (MQ=255)
                                                                                        gttATTCCCGGCCCCCGGGTCGCGGGGCAGTTTAGGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                                  <  2:1013311/145‑1 (MQ=255)
                                                                                            tcccccGCCCCCGGGTCGCGGGGCAGGTTATAAATATTCTGTTCATG‑AGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg                                                            <  1:603725/146‑1 (MQ=255)
                                                                                             tACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                         >  2:878023/1‑151 (MQ=255)
                                                                                               cccGGCCCCTGGTCGCGCGGAAGTTTATGAATTTTCTGTCCAGGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                                       <  1:515249/151‑1 (MQ=255)
                                                                                                              cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt                                        >  1:894650/1‑151 (MQ=255)
                                                                                                               gcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                       >  2:134258/1‑151 (MQ=255)
                                                                                                                cgGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                       <  1:718035/150‑1 (MQ=255)
                                                                                                                cgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                       >  1:519077/1‑150 (MQ=255)
                                                                                                                 ggaagtttATAAATTTTCTGTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                     <  1:804101/145‑1 (MQ=255)
                                                                                                                 ggCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                     <  2:52851/151‑1 (MQ=255)
                                                                                                                 ggCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                                      >  1:122788/1‑150 (MQ=255)
                                                                                                                 gcgaGGTAATAAATTTTCTGTTCGGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                             <  1:1039073/140‑1 (MQ=255)
                                                                                                                            aaTATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                          >  1:437677/1‑151 (MQ=255)
                                                                                                                              tatTCTGTCCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                         <  1:721262/150‑1 (MQ=255)
                                                                                                                                ttCTGTTCATGGAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg                        <  1:617539/149‑1 (MQ=255)
                                                                                                                                      tCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt                >  2:920011/1‑151 (MQ=255)
                                                                                                                                          gaaGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc               <  2:298936/146‑1 (MQ=255)
                                                                                                                                             ggaaaaac‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg             <  2:831869/144‑1 (MQ=255)
                                                                                                                                                aaaaa‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccac      >  1:733956/1‑151 (MQ=255)
                                                                                                                                                      cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACcaca  <  2:367718/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGA‑TATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA  >  NC_000913/3556059‑3556355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑