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	100.0%	66.5 / NA	23	intergenic (+83/‑81)	trpS/zapE	tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (14/9);  total (14/9)

AACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGG  >  NC_002947/1499394‑1499647                                                                                                                    |                                                                                                                                               aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                                                                >  5:214545/1‑147 (MQ=21)  aCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGGGTGGCCGTTGT‑TCGCGGGCGCGCCCGGGCCCCCGGACTGTATGCGATTAAGCgggg                                                                                                               >  6:57985/1‑146 (MQ=16)                     gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCTTGGCACCGGCGTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACg                                                                                              <  5:237849/145‑1 (MQ=21)                        ggACTGAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCg                                                                                         <  3:228914/147‑1 (MQ=21)                                gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCCCGGGGCATGGCACCCGCGTTGCCGGGGT‑TCGCCGGCACGCCCGCACCCACCGACCGGTTGCCCTTTAGCCGGGCTATCGCTTCAAGGTCGGCCCCGGTTttgt                                                                                 >  6:234599/1‑147 (MQ=9)                                          ggTGGGAGCGGGGGTGCCCGGAAAACTGGCGCGGGGGTGCTTGGAACGGGCTTTGCCGGTGTCC‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTa                                                                      <  5:165867/147‑1 (MQ=25)                                           gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCTGTGCATGGCACGGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTaa                                                                     <  1:37185/147‑1 (MQ=21)                                                     gCGTGCCCGGGTAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCGGGTGTGC‑GCGGGCACGCCCGCTCGCACAGACTGATTGCGATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGa                                                           <  2:74627/147‑1 (MQ=25)                                                      cGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGAc                                                          >  8:239299/1‑147 (MQ=25)                                                         gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAATCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACtac                                                       >  8:142525/1‑145 (MQ=25)                                                         gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTcc                                                       >  4:106497/1‑147 (MQ=37)                                                            cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTccc                                                      >  5:113856/1‑145 (MQ=37)                                                            cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTccc                                                      >  5:18978/1‑145 (MQ=37)                                                            cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTccc                                                      >  2:162662/1‑145 (MQ=37)                                                                       cgccgcGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACCGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTAt                                          >  6:13352/1‑146 (MQ=37)                                                                                      caccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAAc                          <  1:163909/144‑1 (MQ=38)                                                                                       accGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACTATACATCACTCCCTTAGAACGCTATCAAGCTGATCTGAAACg                         >  2:116843/3‑147 (MQ=37)                                                                                        ccGGCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGt                        >  7:9432/2‑147 (MQ=39)                                                                                          ggCTTTGCCGGTGT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTcc                      >  1:252510/1‑147 (MQ=39)                                                                                             tttGCCGGTGT‑TCGCGGGCAAGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCTGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCg                    <  1:74893/147‑1 (MQ=38)                                                                                                     gtgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACttcttc            <  2:252510/140‑1 (MQ=255)                                                                                                      tgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTcc           >  8:21250/7‑147 (MQ=255)                                                                                                               gCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcgg  <  2:38253/147‑1 (MQ=255)                                                                                                                    |                                                                                                                                               AACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTC‑GCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGG  >  NC_002947/1499394‑1499647

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