Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,697,352	A→G	N518D (AAT→GAT)	mltF →	membrane‑bound lytic transglycosylase F, murein hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,697,352	0	A	G	100.0%	149.1 / NA	39	N518D (AAT→GAT)	mltF	membrane‑bound lytic transglycosylase F, murein hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/18);  total (21/18)

ATCCGGCGGTATCGCCTACGGAGTTGGGCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTATTCCTTCCGTAATTTCCACTCGGTGATTCATACCCG  >  NC_000913/2697218‑2697487                                                                                                                                       |                                                                                                                                        aTCCGGCGGTATCGCCTACGGAGTTGGGCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTa                                                                                                                                     <  1:859075/139‑1 (MQ=255) aTCCGGCGGTATCGCCTACGGAGTTGGGCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTa                                                                                                                                     >  2:859411/1‑139 (MQ=255)       cggTATCGCCTACGGAGTTGGGCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGt                                                                                                                               >  2:599918/1‑139 (MQ=255)                      aGTTGGGCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCttt                                                                                                                >  2:426698/1‑139 (MQ=255)                        ttGGGCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCttttt                                                                                                              >  1:823062/1‑139 (MQ=255)                         tGGGCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCtttttt                                                                                                             <  2:8811/139‑1 (MQ=255)                           ggCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCt                                                                                                           >  2:503999/1‑139 (MQ=255)                                              ttCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTggcg                                                                                        >  1:107695/1‑139 (MQ=255)                                                                 gtcgtcATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCa                                                                     >  1:138172/1‑139 (MQ=255)                                                                       atcaAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGc                                                               >  1:594925/1‑139 (MQ=255)                                                                                     cccATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCttt                                                                                                    >  2:792194/1‑88 (MQ=255)                                                                                     cccATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCttt                                                                                                    <  1:792194/88‑1 (MQ=255)                                                                                      ccATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCa                                                <  2:594925/139‑1 (MQ=255)                                                                                       cATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCAt                                               <  2:138172/139‑1 (MQ=255)                                                                                          tctcCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTgcg                                                                                                        >  2:869403/1‑79 (MQ=255)                                                                                          tctcCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTgcg                                                                                                        <  1:869403/79‑1 (MQ=255)                                                                                                tctctGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGAtt                                                                                                                       <  1:486748/58‑1 (MQ=255)                                                                                                tctctGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGAtt                                                                                                                       >  2:486748/1‑58 (MQ=255)                                                                                                   ctgctgTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTAttc                                   >  1:212784/1‑139 (MQ=255)                                                                                                        gTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                                                                                     >  2:856018/1‑84 (MQ=255)                                                                                                        gTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGgcggcg                                                                                     <  1:856018/84‑1 (MQ=255)                                                                                                        gTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTAttccttcc                              >  2:777665/1‑139 (MQ=255)                                                                                                        gTTTTCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTAttccttcc                              <  2:908290/139‑1 (MQ=255)                                                                                                         ttttCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCa                                                                                   <  2:579460/85‑1 (MQ=255)                                                                                                         ttttCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCa                                                                                   >  1:579460/1‑85 (MQ=255)                                                                                                          tttCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCg                                                                               <  2:264214/88‑1 (MQ=255)                                                                                                          tttCCAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCg                                                                               >  1:264214/1‑88 (MQ=255)                                                                                                              cAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAAc                                                            >  2:135349/1‑103 (MQ=255)                                                                                                              cAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAAc                                                            <  1:135349/103‑1 (MQ=255)                                                                                                              cAGGAAAGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTATTCCTTCCGTAAtt                        >  1:402866/1‑139 (MQ=255)                                                                                                                  aaaGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTATTCCTTCCGTAATTTCCa                    <  2:402866/139‑1 (MQ=255)                                                                                                                   aaGGGAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTATTCCTTCCGTAATTTCCAc                   <  1:503999/139‑1 (MQ=255)                                                                                                                     gggAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCg                                                                               >  2:776675/1‑77 (MQ=255)                                                                                                                     gggAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCg                                                                               <  1:776675/77‑1 (MQ=255)                                                                                                                      ggAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTATTCCTTCCGTAATTTCCACTca                <  1:426698/139‑2 (MQ=255)                                                                                                                       gAGTGAAGAGAAACAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTATTCCTTCCGTAATTTCCACTCgg               <  2:823062/139‑1 (MQ=255)                                                                                                                                  aaCAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATc                                              >  1:819982/1‑97 (MQ=255)                                                                                                                                  aaCAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATc                                              <  2:819982/97‑1 (MQ=255)                                                                                                                                    cAAGATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTATTCCTTCCGTAATTTCCACTCGGTGATTCATACCCg  >  2:452500/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                        ATCCGGCGGTATCGCCTACGGAGTTGGGCAAAGAGAAATTTCCTTTTCTCTCGTTTCTTTCCCAGTCGTCATCAAACTATTTGACCCATTCTCCCTCTCTGCTGTTTTCCAGGAAAGGGAGTGAAGAGAAACAAAATTAATCCGTCGAGGATTGCGCTTTTTTCTGCGCTTTAATTTCCTGGCGGCGCATGCGAAAGAAGTCACTGAGCAACGCCGCGCACTCATCCGCCAGTATTCCTTCCGTAATTTCCACTCGGTGATTCATACCCG  >  NC_000913/2697218‑2697487

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