Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	97,137	C→T	R390C (CGC→TGC)	murE →	UDP‑N‑acetylmuramoylalanyl‑D‑glutamate‑2, 6‑diaminopimelate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	97,137	0	C	T	100.0%	133.4 / NA	38	R390C (CGC→TGC)	murE	UDP‑N‑acetylmuramoylalanyl‑D‑glutamate‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/19);  total (19/19)

AAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGAC  >  REL606/97068‑97208                                                                      |                                                                        aaaGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATa                                                                   >  1:1068251/1‑76 (MQ=255)  aaGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATaa                                                                  <  2:302230/76‑1 (MQ=255)  aaGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATaa                                                                  <  1:582951/76‑1 (MQ=255)      cTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGt                                                              >  2:1056334/1‑76 (MQ=255)        tACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTaa                                                            <  2:877347/76‑1 (MQ=255)              cggcgCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTcc                                                      <  1:1054950/76‑1 (MQ=255)              cggcgCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTcc                                                      >  2:697142/1‑76 (MQ=255)                  gcgTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTg                                                  <  1:865217/76‑1 (MQ=255)                       tGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATggg                                             >  2:306509/1‑76 (MQ=255)                        gCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGgc                                            >  1:1261374/1‑76 (MQ=255)                        gCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGgc                                            >  2:1029605/1‑76 (MQ=255)                        gCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGgc                                            >  1:884820/1‑76 (MQ=255)                        gCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGgc                                            >  1:481938/1‑76 (MQ=255)                         cACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGgcg                                           <  2:1238255/76‑1 (MQ=255)                            tgtgCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCaa                                        <  2:554960/76‑1 (MQ=255)                                cGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGc                                    <  1:495926/76‑1 (MQ=255)                                cGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGc                                    <  2:191810/76‑1 (MQ=255)                                  ggCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCg                                  <  2:408805/76‑1 (MQ=255)                                  ggCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCg                                  <  2:324945/76‑1 (MQ=255)                                   gCAAGCTGTGGTGTNTCTTTGGCTGTGGTGGCGATTGCGATATAGGTAAGCGTCCACTGAGGGCCGCAATTGCCGa                                 <  1:1464583/76‑1 (MQ=255)                                     aaGGTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAANTGCCGAag                               <  1:1653403/76‑1 (MQ=255)                                     aaGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAag                               <  2:966193/76‑1 (MQ=255)                                         tgtgGTGTGTTTTTGGCTTAGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGtt                           >  1:1616970/1‑76 (MQ=255)                                         tgtgGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAACGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGtt                           >  1:1440572/1‑76 (MQ=255)                                         tgtgGTGTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGtt                           <  2:960699/76‑1 (MQ=255)                                         tgtgGTGTGTCTTTGGCTGTGGTGGCGACTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGtt                           >  1:360089/1‑76 (MQ=255)                                          gtggtgTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGttt                          >  1:406248/1‑76 (MQ=255)                                           tggtgTGTCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTg                         >  1:1240202/1‑76 (MQ=255)                                              tgtgtCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTg                      >  1:1112105/1‑76 (MQ=255)                                              tgtgtCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTg                      <  1:1056336/76‑1 (MQ=255)                                                  tCTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGt                  <  1:1152246/76‑1 (MQ=255)                                                   cTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTg                 <  2:1529069/76‑1 (MQ=255)                                                   cTTTGGCTGTGGTGGCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTg                 <  2:769838/76‑1 (MQ=255)                                                            tggtggCGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGgtggtg        >  1:368812/1‑76 (MQ=255)                                                            tggtggCGATTGCGATAAACGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCggtgggg        >  1:1446390/1‑74 (MQ=255)                                                              gtggCCATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTgac      >  1:105169/1‑76 (MQ=255)                                                                  cGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTgacggac  >  1:442496/1‑76 (MQ=255)                                                                  cGATTGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTgacggac  >  2:996394/1‑76 (MQ=255)                                                                      |                                                                        AAAGCCTTACAGGCGGCGCGTCTGCACTGTGCGGGCAAGCTGTGGTGTGTCTTTGGCTGTGGTGGCGATCGCGATAAAGGTAAGCGTCCACTGATGGGCGCAATTGCCGAAGAGTTTGCTGACGTGGCGGTGGTGACGGAC  >  REL606/97068‑97208

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