Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,716,177 G→C 13.7% T52S (ACC→AGC) ‡ paaI ← hydroxyphenylacetyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,716,1770GC13.7% 135.5 / 14.5 51T52S (ACC→AGC) ‡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 = 7.55e-01

GGGCGGTGAGCACGTCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGT  >  NC_002947/3716037‑3716317
                                                                                                                                            |                                                                                                                                            
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                                                                                                                                      <  1:1200363/130‑2 (MQ=255)
                   gcgcAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCACCTGCCATgg                                                                                                                                      <  1:1200521/130‑2 (MQ=255)
                   gcgcAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCACCTGCCATgg                                                                                                                                      >  2:1200363/1‑129 (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                                                  <  1:681425/150‑1 (MQ=255)
                                                                                   ttGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTgg                                                  <  2:23126/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)
                                                                                                                                            |                                                                                                                                            
GGGCGGTGAGCACGTCATCGCGCAGCGCCGGGGCCAGGTAGTCGATGGTGCAGCCCAGCGCCACCGTGGCCTGGTCGTAGCTATTGCAGGCGAAGGCGAACGCCGAATCGGCCAACGCGAACAGGAAACCGCCATGGCAGGTGCCATGCCCCTGGATCATGTCGGCACGTACCGCCATGCGCAGGCTGGCCCGGCCCGGGCCGGCATCCAGCAGGCTGATACCCAGGCCCTGGGTGGCCGGGTCGCGGGCGTACATGGCGTCGGCACAGGCTTGTGCCAGT  >  NC_002947/3716037‑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: