Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,994,961	G→T	G67V (GGC→GTC)	dapF →	diaminopimelate epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,994,961	0	G	T	100.0%	148.9 / NA	44	G67V (GGC→GTC)	dapF	diaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/22);  total (22/22)

AATGTCTTTTTTTCACCGGAGCTGATTCGTCGCCTGGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGGCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCACCGATG  >  NC_000913/3994825‑3995092                                                                                                                                         |                                                                                                                                    aaTGTCTTTTTTTCACCGGAGCGGATTCGTCCCCTCGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAg                                                                                                                                  <  2:294070/140‑1 (MQ=255)       tttttttCACCGGAGCTGATTCGTCGCCTGGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGt                                                                                                                            >  1:719761/1‑140 (MQ=255)           tttCACCGGAGCTGATTCGTCGCCTGGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGgcg                                                                                                                        <  2:504774/140‑1 (MQ=255)               aCCGGAGCTGATTCGTCGCCTGGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGt                                                                                                                    >  2:125146/1‑140 (MQ=255)                ccGGAGCTGATTCGTCGCCTGGCTGATCGTCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTg                                                                                                                   <  1:256453/140‑1 (MQ=255)                      cTGATTCGTCGCC‑GGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCAGTATCGCATATTCCATGGTGATGTCAGTGAAGTGGCGCAGTGCGGCAAc                                                                                                            >  2:295612/1‑140 (MQ=255)                              tcgCCTGGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTgcgc                                                                                                     >  1:607945/1‑140 (MQ=255)                              tcgCCTGGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTgcgc                                                                                                     >  1:559163/1‑140 (MQ=255)                               cgCCTGGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTgcgcg                                                                                                    <  2:738804/140‑1 (MQ=255)                                 ccTGGCTGATCGGCACCCGGGGGTAGGGTTTAACCGACTCCTGGTGGTAGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCgct                                                                                                  <  2:215292/140‑1 (MQ=255)                                          tCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTAATGTCAGTGAAGAGGCCGAGTGCGGGAACGGAGCGCGCTACTTTGccc                                                                                         >  2:351012/1‑140 (MQ=255)                                                  tGGGGGTAGGGTTTGTCCAACTTCTTGTGGTTGAGCCCCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTtgtg                                                                                 <  1:295612/140‑1 (MQ=255)                                                   gggggTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGc                                                                                >  1:483528/1‑140 (MQ=255)                                                           ggTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGaaa                                                                        >  2:180287/1‑140 (MQ=255)                                                               tGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGAc                                                                    <  2:122361/140‑1 (MQ=255)                                                                gACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACt                                                                   >  1:403554/1‑140 (MQ=255)                                                                 aCCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTg                                                                  <  2:442161/140‑1 (MQ=255)                                                                   cAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGAc                                                                <  1:420119/140‑1 (MQ=255)                                                                      ctgctgGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGa                                                                          <  2:451345/127‑1 (MQ=255)                                                                      ctgctgGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGa                                                                          >  1:451345/1‑127 (MQ=255)                                                                        gctgGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATa                                                           <  1:180287/140‑1 (MQ=255)                                                                                ttGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGat                                                   >  2:11569/1‑140 (MQ=255)                                                                                 tGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGata                                                  >  1:176131/1‑140 (MQ=255)                                                                                   aGCCGCCGTATGATCCTGAACTGGATTTTCATTATTGCATTTTCAATGTTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATc                                                <  2:280548/140‑1 (MQ=255)                                                                                   aGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATc                                                <  2:559163/140‑1 (MQ=255)                                                                                                   tGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCa                                >  1:145286/1‑140 (MQ=255)                                                                                                   tGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCa                                >  2:31316/1‑140 (MQ=255)                                                                                                     aaCTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAAc                              >  2:581152/1‑140 (MQ=255)                                                                                                       cTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGTGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGGCCAATAAACGTGATCTCCGCGACAGCACCGCCAACgg                            >  1:208009/1‑140 (MQ=255)                                                                                                       cTGGATTTTCACTATCGCATTTTCAATGCTGATGTCAGTCAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACgg                            <  2:232665/140‑1 (MQ=255)                                                                                                          gATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCg                         >  2:180395/1‑140 (MQ=255)                                                                                                          gATTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCg                         >  2:119430/1‑140 (MQ=255)                                                                                                           aTTTTCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCgg                        >  2:653415/1‑140 (MQ=255)                                                                                                            ttttCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGa                       <  1:31316/140‑1 (MQ=255)                                                                                                            ttttCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGa                       <  2:419714/140‑1 (MQ=255)                                                                                                               tCACTATCGCATTTTCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGAGTGACCAATGAGCGTGATATCCGCGTAAGCACCGCCAACGGGCGGATgg                    >  2:259447/1‑140 (MQ=255)                                                                                                                        cATTTTCAATGCTGATGTCAGTTAAGTGGCGTAGTTCGGTAACGGTGCCCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCg           <  2:315863/140‑1 (MQ=255)                                                                                                                            ttCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCAc       <  1:104643/140‑1 (MQ=255)                                                                                                                            ttCAATGCTGATGTCAGTGAAGCGGCGCAGTATGGTAACTGTGCGCGCGGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCAc       <  2:208009/140‑1 (MQ=255)                                                                                                                             tCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCAcc      <  1:125146/140‑1 (MQ=255)                                                                                                                             tCAATGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCAcc      >  1:82237/1‑140 (MQ=255)                                                                                                                               aaTGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCACCga    <  2:176131/140‑1 (MQ=255)                                                                                                                                 tGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCACCgatg  <  2:403554/140‑1 (MQ=255)                                                                                                                                 tGCTGATGTCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCACCgatg  >  1:529404/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                    AATGTCTTTTTTTCACCGGAGCTGATTCGTCGCCTGGCTGATCGGCACCTGGGGGTAGGGTTTGACCAACTGCTGGTGGTTGAGCCGCCGTATGATCCTGAACTGGATTTTCACTATCGCATTTTCAATGCTGATGGCAGTGAAGTGGCGCAGTGCGGCAACGGTGCGCGCTGCTTTGCCCGTTTTGTGCGTCTGAAAGGACTGACCAATAAGCGTGATATCCGCGTCAGCACCGCCAACGGGCGGATGGTTCTGACCGTCACCGATG  >  NC_000913/3994825‑3995092

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