Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,704,886 A→G 7.3% L387L (TTG→CTG)  paaZ ← bifunctional oxepin‑CoA hydrolase/3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,704,8860AG 7.3% 156.8 / 2.8 55L387L (TTG→CTG) paaZbifunctional oxepin‑CoA hydrolase/3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base G (2/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

GCCACCAGGCTGCCTTTGCCACGGGCAGCCAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCGCGCTGCTGGTCATGGGAGGCCAGGGCG  >  NC_002947/3704750‑3705034
                                                                                                                                        |                                                                                                                                                    
gCCACCAGGCTGCCTTTGCCACGGGCAGCCAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCa                                                                                                                                         >  2:492348/1‑150 (MQ=255)
  caccaGGCTGCCTTTGCCACGGGCAGCCAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCaaa                                                                                                                                       <  2:434095/150‑1 (MQ=255)
   accaGGCTGCCTTTGCCACGGGCAGCCAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAg                                                                                                                                      >  1:963525/1‑150 (MQ=255)
        gCTGCCTTTGCCACGGGCAGCCAGTGCCAAGGCTTAGTCGAGGTCGTCATAGGCCGTCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGCTCGCGGGCCTACAACAAGGTTGGCCCaacaaccac                                                                                                                                 >  1:454888/1‑143 (MQ=255)
             cTTTGCCACGGGCAGCCAGTGCCAAGGCTACGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCtt                                                                                                                            >  1:929331/1‑150 (MQ=255)
              tttGCCACGGGCAGCCAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGAACAGACACCTGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTc                                                                                                                           >  1:806572/1‑150 (MQ=255)
                        gCAGCCAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCt                                                                                                                 >  1:513086/1‑150 (MQ=255)
                            ccAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCtt                                                                                                                            <  1:1016406/135‑1 (MQ=255)
                            ccAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCtt                                                                                                                            >  2:1016406/1‑135 (MQ=255)
                                  ccAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCcgcg                                                                                                       >  1:635051/1‑150 (MQ=255)
                                   cAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCgg                                                                                                      <  2:963525/150‑1 (MQ=255)
                                   cAAGGATTCGTCGGCTTTCTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCgg                                                                                                      <  2:1050420/150‑1 (MQ=255)
                                    aaGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCggg                                                                                                     >  1:505756/1‑150 (MQ=255)
                                    aaGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCggg                                                                                                     <  2:1039481/150‑1 (MQ=255)
                                    aaGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCggg                                                                                                     <  2:95333/150‑1 (MQ=255)
                                         ttCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCaaa                                                                                                <  2:284912/150‑1 (MQ=255)
                                         ttCGTCGAGGTCGTCAGAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCaaa                                                                                                <  1:816765/150‑1 (MQ=255)
                                            gtcgAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGcc                                                                                             >  2:632735/1‑150 (MQ=255)
                                                  gtcgtcATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCt                                                                                       >  2:901581/1‑150 (MQ=255)
                                                     gtcATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGgc                                                                                    <  1:901581/150‑1 (MQ=255)
                                                      tcATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGgcg                                                                                   >  2:620495/1‑150 (MQ=255)
                                                      tcATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGgcg                                                                                   >  1:123242/1‑150 (MQ=255)
                                                          aGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGa                                                                               >  1:1046990/1‑150 (MQ=255)
                                                           ggCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGaa                                                                              <  1:739795/150‑1 (MQ=255)
                                                            gCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAGCAAGGTTGgc                                                                                                                                            <  1:674435/87‑1 (MQ=255)
                                                            gCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAGCAAGGTTGgc                                                                                                                                            >  2:674435/1‑87 (MQ=255)
                                                            gCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAGCAAGGTTGgc                                                                                                                                            >  2:670778/1‑87 (MQ=255)
                                                            gCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAGCAAGGTTGgc                                                                                                                                            <  1:670778/87‑1 (MQ=255)
                                                             ccATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACa                                                                            >  1:6770/1‑150 (MQ=255)
                                                               aTCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCaaa                                                                                                >  2:662728/1‑128 (MQ=255)
                                                               aTCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCaaa                                                                                                <  1:662728/128‑1 (MQ=255)
                                                                  aGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTgg                                                                       >  1:529997/1‑150 (MQ=255)
                                                                        cTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACaa                                                                 <  2:961979/150‑1 (MQ=255)
                                                                                   gAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGcagcag                                                      >  1:97437/1‑150 (MQ=255)
                                                                                    aaCGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACaa                                                                 >  2:437802/1‑138 (MQ=255)
                                                                                    aaCGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACaa                                                                 <  1:437802/138‑1 (MQ=255)
                                                                                         ttCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTgcg                                                >  1:79025/1‑150 (MQ=255)
                                                                                         ttCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTgcg                                                >  1:1067995/1‑150 (MQ=255)
                                                                                            gATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCAc                                             >  1:301863/1‑150 (MQ=255)
                                                                                             atatCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCAcc                                            <  2:635051/150‑1 (MQ=255)
                                                                                                        gccgccCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTc                                                                     >  1:485520/1‑114 (MQ=255)
                                                                                                        gccgccCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTc                                                                     <  2:485520/114‑1 (MQ=255)
                                                                                                        gccgccCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGAc                                 >  2:56619/1‑150 (MQ=255)
                                                                                                         ccgccCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACa                                <  2:756201/150‑1 (MQ=255)
                                                                                                            cccTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATcg                             <  2:1046990/150‑1 (MQ=255)
                                                                                                             ccTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATcgc                            >  2:416619/1‑150 (MQ=255)
                                                                                                               tCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCgcgc                          >  2:610866/1‑150 (MQ=255)
                                                                                                                cGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCGCgct                         <  2:6770/150‑1 (MQ=255)
                                                                                                                 ggCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCGCgctg                        <  2:629456/150‑1 (MQ=255)
                                                                                                                    aTGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAAc                                                                             >  1:853091/1‑94 (MQ=255)
                                                                                                                    aTGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAAc                                                                             <  2:853091/94‑1 (MQ=255)
                                                                                                                       cGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCGCGCTGCTGGTc                  <  2:187459/150‑1 (MQ=255)
                                                                                                                       cGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCGCGCTGCTGGTc                  >  1:187459/1‑150 (MQ=255)
                                                                                                                                  ccTTCAACAATGTTGGCGCAAAGAACGCGCCTTCCGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCGCGCTGCTGGTCATGGGAGGCCa       <  2:505756/150‑1 (MQ=255)
                                                                                                                                     gCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCTGGCTGCGCACCCGCTCGCCGACATCGCGCTGCTGGTCATGGGAGGCCAggg    <  2:473050/150‑1 (MQ=255)
                                                                                                                                       aacaaGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCGCGCTGCTGGTCATGGGAGGCCAGGgcg  <  2:123242/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
GCCACCAGGCTGCCTTTGCCACGGGCAGCCAGTGCCAAGGCTTCGTCGAGGTCGTCATAGGCCATCAGCGTACTGACCGGGCCGAACGCTTCGATATCGTGGGCGCCGCCCTCGGCATGCGGATCGCGGGCCTGCAACAAGGTTGGCGCAAAGAACGCGCCTTCAGCCACGCCTTCGCCCCGCGGGGCAAAGCCATCGCTGGCGCCGAACAGCTGGTCACAACTCTGCAGCAGGCTGCGCACCCGCTCGCCGACATCGCGCTGCTGGTCATGGGAGGCCAGGGCG  >  NC_002947/3704750‑3705034

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