Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,711,792 A→G 7.3% Y78H (TAT→CAT)  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,7920AG 7.3% 157.5 / 2.7 55Y78H (TAT→CAT) paaBring 1,2‑phenylacetyl‑CoA epoxidase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base G (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GGCCACACCATTCGCAGAGGCGCTGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATC  >  NC_002947/3711646‑3711939
                                                                                                                                                  |                                                                                                                                                   
ggCCACACCATTCGCAGAGGCGCTGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGAc                                                                                                                                                  <  1:512925/150‑1 (MQ=255)
  ccACACCATTCGCAGAGGCGCTGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCt                                                                                                                                                <  2:1104820/150‑1 (MQ=255)
        cATTCGCAGAGGCGCTGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTgtcgt                                                                                                                                          <  1:100493/150‑1 (MQ=255)
         aTTCGCAGAGGCGCTGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTgtcgtc                                                                                                                                         <  1:360753/150‑1 (MQ=255)
          ttCGCAGAGGCGCTGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTcc                                                                                                                                        >  2:78694/1‑150 (MQ=255)
                     gcTGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGaa                                                                                                                             <  2:955917/150‑1 (MQ=255)
                       tGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACa                                                                                                                           <  1:1199598/150‑1 (MQ=255)
                         gCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGc                                                                                                                         <  1:931202/150‑1 (MQ=255)
                           ccTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCgg                                                                                                                       <  1:692738/150‑1 (MQ=255)
                            cTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCggg                                                                                                                      <  1:794775/150‑1 (MQ=255)
                              ggACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTc                                                                                                                    >  2:530535/1‑150 (MQ=255)
                                 ccAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCttt                                                                                                                 <  2:623596/150‑1 (MQ=255)
                                       cgcTGTCGCCGGGTAACAGCAGGGAGGGGAGAAGGCCTTCTTTATGCACGGCCCACCTCACATGTGTCCGACTTCGTCCGGCAGCTCGTAGAACCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCtcgtc                                                                                                           <  1:557202/150‑1 (MQ=255)
                                         cTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGCCgg                                                                                                         >  2:275052/1‑150 (MQ=255)
                                                      aCAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCgg                                                                                                         <  2:228388/137‑1 (MQ=255)
                                                      aCAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCgg                                                                                                         >  1:228388/1‑137 (MQ=255)
                                                      aCAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCAAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGt                                                                                            >  2:197203/1‑150 (MQ=255)
                                                         gcagGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCg                                                                                                                        >  1:178698/1‑119 (MQ=255)
                                                         gcagGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCg                                                                                                                        <  2:178698/119‑1 (MQ=255)
                                                             gTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAgc                                                                                     >  2:978719/1‑150 (MQ=255)
                                                             gTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAgc                                                                                     >  2:471062/1‑150 (MQ=255)
                                                                   gATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGaa                                                                               <  1:530535/150‑1 (MQ=255)
                                                                         ggCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGcacc                                                                         >  2:839977/1‑150 (MQ=255)
                                                                            ttcttTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCaccacc                                                                      >  2:334833/1‑150 (MQ=255)
                                                                                ttATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCcaca                                                                  <  1:839977/150‑1 (MQ=255)
                                                                                  aTGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCAAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGa                                                                <  1:197203/150‑1 (MQ=255)
                                                                                         gCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCacac                                                         >  1:173584/1‑150 (MQ=255)
                                                                                           ccACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGTGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCt                                                       >  1:686775/1‑150 (MQ=255)
                                                                                                 cacaTGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTc                                                                                                                    <  1:1091305/83‑1 (MQ=255)
                                                                                                 cacaTGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTc                                                                                                                    >  2:1091305/1‑83 (MQ=255)
                                                                                                  acaTGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATGGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTAc                                                >  2:271368/1‑150 (MQ=255)
                                                                                                  acaTGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATGGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTAc                                                >  2:271736/1‑150 (MQ=255)
                                                                                                    aTGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGTGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTAcgg                                              >  2:580924/1‑150 (MQ=255)
                                                                                                       tgtCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTAcggcgg                                           <  2:173584/150‑1 (MQ=255)
                                                                                                        gtCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCGCACATTCGGTAcgacagc                                          >  2:858583/1‑145 (MQ=255)
                                                                                                         tCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGc                                                                                                                         <  1:355258/70‑1 (MQ=255)
                                                                                                         tCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGc                                                                                                                         >  2:355258/1‑70 (MQ=255)
                                                                                                         tCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGaa                                                                               <  1:536018/112‑1 (MQ=255)
                                                                                                         tCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGaa                                                                               >  2:536018/1‑112 (MQ=255)
                                                                                                         tCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGgt                                         <  1:517290/150‑1 (MQ=255)
                                                                                                         tCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGgt                                         >  1:751339/1‑150 (MQ=255)
                                                                                                            gACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCTAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCa                                                                                       <  1:626899/101‑1 (MQ=255)
                                                                                                            gACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCa                                                                                       >  2:626899/1‑101 (MQ=255)
                                                                                                                      ggCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGTGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCgg                            <  1:580924/150‑1 (MQ=255)
                                                                                                                          gCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGt                        >  1:761800/1‑150 (MQ=255)
                                                                                                                                aGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGa                  >  1:999439/1‑150 (MQ=255)
                                                                                                                                  aaGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATg                <  2:761800/150‑1 (MQ=255)
                                                                                                                                  aaGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATg                <  1:847122/150‑1 (MQ=255)
                                                                                                                                        gCATGGCGATGGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATg          <  1:271736/150‑1 (MQ=255)
                                                                                                                                        gCATGGCGATGGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATg          <  1:271368/150‑1 (MQ=255)
                                                                                                                                        gCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATg          >  2:543671/1‑150 (MQ=255)
                                                                                                                                          aTGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATg                <  1:849552/142‑1 (MQ=255)
                                                                                                                                          aTGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATg                >  2:849552/1‑142 (MQ=255)
                                                                                                                                           tGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATggcg       <  1:543671/150‑1 (MQ=255)
                                                                                                                                                aTAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATc  >  1:762149/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCCACACCATTCGCAGAGGCGCTGGCCCTGGACCAGGGCGCTGTCGCCGAGCAACAGCAGGTAGGGGATAAGGGCTTCTTTATGCATGGCCCACCTCACATGTGTCCGACTTCGTCGGGCAGCTCGTAGAAGCTGGCATGGCGATAGACCTTGTCGTCCGAGGGGGCGAACAGCGGGTCTTTCTCGTCGGGAGAGGAGGCGGTGATCAGCGCCGAAGGCACCACCCACAGACTCACACCTTCGCTACGGCGGGTGTAGAGCTCGCGGGCGTTCTCGATGGCCATGGCGGCATC  >  NC_002947/3711646‑3711939

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