Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,304,640	G→A	11.8%	P887L (CCC→CTC)	recC ←	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,304,640	0	G	A	11.8%	82.6 / 3.6	34	P887L (CCC→CTC)	recC	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit
Reads supporting (aligned to +/- strand):  ref base G (12/18);  new base A (2/2);  total (14/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 = 8.22e-01

GCCGAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGGC  >  NC_002947/5304492‑5304786                                                                                                                                                     |                                                                                                                                                   gCCGAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAggg                                                                                                                                                   >  2:510447/1‑150 (MQ=255)             tAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGAGCCCTTCGTcctc                                                                                                                                       <  2:206222/150‑3 (MQ=255)               gCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAggt                                                                                                                                                   >  2:622401/1‑135 (MQ=255)               gCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAggg                                                                                                                                                   <  1:622401/136‑1 (MQ=255)               gCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCgggg                                                                                                                                     >  1:41772/1‑150 (MQ=255)                       agACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCt                                                                                                                             >  2:166729/1‑150 (MQ=255)                             cAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAgc                                                                                                                       <  1:166729/150‑1 (MQ=255)                               ggCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAgcgc                                                                                                                     <  1:516399/150‑1 (MQ=255)                               ggCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAgcgc                                                                                                                     <  2:356483/150‑1 (MQ=255)                                           tCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTacac                                                                                                         <  2:175075/150‑1 (MQ=255)                                                          cAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTc                                                                                                                                             <  1:334879/99‑1 (MQ=255)                                                          cAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTc                                                                                                                                             >  2:334879/1‑99 (MQ=255)                                                           aaGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGCATCGCGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggctg                                                                                         >  2:244488/1‑150 (MQ=255)                                                                        cGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCa                                                                                                                         >  2:94355/1‑105 (MQ=255)                                                                        cGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCa                                                                                                                         <  1:94355/105‑1 (MQ=255)                                                                                   tCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGgcgc                                                                 <  1:761008/150‑1 (MQ=255)                                                                                       ccAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaaga                                                             >  2:729708/1‑150 (MQ=255)                                                                                         aaGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaagaag                                                           <  2:159447/150‑1 (MQ=255)                                                                                                   cAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCt                                                 <  1:510447/150‑1 (MQ=255)                                                                                                      gCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGgc                                                                   >  2:350007/1‑129 (MQ=255)                                                                                                      gCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGgc                                                                   <  1:350007/129‑1 (MQ=255)                                                                                                         ggCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCa                                           <  2:41772/150‑1 (MQ=255)                                                                                                            tttCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGAc                                        >  1:396655/1‑150 (MQ=255)                                                                                                                aCTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAg                                    <  2:724287/150‑1 (MQ=255)                                                                                                                         gTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGc                           <  2:785497/150‑1 (MQ=255)                                                                                                                          tAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGa                                         <  1:357858/135‑1 (MQ=255)                                                                                                                          tAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGa                                         >  2:357858/1‑135 (MQ=255)                                                                                                                                    gggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGTAGACTCAGGGCTTCAGCGTCCTGATAAg                >  1:563728/1‑150 (MQ=255)                                                                                                                                     ggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggctg                                                                                         <  2:443527/76‑1 (MQ=255)                                                                                                                                     ggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggctg                                                                                         >  1:443527/1‑76 (MQ=255)                                                                                                                                       tATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGTAGACTCAGGGCTTCAGCGTCCTGATAAggcg             <  2:563728/150‑1 (MQ=255)                                                                                                                                         tCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAggcggc           <  2:396655/150‑1 (MQ=255)                                                                                                                                          cccggACGAAGAGCCCTTCGTCCTGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggc                                                                                           <  2:143110/66‑1 (MQ=25)                                                                                                                                          cccggACGAAGAGCCCTTCGTCCTGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggc                                                                                           >  1:143110/4‑69 (MQ=25)                                                                                                                                            aGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAg        >  1:407203/1‑150 (MQ=255)                                                                                                                                             ggACGAAGAGCCCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTTCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGc       >  2:328679/1‑150 (MQ=255)                                                                                                                                                  aaGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGgc  <  1:244488/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCCGAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGGC  >  NC_002947/5304492‑5304786

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: ‑