Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,453,278 G→A 8.7% T122M (ACG→ATG)  paaZ ← oxepin‑CoA hydrolase and 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_0009131,453,2780GA 8.7% 135.9 / 3.3 46T122M (ACG→ATG) paaZoxepin‑CoA hydrolase and 3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base A (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGAAACGCTCTTTTTCACTCAGCAGATGTTT  >  NC_000913/1453130‑1453426
                                                                                                                                                    |                                                                                                                                                    
cGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCg                                                                                                                                                      >  5:314772/1‑149 (MQ=255)
        ccAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTc                                                                                                                                               <  6:324573/148‑1 (MQ=255)
         cAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAg                                                                                                                                             >  6:261940/1‑149 (MQ=255)
              tttCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGc                                                                                                                                        >  1:131015/1‑149 (MQ=255)
                 ccAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTc                                                                                                                                      <  6:35704/148‑1 (MQ=255)
                       ttCCCACGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTa                                                                                                                               <  8:282255/149‑1 (MQ=255)
                          cccAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTAcc                                                                                                                             >  3:305222/1‑148 (MQ=255)
                           ccAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCATATCGTCAGGCAGCTCCCGGCTAcc                                                                                                                             >  8:267384/1‑147 (MQ=255)
                           ccAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCATATCGTCAGGCAGCTCCCGGCTAcc                                                                                                                             >  8:272050/1‑147 (MQ=255)
                                agGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTg                                                                                                                      <  2:239821/149‑1 (MQ=255)
                                    gAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCgt                                                                                                                  <  8:250855/149‑1 (MQ=255)
                                    gAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCg                                                                                                                   >  2:72561/1‑148 (MQ=255)
                                     aaGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCgt                                                                                                                  >  5:328444/1‑148 (MQ=255)
                                       gTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCgtaag                                                                                                               >  7:230280/1‑149 (MQ=255)
                                                gTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCATATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACg                                                                                                      <  7:272050/149‑1 (MQ=255)
                                                gTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCATATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACg                                                                                                      <  7:267384/149‑1 (MQ=255)
                                                    atatGCACTGCCAGGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCGACCGAGGCTGGCGTAAGTAAATAACGTcc                                                                                                   <  6:197091/148‑1 (MQ=255)
                                                       tGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCaa                                                                                                <  6:314772/148‑1 (MQ=255)
                                                       tGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAAt                                                                                               <  2:131015/149‑1 (MQ=255)
                                                       tGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAAt                                                                                               >  3:153212/1‑149 (MQ=255)
                                                         cACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGc                                                                                             >  4:16289/1‑149 (MQ=255)
                                                         aaCTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATg                                                                                              <  5:215807/147‑1 (MQ=255)
                                                               cACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCAACTTCTTTCGATAAGGGGCTCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAAAAACGTCCCAATGCCACCtt                                                                                       <  5:4717/149‑1 (MQ=255)
                                                                  gCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCaa                                                                                    <  5:261940/149‑1 (MQ=255)
                                                                            gAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAAccc                                                                           >  3:50127/1‑148 (MQ=255)
                                                                               gTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACt                                                                       >  2:116518/1‑149 (MQ=255)
                                                                                 cAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCAACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTg                                                                      >  7:66413/1‑148 (MQ=255)
                                                                                              gcgGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTCTCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTgc                                                        <  6:328444/149‑1 (MQ=255)
                                                                                              gcgGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTgc                                                        <  4:139454/149‑1 (MQ=255)
                                                                                               cgGCAAATCCACCTTCTTTAGATAAGGGGATCAATTCATCTTCAGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGATgcg                                                       >  5:69186/1‑149 (MQ=255)
                                                                                                     aTCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGtt                                                 >  5:159579/1‑149 (MQ=255)
                                                                                                      tCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTTTCGTCAGGCAGCTCCCGGCTAGCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGttt                                                <  3:59979/149‑1 (MQ=255)
                                                                                                      tCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGttt                                                >  2:183219/1‑149 (MQ=255)
                                                                                                        cACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTgc                                              >  8:282275/1‑149 (MQ=255)
                                                                                                        cACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTCCCCGCGTTGCGCCTGTTTgc                                              >  8:60008/1‑149 (MQ=255)
                                                                                                           cttcttTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCa                                           <  4:50127/149‑1 (MQ=255)
                                                                                                            ttcttTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTgcg                                             >  7:191348/1‑146 (MQ=255)
                                                                                                              ctatCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGa                                         <  2:265572/145‑1 (MQ=255)
                                                                                                                 tCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGaaa                                       >  5:284700/1‑147 (MQ=255)
                                                                                                                   gATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAagag                                    <  4:305222/148‑1 (MQ=255)
                                                                                                                      aaGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATa                                <  7:282275/149‑1 (MQ=255)
                                                                                                                         gggATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGa                              >  6:79645/1‑148 (MQ=255)
                                                                                                                         gggATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGa                              >  1:121935/1‑148 (MQ=255)
                                                                                                                          ggATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGaaa                            <  8:230280/149‑1 (MQ=255)
                                                                                                                                   cATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGAAACGCTCttt                    <  2:121935/148‑1 (MQ=255)
                                                                                                                                       ttCCGGCCACAGCGTATCGTCAGGCCGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAACATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGAAACGCTCTTTTTccc               >  6:68467/1‑147 (MQ=255)
                                                                                                                                                caGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGAAACGCTCTTTTTCACTCAGCAGAt      <  3:16060/149‑1 (MQ=255)
                                                                                                                                                    gTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGAAACGCTCTTTTTCACTCAGCAGATGttt  <  6:156007/149‑1 (MQ=255)
                                                                                                                                                    gTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGAAACGCTCTTTTTCACTCAGCAGATGttt  <  1:72561/149‑1 (MQ=255)
                                                                                                                                                    gTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGAAACGCGCTTTTTCACTCAGCAGATGtt   >  8:281604/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTTGGTGCCAGCTTTTCCAGCATTCCCCAGCAGGGGAAGTTAAAGGCGTTAATATGCACTGCCACGCCTGACTTTGAGGTCAGTAAATGGCGCGCGGCAAATCCACCTTCTTTCGATAAGGGGATCAATTCATCTTCCGGCCACAGCGTATCGTCAGGCAGCTCCCGGCTACCGAGGCTGGCGTAAGTAAATAACGTCCCAATGCCACCTTCAATATCAACCCAACTGTCTGCCCGCGTTGCGCCTGTTTGCGCAGAAAGAGCATAGAAACGCTCTTTTTCACTCAGCAGATGTTT  >  NC_000913/1453130‑1453426

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