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%	89.2 / NA	33	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/19);  total (14/19)

GCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑‑TCGGGCACG‑CCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCA  >  NC_002947/1499367‑1499639                                                                                                                                                 |                                                                                                                                       gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCATCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTcc                                                                                                                                 >  1:74017/1‑149 (MQ=25) gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACCGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCCCAA‑CAGTCGCCGCGGTGCAAGGCCCCCGCGTTGCCGGTGT‑TC‑GCGGGCAcc‑cccccgcc                                                                                                                                >  5:289016/1‑143 (MQ=255)    cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TA‑GCGGGCACGCCCGCTCCca                                                                                                                               >  2:100375/1‑148 (MQ=25)     gTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCcaca                                                                                                                             <  5:1435/149‑1 (MQ=21)         tGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGAc                                                                                                                          <  1:308303/149‑1 (MQ=21)               aaTTCATAGCCTTAACCAAGCCTTGGAACCACACCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCTAAA‑AAGGCGCCGCGGTGCAGGGCACCGGCTTGGCCGGTTTGC‑‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGc                                                                                                                  <  2:11437/148‑1 (MQ=18)                        ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCg                                                                                                         >  8:89542/1‑149 (MQ=21)                            aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGc                                                                                                      <  7:89542/148‑1 (MQ=21)                               gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCCCAA‑CACGCGCCCCCGTGCAATGCACCCGCGTTGGCGGTGT‑TA‑GCGGGCCCGCCCCCACCCCCAGACAGATTTCCCATTAACCGGGcta                                                                                                    >  8:15027/1‑146 (MQ=255)                                aaGCCTTGCCACCACAGCTGGACTTCGGCTCCCAAACGGTGGGAGCGGGCGTCCCCGCGAAA‑CAGGCGCCGGGGTGCATGGCACGGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATc                                                                                                 <  7:324099/149‑1 (MQ=21)                                     ttGCCACCACCGCCGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGAc                                                                                                                          <  3:69026/120‑1 (MQ=14)                                        gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TCGTCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACa                                                                                           <  1:100375/148‑1 (MQ=25)                                          accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCAGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                                       >  7:204099/1‑149 (MQ=21)                                                         aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                          <  7:359346/147‑1 (MQ=21)                                                           gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGTCTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                          <  1:91240/145‑1 (MQ=21)                                                           gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                          >  2:91240/1‑145 (MQ=21)                                                            cTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAA‑CCCGCGCCGCGGTGCCTGGCACCGGGGTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc                                                                     <  5:367781/149‑1 (MQ=21)                                                                  aaCGGTGGGAGCGGGCGGGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGc                                                                                   <  3:353454/129‑1 (MQ=18)                                                                   aCGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACGGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACa                                                                                           <  5:71002/120‑1 (MQ=17)                                                                      gTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TA‑ACGGGCACGCCCGCGCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc                                                            >  7:347154/1‑148 (MQ=25)                                                                                     cccGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGAACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                             <  2:312896/148‑1 (MQ=37)                                                                                      ccGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTAt                                                                                                  >  6:214450/1‑94 (MQ=17)                                                                                      ccGCGAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTAt                                                                                                  <  5:214450/94‑1 (MQ=17)                                                                                       cgcgAAA‑CAGGCGCCGCGGTGCCTGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                                                            >  7:246076/1‑115 (MQ=21)                                                                                       cgcgAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                                                            <  8:246076/115‑1 (MQ=21)                                                                                        gcgAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                       <  4:66978/119‑1 (MQ=21)                                                                                        gcgAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                       >  3:66978/1‑119 (MQ=21)                                                                                         cgAAA‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                        <  7:15027/149‑1 (MQ=37)                                                                                            aa‑CAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc                                     >  6:303859/1‑149 (MQ=37)                                                                                                       gcgGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                 >  7:150151/1‑143 (MQ=37)                                                                                                       gcgGTGCATGGCACCGGCTTTGCCGGTGT‑TC‑GCGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                 <  8:150151/143‑1 (MQ=37)                                                                                                                                  gtgttcgcGGGCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa  >  8:101755/8‑149 (MQ=255)                                                                                                                                          gggCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                             <  4:71484/114‑1 (MQ=255)                                                                                                                                          gggCACGCCCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc                             >  3:71484/1‑114 (MQ=255)                                                                                                                                                 |                                                                                                                                       GCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAA‑CAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGT‑T‑‑TCGGGCACG‑CCGCTCCCACAGACTGA‑TTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCA  >  NC_002947/1499367‑1499639

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