Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,593,193 T→C 7.3% D484D (GAT→GAC)  mccA → methylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,593,1930TC 7.3% 157.1 / 2.7 55D484D (GAT→GAC) mccAmethylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base C (2/2);  total (27/28)
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.82e-01

GCTGCAGCGGAGTTGGATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAGCTGCAAGGTGAGCAG  >  NC_002947/4593050‑4593337
                                                                                                                                               |                                                                                                                                                
gctgcAGCGGAGTTGGATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGcc                                                                                                                                              <  8:291017/148‑1 (MQ=255)
 atgcAGCGGAGTTGGATACTGGCTACATTCCACGACACCAGGATGTGCTGCTGCCGGCGCAACGAGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCa                                                                                                                                             >  3:283037/2‑148 (MQ=255)
      gCGGAGTTGGATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGcc                                                                                                                                       <  6:120492/149‑1 (MQ=255)
       cGGAGTTGGATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCg                                                                                                                                      >  5:182422/1‑149 (MQ=255)
       cGGAGTTGGATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCg                                                                                                                                      >  5:177904/1‑149 (MQ=255)
         gAGTTGGATACTGGCTTCATTCCCCGCCAAGAGGTTGTGCTGCTGCCGGCGCCAGGTGCGCTGCCTGCTGCCTTTTGGGAGGCGGCAGCCGATGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTg                                                                                                                                    <  4:283037/149‑1 (MQ=255)
              ggATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGACGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCg                                                                                                                               >  8:155937/1‑149 (MQ=255)
              ggATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCg                                                                                                                               >  8:159764/1‑149 (MQ=255)
               gATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCgg                                                                                                                              >  1:21566/1‑149 (MQ=255)
               gATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCgg                                                                                                                              <  5:163926/149‑1 (MQ=255)
                 tACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGa                                                                                                                            <  2:21566/149‑1 (MQ=255)
                       cTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGTGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATg                                                                                                                      <  1:79182/149‑1 (MQ=255)
                           aTTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCt                                                                                                                  >  2:218187/1‑149 (MQ=255)
                                ccGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCt                                                                                                             >  4:174076/1‑149 (MQ=255)
                                   ccaccaGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTgg                                                                                                          <  1:218187/149‑1 (MQ=255)
                                   acaccaGGATGTGCTGCTGCCGGCGCCACGTGCTCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTgg                                                                                                          >  5:155111/2‑149 (MQ=255)
                                     acccGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGc                                                                                                        <  6:255806/149‑1 (MQ=255)
                                                 tgctgcCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATg                                                                                                                      <  3:210163/123‑1 (MQ=255)
                                                 tgctgcCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATg                                                                                                                      >  4:210163/1‑123 (MQ=255)
                                                     gcCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCAGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACgcag                                                                                          >  6:40425/1‑147 (MQ=255)
                                                      ccGGTGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGcagcag                                                                                       >  1:120416/1‑149 (MQ=255)
                                                      ccGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACgcagc                                                                                         <  8:150174/147‑1 (MQ=255)
                                                       cGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCagcagc                                                                                      <  6:148098/149‑1 (MQ=255)
                                                           gcCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTg                                                                                   >  6:136865/1‑148 (MQ=255)
                                                            ccACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCa                                                                                 >  6:101256/1‑149 (MQ=255)
                                                                          ggCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTg                                                                   >  5:75270/1‑149 (MQ=255)
                                                                          ggCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTg                                                                   <  7:159764/149‑1 (MQ=255)
                                                                                cTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGTGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGc                                                                >  8:226230/1‑146 (MQ=255)
                                                                                   ttGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTg                                                                                                                                    >  7:210160/1‑75 (MQ=255)
                                                                                   ttGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTg                                                                                                                                    <  8:210160/75‑1 (MQ=255)
                                                                                   ttGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGa                                                            <  2:150906/147‑1 (MQ=255)
                                                                                      ggAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTggg                                                                 >  6:198215/1‑139 (MQ=255)
                                                                                      ggAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTggg                                                                 <  5:198215/139‑1 (MQ=255)
                                                                                               aGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTg                                                  >  1:135512/1‑145 (MQ=255)
                                                                                               aGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTg                                                  <  2:135512/145‑1 (MQ=255)
                                                                                                 ccGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTg                                            <  7:226230/149‑1 (MQ=255)
                                                                                                        cTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACg                                       >  4:12607/1‑147 (MQ=255)
                                                                                                         tGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGACCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATc                                                          <  5:52703/127‑1 (MQ=255)
                                                                                                         tGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGACCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATc                                                          >  6:52703/1‑127 (MQ=255)
                                                                                                         tGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGACCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATc                                                          <  5:19325/127‑1 (MQ=255)
                                                                                                         tGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGACCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATc                                                          >  6:19325/1‑127 (MQ=255)
                                                                                                           gCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAg                                    <  6:75270/147‑1 (MQ=255)
                                                                                                                   gCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTgg                                           <  6:2980/132‑1 (MQ=255)
                                                                                                                   gCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTgg                                           >  5:2980/1‑132 (MQ=255)
                                                                                                                     cAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTAc                        >  4:252220/1‑149 (MQ=255)
                                                                                                                     cAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTAc                        <  6:155111/149‑1 (MQ=255)
                                                                                                                               aTCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAGCTg              >  3:105076/1‑149 (MQ=255)
                                                                                                                                  aGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAGCTGCaa           <  3:252220/149‑1 (MQ=255)
                                                                                                                                   gcgcGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAg                 >  6:57120/1‑142 (MQ=255)
                                                                                                                                   gcgcGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAg                 <  5:57120/142‑1 (MQ=255)
                                                                                                                                   gcgcGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAGCTGCAAg          <  4:105076/149‑1 (MQ=255)
                                                                                                                                     gcgAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAGCTGCAAGGt        <  5:101256/149‑1 (MQ=255)
                                                                                                                                       gAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAGCTGCAAGGt        >  2:272433/1‑147 (MQ=255)
                                                                                                                                       gAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAGCTGCAAGGTGa      <  1:272433/149‑1 (MQ=255)
                                                                                                                                           aCGATCGCCATTCGCCGTGGAGCAGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTTGCAGCTGCCAGGtagaagc  <  5:40425/149‑7 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCTGCAGCGGAGTTGGATACTGGCTTCATTCCCCGCCACCAGGATGTGCTGCTGCCGGCGCCACGTGCGCTGCCGGCTGCCTTTTGGGAGGCGGCAGCCGAAGCCTGGCTGCAGGGCCAGGCGGGGCATCAGCGCGAAGACGATCGCCATTCGCCGTGGAGCGGGAGTAATGGCCTGCGCTTGGGCCTACCGGCACGCAGCAGCTTGCATCTGGTGAGCGCTGGGCAAGATCAGGCAGTGGCCCTGGAACGCAGCGCTGCTTCTACCTGGCAGCTGCAAGGTGAGCAG  >  NC_002947/4593050‑4593337

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