Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,716,176 G→C 13.8% T52T (ACC→ACG) ‡ paaI ← hydroxyphenylacetyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,716,1760GC13.8% 135.3 / 15.3 51T52T (ACC→ACG) ‡paaIhydroxyphenylacetyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base C (4/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

TCGCTGGCGTGGGCGGTGAGCACGTCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGT  >  NC_002947/3716027‑3716317
                                                                                                                                                     |                                                                                                                                             
tCGCTGGCGTGGGCGGTGAGCACGTCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAg                                                                                                                                               <  2:1052650/150‑1 (MQ=255)
          gggCGGTGAGCACGTCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAACAATGTGCTACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGcc                                                                                                                                     <  2:766570/150‑1 (MQ=255)
             cGGTGAGCACGTCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCt                                                                                                                                  >  1:5363/1‑150 (MQ=255)
              ggTGAGCACGTCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTg                                                                                                                                 <  1:622508/150‑1 (MQ=255)
                          aTCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTTGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGc                                                                                                                     >  2:587629/1‑150 (MQ=255)
                            cgcgCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGACGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCAc                                                                                                                   >  1:878866/1‑150 (MQ=255)
                             gcgcAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCACCTGCCATgg                                                                                                                                      >  2:1200363/1‑129 (MQ=255)
                             gcgcAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCACCTGCCATgg                                                                                                                                      <  1:1200363/130‑2 (MQ=255)
                             gcgcAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCACCTGCCATgg                                                                                                                                      <  1:1200521/130‑2 (MQ=255)
                             gcgcAGCGCCGGGGCAAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCACCTGCCATgg                                                                                                                                      >  2:1200521/1‑129 (MQ=255)
                                 agggacGGGGCCAGGTAGTCGATGGTTCAGCACAGCGCGACTGTGGCTTGGTCGTAGCTATTGCAGGCGAAGGCGAGGGCCGAATCGGTCGCCGCGGATCGGAAGCCGCCATGCCAGGTGCCATTCTCGTGGAACATTTCGGCACGTAcc                                                                                                              <  2:104167/145‑1 (MQ=255)
                                      cGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCAt                                                                                                         >  2:1153606/1‑150 (MQ=255)
                                            cAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAg                                                                                                   >  2:681425/1‑150 (MQ=255)
                                                   tCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTggcc                                                                                            >  2:701610/1‑150 (MQ=255)
                                                       gggTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTggcccggc                                                                                        <  2:803525/149‑1 (MQ=255)
                                                          tGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATTATTTCGGCACGTACCGCCATGCGCAGGCTGgcccggcccg                                                                                     >  1:615935/1‑150 (MQ=255)
                                                           gCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGcccggcccgg                                                                                    >  1:307076/1‑150 (MQ=255)
                                                                 cAGCGCCACTGTGGCCTGGTGGTAGCTATTGCAGGCGAAGGCTGACGCCGAATCGGCTAACGCGAAAAGGAAACCGCCATGGAAGGTGCCATCCCCGCGTATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGGCCGGGCCAGc                                                                              <  2:203784/150‑1 (MQ=255)
                                                                   gcgcCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCAt                                                                            >  2:249979/1‑150 (MQ=255)
                                                                      ccACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCca                                                                         <  2:844526/150‑1 (MQ=255)
                                                                             ggCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCt                                                                  <  2:1910/150‑1 (MQ=255)
                                                                                cttgTCGTAGCTCTTGCAGGCGAGGGAGGAGGTCTTGTCGATGGACGGGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCTGCATCCAGCAGGCTGAt                                                               <  2:615935/147‑1 (MQ=255)
                                                                                cTGGTCGTAGCTATTGCAGTCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGACTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGAt                                                               >  1:553762/1‑150 (MQ=255)
                                                                                cTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCAt                                                                                                         >  1:56477/1‑108 (MQ=255)
                                                                                cTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCAt                                                                                                         <  2:56477/108‑1 (MQ=255)
                                                                                cTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGTGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGAt                                                               >  1:35849/1‑150 (MQ=255)
                                                                                cTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGAt                                                               >  1:912425/1‑150 (MQ=255)
                                                                                      gTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCca                                                                         >  2:485868/1‑134 (MQ=255)
                                                                                      gTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCca                                                                         <  1:485868/134‑1 (MQ=255)
                                                                                             ttGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTgg                                                  <  2:23126/150‑1 (MQ=255)
                                                                                             ttGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTgg                                                  <  1:681425/150‑1 (MQ=255)
                                                                                               gCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGTGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGt                                                <  2:35849/150‑1 (MQ=255)
                                                                                                          cgaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTcgc                                     <  1:609051/150‑1 (MQ=255)
                                                                                                           gaaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTcgcg                                    >  1:1029974/1‑150 (MQ=255)
                                                                                                            aaCGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATACAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCgg                                   <  1:1153606/150‑1 (MQ=255)
                                                                                                                     tCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATATCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATg                          <  1:1150248/150‑1 (MQ=255)
                                                                                                                          cAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTc                     >  2:201197/1‑150 (MQ=255)
                                                                                                                              gcgAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGca                 >  1:307320/1‑150 (MQ=255)
                                                                                                                                     ggAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCtt          <  2:878866/150‑1 (MQ=255)
                                                                                                                                      gAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCCGCATCCAGCATGCCGATACCCCGGCCCTGGGTGGCCTGGTCGCGGGCGTACATGGCGTCGGCACAGGCTtg         >  2:540207/1‑150 (MQ=255)
                                                                                                                                       aaaCCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTtgt        >  2:376188/1‑150 (MQ=255)
                                                                                                                                       aaaCCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTtgt        >  2:376092/1‑150 (MQ=255)
                                                                                                                                           gggcATGGCACCTGCCATGGCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCggg                                         <  1:772360/110‑1 (MQ=255)
                                                                                                                                           gggcATGGCACCTGCCATGGCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCggg                                         >  2:772360/4‑113 (MQ=255)
                                                                                                                                            ggcATGGCACCTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATCAGGAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAg   >  1:12429/3‑150 (MQ=255)
                                                                                                                                            gCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCCGGCTGATACCCAGGCCCTGGGTGGCCGGGTcgc                                     >  1:85165/1‑116 (MQ=255)
                                                                                                                                            gCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTcgc                                     <  2:85165/116‑1 (MQ=255)
                                                                                                                                             ccATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGACGGCATAATGCAGGCTGCTACCCAGGCCCGGGGTGGGCGGGTGGTGGTCGTATATGGCTTCGGCACAGGCCTGTGccgac  <  2:740526/150‑4 (MQ=255)
                                                                                                                                               aTGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCg                                                                                >  2:602963/1‑70 (MQ=255)
                                                                                                                                               aTGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCg                                                                                <  1:602963/70‑1 (MQ=255)
                                                                                                                                                tGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTggcccggc                                                                                        <  2:138995/61‑1 (MQ=255)
                                                                                                                                                tGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTggcccggc                                                                                        >  1:138995/1‑61 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
TCGCTGGCGTGGGCGGTGAGCACGTCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGT  >  NC_002947/3716027‑3716317

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: ‑