Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,992,076 A→T 25.0% intergenic (‑183/‑16) bkdR ← / → bkdAA Bkd operon transcriptional regulator/branched‑chain alpha‑keto acid dehydrogenase complex subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,992,0760AT25.0% 19.8 / 6.4 16intergenic (‑183/‑16)bkdR/bkdAABkd operon transcriptional regulator/branched‑chain alpha‑keto acid dehydrogenase complex subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base T (2/2);  total (9/7)
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.13e-01

CTCAGCAAAGGCAATCAACTTGAGAGAAAAATTCTCCTGCCGTGCCACTAAGATGTAGGTGACGCTGACTTACCAGTCACAAGCCGGTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCAGGCTGCCAGACCGATTTTTCCTACCTGCGCCTCAACGATGCAGGTCAAGCCCGTAAACCCGCGATCGA  >  NC_002947/4991945‑4992213
                                                                                                                                   |                                                                                                                                         
cTCAGCAAAGGCAATCAACTTGAGAGAAAAATTCTCCTGCCGTGCCACTAAGATGTAGGTGACGCTGACTTACCAGTCACAAGCCGGTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATg                                                                                                                         <  2:285334/150‑1 (MQ=255)
   aGCAAAGGCAATCAACTTGAGAGAAAAATTCTCCTGCCGTGCCACTAAGATGTAGGTGACGCTGACTTACCAGTCACAAGCCGGTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAAc                                                                                                                      >  1:549893/1‑150 (MQ=255)
                     gagagaAAAATTCTCCTGCCGTGCCACTAAGATGTAGGTGACGCTGACTTACCAGTCACAAGCCGGTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGt                                                                                                    <  2:549893/150‑1 (MQ=255)
                                  tcCTGCCGTGCCACTAAGATGTAGGTGACGCTGACTTACCAGTCACAAGCCGGTACTCAGCGGCGGCCGCTGCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAACGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCg                                                                                       >  1:140828/1‑150 (MQ=255)
                                                       tAGGTGACGCTGACTTACCAGTCACAAGCCGGTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCAGGCt                                                                  <  2:140828/150‑1 (MQ=255)
                                                                              acaAGCCGGTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCAGGCTGCCAGACCGATTTTTCCTACCTg                                           >  1:364475/1‑150 (MQ=255)
                                                                               caAGCCGGTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCCGGCTGCCAGACCGATTTTTCCTACCTgc                                          >  2:111034/1‑150 (MQ=255)
                                                                                      gTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCAGGCTGCCAGACCGATTTTTCCTACCTGCGCCTCaa                                   <  2:515674/150‑1 (MQ=255)
                                                                                       tACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGTGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCAtgt                                                                                            >  1:158792/1‑92 (MQ=255)
                                                                                       tACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGTGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCAtgt                                                                                            >  1:158793/1‑92 (MQ=255)
                                                                                       tACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGTGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCAtgt                                                                                            <  2:158792/92‑1 (MQ=255)
                                                                                       tACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGTGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCAtgt                                                                                            <  2:158793/92‑1 (MQ=255)
                                                                                            aGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCAGGCTGCCAGACCGATTTTTCCTACCTGCGCCTCAACGATGc                             >  1:489889/1‑150 (MQ=255)
                                                                                             gcggagGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTCCGTTTGCATGTGCCCGCGCCCACCGGCCGGCCAGGCTCCCAGACCCATTTTTCCTACCTCCGCCTCAACGATGCa                            <  2:23578/150‑1 (MQ=255)
                                                                                                                      aaaCAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCAGGCTGCCAGACCGATTTTTCCTACCTGCGCCTCAACGATGCAGGTCAAGCCCGTAAACCCGcgatcg   >  1:509932/1‑150 (MQ=255)
                                                                                                                       aaCAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCAGGCTGCCAGACCGATTTTTCCTACCTGCGCCTCAACGATGCAGGTCAAGCCCGTAAACCCGcgatcga  >  1:330410/1‑150 (MQ=255)
                                                                                                                                   |                                                                                                                                         
CTCAGCAAAGGCAATCAACTTGAGAGAAAAATTCTCCTGCCGTGCCACTAAGATGTAGGTGACGCTGACTTACCAGTCACAAGCCGGTACTCAGCGGCGGCCGCTTCAGAGCTCACAAAAACAAATACCCGAGCGAGCGTAAAAAGCATGAACGAGTACGCCCCCCTGCGTTTGCATGTGCCCGAGCCCACCGGCCGGCCAGGCTGCCAGACCGATTTTTCCTACCTGCGCCTCAACGATGCAGGTCAAGCCCGTAAACCCGCGATCGA  >  NC_002947/4991945‑4992213

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