Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,592,341 T→G 10.3% L185R (CTG→CGG) ‡ 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,3410TG10.3% 88.1 / 5.3 39L185R (CTG→CGG) ‡mccAmethylcrotonyl‑CoA carboxylase biotin‑containing subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (17/18);  new base G (2/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

TACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCAACGAACGCGACTGT  >  NC_002947/4592199‑4592486
                                                                                                                                              |                                                                                                                                                 
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                                                                                                              <  1:564551/111‑1 (MQ=255)
                                                                     ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGc                                                                                                              >  2:564551/1‑111 (MQ=255)
                                                                     ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATg                                                                                                  >  2:340745/1‑123 (MQ=255)
                                                                     ctCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATg                                                                                                  <  1:340745/123‑1 (MQ=255)
                                                                      tCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGccagtt                                                                                                                                     <  1:556920/87‑5 (MQ=37)
                                                                      tCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGAGGCCAGGGCGTCGGccagtt                                                                                                                                     >  2:556920/1‑83 (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                          >  2:555287/1‑150 (MQ=255)
                                                                                                                  gAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGc                          >  1:129626/1‑150 (MQ=255)
                                                                                                                    ggAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg                                                     <  1:516953/121‑1 (MQ=255)
                                                                                                                    ggAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAg                                                     >  2:516953/1‑121 (MQ=255)
                                                                                                                        gaaAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTAc                    <  1:555287/150‑1 (MQ=255)
                                                                                                                                  tGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCaacgaac          <  2:491639/150‑1 (MQ=255)
                                                                                                                                          gCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCAACGAACGCGACTGt  <  2:572896/150‑1 (MQ=255)
                                                                                                                                           cccTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGgaga                                                          >  1:331670/1‑93 (MQ=255)
                                                                                                                                           cccTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGgaga                                                          <  2:331670/93‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TACCACGGCGAGGCCCAGGACCTGGACACCTTCCGCGCTGCGGCCGAACGCATCGGCTACCCGGTGCTGCTCAAGGCCAGCGCGGGAGGCGGTGGCAAGGGCATGAAAGTGGTCGAGGAAGAAAGCCAACTGGCCGACGCCCTGGCCTCGGCCCAGCGTGAGGCACAGTCGTCGTTCGGCGATGCGCGCATGCTGGTGGAAAAGTACGTGCTCAAGCCACGCCACGTGGAGATTCAGGTGTTCGCCGATCAGCACGGCAACTGCCTGTACCTCAACGAACGCGACTGT  >  NC_002947/4592199‑4592486

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: