Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,499,508	+C	intergenic (+85/‑79)	trpS → / → zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,499,506	1	.	C	94.4%	48.8 / ‑3.3	18	intergenic (+83/‑81)	trpS/zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (8/9);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCG  >  NC_002947/1499369‑1499646                                                                                                                                              |                                                                                                                                               ccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCGGCCCCCGCAACACGCGCCGCCGGGCAAGGCACCCGCGTTTGCCGGGTT‑G‑GCCGGCAcg‑c‑ccgcgac                                                                                                                                       >  8:12122/1‑141 (MQ=255)              aTTCAGAGCCTGAACGAAGCCTTGGCACCACCGCAGGCCTGCGGCGGCCAAACGGGGGGCGCGGGCGGGCCCGCGAAACAGGCGCCGCGGGGCAGGGC‑ACCGGCTTTGCCGGTGTGC‑GCGGGCACGCC‑CGCTCCCACAGACTGATTGc                                                                                                                          <  3:77382/148‑1 (MQ=21)                        tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTACCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                                <  4:141722/148‑1 (MQ=21)                                        accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTg                                                                                                <  6:28800/148‑1 (MQ=21)                                                       aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTttg                                                                                >  3:131519/1‑148 (MQ=21)                                                        ggCTGCCAAACGGGGGGCGCGGGCGGGCCCGCTAAACTGGCGCCGCGGGGCAGGGC‑ACCGGGTTTGCCGGTGTGC‑GCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTttgt                                                                               <  5:121023/148‑1 (MQ=25)                                                           tGCCAAACGGTGGGAGCGGGCGTGGCCCCGCAACCCGCGCCCCGGTGGATGGC‑ACCCGCGTTGCCCGTGT‑TAGCGGGCCCGCC‑CCCGCCCCCACACGGATTGGCAATAAGCCGGCCGATAGCGACAGGGCCCGCCCCGGTTTttgtgg                                                                             >  2:20509/1‑146 (MQ=255)                                                                          gcgggcggGCCCGCGAAACAGGCGCCGCGGTGCTTGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGAACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                             <  2:118341/148‑1 (MQ=25)                                                                             ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                                                          <  2:11781/148‑1 (MQ=37)                                                                                    ccGCGAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                     >  2:155036/1‑146 (MQ=37)                                                                                    ccGCGAAACAGGCGCCGCGGTGCATGGC‑ACAGGCTTTGCGGGTGT‑TCGCGGGCACGCC‑CGATCCCACAGACTGATTGCCATTAAGCTGGCCTATCGCTACAGTT‑ACGGCCCGTTTTTTGTTGCCTCTCTAACGAATCATGACTCCCt                                                    <  4:135919/147‑1 (MQ=25)                                                                                     cgcgAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                     >  7:31576/1‑145 (MQ=37)                                                                                     cgcgAAACAGGCGCCGCGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                     >  5:89203/1‑145 (MQ=37)                                                                                                 gccgcGGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                                      >  1:37324/1‑148 (MQ=37)                                                                                                   cgcgGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                                    <  4:131519/148‑1 (MQ=37)                                                                                                    gcgGTGCATGGC‑ACCGGCTTTGCCGGTGT‑TCGCGGGCACG‑CACGCTACCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTAGCTAACGAATCATGACTCCCTTAGAACGCTATCaa                                      >  2:104701/1‑145 (MQ=37)                                                                                                                           ccGGTGT‑TCGCGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct             >  4:83942/1‑148 (MQ=255)                                                                                                                                     gcGGGCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcg    >  3:172901/2‑148 (MQ=255)                                                                                                                                       tggCACGCC‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcgcg  <  5:4016/147‑1 (MQ=255)                                                                                                                                              |                                                                                                                                               CCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGC‑CACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑C‑CGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTG‑ACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCG  >  NC_002947/1499369‑1499646

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