Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,456,162	A→G	E126E (GAA→GAG)	mpl →	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,456,162	0	A	G	100.0%	207.0 / NA	58	E126E (GAA→GAG)	mpl	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/24);  total (34/24)

ACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGC  >  NC_000913/4456027‑4456294                                                                                                                                        |                                                                                                                                     aCTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAgt                                                                                                                                  <  1:521209/140‑1 (MQ=255)    ggAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGtgtg                                                                                                                               <  1:541657/140‑1 (MQ=255)     gAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTgg                                                                                                                              >  1:274703/1‑140 (MQ=255)          aaacatCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACa                                                                                                                         >  2:919130/5‑140 (MQ=255)               tCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCg                                                                                                                    >  2:612166/1‑140 (MQ=255)                  cTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGc                                                                                                                 <  2:274703/140‑1 (MQ=255)                   ttATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCt                                                                                                                <  1:843552/140‑1 (MQ=255)                    tataTGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCtt                                                                                                               >  2:760798/1‑140 (MQ=255)                             ggTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATcgg                                                                                                      >  2:929867/1‑140 (MQ=255)                             ggTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATcgg                                                                                                      >  2:498012/1‑140 (MQ=255)                             ggTCCACAGTAGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATcgg                                                                                                      >  2:490628/1‑140 (MQ=255)                              gTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATcggc                                                                                                     >  1:917647/1‑140 (MQ=255)                              gTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATcggc                                                                                                     >  2:339040/1‑140 (MQ=255)                               tCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATcggcg                                                                                                    <  2:917647/140‑1 (MQ=255)                                        gCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCgg                                                                                           >  2:583876/1‑140 (MQ=255)                                        gCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCgg                                                                                           >  2:744709/1‑140 (MQ=255)                                         cTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGTATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCggg                                                                                          >  1:898088/1‑140 (MQ=255)                                          tGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCgggg                                                                                         >  2:355277/1‑140 (MQ=255)                                           gCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGTATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGa                                                                                        <  2:898088/140‑1 (MQ=255)                                           gCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGa                                                                                        >  1:318120/1‑140 (MQ=255)                                           gCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAACGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGa                                                                                        <  2:846268/140‑1 (MQ=255)                                            cACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGaa                                                                                       >  1:769865/1‑140 (MQ=255)                                               gATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTATAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACtt                                                                                    <  1:243228/140‑1 (MQ=255)                                                     gtgCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGt                                                                              <  2:116695/140‑1 (MQ=255)                                                     gtgCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGt                                                                              <  1:339040/140‑1 (MQ=255)                                                      tgctgcGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGtt                                                                             >  2:608129/1‑140 (MQ=255)                                                      tgctgcGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGtt                                                                             >  1:758208/1‑140 (MQ=255)                                                          gcgcGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCgg                                                                         <  1:919130/140‑1 (MQ=255)                                                          gcgcGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCgg                                                                         <  1:583876/140‑1 (MQ=255)                                                           cgcgACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGc                                                                        >  1:546302/1‑140 (MQ=255)                                                           cgcgACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGc                                                                        >  2:283511/1‑140 (MQ=255)                                                             cgACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTc                                                                      >  2:58275/1‑140 (MQ=255)                                                                      ggTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCg                                                             <  1:760798/140‑1 (MQ=255)                                                                       gTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGa                                                            >  1:908683/1‑140 (MQ=255)                                                                           tGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGc                                                        <  2:6810/140‑1 (MQ=255)                                                                            ggCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCg                                                       >  2:900054/1‑140 (MQ=255)                                                                                  tGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGActtct                                                 >  1:950612/1‑140 (MQ=255)                                                                                            cacGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGa                                       <  1:608129/140‑1 (MQ=255)                                                                                             acGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGaa                                      <  1:58275/140‑1 (MQ=255)                                                                                             acGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGaa                                      <  2:546302/140‑1 (MQ=255)                                                                                               ggCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGc                                    >  1:89971/1‑140 (MQ=255)                                                                                                 cAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCgg                                  <  2:769865/140‑1 (MQ=255)                                                                                                 cAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCgg                                  <  1:744709/140‑1 (MQ=255)                                                                                                  aaaaCCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGa                                 >  2:100638/1‑140 (MQ=255)                                                                                                  aaaaCCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGa                                 >  1:571098/1‑140 (MQ=255)                                                                                                   aaaCCACCACCGCGGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGAt                                <  2:571098/140‑1 (MQ=255)                                                                                                               cgGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTgc                    >  1:845888/1‑140 (MQ=255)                                                                                                               cgGGAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTgc                    >  1:893442/1‑140 (MQ=255)                                                                                                                 ggAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTgcg                   >  1:962167/1‑139 (MQ=255)                                                                                                                 ggAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTgcg                   <  2:962167/139‑1 (MQ=255)                                                                                                                  gAATGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGcc                 <  1:498012/140‑1 (MQ=255)                                                                                                                    aTGGCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCctt               >  2:223470/1‑140 (MQ=255)                                                                                                                      ggCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCcttct             >  2:548356/1‑140 (MQ=255)                                                                                                                       gCGACCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCcttctt            >  2:312111/1‑140 (MQ=255)                                                                                                                          aCCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTgcg                   <  1:42024/130‑1 (MQ=255)                                                                                                                          aCCTGGATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCAGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTgcg                   >  2:42024/1‑130 (MQ=255)                                                                                                                              ggATTCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACaaa     <  2:950612/140‑1 (MQ=255)                                                                                                                                 ttCTGGAGCAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGc  >  1:275091/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                     ACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTCTGGAACAGTGTGGTTACAAACCGGGCTTTGTAATCGGCGGTGTGCCGGGGAACTTTGAGGTTTCGGCTCATCTGGGCGAAAGCGACTTCTTTGTTATCGAAGCGGATGAGTATGACTGCGCCTTCTTCGACAAACGC  >  NC_000913/4456027‑4456294

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