Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,669,936 G→A R302H (CGC→CAC)  hcaE → 3‑phenylpropionate dioxygenase, large (alpha) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,669,9360GA97.4% 244.6 / NA 77R302H (CGC→CAC) hcaE3‑phenylpropionate dioxygenase, large (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (42/32);  minor base C (0/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATGCTAATGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCGCCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGTTG  >  NC_000913/2669841‑2670031
                                                                                               |                                                                                               
ggATGCTAATGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAGCAACGTTTAGGTGAAGTTCGCGCCCTGCACCTg                                                                                             >  2:271323/1‑100 (MQ=255)
    gCTAATGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCgg                                                                                         >  1:2150506/1‑100 (MQ=255)
     cTAATGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCggg                                                                                        <  1:1250270/100‑1 (MQ=255)
     cTAATGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCggg                                                                                        <  2:1895888/100‑1 (MQ=255)
     cTAATGTCTGGGTCGAGTGCGCAGTTTCAAGCTATTACCGCTAAACCTATGCCGAAGCAGAACAACGTTTATGTTAAGTTCGCGCCCTGCACCTGGCggg                                                                                        <  2:1557960/100‑1 (MQ=255)
      tAATGTTTGGGTCGGTGGCGCAGGTGAACGCTATTACCGCGAAACCTAGGCCGAAGCAGAACAACGGTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGt                                                                                       <  2:1668966/100‑1 (MQ=255)
      tAATGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACACCGTTTAGGTGACGTTCGCGCCCTGCACCTCGCGGGt                                                                                       >  2:2061799/1‑100 (MQ=255)
       aaTGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTc                                                                                      >  2:964380/1‑100 (MQ=255)
       aaTGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTc                                                                                      <  1:2028245/100‑1 (MQ=255)
        aTGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCa                                                                                     <  1:969563/100‑1 (MQ=255)
        aTGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCa                                                                                     <  2:263920/100‑1 (MQ=255)
         tGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCAt                                                                                    <  2:2131457/100‑1 (MQ=255)
         tGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCAt                                                                                    >  1:248475/1‑100 (MQ=255)
          gTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATa                                                                                   <  2:742572/100‑1 (MQ=255)
          gTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATa                                                                                   >  2:1780174/1‑100 (MQ=255)
           tCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGCCGGAGAACAAAGTTTAGGTGGAGTTCGCGACCTGCCTCTGGCGGGTAATac                                                                                  <  2:1296233/100‑2 (MQ=255)
           tCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATaa                                                                                  <  2:2372586/100‑1 (MQ=255)
            cTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAAc                                                                                 <  1:1136460/100‑1 (MQ=255)
               ggTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAat                                                                              >  2:1866350/1‑100 (MQ=255)
               ggTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCAGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAat                                                                              >  2:355180/1‑100 (MQ=255)
                      ggCGCAGTTTCAAGCTATTACCGCGAAACCTCTGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTcc                                                                       <  1:2483490/100‑1 (MQ=255)
                      ggCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTcc                                                                       <  1:744328/100‑1 (MQ=255)
                       tcgcAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTccc                                                                      <  1:750593/99‑1 (MQ=255)
                        cgcAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTccca                                                                     >  2:736580/1‑99 (MQ=255)
                           aGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACg                                                                  <  1:172672/100‑1 (MQ=255)
                            gTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGc                                                                 <  1:2447843/100‑1 (MQ=255)
                              ttCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCtt                                                               >  1:1078044/1‑100 (MQ=255)
                              ttCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCtt                                                               >  1:5304/1‑100 (MQ=255)
                                 aaGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCa                                                            >  1:234321/1‑100 (MQ=255)
                                   gCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATg                                                          >  1:1475497/1‑100 (MQ=255)
                                   gCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATg                                                          <  1:1657421/100‑1 (MQ=255)
                                   gCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATg                                                          >  1:376162/1‑100 (MQ=255)
                                    ttATTACCTCTAAACCTATGCCGAAGCAGAACAAAGTTTAGGTGTAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATgg                                                         <  2:64512/99‑1 (MQ=255)
                                      aTTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCt                                                       >  2:78413/1‑100 (MQ=255)
                                       ttACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTa                                                      <  2:58015/100‑1 (MQ=255)
                                       ttACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTa                                                      >  1:579200/1‑100 (MQ=255)
                                       ttACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTa                                                      <  1:327343/100‑1 (MQ=255)
                                       ttACCGCGAAACCTATGCCGAAGCAGAACAACGTATAGGTGAAGTACCCGCCCTGCACCTGGCGTGTAATACCAATATTTGCCCCACGCTTTCATGGCTa                                                      >  2:2170584/1‑100 (MQ=255)
                                        tACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTaa                                                     >  2:388963/1‑100 (MQ=255)
                                           cgcgAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACg                                                  >  1:182864/1‑100 (MQ=255)
                                                  ccTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCAAGCTTTCATGGCTAAACGGCACTGc                                           >  2:1141227/1‑100 (MQ=255)
                                                   cTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGcc                                          >  1:1425620/1‑100 (MQ=255)
                                                   cTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGcc                                          <  1:1736754/100‑1 (MQ=255)
                                                   cTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGACTAAACGGCACTGcc                                          <  2:1383240/100‑1 (MQ=255)
                                                    tATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCca                                         >  2:110155/1‑100 (MQ=255)
                                                    tATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGACATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGccc                                         >  1:1287815/1‑99 (MQ=255)
                                                     aTGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCcac                                        <  1:190977/100‑1 (MQ=255)
                                                       gCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCacac                                      >  2:1066074/1‑100 (MQ=255)
                                                       gCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCacac                                      >  2:741749/1‑100 (MQ=255)
                                                        ccGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCGCTGCACCTGGCGGGTCATGACAATATTCTCCCCACGCTTGCATGGCTAAACGGCACTGCCACACt                                     >  1:2457889/1‑100 (MQ=255)
                                                           aaGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCcg                                  >  1:793835/1‑100 (MQ=255)
                                                           aaGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCcg                                  <  1:1724012/100‑1 (MQ=255)
                                                             gCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCcgcg                                <  1:2584027/100‑1 (MQ=255)
                                                             gCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCcgcg                                >  2:452217/1‑100 (MQ=255)
                                                             gCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCcgcg                                >  1:48218/1‑100 (MQ=255)
                                                              cAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGt                               >  1:2607086/1‑100 (MQ=255)
                                                              cAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGt                               >  1:547709/1‑100 (MQ=255)
                                                               aGAACAACGTTTAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTc                              >  2:1852145/1‑100 (MQ=255)
                                                                  acaacGTTTAGGTGAAGTTCGCGCCCTGCGCCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTgg                           <  2:794834/100‑1 (MQ=255)
                                                                          tAGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCcgcg                   <  1:17653/100‑1 (MQ=255)
                                                                          tAGGTGAAGTGCGCGCCCTGTATCGGGCGGGTCATAAGAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCAGCcgcg                   <  2:888773/100‑1 (MQ=255)
                                                                           aGGTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCcgcgc                  >  2:1220485/1‑100 (MQ=255)
                                                                            ggTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCcgcgcg                 >  2:662865/1‑100 (MQ=255)
                                                                             gTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCgg                >  2:265940/1‑100 (MQ=255)
                                                                             gTGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTACCAACGCTTTCATGGCTCAACGGCCCTGCCACACTCCGCGTATGGCATCCGAGCgg                >  1:2515523/1‑100 (MQ=255)
                                                                              tGAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGc               >  1:326669/1‑100 (MQ=255)
                                                                               gAAGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGcc              <  1:901365/100‑1 (MQ=255)
                                                                               gAAGTTCGCGACCTGCACCTGGCGGGTCATAACAATATTTTCCACACGCTATCATGACTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGcc              >  2:1141602/1‑100 (MQ=255)
                                                                                 aGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCt            >  2:727333/1‑100 (MQ=255)
                                                                                 aGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCt            >  2:1191811/1‑100 (MQ=255)
                                                                                 aGTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCt            >  1:1700537/1‑100 (MQ=255)
                                                                                  gTTCGCGCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTg           <  1:1499399/100‑1 (MQ=255)
                                                                                     cgcgCCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATc        <  1:196193/100‑1 (MQ=255)
                                                                                      gcgcCCTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCa       <  1:1093124/100‑1 (MQ=255)
                                                                                         cccTTCCCCGGGCGGCGCAAAACACTAGTATCCCCACGCTTTCATGGCTAAGGGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGt    <  2:2119047/100‑1 (MQ=255)
                                                                                          ccTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGtt   >  2:2409717/1‑100 (MQ=255)
                                                                                           cTGCACCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGTTg  <  1:332122/100‑1 (MQ=255)
                                                                                               |                                                                                               
GGATGCTAATGTCTGGGTCGATGGCGCAGTTTCAAGCTATTACCGCGAAACCTATGCCGAAGCAGAACAACGTTTAGGTGAAGTTCGCGCCCTGCGCCTGGCGGGTCATAACAATATTTTCCCCACGCTTTCATGGCTAAACGGCACTGCCACACTCCGCGTCTGGCATCCGCGCGGCCCTGATCAAGTTG  >  NC_000913/2669841‑2670031

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