Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,456,022	G→T	A80S (GCG→TCG)	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,022	0	G	T	100.0%	181.8 / NA	52	A80S (GCG→TCG)	mpl	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (27/25);  total (27/25)

GCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTC  >  NC_000913/4455892‑4456157                                                                                                                                   |                                                                                                                                        gccgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTgg                                                                                                                                >  1:454556/1‑140 (MQ=255)  ccgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTATGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGa                                                                                                                               >  1:802323/1‑140 (MQ=255)   cgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGaa                                                                                                                              >  1:426216/1‑140 (MQ=255)    gATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGaaa                                                                                                                             >  2:685500/1‑140 (MQ=255)      tGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGaaaaa                                                                                                                           >  1:738151/1‑140 (MQ=255)        aGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGaaaaaaa                                                                                                                         <  2:176990/140‑1 (MQ=255)          cACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAata                                                                                                                       <  2:426216/140‑1 (MQ=255)           aCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAatat                                                                                                                      <  2:426239/140‑1 (MQ=255)            ccTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATc                                                                                                                     <  1:176934/140‑1 (MQ=255)            ccTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATc                                                                                                                     <  2:146033/140‑1 (MQ=255)                aCTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCtt                                                                                                                 >  1:497502/1‑140 (MQ=255)                      gaAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                                           <  2:738151/140‑1 (MQ=255)                      gaAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                                           <  2:454556/140‑1 (MQ=255)                       aaGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTc                                                                                                          <  2:96833/140‑1 (MQ=255)                          cAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAgg                                                                                                       <  1:685500/140‑1 (MQ=255)                          cAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAACCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAgg                                                                                                       >  2:765283/1‑140 (MQ=255)                                         aTTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCa                                                                                        >  2:725370/1‑140 (MQ=255)                                              ggggTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGAtt                                                                                   >  2:255733/1‑140 (MQ=255)                                               gggTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGAttt                                                                                  >  1:6568/1‑140 (MQ=255)                                               gggTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGAttt                                                                                  >  2:184674/1‑140 (MQ=255)                                                   aCGATGCCAGTCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCa                                                                                        >  1:332577/1‑130 (MQ=255)                                                   aCGATGCCAGTCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCa                                                                                        <  2:332577/130‑1 (MQ=255)                                                       tGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGtgctgc                                                                          >  2:245477/1‑140 (MQ=255)                                                        gccagccaGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTgcg                                                                         <  2:523251/140‑1 (MQ=255)                                                          cagccagcTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGcgcg                                                                       >  1:516331/1‑140 (MQ=255)                                                                gcTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCt                                                                 <  1:47965/140‑1 (MQ=255)                                                                gcTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCt                                                                 <  1:725370/140‑1 (MQ=255)                                                                      gCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGc                                                           >  2:515745/1‑140 (MQ=255)                                                                           aGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGcc                                                      <  1:282758/140‑1 (MQ=255)                                                                                gATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGc                                                 >  2:143150/1‑140 (MQ=255)                                                                                  tCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCg                                               <  2:497502/140‑1 (MQ=255)                                                                                       tGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTaca                                          <  2:6568/140‑1 (MQ=255)                                                                                                gCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCaaa                                 >  1:75684/1‑140 (MQ=255)                                                                                                gCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCaaa                                 >  2:368991/1‑140 (MQ=255)                                                                                                gCAACGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCACTGGGTGCTGGCCGTTGCCGGTACACACGGCaaa                                 >  1:751251/1‑140 (MQ=255)                                                                                                   aCGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAacc                              >  1:651548/1‑140 (MQ=255)                                                                                                    cGCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAacca                             <  1:515745/140‑1 (MQ=255)                                                                                                     gCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAaccac                            <  2:516331/140‑1 (MQ=255)                                                                                                     gCCATGACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAaccac                            <  1:255733/140‑1 (MQ=255)                                                                                                          gACCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACcg                       >  1:356811/1‑140 (MQ=255)                                                                                                           aCCCGTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCACTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACcgc                      <  2:751251/140‑1 (MQ=255)                                                                                                               gTGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGa                  >  1:766003/1‑140 (MQ=255)                                                                                                                tGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGaa                 >  2:510408/1‑140 (MQ=255)                                                                                                                tGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGaa                 >  1:447427/1‑140 (MQ=255)                                                                                                                tGGAAATCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGaa                 >  1:770737/1‑140 (MQ=255)                                                                                                                    aaTCCGTGTGTGGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGc             <  2:447427/140‑1 (MQ=255)                                                                                                                          tgtgtgGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTg       <  1:368991/140‑1 (MQ=255)                                                                                                                          tgtgtgGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTg       <  1:245477/140‑1 (MQ=255)                                                                                                                              tgGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTTGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGAtt   >  2:294449/1‑140 (MQ=255)                                                                                                                              tgGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGAt    <  2:240677/139‑1 (MQ=255)                                                                                                                              tgGAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGAt    >  1:240677/1‑139 (MQ=255)                                                                                                                               ggAATCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTc  <  2:770737/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                        GCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCACCGCGGGAATGGCGACCTGGATTC  >  NC_000913/4455892‑4456157

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: ‑