Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	665,850	G→T	A155D (GCT→GAT)	mrdB ←	cell wall shape‑determining protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	665,850	0	G	T	100.0%	204.3 / NA	53	A155D (GCT→GAT)	mrdB	cell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (24/29);  total (24/29)

ATCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCAGCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAAC  >  NC_000913/665706‑665992                                                                                                                                                 |                                                                                                                                               aTCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATccacca                                                                                                                                          <  1:870993/151‑1 (MQ=255) aTCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATccacca                                                                                                                                          <  2:11201/151‑1 (MQ=255)   cAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACcagc                                                                                                                                        >  2:77739/1‑151 (MQ=255)   cAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACcag                                                                                                                                         <  1:684205/150‑1 (MQ=255)          aaCCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTgg                                                                                                                                 <  2:778104/151‑1 (MQ=255)              acaGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCAt                                                                                                                             <  1:134735/151‑1 (MQ=255)              acaGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCAt                                                                                                                             >  1:280310/1‑151 (MQ=255)                  aaTCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAat                                                                                                                         >  1:143963/1‑151 (MQ=255)                      ggAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCa                                                                                                                     <  1:43428/151‑1 (MQ=255)                           gAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCAcc                                                                                                                <  2:731422/151‑1 (MQ=255)                           gAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCAcc                                                                                                                <  2:64434/151‑1 (MQ=255)                            aaCGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCTGCAGCGTGGGCATAAATATCAGCAcc                                                                                                                >  1:734703/1‑150 (MQ=255)                            aaCGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCa                                                                                                               <  2:746216/151‑1 (MQ=255)                              cgccgcTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAgc                                                                                                             <  1:66005/151‑1 (MQ=255)                               gccgcTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAgcg                                                                                                            <  2:223874/151‑1 (MQ=255)                                ccgcTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAgcgc                                                                                                           <  1:444126/151‑1 (MQ=255)                                        aGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGcc                                                                                                     >  1:326892/1‑149 (MQ=255)                                            ctactGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAat                                                                                                                         >  2:249750/1‑125 (MQ=255)                                            ctactGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAat                                                                                                                         <  1:249709/125‑1 (MQ=255)                                             tactGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGttc                                                                                              >  2:767003/1‑151 (MQ=255)                                               ctGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCGACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGttctt                                                                                            >  2:108264/1‑151 (MQ=255)                                                   gACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAAc                                                                                        >  2:615805/1‑151 (MQ=255)                                                   gACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAAc                                                                                        >  1:8255/1‑151 (MQ=255)                                                              gACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGc                                                                              >  2:878797/1‑150 (MQ=255)                                                              gACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCa                                                                             <  1:466119/151‑1 (MQ=255)                                                                cGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCaa                                                                            >  1:366361/1‑150 (MQ=255)                                                                   cAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGt                                                                        <  1:358233/151‑1 (MQ=255)                                                                    aGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGTATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTc                                                                       <  1:480231/151‑1 (MQ=255)                                                                          ggCCAGAGCGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGt                                                                  >  1:433562/1‑150 (MQ=255)                                                                          ggCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGtt                                                                 >  2:339379/1‑151 (MQ=255)                                                                             cAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGCATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGAt                                                              <  2:280331/151‑1 (MQ=255)                                                                              agagagGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGAt                                                              <  2:205988/150‑1 (MQ=255)                                                                                   ggAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAgcg                                                        >  2:280248/1‑151 (MQ=255)                                                                                     aaCAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAgcgcg                                                      <  1:594852/151‑1 (MQ=255)                                                                                      aCAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAgcgcgc                                                     >  2:684678/1‑151 (MQ=255)                                                                                       cAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCa                                                    <  2:280355/151‑1 (MQ=255)                                                                                        aGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCa                                                    >  2:692911/1‑150 (MQ=255)                                                                                        aGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGCCGCGGTTGATAAAGCGCGCaa                                                   <  2:851814/151‑1 (MQ=255)                                                                                         gTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAAc                                                  >  2:374113/1‑151 (MQ=255)                                                                                            cAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCAt                                               <  1:689198/151‑1 (MQ=255)                                                                                                cAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGt                                           <  1:258773/151‑1 (MQ=255)                                                                                                  gACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTgg                                         >  1:682906/1‑151 (MQ=255)                                                                                                        aaaGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGc                                   >  2:296334/1‑151 (MQ=255)                                                                                                              cAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTAtttt                             <  2:432341/151‑1 (MQ=255)                                                                                                               aaCGAGGATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTg                            >  2:504026/1‑151 (MQ=255)                                                                                                                    ggATTGATGTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAAt                       >  2:70337/1‑151 (MQ=255)                                                                                                                            gTTCCCAGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGAc                <  2:869843/150‑1 (MQ=255)                                                                                                                                  aGGTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGaa         <  1:266703/151‑1 (MQ=255)                                                                                                                                   ggTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGaaa        <  1:266097/151‑1 (MQ=255)                                                                                                                                   ggTCAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGaaa        <  2:267448/151‑1 (MQ=255)                                                                                                                                      cAGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAAc      <  2:639803/150‑1 (MQ=255)                                                                                                                                       aGGCTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAAc      >  1:444651/1‑149 (MQ=255)                                                                                                                                          cTGTGCATCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAAc  >  2:618496/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                               ATCAGGAAGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCAGAGAGGAACAGTACAAACAGACCGGAAAGCGCAACGAGGATTGATGTTCCCAGGTCAGGCTGTGCAGCCACCAGCAGCGTGGGCATAAATATCAGCACCAGCGCGATGCCAGTGTTCTTCAACGATGGCGGGCAAACGTCGCGGTTGATAAAGCGCGCAACCATCAGTGGTACGGCTATTTTGGCAATTTCCGACGGCTGAAAACGAAC  >  NC_000913/665706‑665992

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