Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	732,498	G→A	19.2%	V679M (GTG→ATG) ‡	clpB →	chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	732,498	0	G	A	19.2%	54.9 / 7.4	26	V679M (GTG→ATG) ‡	clpB	chaperone protein ClpB
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base A (2/3);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGC  >  NC_002947/732380‑732631                                                                                                                       |                                                                                                                                      gaAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                                                      <  1:348743/136‑1 (MQ=255) gaAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                                                      <  2:15544/136‑1 (MQ=255)                       attggTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGtgctgctgc                                                                                                <  2:508314/133‑1 (MQ=255)                               accaccAGGCTATGTAGGGTATGAAGAGGGTGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTg                                                                                        >  2:112659/1‑136 (MQ=255)                               accaccAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTg                                                                                        >  2:394523/1‑136 (MQ=255)                               accaccAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTg                                                                                        >  2:110769/1‑136 (MQ=255)                                         tATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgcc                                                                              >  2:276681/1‑136 (MQ=255)                                             tcggctATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCt                                                                          >  1:374781/6‑136 (MQ=255)                                              aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTg                                                                         >  1:134095/1‑136 (MQ=255)                                              aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTg                                                                         >  2:281694/1‑136 (MQ=255)                                               ggctATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGa                                                                        <  2:502383/133‑1 (MQ=255)                                                ggTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                                                       <  1:110769/136‑1 (MQ=255)                                                ggTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGGTCACCGGGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                                                       <  1:276681/136‑1 (MQ=255)                                                      aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAg                                                                 <  2:511257/136‑1 (MQ=255)                                                      aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGCTGCAGGTGCTGGAGGACGGCCGCCTGACCGACAg                                                                 <  2:335251/136‑1 (MQ=255)                                                              ggTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGgc                                                         >  1:474384/1‑136 (MQ=255)                                                              ggTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGgc                                                         >  1:426130/1‑136 (MQ=255)                                                                               cGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGc                                        >  1:186546/1‑136 (MQ=255)                                                                                    ggcgCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAacac                                   >  1:412922/1‑136 (MQ=255)                                                                                                  tgggtggtGCTGCTGGACGTGGTTGAAAAAGCCTACCCGTATGTGTTCAACGTGCTTTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGtgttga                     <  2:412922/134‑1 (MQ=255)                                                                                                  tCGGTTGTGCTGATGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGtgatga                     <  1:486396/136‑1 (MQ=255)                                                                                                  caggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGa                                                                        <  1:470901/83‑1 (MQ=255)                                                                                                  caggtgatCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGa                                                                        >  2:470901/3‑85 (MQ=255)                                                                                                        gtgCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCa               >  1:447534/1‑136 (MQ=255)                                                                                                         tgctgcTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCaa              >  2:459689/1‑136 (MQ=255)                                                                                                                     aGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGgc  <  2:134095/136‑1 (MQ=255)                                                                                                                       |                                                                                                                                      GAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGACCTCCAACCTGGGCTCGGC  >  NC_002947/732380‑732631

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑