Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,590,047 C→T 7.7% G243G (GGC→GGT)  mccB → methylcrotonyl‑CoA carboxylase biotin‑containing subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,590,0470CT 7.7% 149.1 / 2.5 52G243G (GGC→GGT) mccBmethylcrotonyl‑CoA carboxylase biotin‑containing subunit beta
Reads supporting (aligned to +/- strand):  ref base C (21/27);  new base T (2/2);  total (23/29)
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.72e-01

ACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGCCGATGAGC  >  NC_002947/4589899‑4590183
                                                                                                                                                     |                                                                                                                                         
aCGTGCCGGCCATGGCCGCCGCAGCGATCATGGTGCGCCAGCAGGCGACCAGCTTCCTCGCCGGCCCG‑CCGCGGGTGAAAGCCGCCCCCGGCGAAGTGGTGAGCGGCGAAGACCTCGGCGGCGCCGATGGGCATGGCCGCCCCCgcggc                                                                                                                                           >  6:370746/1‑149 (MQ=255)
aCGTGCCGGCCATGGCCGACGAAGCGATAATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCAACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCCCACgcggc                                                                                                                                           >  7:334276/1‑149 (MQ=255)
 cGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAgcggc                                                                                                                                           <  6:139935/148‑1 (MQ=255)
 cGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAgcggc                                                                                                                                           <  2:33646/148‑1 (MQ=255)
                  acgaAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGACAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGcc                                                                                                                         <  5:370746/149‑1 (MQ=255)
                    gaAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGa                                                                                                                       <  2:62451/149‑1 (MQ=255)
                      aGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCGACCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGa                                                                                                                       >  1:81957/1‑148 (MQ=255)
                         gATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAacg                                                                                                                  <  2:349090/149‑1 (MQ=255)
                          aTCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAacga                                                                                                                 >  2:48061/1‑149 (MQ=255)
                                  gcgcCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACg                                                                                                        >  5:328525/1‑149 (MQ=255)
                                        gcagGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTgg                                                                                                  <  3:216966/149‑1 (MQ=255)
                                              gACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGAcc                                                                                                                                 <  4:293834/113‑1 (MQ=255)
                                              gACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGAcc                                                                                                                                 >  3:293834/1‑113 (MQ=255)
                                               aCCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCg                                                                                                                        <  1:115397/121‑1 (MQ=255)
                                               aCCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCg                                                                                                                        >  2:115397/1‑121 (MQ=255)
                                                ccATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATcgcg                                                                                          <  5:329286/149‑1 (MQ=255)
                                                       ccTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGa                                                                                       <  8:92462/145‑1 (MQ=255)
                                                         tCGCCGGCCCG‑CCGCTGTTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAgt                                                                                 >  8:162792/1‑149 (MQ=255)
                                                         tCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGAcc                                                                                                                                 >  5:382932/1‑102 (MQ=255)
                                                         tCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGAcc                                                                                                                                 <  6:382932/102‑1 (MQ=255)
                                                         tCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAgt                                                                                 >  3:368865/1‑149 (MQ=255)
                                                               gCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGc                                                                            >  3:30460/1‑148 (MQ=255)
                                                                  cg‑ccgcTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAAc                                                                        >  6:180712/1‑149 (MQ=255)
                                                                  cg‑ccgcTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAAc                                                                        >  6:35146/1‑149 (MQ=255)
                                                                   g‑ccgcTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAAcc                                                                       <  7:111322/149‑1 (MQ=255)
                                                                       gcTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTc                                                                     >  2:11333/1‑148 (MQ=255)
                                                                           gTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTg                                                                <  5:320177/149‑1 (MQ=255)
                                                                              aaaGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGACACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTgg                                                               <  1:270207/148‑1 (MQ=255)
                                                                              aaaGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGc                                                              <  4:30460/148‑1 (MQ=255)
                                                                               aaGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGcac                                                            <  5:35146/149‑1 (MQ=255)
                                                                                aGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCGGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCCACCTCAACTGGccc                                                            >  4:255278/1‑146 (MQ=255)
                                                                                  ccgccACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAg                                                         <  6:204333/149‑1 (MQ=255)
                                                                                  ccgccACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAg                                                         <  2:360102/149‑1 (MQ=255)
                                                                                       aCCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCgg                                                    <  6:198555/149‑1 (MQ=255)
                                                                                                  ggtgAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCACCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGc                                         >  6:288789/1‑149 (MQ=255)
                                                                                                      aGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGtt                                     >  7:161192/1‑149 (MQ=255)
                                                                                                      aGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCACCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGtt                                     <  5:288789/149‑1 (MQ=255)
                                                                                                         gcCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGc                                  <  5:180712/149‑1 (MQ=255)
                                                                                                                aCCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTg                           <  6:328525/149‑1 (MQ=255)
                                                                                                                           gcCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACAA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTa                >  1:350618/1‑149 (MQ=255)
                                                                                                                              gATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGc             <  2:320245/149‑1 (MQ=255)
                                                                                                                               aTGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCg            >  5:189764/1‑149 (MQ=255)
                                                                                                                               aTGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCg            <  7:162792/149‑1 (MQ=255)
                                                                                                                               aTGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACAA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCg            <  2:350618/149‑1 (MQ=255)
                                                                                                                                tgtgCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCgg           >  4:317488/1‑149 (MQ=255)
                                                                                                                                  tgCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTTGCACCAGTGGCACCGCTTTATGCGGc          <  1:11333/148‑1 (MQ=255)
                                                                                                                                    cATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGCCg        <  4:368865/148‑1 (MQ=255)
                                                                                                                                    cATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGCCGa       >  4:288606/1‑149 (MQ=255)
                                                                                                                                     aTTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCACCTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCCCCGCTTTATGCGGCCGAt      >  8:161657/1‑149 (MQ=255)
                                                                                                                                         ccGCACCAGCGGTGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGc          <  7:193410/141‑1 (MQ=255)
                                                                                                                                         ccGCACCAGCGGTGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGc          <  7:198698/141‑1 (MQ=255)
                                                                                                                                         ccGCACCAGCGGTGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGc          >  8:193410/1‑141 (MQ=255)
                                                                                                                                         ccGCACCAGCGGTGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGc          >  8:198698/1‑141 (MQ=255)
                                                                                                                                         ccGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGCCGAt      >  7:213055/1‑145 (MQ=255)
                                                                                                                                         ccGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGCCGAt      <  8:213055/145‑1 (MQ=255)
                                                                                                                                         ccGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGCCGATGAGc  <  7:45226/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
ACGTGCCGGCCATGGCCGACGAAGCGATCATGGTGCGCCAGCAGGCGACCATCTTCCTCGCCGGCCCG‑CCGCTGGTGAAAGCCGCCACCGGCGAAGTGGTGAGCGCCGAAGACCTCGGCGGCGCCGATGTGCATTGCCGCACCAGCGGCGTCGCCGACCATTATGCCGACAACGACGA‑ACACGCCCTGGCCATCGCGCGACGCAGTGTGGCCAACCTCAACTGGCACAAGCTCGGCAAGCTGCAGCGGTTGGCACCAGTGGCACCGCTTTATGCGGCCGATGAGC  >  NC_002947/4589899‑4590183

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