Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	697,097	A→G	intergenic (+21/‑730)	mmsA‑I → / → PP_16SD	methylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	697,097	0	A	G	84.6%	24.7 / ‑2.7	13	intergenic (+21/‑730)	mmsA‑I/PP_16SD	methylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/0);  total (12/1)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGAGGCC‑GTCA‑‑GGC‑C‑TCCCCGC‑TTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCG  >  NC_002947/696968‑697219                                                                                                                                   |                                                                                                                                 cTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑Ggggggggggg                                                                                                                            >  1:1896596/1‑132 (MQ=255)      gCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGC‑GGCacgg                                                                                                                      >  1:1987656/1‑133 (MQ=255)              ggTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGC‑GGGA‑‑GGG‑C‑Gcccccg                                                                                                          >  1:106793/1‑134 (MQ=255)                           ataCCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGAGGCC‑GGCA‑‑GGC‑C‑TCCCCGC‑GTTGTTTTTgggg                                                                                            >  2:1235832/1‑136 (MQ=255)                             aCCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCG‑GGCG‑‑GGC‑C‑GCCCCCC‑GGTGTTTTTggggg                                                                                           >  1:146006/1‑136 (MQ=255)                                    gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCG‑GGCAGCGGC‑C‑CCCCCGG‑TGTGTTTTTggggtggtgg                                                                                      >  1:1187819/1‑130 (MQ=255)                                    gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCC‑GGCG‑‑GGC‑G‑CCCCCCC‑GGTGTTTTTTGGgggtggttt                                                                                    >  1:1202372/1‑135 (MQ=255)                                    gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGC‑GGGA‑‑CGC‑C‑GCCCCCCGTTTTTTTTTTGGggttgtttt                                                                                    >  1:2039618/1‑135 (MQ=255)                                    gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGCCGGCCG‑‑GGC‑G‑CCCCCCC‑ATGGTTTTTTGGggtgggttt                                                                                    >  1:1604936/1‑135 (MQ=255)                                    gCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGCCGGCCC‑‑GGCGC‑CCCCCGC‑TGT‑TTTTTGGTggttggttt                                                                                    >  1:2442037/1‑135 (MQ=255)                                               aCCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGCGGGGGGGGCC‑GGCG‑‑GCC‑C‑CCCCCCG‑GTGGGTTTTGGGGGTGGGTT‑TTGGTGGTgggg                                                                        >  1:2739203/1‑136 (MQ=255)                                                cccAGCGTTGGCCGCAGCGGGCCAGCCATGCAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGC‑GGCC‑‑GGC‑CGAGCCCGG‑TTGTTGTTTGGGGGGGGGTTGTTGGGGTGgggg                                                                        >  2:1507288/1‑136 (MQ=255)                                                cccAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGGGGGC‑GGCA‑‑GGG‑C‑GCCCCCC‑GTTGTTTTTTGGGGGTGGTTGTTTGGGTTggggg                                                                       >  2:1585208/1‑135 (MQ=255)                                                                                                                     ataTGAAGC‑GGGGAGGCC‑GTCA‑‑GGC‑C‑TCCCCGC‑TTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCg  <  2:2130889/134‑1 (MQ=255)                                                                                                                                   |                                                                                                                                 CTATGGCCCGGATGGTGTGCGCTTCTATACCCGTCGCAAGGCCATTACCCAGCGTTGGCCGCAGCGGGCCAGCCATGAAGCGTCGCAGTTTGCATTCCCTAGCCTGTAAGGGAAGTAGATGAAGC‑GGGGAGGCC‑GTCA‑‑GGC‑C‑TCCCCGC‑TTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCG  >  NC_002947/696968‑697219

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