Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,570,438 A→C 13.2% intergenic (+46/‑142) pcaF‑I → / → pcaT beta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,570,4380AC13.2% 117.2 / 3.3 56intergenic (+46/‑142)pcaF‑I/pcaTbeta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base C (3/5);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACC  >  NC_002947/1570293‑1570585
                                                                                                                                                   |                                                                                                                                                   
tgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACg                                                                                                                                                   >  8:289393/1‑148 (MQ=255)
        aCCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTAGCGGGGAAGGCCCCGGCAAAAGGGCTCCACGCAATAAccc                                                                                                                                          >  7:125023/1‑148 (MQ=255)
          cAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                                        <  4:283319/149‑1 (MQ=255)
          cAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                                        <  7:289393/149‑1 (MQ=255)
                 agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                                                 <  4:81953/149‑1 (MQ=255)
                  gaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                                                 <  1:76920/149‑1 (MQ=255)
                            ggTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATCCCTGTTGGAGCGGGCTTCCCcgc                                                                                                                      <  7:240456/149‑1 (MQ=18)
                                gCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTACCGGGGGAGCCCCCACCAACAGGGCTCCACGCAATACACGGTGGGAGCGCGGTTGCCcgcg                                                                                                                     >  2:282242/1‑146 (MQ=255)
                                  aaGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGACCGCGGCAACAGGGATCCACGCAATCCCCGGTGGGAGCGGGCGTCCCCGCGAagcg                                                                                                                >  4:67606/1‑147 (MQ=255)
                                    gggCCTGGCGACCATGTGGGGGGGTGTCGGCCAAGGTCGGGCGTGGGCCATCGAGCGGGTTTGCCCGGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCACCAGGGAGCCACGCAATCACCTGGGGGAGCGGGCTTCCCCGCGAAGAGGc                                                                                                              <  6:101454/149‑1 (MQ=255)
                                        cTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGCAGCCCGCTCCAACAGGGGTCCCCGCCATCACCTGTGGGCGCGGGCTTGCCCGCGAAGAGGCcgct                                                                                                          >  3:215439/1‑149 (MQ=255)
                                        cTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGCTCAACGCCATCACCTGTGGGAGCGGGCTTCCCCGCGCAGAGGcagag                                                                                                          >  3:219745/1‑145 (MQ=255)
                                           gCGACCATGTGTGTGGTGTCGGCTAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCCCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCtg                                                                                                      >  2:105310/1‑149 (MQ=255)
                                             gACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGc                                                                                                              <  1:101939/140‑1 (MQ=255)
                                             gACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATAACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGc                                                                                                              >  2:101939/1‑140 (MQ=255)
                                                 aTGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTgaga                                                                                                 >  8:101641/1‑149 (MQ=255)
                                                      tgtgGGTGCCGGCCACGGTCTGGCGTTGCCCATGGCGCGGGTTGCCCCTGTCCCGCCCGCTCCGCGGGCACCCCCGCCCCACCCGGTACCCCCCCAAACCCCTGGGGGCGGGGGCTCCCCCGCGAAGAGGCCTCTCCTGTGAGACATaa                                                                                            <  4:320491/149‑1 (MQ=255)
                                                       gggggTGGCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAg                                                                                           <  2:48965/147‑1 (MQ=255)
                                                        tgGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTcc                                                                                                                          >  6:207282/1‑117 (MQ=255)
                                                        tgGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTcc                                                                                                                          <  5:207282/117‑1 (MQ=255)
                                                          ggTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGcc                                                                                        >  2:56424/1‑149 (MQ=255)
                                                                ggCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGACCGCTACAACAGGGGTGAACGGAATCACCTGTGGGGGCGGGGGTGCCCGCGAAGCGGCCGGGCCGGGGGGACATAAGGGCCATTgg                                                                                  >  7:215170/1‑148 (MQ=255)
                                                                 gCCAAGGTCGGGCGTTGCCCAGCGAGCGGGTATGACCGGTACCGGCCGCTTCGCGGGGAAGCCCGCTCCAACCGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGcc                                                                                 <  8:227144/149‑1 (MQ=255)
                                                                         cTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCGCCAACACGGATGCCCCCAATCACCAGGGGGGGCGGGGTTCCCCGCGAAGAGGCCGCGCCGGTGAGAAATAAGGCCCATTgcccgccccct                                                                         >  8:198466/1‑146 (MQ=255)
                                                                          tGGCGTTGGCCACCGAGCGGGTTGGACCTGTACCGCCCTCTTCGCGGGCACGCCCGCTCCAACAGGTATCCACGAAATCCCCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATa                                                                        <  8:297711/149‑1 (MQ=38)
                                                                            gCGTTGCCCATCCAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGCGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATa                                                                        <  7:198466/147‑1 (MQ=255)
                                                                                    cATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTACAACAGGGATCCACGCAATCACCAGTGGGAGCGGGCTTCCCCCCGAAGAGGCCGATCCGGTGAGACAAAAGGc                                                                                         >  2:261592/1‑122 (MQ=38)
                                                                                     aTCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGTAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgt                                                             <  5:137775/149‑1 (MQ=255)
                                                                                     aTCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgt                                                             >  2:285068/1‑149 (MQ=255)
                                                                                     aTCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgt                                                             <  1:285068/149‑1 (MQ=255)
                                                                                        gAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgtgtg                                                          <  1:105310/149‑1 (MQ=255)
                                                                                         tgCGGGTTTGACCTTACCCGGCGTCTTCGGGGGCACCCCCGCTCCAACTGGGATCCCCGAAATCACCTGTGGGAGGGGGCTGCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgtgtgt                                                         <  1:152654/148‑1 (MQ=255)
                                                                                         aGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgtgt                                                           <  2:234990/147‑1 (MQ=255)
                                                                                         aGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgtgt                                                           >  1:234990/1‑147 (MQ=255)
                                                                                           cGGGTTTTACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTa                                                       <  1:282242/149‑1 (MQ=255)
                                                                                                   aCCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGc                                                                                  <  2:186325/114‑1 (MQ=255)
                                                                                                   aCCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGc                                                                                  >  1:186325/1‑114 (MQ=255)
                                                                                                           cGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGaa                                       >  7:283917/1‑149 (MQ=255)
                                                                                                                tctTCGCGGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATCCCTGTTGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTAcc                                  >  6:90039/1‑149 (MQ=18)
                                                                                                                 ctTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTAcc                                  >  2:253892/1‑148 (MQ=255)
                                                                                                                    cgcgGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAg                               <  1:253892/148‑1 (MQ=255)
                                                                                                                    cgcgGGGAACCCGGCGACAACAGGGATCCACCAAATCACCTGTGGCAGCGGGATTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTa                                <  3:67606/147‑1 (MQ=255)
                                                                                                                     gcgGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATCCCTGTTGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGa                             <  5:90039/149‑1 (MQ=18)
                                                                                                                        ggggAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGAcc                           >  1:3657/1‑149 (MQ=255)
                                                                                                                           gAAGCCCGCTCCACCAGGGCTCCCCCCAATCACCGGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTc                        <  8:283917/149‑1 (MQ=255)
                                                                                                                           gAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTc                        >  8:179580/1‑149 (MQ=255)
                                                                                                                               cccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg                     >  5:317799/1‑148 (MQ=255)
                                                                                                                                      cAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAc             >  1:113434/1‑149 (MQ=255)
                                                                                                                                        aCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAg                                                                       >  1:225933/1‑89 (MQ=255)
                                                                                                                                        aCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAg                                                                       <  2:225933/89‑1 (MQ=255)
                                                                                                                                        aCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgt                                                             >  8:149998/1‑99 (MQ=255)
                                                                                                                                        aCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgt                                                             <  7:149998/99‑1 (MQ=255)
                                                                                                                                        aCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg                     <  2:31692/139‑1 (MQ=255)
                                                                                                                                        aCAGGGATCCACGCAATCACCTGTGGGAGCCGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg                     >  1:31692/1‑139 (MQ=255)
                                                                                                                                                  cACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGAcc  <  7:179580/149‑1 (MQ=255)
                                                                                                                                                  cACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGAc   <  1:59726/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACC  >  NC_002947/1570293‑1570585

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