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

GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA  >  NC_000913/3556059‑3556355
                                                                                                                                                       |                                                                                                                                                    
gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcc                                                                                                                                                    <  1:725141/151‑1 (MQ=255)
  cGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcccc                                                                                                                                                  <  2:21358/151‑1 (MQ=255)
        caccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAgt                                                                                                                                            >  2:909987/1‑151 (MQ=255)
         cccgcGGCACTGACGCGTTTTCCACCGACGAAAAAAAAGCGTTGGGTGATATTTTCTTTTTCGACTTAAAGGGAGACGTTATTCCCCGCCCCCGGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAgtg                                                                                                                                           <  2:291117/147‑1 (MQ=255)
                      cgcgTTTTCCCCCGACGAAAAAAAGGCGTTCGGTGATTTTTCCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATAAATTTTTTGTTCATGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTc                                                                                                                                 <  1:541084/148‑1 (MQ=255)
                      ccctcTTTCCCCCGACAAAAAAAAGGCGTCCGGTAATTTTTTCCCTTCCGACTCCAAGCGAGAGGTTATCCCCCGCCCCCTGGCCCCGCGGCGGTTTATAAATTTTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTc                                                                                                                                 <  2:931736/144‑1 (MQ=255)
                       gctcTTTCCACGGACGAAAAAAAAGCGTCGGTTATTTTTTTCCCTTTCGACTAAAAGCGAGCCGTTATTCCCCGCCCCCTGGCCGCGGGGAAGTTTAGGAATATTCTTTTCAAGCAGAAAAAACCCCCCCGCCAGTGTAGCGATTTc                                                                                                                                 <  2:552723/144‑1 (MQ=255)
                       gcgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGc                                                                                                                             >  1:877924/1‑151 (MQ=255)
                         gcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGAt                                                                                                                                    >  1:874032/1‑142 (MQ=255)
                              tCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCa                                                                                                                       >  1:700867/1‑150 (MQ=255)
                                  cgacgaCAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                  >  2:173524/1‑151 (MQ=255)
                                  cgacgaCAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                  >  1:570266/1‑151 (MQ=255)
                                     cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                                               >  1:367739/1‑151 (MQ=255)
                                       acAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                             >  1:446185/1‑151 (MQ=255)
                                       aaaaaaaaGCGTTCGGTGATTTTTCCCCTTTCGACTTAAAGCGAGACGTTATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                              <  1:105634/148‑1 (MQ=255)
                                         aaaaaaGGGGTCGGGTGATTTTTTCCTTTTCGACTAAAAGCGGGACTTTTTCCCCCGCCCCTTGGCCGCGCGGCGGTTTATGAATTTTCTGTTCAGGAAGAAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                   <  2:553473/144‑1 (MQ=255)
                                          aaaaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                           >  1:916123/1‑150 (MQ=255)
                                                   tCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                  >  2:593195/1‑150 (MQ=255)
                                                       ggtattttttccccTTACGCCTTCAAGCGAGACGTTTTTCCCCGCCCCCGGGTCGCGGGGGAGTTTATGAATTTTTTTTTCATGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                  <  2:726285/135‑1 (MQ=255)
                                                            tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                          >  1:798944/1‑151 (MQ=255)
                                                                ttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                                      >  2:1039024/1‑151 (MQ=255)
                                                                     ttccGACTTCAGGGGAGACGGTTTTCCCCGGCCCCAGGCCGCGGGGCAGTTAAGGAATTTTCTGTTCAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                    <  2:567272/145‑1 (MQ=255)
                                                                      tcggcctaaaggggagaGGGTATCCCCCGGCCCCGGGGCCGGCGGAAGTTTATGATTTTTTTGTTCAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                 <  1:155110/137‑1 (MQ=255)
                                                                            ttaaaGCGAGAGGTTATACCCCGCCCCCGGGCCGCGCGGAGGGTTATAAATTTTCTGTTGATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                           <  1:910028/147‑1 (MQ=255)
                                                                               aaGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                        >  2:932443/1‑150 (MQ=255)
                                                                                   gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                   >  2:316510/1‑151 (MQ=255)
                                                                                     gaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCGCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTc                                                                 <  1:158158/151‑1 (MQ=255)
                                                                                      acgttatccccgggCCCCGGGGCGCGCGCCTGTTTATAAATTTTTTTTTCTTGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                   <  1:587039/136‑1 (MQ=255)
                                                                                        gttATTCCCGGCCCCCGGGTCGCGGGGCAGTTTAGGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                                  <  2:1013311/145‑1 (MQ=255)
                                                                                            tcccccGCCCCCGGGTCGCGGGGCAGGTTATAAATATTCTGTTCATGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgcg                                                            <  1:603725/146‑1 (MQ=255)
                                                                                             tACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                         >  2:878023/1‑151 (MQ=255)
                                                                                               cccGGCCCCTGGTCGCGCGGAAGTTTATGAATTTTCTGTCCAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                                       <  1:515249/151‑1 (MQ=255)
                                                                                                              cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtgt                                        >  1:894650/1‑151 (MQ=255)
                                                                                                               gcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                       >  2:134258/1‑151 (MQ=255)
                                                                                                                cgGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                       <  1:718035/150‑1 (MQ=255)
                                                                                                                cgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                                       >  1:519077/1‑150 (MQ=255)
                                                                                                                 ggaagtttATAAATTTTCTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                     <  1:804101/145‑1 (MQ=255)
                                                                                                                 ggCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                     <  2:52851/151‑1 (MQ=255)
                                                                                                                 ggCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                                      >  1:122788/1‑150 (MQ=255)
                                                                                                                 gcgaGGTAATAAATTTTCTGTTCGGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                             <  1:1039073/140‑1 (MQ=255)
                                                                                                                            aaTATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                          >  1:437677/1‑151 (MQ=255)
                                                                                                                              tatTCTGTCCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                         <  1:721262/150‑1 (MQ=255)
                                                                                                                                ttCTGTTCATGGAGAAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACggg                        <  1:617539/149‑1 (MQ=255)
                                                                                                                                      tCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt                >  2:920011/1‑151 (MQ=255)
                                                                                                                                          gaaGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc               <  2:298936/146‑1 (MQ=255)
                                                                                                                                             ggaaaaacCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg             <  2:831869/144‑1 (MQ=255)
                                                                                                                                                aaaaaCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccac      >  1:733956/1‑151 (MQ=255)
                                                                                                                                                      cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACcaca  <  2:367718/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACA  >  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: