Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,995,572:1 +G 9.8% coding (1188/1272 nt) bkdB → branched‑chain alpha‑keto acid dehydrogenase complex lipoamide acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,995,5721.G 9.8% 150.4 / 11.3 41coding (1188/1272 nt)bkdBbranched‑chain alpha‑keto acid dehydrogenase complex lipoamide acyltransferase
Reads supporting (aligned to +/- strand):  ref base . (20/17);  new base G (2/2);  total (22/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 = 9.54e-01

CGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTAT  >  NC_002947/4995445‑4995715
                                                                                                                                |                                                                                                                                               
cGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATgccg                                                                                                                           <  1:565725/150‑1 (MQ=255)
        gtCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTc                                                                                                                   >  1:288457/1‑150 (MQ=255)
           aaCACCCCGGAAGTGGCAATCGTCTGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATc                                                                                                                >  2:316659/1‑150 (MQ=255)
            ccacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATcc                                                                                                               >  2:546929/2‑150 (MQ=255)
             cacCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCa                                                                                                              <  2:570653/150‑1 (MQ=255)
              acCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAg                                                                                                             >  2:797145/1‑150 (MQ=255)
              acCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAg                                                                                                             <  1:300773/150‑1 (MQ=255)
                  cGGCAGTGGCACTAGTCGGGGTAAACCGCATGGTCGAACGGCCAGGGGTGATGTACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCg                                                                                                         >  2:729925/1‑150 (MQ=255)
                            aaTCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCt                                                                                               >  2:200083/1‑150 (MQ=255)
                            aaTCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCt                                                                                               >  2:199718/1‑150 (MQ=255)
                                gtcgGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTc                                                                                           <  2:625020/150‑1 (MQ=255)
                                           ccGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCgc                                                                                >  2:434111/1‑150 (MQ=255)
                                                       aCGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTgg                                                                                                  <  2:519948/120‑1 (MQ=255)
                                                       aCGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTgg                                                                                                  >  1:519948/1‑120 (MQ=255)
                                                           ccAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAg                                                                <  1:199718/150‑1 (MQ=255)
                                                           ccAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAg                                                                <  1:200083/150‑1 (MQ=255)
                                                                  tgATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGt                                                                    >  1:160066/1‑139 (MQ=255)
                                                                  tgATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGt                                                                    <  2:160066/139‑1 (MQ=255)
                                                                    aTCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGc                                                       >  2:457555/1‑150 (MQ=255)
                                                                                        gcgcAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTccc                                   <  1:797145/150‑1 (MQ=255)
                                                                                            aaGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGcc                                                                                                >  2:212691/1‑85 (MQ=255)
                                                                                            aaGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGcc                                                                                                <  1:212691/85‑1 (MQ=255)
                                                                                             aGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGtt                              <  2:623923/150‑1 (MQ=255)
                                                                                                 gatgaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatc                          >  1:66651/1‑150 (MQ=255)
                                                                                                 gatgaACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatc                          <  2:66651/150‑1 (MQ=255)
                                                                                                          gTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaacaac                                                                                    <  1:550668/83‑1 (MQ=255)
                                                                                                          gTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGaacaac                                                                                    >  2:550668/1‑83 (MQ=255)
                                                                                                             cAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGATCATGCAACAGATTATCCAGACTACCCTGTTGATCATCggcggcggc              >  2:417940/1‑150 (MQ=255)
                                                                                                             cAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCggcggcggc              >  2:633328/1‑150 (MQ=255)
                                                                                                                 tCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                   <  1:335385/93‑1 (MQ=255)
                                                                                                                 tCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                   >  2:335385/1‑93 (MQ=255)
                                                                                                                   gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                   <  1:673271/91‑1 (MQ=255)
                                                                                                                   gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                   >  2:673271/1‑91 (MQ=255)
                                                                                                                   gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggc        >  2:319446/1‑150 (MQ=255)
                                                                                                                   catCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatcatc                       <  1:437345/133‑1 (MQ=255)
                                                                                                                   catCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatcatc                       >  2:437345/3‑135 (MQ=255)
                                                                                                                     tCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCAt                                                         >  1:200837/1‑99 (MQ=255)
                                                                                                                     tCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCAt                                                         <  2:200837/99‑1 (MQ=255)
                                                                                                                         ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCAt                                                         <  2:704238/95‑1 (MQ=255)
                                                                                                                         ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCAt                                                         >  1:704238/1‑95 (MQ=255)
                                                                                                                         ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCACCGAGGTCGTACGTGACCAGTGTGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTAt  <  1:729925/150‑1 (MQ=255)
                                                                                                                                |                                                                                                                                               
CGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTAT  >  NC_002947/4995445‑4995715

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