Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,375,801	T→G	E54D (GAA→GAC)	menC ←	O‑succinylbenzoyl‑CoA synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,375,801	0	T	G	100.0%	103.6 / NA	31	E54D (GAA→GAC)	menC	O‑succinylbenzoyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/14);  total (17/14)

CGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGTTTGAATTTGCTGAAGTCAGGCTGACGTTTCTGG  >  NC_000913/2375573‑2376007                                                                                                                                                                                                                                      |                                                                                                                                                                                                               cGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCctgactga                                                                                                                                                                                                      >  2:511071/1‑239 (MQ=255)     ccTTTCTCGCACTGCATCTCTG‑GATTTTGACGGATCAGCTCGACGGGATCGCCATTCCCCAGCGGTGCCGCACGGGCTTTGGCTGCTTGCGGCAACGTATCTGACAGTTCTGCCAATGCCCAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTGATTTACCCAGGCAAGCAGCACGCGTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGcc                                                                                                                                                                                                  <  2:670393/238‑1 (MQ=255)        ttCTCGCCTGGCATATCTGCAAGTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCgg                                                                                                                                                                                               <  1:879900/239‑1 (MQ=255)          ctcCCCTGGCATATCTGCAAGTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCa                                                                                                                                                                                             <  2:1045353/239‑1 (MQ=255)           tcGCCTGGCATATCTGCAAGTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCa                                                                                                                                                                                             >  2:828043/1‑238 (MQ=255)           tcGCCTGGCATATCTGCAAGTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCa                                                                                                                                                                                             <  1:828043/238‑1 (MQ=255)                   cATATCTGCAAGTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACATCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACATTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTAGGgaca                                                                                                                                                                                    >  1:815100/1‑237 (MQ=255)                             cgTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCTCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCa                                                                                                                                                                          <  2:815100/238‑1 (MQ=255)                                                      tCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCa                                                                                                                                                   <  2:267261/238‑1 (MQ=255)                                                                     ggTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACAt                                                                                                                                          <  2:456557/232‑1 (MQ=255)                                                                       tGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCa                                                                                                                                                                                             <  1:252146/179‑1 (MQ=255)                                                                       tGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCa                                                                                                                                                                                             >  2:252146/1‑179 (MQ=255)                                                                       tGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGcc                                                                                                                                                   >  1:547991/1‑220 (MQ=255)                                                                       tGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCa                                                                                                                                                   <  2:547991/221‑1 (MQ=255)                                                                       tGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCCCTCGCCTTCCAGCAGGCAAACATACAg                                                                                                                                      >  1:195596/1‑234 (MQ=255)                                                                           gCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTcgc                                                                                                                              <  1:511071/238‑1 (MQ=255)                                                                            cACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCCGGCAAAAATAAAGCCCGTAGCgg                                                                                                                            >  1:323466/1‑239 (MQ=255)                                                                                           ctTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTggg                                                                                                                                                                                        >  1:815565/1‑164 (MQ=255)                                                                                           ctTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTggg                                                                                                                                                                                        <  2:815565/164‑1 (MQ=255)                                                                                           ctTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGAt                                                                                                                                                                                   >  2:255557/1‑169 (MQ=255)                                                                                           ctTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGAt                                                                                                                                                                                   <  1:255557/169‑1 (MQ=255)                                                                                                            gTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGGTTTTAACCGACTGCCGCCCAGAACCAccc                                                                                             >  1:44571/1‑238 (MQ=255)                                                                                                                                 cTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCa                                                                       >  1:760143/1‑239 (MQ=255)                                                                                                                                  ttACGCCAAAGGCCACTGAAGGCATCTGCGGTAGCTCGCACTCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGTTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGc                                                                         <  1:842047/236‑1 (MQ=255)                                                                                                                                                       gCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTTTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCCGTATACCTGCCCTCTACGc                                                 >  2:1055287/1‑239 (MQ=255)                                                                                                                                                                                   cagccagTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGTTTGAATTTGCTGaa                     >  1:848088/1‑239 (MQ=255)                                                                                                                                                                                                     ggCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGCCGCGGGGTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCCACGCATTACGGATTCCGTTTGAATTTGCTGAAGTCAGGCTGACggt       >  2:761098/1‑233 (MQ=255)                                                                                                                                                                                                      gCAAGCAGCACACTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGGCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACCCATTACGGACTCCGTTTGAATTTGCTGAAGTCAGGCCGACGTTTCTgg  >  1:121335/1‑239 (MQ=255)                                                                                                                                                                                                              cacaCTTTGCGCCTCTTCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGCATACCTGCGCGCTACGCATTACGGAtt                                       >  1:401110/1‑194 (MQ=255)                                                                                                                                                                                                                   tttGCCCCTCTTCCCAGGTGTCCTGCCGGAACCCCGGCAGGGGGGAGACCCCCCCCCCCCCTTCGCGCCCGCCTCCACGCAGGAAAACATCCAGCCCGTCGGGGTTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGcc                                                                        <  2:456112/156‑1 (MQ=255)                                                                                                                                                                                                                             ttCCCAGGTGTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCCCCCTTCCCCCTCCCCTTCACGCAGGCAAACATACAGCCCGCCCCGGGTTTTTAACCGCCTGTCCCGCAGAACCACCCCCCCCTCCATGGGGATCTGCCCGCGGTATAc                                                              >  1:456112/1‑156 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                               CGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGA‑GGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCCAGTTATTTACCCAGGCAAGCAGCACACTTTGCGCCTCTTCCCAGGTTTCCTGACTGAAGCCCGGCAGTGGGGAGATCTCCCCCCACCCTTCGCGCTCGCCTTCACGCAGGCAAACATACAGCCCGTCGCGGGTTTTTAACCGCCTGTCGCGCAGAACCACCCCCGCGTCCATGGGGATCTGCCAGCGGTATACCTGCGCGCTACGCATTACGGATTCCGTTTGAATTTGCTGAAGTCAGGCTGACGTTTCTGG  >  NC_000913/2375573‑2376007

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑