Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,714,216 C→G 7.0% G100R (GGC→CGC)  paaK ← phenylacetate‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,714,2160CG 7.0% 171.7 / 2.3 55G100R (GGC→CGC) paaKphenylacetate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base G (2/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GTGAACAGCCCGTAGCCGTAGGACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAGGCACTTGAGGTCATCAGGGTGGATGCCGTA  >  NC_002947/3714074‑3714360
                                                                                                                                              |                                                                                                                                                
gTGAACAGCCCGTAGCCGTAGGACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTgg                                                                                                                                             <  7:186504/148‑1 (MQ=255)
gTGAACAGCCCGTAGCCGTAGGACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTgg                                                                                                                                             >  8:186504/1‑148 (MQ=255)
      aGCCCGTAGCCGTAGGACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATg                                                                                                                                      >  8:94096/1‑149 (MQ=255)
        cccGTAGCCGTAGGACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATg                                                                                                                                      >  8:132397/1‑147 (MQ=255)
          cgtagccgtagGACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAgg                                                                                                                                  <  7:94096/149‑1 (MQ=255)
            tagccgtagGACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACCGGAGGCATGCAGgcg                                                                                                                                >  5:171088/1‑149 (MQ=255)
                    ggACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCGCCTGCATGCCTCCAGTGGcaccacc                                                                                                                         <  3:54456/148‑2 (MQ=255)
                    ggACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCAct                                                                                                                         <  8:153926/148‑1 (MQ=255)
                     gACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCGCCTGCATGCCTCCAGTGGcaccaccgg                                                                                                                       <  5:246564/149‑4 (MQ=38)
                     gACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCActtc                                                                                                                       <  2:284582/149‑1 (MQ=255)
                      acacAGGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCActtct                                                                                                                      <  5:268993/149‑1 (MQ=255)
                        acaTGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACttcttc                                                                                                                    <  7:68311/149‑1 (MQ=255)
                        acaTGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACttcttc                                                                                                                    <  7:132397/149‑1 (MQ=255)
                         caTGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTcc                                                                                                                   <  5:55976/149‑1 (MQ=255)
                           tGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACCGGAGGCATGCAGGCGCACCACTTCTTCCTg                                                                                                                 <  6:171088/149‑1 (MQ=255)
                                      aCCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCa                                                                                                       >  4:126058/1‑148 (MQ=255)
                                      aCCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCa                                                                                                       >  6:3273/1‑148 (MQ=255)
                                       ccTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCaaa                                                                                                     <  2:115720/149‑1 (MQ=255)
                                       ccTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCaa                                                                                                      >  7:169905/1‑148 (MQ=255)
                                        cTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAAc                                                                                                    <  5:3273/149‑1 (MQ=255)
                                            gcgcCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGTGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATg                                                                                                 <  5:165429/148‑1 (MQ=255)
                                             cgcccgccGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCGAGTGGAGGCATGCAGGCGCAGCACTTCGTCCTGGGGCACGGCAAACAGGcc                                                                                               >  8:261474/1‑149 (MQ=255)
                                             cgcccgccGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGc                                                                                                <  3:126058/148‑1 (MQ=255)
                                                    cggccgCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAgggg                                                                                        >  8:17869/1‑149 (MQ=255)
                                                    cggccgCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGCGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAgggg                                                                                        >  5:95234/1‑149 (MQ=255)
                                                       ccgCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAg                                                                                     <  1:151119/149‑1 (MQ=255)
                                                        cgCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGt                                                                                    >  2:4704/1‑149 (MQ=255)
                                                            cGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTc                                                                                >  3:289060/1‑149 (MQ=255)
                                                                  gAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAg                                                                          >  8:138130/1‑149 (MQ=255)
                                                                         cAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTc                                                                   <  7:138130/149‑1 (MQ=255)
                                                                            ccACATTGGCCCAGGTGTCGATATCGTTCTGTGAGTAGACGACCACCGTCGGCTTTCCGGTGGTGCCACTGGATGCATGCAGGCTCACCACTTATTCCTGCCGCACGGCAAACATGCAGTAGGGGTACTTGACGCGCAGGTCGTTCTTg                                                                >  5:65830/1‑149 (MQ=255)
                                                                            ccACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTc                                                                                <  4:267964/133‑1 (MQ=255)
                                                                            ccACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTc                                                                                >  3:267964/1‑133 (MQ=255)
                                                                                aTTGGCCCAGGTGTCGATATCGTTCTGGGTGTACCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCATCCTGGGGCACCGCAAACATTCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCgg                                                            >  7:280734/1‑149 (MQ=255)
                                                                                    gCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTg                                                          >  6:144676/1‑147 (MQ=255)
                                                                                    gCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACAGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTATTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGtt                                                          >  2:139997/1‑146 (MQ=255)
                                                                                     cccAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAg                                                       <  7:17869/149‑1 (MQ=255)
                                                                                      ccAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGTGTAGTTGTCGCGCAGGTAGTTCTTGCCGGTGAAgg                                                      <  8:169905/149‑1 (MQ=255)
                                                                                       cAGGTGTCGATATCGTTCTGGGCGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAggg                                                     <  6:95234/149‑1 (MQ=255)
                                                                                               gATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCg                                             >  7:92681/1‑149 (MQ=255)
                                                                                                atatGGTTCTGGGTGTAGCCGACGACCGTCGGCTTTCCGGTGGTGCCTCTGGAGGCATGGTGGCGGACCAGTTGTTGCTGGGGCACGGCAAAGATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGc                                            <  8:92681/149‑1 (MQ=255)
                                                                                                atatCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGc                                            >  2:138606/1‑149 (MQ=255)
                                                                                                  atCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGc                                            >  4:48860/1‑147 (MQ=255)
                                                                                                   tCGTTCTGTGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGTAACTTCGCCAg                                         <  2:166597/149‑1 (MQ=255)
                                                                                                     gTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAgg                                        <  5:144676/148‑1 (MQ=255)
                                                                                                           ggTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTc                                   <  1:138606/147‑1 (MQ=255)
                                                                                                             tgtAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAgg                               >  6:184555/1‑149 (MQ=255)
                                                                                                                    gACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAGGCACtt                          >  4:210602/1‑147 (MQ=255)
                                                                                                                         ccGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGAAGGTCGTTCTTGCCGGTGAAGGGGAAATTCGACAGGTCTTCCAGGCACTTGAGGtca                   >  4:36955/1‑149 (MQ=255)
                                                                                                                              ggCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAGGCACTTGAGGTCATCAgg              >  8:107200/1‑149 (MQ=255)
                                                                                                                                 tttCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAGGCACTTGAGGTCATCAGGGTg           >  5:113540/1‑149 (MQ=255)
                                                                                                                                  ttCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAGGCACTTGAGGTCATCAGGGTg           >  7:229222/1‑148 (MQ=255)
                                                                                                                                    gaagtggtgCGCCTGCATGCCTCCAGTCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGTAACTTCGCCAGGTCTTCCAGGCACTTGAGGTCATCAGGGTGGAt        >  6:131346/4‑149 (MQ=255)
                                                                                                                                    ccGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAGGCACTTGAGGTCATCTGGGTGGAt        >  6:300061/1‑149 (MQ=255)
                                                                                                                                          gtgCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAGGCACTTGAGGTCATCAGGGTGGATGCCGTa  <  5:184555/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GTGAACAGCCCGTAGCCGTAGGACACATGCACCTTGTCACCTTTGCGCCCGCCGGCCGCACGAATCGAACGGGCAACCACATTGGCCCAGGTGTCGATATCGTTCTGGGTGTAGCCGACCACCGTCGGCTTTCCGGTGGTGCCACTGGAGGCATGCAGGCGCACCACTTCTTCCTGGGGCACGGCAAACATGCCGTAGGGGTAGTTGTCGCGCAGGTCGTTCTTGCCGGTGAAGGGGAACTTCGCCAGGTCTTCCAGGCACTTGAGGTCATCAGGGTGGATGCCGTA  >  NC_002947/3714074‑3714360

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