Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,591,633	A→G	T9T (ACA→ACG)	dapE →	succinyl‑diaminopimelate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,591,633	0	A	G	100.0%	88.7 / NA	26	T9T (ACA→ACG)	dapE	succinyl‑diaminopimelate desuccinylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/12);  total (14/12)

AATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGC  >  NC_000913/2591501‑2591769                                                                                                                                     |                                                                                                                                         aaTGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCGATGAGGTGGAGTATATGTCGTGCCCGGTTAATGAGCTGACGCAACAGc                                                                                                                                   >  2:168649/1‑140 (MQ=255)            ttATCAAACGTCCATTGGTCTGAGTGCCCGGGAAGGCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTcgccgc                                                                                                                        <  2:240713/140‑1 (MQ=255)               tCAAACGTCCATTGCTATGCGTCCCCGGTATGCGTTTGTGGCTGGGTTGCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGTTGCCGTAACAGCTTATTCGCCGcccg                                                                                                                     <  2:327467/140‑2 (MQ=255)                   aCGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCt                                                                                                                     <  1:270303/136‑1 (MQ=255)                   aCGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCt                                                                                                                     >  2:270303/1‑136 (MQ=255)                    cGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCt                                                                                                                >  2:374264/1‑140 (MQ=255)                      tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGa                                                                                                              <  2:606571/140‑1 (MQ=255)                      tCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGa                                                                                                              >  1:490519/1‑140 (MQ=255)                            gCTCTGCGTGCCCGGTAAGCCTTTGCTGCCGGGTTTCGGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtg                                                                                                        <  1:215471/140‑1 (MQ=255)                             ctctGCGTGCCCGGTTAGCCTAAGCTGCTTGGTTTCAATGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtga                                                                                                       <  1:244828/140‑1 (MQ=255)                             ctctGCGTGCCCGGTAAGCCTACGCTGCTTGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCtga                                                                                                       <  1:290617/140‑1 (MQ=255)                               ctGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCt                                                                                                                >  1:585574/1‑129 (MQ=255)                               ctGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCt                                                                                                                <  2:585574/129‑1 (MQ=255)                                    tGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCa                                                                                                <  1:374264/140‑1 (MQ=255)                                                   tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAg                                                                                       >  1:103319/1‑134 (MQ=255)                                                   tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAg                                                                                       <  2:103319/134‑1 (MQ=255)                                                   tgctgcTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTttg                                                                                 >  1:357482/1‑140 (MQ=255)                                                        ctgGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGttgat                                                                            >  2:415203/1‑140 (MQ=255)                                                                                              tGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCa                                      <  1:330797/140‑1 (MQ=255)                                                                                                  gtgtAGTCTATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGa                                  >  2:441102/1‑140 (MQ=255)                                                                                                          tATGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCg                          >  1:88499/1‑140 (MQ=255)                                                                                                            tGTCGTGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGAt                        >  1:425685/1‑140 (MQ=255)                                                                                                                 tGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGTCGATACGCa                   >  1:396822/1‑140 (MQ=255)                                                                                                                 tGCCCGGTTATTGAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCa                   >  2:67420/1‑140 (MQ=255)                                                                                                                             gAGCTGACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGc       <  1:441102/140‑1 (MQ=255)                                                                                                                             acagtGACGCAACAGCTTATTCCCCTCCCTTCCCTGACTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGc       <  1:183226/136‑1 (MQ=255)                                                                                                                                  gACGCAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGc  >  1:63231/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                         AATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGGTTTCAGTGATTCCAGTTATCAGCAATTTTTCCATGAGGTGTAGTCTATGTCGTGCCCGGTTATTGAGCTGACACAACAGCTTATTCGCCGCCCTTCCCTGAGTCCTGATGATGCAGGATGCCAGGCTTTGTTGATTGAACGTTTGCAGGCGATCGGTTTTACCGTTGAACGCATGGACTTTGCCGATACGCAGAATTTTTGGGCATGGC  >  NC_000913/2591501‑2591769

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