Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,995,569 Δ1 bp 10.5% coding (1185/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,5690C.10.5% 113.1 / 8.5 38coding (1185/1272 nt)bkdBbranched‑chain alpha‑keto acid dehydrogenase complex lipoamide acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (19/15);  new base . (2/2);  total (21/17)
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.46e-01

GCGCCCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTAT  >  NC_002947/4995421‑4995715
                                                                                                                                                    |                                                                                                                                                   
gcgcCCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGaccacc                                                                                                                                                    <  1:134734/150‑1 (MQ=255)
                        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                                                                                                             <  1:300773/150‑1 (MQ=255)
                                      acCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAg                                                                                                             >  2:797145/1‑150 (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                                                                                                  >  1:519948/1‑120 (MQ=255)
                                                                               aCGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTgg                                                                                                  <  2:519948/120‑1 (MQ=255)
                                                                                   ccAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAg                                                                <  1:200083/150‑1 (MQ=255)
                                                                                   ccAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAg                                                                <  1:199718/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                                                                                                <  1:212691/85‑1 (MQ=255)
                                                                                                                    aaGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGcc                                                                                                >  2:212691/1‑85 (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                                                                   >  2:335385/1‑93 (MQ=255)
                                                                                                                                         tCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                   <  1:335385/93‑1 (MQ=255)
                                                                                                                                           gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                   <  1:673271/91‑1 (MQ=255)
                                                                                                                                           gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTg                                                                   >  2:673271/1‑91 (MQ=255)
                                                                                                                                           gTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTggc        >  2:319446/1‑150 (MQ=255)
                                                                                                                                           catCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatcatc                       >  2:437345/3‑135 (MQ=255)
                                                                                                                                           catCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGatcatc                       <  1:437345/133‑1 (MQ=255)
                                                                                                                                             tCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCAt                                                         <  2:200837/99‑1 (MQ=255)
                                                                                                                                             tCGACCACGCGGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCAt                                                         >  1:200837/1‑99 (MQ=255)
                                                                                                                                                 ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCAt                                                         >  1:704238/1‑95 (MQ=255)
                                                                                                                                                 ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCAt                                                         <  2:704238/95‑1 (MQ=255)
                                                                                                                                                 ccaccGCGTGGTCGATGGCATGGATGCCGCCCTGTTCACCGAGGTCGTACGTGACCAGTGTGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTAt  <  1:729925/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCCCTGGGTGGCATTGTCAGCACGCCGGTGGTCAACACCCCGGAAGTGGCAATCGTCGGGGTCAACCGCATGGTCGAACGGCCAGTGGTGATCGACGGCCAGATCGTCGTGCGCAAGATGATGAACCTGTCCAGCTCGTTCGACCACCGCGTGGTCGATGGCATGGATGCCGCCCTGTTCATCCAGGCCGTACGTGGCCTGCTCGAACAACCCGCCTGCCTGTTCGTGGAGTGAGCATGCAACAGATTATCCAGACTACCCTGTTGATCATCGGCGGCGGCCCTGGCGGCTAT  >  NC_002947/4995421‑4995715

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: