Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,711,911 T→C 14.3% S24G (AGC→GGC)  paaB ← ring 1,2‑phenylacetyl‑CoA epoxidase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,711,9110TC14.3% 62.2 / 4.2 28S24G (AGC→GGC) paaBring 1,2‑phenylacetyl‑CoA epoxidase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (12/12);  new base C (2/2);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACTGCGG  >  NC_002947/3711768‑3712052
                                                                                                                                               |                                                                                                                                             
ggCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGa                                                                                                                                          <  8:88944/149‑1 (MQ=255)
            gTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGt                                                                                                                              >  7:122787/1‑149 (MQ=255)
                      tcgGGAGAGGTGGCGGTGATCAGCGCCGAAGGCACCACCCTCAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGa                                                                                                                                          <  3:7022/127‑1 (MQ=255)
                      tcgGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGATCTCGCTGGACTACTAAAAATAAATAACAACATCGGCGTCGTGCACGCTGCCGa                                                                                                                                          >  4:7022/1‑127 (MQ=255)
                                      ttaTCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAac                                                                                                    <  5:264426/147‑1 (MQ=255)
                                            gcgcCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCACGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCg                                                                                              <  1:264318/149‑1 (MQ=255)
                                                                agccTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGGACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGCCAt                                                                          >  5:82036/4‑149 (MQ=255)
                                                                    tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGTCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTa                                                                                            <  6:82880/127‑1 (MQ=255)
                                                                    tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGTCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTa                                                                                            >  5:82880/1‑127 (MQ=255)
                                                                    tCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                      >  1:185741/1‑149 (MQ=255)
                                                                      acacCTTCGCTACGGCGGGTGTAGAGCTCGCGTGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                      <  4:288410/147‑1 (MQ=255)
                                                                        acCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCTGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGAtc                                                                  <  6:285642/149‑1 (MQ=255)
                                                                                                     gggCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCATGCGGCGt                                      >  7:291721/1‑148 (MQ=255)
                                                                                                             tcGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTGGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCt                             >  5:281751/1‑149 (MQ=255)
                                                                                                                   gCCATGGCGGCATCGGCTGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACAGTTCGTAGAGGTTCAAGAGAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGttgc                        <  1:98441/148‑1 (MQ=255)
                                                                                                                      aTGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcg                                          >  1:172799/1‑127 (MQ=255)
                                                                                                                      aTGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcg                                          <  2:172799/127‑1 (MQ=255)
                                                                                                                                 ggcggcGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcg                                          >  1:34519/1‑116 (MQ=255)
                                                                                                                                 ggcggcGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAggcg                                          <  2:34519/116‑1 (MQ=255)
                                                                                                                                 ggcagcGTGCACGCCGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCtt                                  >  3:292685/1‑124 (MQ=255)
                                                                                                                                 ggcagcGTGCACGCCGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCtt                                  <  4:292685/124‑1 (MQ=255)
                                                                                                                                 ggcagcGTGCACGCCGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCtt                                  >  3:297311/1‑124 (MQ=255)
                                                                                                                                 ggcagcGTGCACGCCGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTAtt                                  <  4:297311/124‑1 (MQ=255)
                                                                                                                                       gTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACTGCg   >  4:88245/1‑149 (MQ=255)
                                                                                                                                       gTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACTGCg   <  8:122787/149‑1 (MQ=255)
                                                                                                                                        tGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACTGCgg  >  5:56846/1‑149 (MQ=255)
                                                                                                                                          cACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                      >  2:282282/1‑79 (MQ=255)
                                                                                                                                          cACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTc                                                                      <  1:282282/79‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATCGGCGGCGTGCACGCTGCCGACATGTTTGTGGTTAAGGCCGTGCTTGCTGCGCACGAACACTTCGTAGAGGGTCCAGACAGACATTTTCGATCTCCGCATCTATCGCCGCTCAGGCGGCGTTCTTGTTCTGTTGCTTGCGCGCATAGGCCACTGCGG  >  NC_002947/3711768‑3712052

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: