Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,447,483 G→T 67.4% intergenic (‑10/+25) aroC ← / ← prmB chorismate synthase/50S ribosomal subunit protein L3 N(5)‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,447,4830GT67.4% 31.6 / 21.7 35intergenic (‑10/+25)aroC/prmBchorismate synthase/50S ribosomal subunit protein L3 N(5)‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base T (16/7);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

TTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCC  >  NC_000913/2447348‑2447619
                                                                                                                                       |                                                                                                                                         
ttGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTGAttgt                                                                                                                                    >  1:714960/1‑143 (MQ=255)
      gTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTGGTCCAATTGTGGTTCCAGCCATCACGGCTCCTTTATTGTTGTGtt                                                                                                                              >  2:39042/1‑143 (MQ=255)
              ccGGCAGCGAGATGCCTGGCGGAACACCATCGACGAGGCAGCCGGGCGCCAGCCCGTGCGATTCGAGCCAGGTGGTGACGCTATAGCGTTGTCCAATTGGGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGtt                                                                                                                     <  2:290447/143‑1 (MQ=255)
                 cccgCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTAc                                                                                                                  <  1:105759/140‑1 (MQ=255)
                         aggCCGGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGTAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCtt                                                                                                          <  2:530794/141‑1 (MQ=255)
                          tGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCttt                                                                                                         <  1:354705/143‑1 (MQ=255)
                              tGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATaa                                                                                                     >  1:305754/1‑143 (MQ=255)
                               ggCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGGGTTTACTTCATCTTTATaaa                                                                                                    >  2:451275/1‑143 (MQ=255)
                                              acgaTGCAGCCGACCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtctc                                                                                     <  2:821457/143‑1 (MQ=255)
                                                     aGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTCATCAATCGCGAAATGTTCTCGTGCGGc                                                                              >  2:792621/1‑143 (MQ=255)
                                                        ggggcgcCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCGATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAAt                                                                           <  2:699062/140‑1 (MQ=255)
                                                        ccaGCGCCCGCCCGTCCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCGCCTTTATTGTTGTGTTTGCGTGTTGACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAAt                                                                           <  1:402695/141‑1 (MQ=255)
                                                               cAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAgctgc                                                                    >  1:632385/1‑143 (MQ=255)
                                                                aGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAATGAGTTGTCCATTTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAgctgct                                                                   <  2:646647/143‑1 (MQ=255)
                                                                 gCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTGCTTAATCTTTATAACTCGCGACATGGTCTCGTGCGGCAATCAgctgccc                                                                  >  2:323385/1‑141 (MQ=255)
                                                                       gCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTTTGTTATTCTCGAGAACGCGTGTATCACATTGCATTAGTAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGt                                                            <  1:604267/143‑1 (MQ=255)
                                                                           ttCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGc                                                        >  2:160984/1‑143 (MQ=255)
                                                                               ccGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGCAATGTTCTCGTGCGGCAATAAGCTGCTCTCTGGTGAGcggca                                                    >  2:810368/1‑139 (MQ=255)
                                                                                 gAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGGTTACTTAATCTTTATACATCGCGAAATGTTCTGGTGCGGCAATAAGGTGATGTATGGTGGACATaaaa                                                  >  2:237941/1‑142 (MQ=255)
                                                                                          taggcgAAGGAGTTGTCCGATTGTGTTTCCAGCCATCGCGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATcg                                         <  2:638648/140‑1 (MQ=255)
                                                                                          tACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATcg                                         <  2:632385/143‑1 (MQ=255)
                                                                                             gcgAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATcgccg                                      >  1:816770/1‑143 (MQ=255)
                                                                                                       tGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCaa                            >  1:609213/1‑143 (MQ=255)
                                                                                                        gTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCaaa                           <  1:160984/143‑1 (MQ=255)
                                                                                                         tCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATAGCCGCCGTTATCAAAc                          >  2:393165/1‑143 (MQ=255)
                                                                                                           cAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTGACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGTGCATAAACACACCATCGCCGCCGTTATCGAACTc                        >  1:806310/1‑143 (MQ=255)
                                                                                                                gtgtTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAACTGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCACCGCCGCCGTCACCAAACTccggcc                   >  2:80326/1‑143 (MQ=255)
                                                                                                                gtgtTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACCCACCATCGCCGCCGCTATCAAACTccagcc                   >  2:779127/1‑143 (MQ=255)
                                                                                                                  gtTTCCAGCCATCACGGCGCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTccagccag                 <  1:171292/143‑1 (MQ=255)
                                                                                                                    ttCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGt               >  1:693710/1‑143 (MQ=255)
                                                                                                                    ttCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGt               >  1:454988/1‑143 (MQ=255)
                                                                                                                     tCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGAAATAAGCTGCTCTTAGCTGAACATACACCCACCATCGCCGCAGTTATCAAACTccagcgcggcg              >  2:665458/1‑137 (MQ=255)
                                                                                                                      ccAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGa             >  1:280438/1‑143 (MQ=255)
                                                                                                                              cACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACATACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACa     >  1:212612/1‑143 (MQ=255)
                                                                                                                                 ggCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCCAACTCCAGCCAGGTGAACGGAACATcc  >  1:52784/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                         
TTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCC  >  NC_000913/2447348‑2447619

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: