Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	102,228	Δ1 bp	coding (1464/1476 nt)	murC →	UDP‑N‑acetylmuramate:L‑alanine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	102,227	0	A	.	100.0%	159.4 / NA	33	coding (1463/1476 nt)	murC	UDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (19/14);  total (19/14)

TGGTGCCGGATCCGGCGCGGGTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAGAACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTCGACCCGAAAGAAGTCG  >  NC_000913/102080‑102368                                                                                                                                                    |                                                                                                                                              tGGTGCCGGATCCGGCGCGGGTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGgaagaca                                                                                                                                             >  1:163649/1‑148 (MQ=255)  ggTGCCGGATCCGGCGCGGGGAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTTAAGCCGCAAACTCCGGAGgaagacaac                                                                                                                                           <  2:37378/151‑5 (MQ=255)   gTGCCGGATCCGGCGCGGGTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGACTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑acaaca                                                                                                                                         <  1:288345/151‑1 (MQ=255)       cGGATCCGGCGCGGGTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGAc                                                                                                                                     >  1:573536/1‑151 (MQ=255)          aTCCGGCGCGGGTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTg                                                                                                                                   >  2:54206/1‑150 (MQ=255)             cGGCGCGGGTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAGCTCCGGAGGAAG‑ACAACATGACTGATaa                                                                                                                               <  2:784438/151‑1 (MQ=255)                 gcgGGTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATc                                                                                                                           >  1:715622/1‑151 (MQ=255)                    ggTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATcgc                                                                                                                         <  2:552262/150‑1 (MQ=255)                              aTGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTtgg                                                                                                              >  2:367588/1‑151 (MQ=255)                                gCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTtgggt                                                                                                            >  1:685381/1‑151 (MQ=255)                                          gTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCg                                                                                                  <  2:790998/151‑1 (MQ=255)                                          gTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCg                                                                                                  >  2:412006/1‑151 (MQ=255)                                                         gACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGatt                                                                                   >  2:247195/1‑151 (MQ=255)                                                             tGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTctct                                                                               >  2:744729/1‑151 (MQ=255)                                                                   tcGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAAtt                                                                          >  2:390151/1‑150 (MQ=255)                                                                    cGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCt                                                                        <  2:214404/151‑1 (MQ=255)                                                                     gTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTg                                                                       <  1:159434/151‑1 (MQ=255)                                                                                  gTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGtt                                                          <  1:459189/151‑1 (MQ=255)                                                                                   tAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGgtgt                                                           >  2:393596/1‑149 (MQ=255)                                                                                       atTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCg                                                     >  1:252333/1‑151 (MQ=255)                                                                                            aaaaaTTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTg                                                >  2:835682/1‑151 (MQ=255)                                                                                                      cGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGc                                       >  2:491066/1‑150 (MQ=255)                                                                                                       gTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCgg                                     <  2:421093/151‑1 (MQ=255)                                                                                                        ttCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCgg                                     >  2:486755/1‑150 (MQ=255)                                                                                                            ttAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTAtt                                 <  1:25356/150‑1 (MQ=255)                                                                                                                 tGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGAcgcg                           >  2:622179/1‑151 (MQ=255)                                                                                                                  gAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGAcgc                            <  1:177047/149‑1 (MQ=255)                                                                                                                  gAAATCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGt                          >  1:826986/1‑151 (MQ=255)                                                                                                                     aTCAAACTGAAGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATc                       <  1:367680/151‑1 (MQ=255)                                                                                                                          aCTGAAGCCGCAAACGCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTc                  <  2:230201/151‑1 (MQ=255)                                                                                                                              aaGCCGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTCGAcc              >  2:272109/1‑151 (MQ=255)                                                                                                                                  cGCAAACTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTCGACCCGa           >  1:537719/1‑150 (MQ=255)                                                                                                                                        cTCCGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTCGACCCGAaagaag     <  2:104218/150‑1 (MQ=255)                                                                                                                                          ccGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTCGACCCGAAAGAAGt    <  1:577079/149‑1 (MQ=255)                                                                                                                                          ccGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTCGACCCGAAAGAAGTc   <  2:801287/150‑1 (MQ=255)                                                                                                                                          ccGGAGGAAG‑ACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTCGACCCGAAAGAAGTCg  >  1:661147/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                              TGGTGCCGGATCCGGCGCGGGTAGCCGAGATGCTGGCACCGGTATTAACCGGTAACGACCTGATTCTCGTTCAGGGGGCTGGTAATATTGGAAAAATTGCCCGTTCTTTAGCTGAAATCAAACTGAAGCCGCAAACTCCGGAGGAAGAACAACATGACTGATAAAATCGCGGTCCTGTTGGGTGGGACCTCCGCTGAGCGGGAAGTTTCTCTGAATTCTGGCGCAGCGGTGTTAGCCGGACTGCGTGAAGGCGGTATTGACGCGTATCCTGTCGACCCGAAAGAAGTCG  >  NC_000913/102080‑102368

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