Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,592,336 C→G 13.6% D183E (GAC→GAG 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,592,3360CG13.6% 94.9 / 9.2 37D183E (GAC→GAGmccAmethylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base G (2/3);  total (18/19)
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.82e-01

GCCGGGCTACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCAACGAAC  >  NC_002947/4592192‑4592478
                                                                                                                                                |                                                                                                                                              
gCCGGGCTACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGccag                                                                                                                                           <  2:310411/150‑3 (MQ=255)
       tACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGTGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCg                                                                                                                                    >  1:589454/1‑150 (MQ=255)
                gAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGt                                                                                                                           <  1:277225/150‑1 (MQ=255)
                        ggACCGGGACACCGTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGcaca                                                                                                                   <  2:229313/150‑1 (MQ=255)
                            cTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAgtcg                                                                                                               >  1:314315/1‑150 (MQ=255)
                                  acCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTg                                                                                                                                          <  2:8246/117‑1 (MQ=255)
                                  acCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTg                                                                                                                                          >  1:8246/1‑117 (MQ=255)
                                                  gcGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGc                                                                                         <  2:213403/149‑1 (MQ=255)
                                                    gAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCtg                                                                                       >  2:92697/1‑150 (MQ=255)
                                                        gCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCtggtgg                                                                                   <  2:589454/150‑1 (MQ=255)
                                                          aTCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGaa                                                                                 >  1:443199/1‑150 (MQ=255)
                                                                tACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATAAAAGTGGTCGAGGAACAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTAc                                                                           <  2:531680/150‑1 (MQ=255)
                                                                 aCCCGGTGCTGCTCGAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACg                                                                          >  1:228677/1‑150 (MQ=255)
                                                                  cccGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGt                                                                         >  1:572896/1‑150 (MQ=255)
                                                                            ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGc                                                                                                      >  2:564551/1‑111 (MQ=255)
                                                                            ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGc                                                                                                      <  1:564551/111‑1 (MQ=255)
                                                                            ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATg                                                                                          <  1:340745/123‑1 (MQ=255)
                                                                            ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATg                                                                                          >  2:340745/1‑123 (MQ=255)
                                                                             tCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGccagtt                                                                                                                             >  2:556920/1‑83 (MQ=37)
                                                                             tCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGccagtt                                                                                                                             <  1:556920/87‑5 (MQ=37)
                                                                             tCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAgccac                                                              <  1:92697/150‑1 (MQ=255)
                                                                                    cAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGCGGGCTAGAGCCAAATGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGt                                                       >  2:21578/1‑150 (MQ=255)
                                                                                     aGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGcca                                                                                                                                >  2:190391/1‑75 (MQ=37)
                                                                                     aGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGcca                                                                                                                                <  1:190391/76‑2 (MQ=37)
                                                                                     aGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTg                                                      <  2:228677/150‑1 (MQ=255)
                                                                                              ggCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg                                             <  1:20032/150‑1 (MQ=255)
                                                                                                  gTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGgtgt                                         >  2:75045/1‑150 (MQ=255)
                                                                                                           gCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCtggt                                                                                     >  2:272500/1‑97 (MQ=255)
                                                                                                           gCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCtggt                                                                                     <  1:272500/97‑1 (MQ=255)
                                                                                                           gCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg                                             >  1:591102/1‑137 (MQ=255)
                                                                                                           gCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg                                             <  2:591102/137‑1 (MQ=255)
                                                                                                                         gAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGc                  >  1:129626/1‑150 (MQ=255)
                                                                                                                         gAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGc                  >  2:555287/1‑150 (MQ=255)
                                                                                                                           ggAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg                                             >  2:516953/1‑121 (MQ=255)
                                                                                                                           ggAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg                                             <  1:516953/121‑1 (MQ=255)
                                                                                                                               gaaAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTAc            <  1:555287/150‑1 (MQ=255)
                                                                                                                                         tGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCaacgaac  <  2:491639/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GCCGGGCTACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCAACGAAC  >  NC_002947/4592192‑4592478

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