Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	97,947	T→C	R287R (CGT→CGC)	murD →	UDP‑N‑acetylmuramoyl‑L‑alanine:D‑glutamate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	97,947	0	T	C	100.0%	118.0 / NA	36	R287R (CGT→CGC)	murD	UDP‑N‑acetylmuramoyl‑L‑alanine:D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/19);  total (17/19)

CACCTGAATCATCAGCAGGGCGAAACCTGGCTGCGGGTTAAAGGCGAGAAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGTGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAGCGGCGCTGAATGGCCTGCACGT  >  NC_000913/97801‑98093                                                                                                                                                   |                                                                                                                                                   cACCTGAATCATCAGCAGGGCGAAACCTGGCTGCGGGTTAAAGGCGAGAAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGcc                                                                                                                                                 >  1:110202/1‑150 (MQ=255)              agcagGGCGAAACCTGGCTGCGGGTTAAAGGCGAGAAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGaaa                                                                                                                                     >  1:103816/1‑149 (MQ=255)              agcagGGCGAAACCTGGCTGCGGGTTAAAGGCGAGAAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAg                                                                                                                                    >  1:265657/1‑150 (MQ=255)                     cGAAACCTGGCTGCGGGTTAAAGGCGAGAAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAAcc                                                                                                                            >  1:118863/1‑151 (MQ=255)                      gAAACCTGGCTGCGGGTTAAAGGCGAGAAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAAcc                                                                                                                            >  2:767877/1‑150 (MQ=255)                              gCTGCGGGTTAAAGGCGAGAAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACt                                                                                                                   >  1:475493/1‑151 (MQ=255)                                              gagaAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCt                                                                                                   >  1:696972/1‑151 (MQ=255)                                                         gAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTtgtg                                                                                        <  1:636565/151‑1 (MQ=255)                                                         gAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTtgtg                                                                                        <  2:671242/151‑1 (MQ=255)                                                               gAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGa                                                                                   <  2:406232/150‑1 (MQ=255)                                                                aaaGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAg                                                                                  <  1:740464/150‑1 (MQ=255)                                                                  agagATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTgg                                                                                    >  1:177503/1‑146 (MQ=255)                                                                      aTGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACg                                                                           <  1:335219/151‑1 (MQ=255)                                                                      aTGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACg                                                                           <  1:599960/151‑1 (MQ=255)                                                                          aaCTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACgg                                                                          >  1:526877/1‑148 (MQ=255)                                                                          aaCTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACgg                                                                          <  2:526929/148‑1 (MQ=255)                                                                          aaCTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGt                                                                       >  2:392536/1‑151 (MQ=255)                                                                                ccGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGtt                                                                  <  1:262142/150‑1 (MQ=255)                                                                                           aaCTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGAtt                                                      >  1:223041/1‑151 (MQ=255)                                                                                             cTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCg                                                    >  1:654635/1‑151 (MQ=255)                                                                                                cacCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGaa                                                  >  2:444254/1‑150 (MQ=255)                                                                                                    aaTGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCg                                              >  2:307390/1‑150 (MQ=255)                                                                                                       gcgctggcggcgctggcgCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGAcc                                           <  1:721910/150‑1 (MQ=255)                                                                                                       gcgctggcggcgctggcgCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGAcc                                           >  2:337613/1‑150 (MQ=255)                                                                                                       gcgctggcggcgctggcgCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCa                                          <  1:767802/151‑1 (MQ=255)                                                                                                        cgctggcggcgctggcgCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCaa                                         <  1:135814/151‑1 (MQ=255)                                                                                                         gctggcggcgctggcgCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCaa                                         <  1:756747/150‑1 (MQ=255)                                                                                                           tggcggcgctggcgCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGt                                      >  2:713037/1‑151 (MQ=255)                                                                                                            ggcggcgctggcgCTGGCAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACg                                       <  1:90785/149‑1 (MQ=255)                                                                                                                         ctggcAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGa                          >  1:51884/1‑149 (MQ=255)                                                                                                                           ggcAGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAg                        <  1:712963/149‑1 (MQ=255)                                                                                                                              aGATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAgcggcg                   <  1:741488/151‑1 (MQ=255)                                                                                                                               gATGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAGCGgcgc                  <  2:289838/151‑1 (MQ=255)                                                                                                                                aTGCTGCAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAGCGGCGCt                 <  2:449575/151‑1 (MQ=255)                                                                                                                                     gcAGGGTTACCGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAGCGGCGCTGAAt             <  1:644086/150‑1 (MQ=255)                                                                                                                                              ccGCGCGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAGCGGCGCTGAATGGCCTGc      <  1:75924/148‑1 (MQ=255)                                                                                                                                               cgcgcgCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAGCGGCGCTGAATGGCCTGCACGt  <  2:417673/151‑1 (MQ=255)                                                                                                                                               cgcgcgCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAGCGGCGCTGAATGGCCTGCACGt  <  2:736760/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CACCTGAATCATCAGCAGGGCGAAACCTGGCTGCGGGTTAAAGGCGAGAAAGTGCTGAATGTGAAAGAGATGAAACTTTCCGGGCAGCATAACTACACCAATGCGCTGGCGGCGCTGGCGCTGGCAGATGCTGCAGGGTTACCGCGTGCCAGCAGCCTGAAAGCGTTAACCACATTCACTGGTCTGCCGCATCGCTTTGAAGTTGTGCTGGAGCATAACGGCGTACGTTGGATTAACGATTCGAAAGCGACCAACGTCGGCAGTACGGAAGCGGCGCTGAATGGCCTGCACGT  >  NC_000913/97801‑98093

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