Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,941,225	T→C	Q137Q (CAA→CAG)	mepM ←	peptidoglycan DD‑endopeptidase MepM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,941,225	0	T	C	100.0%	189.7 / NA	53	Q137Q (CAA→CAG)	mepM	peptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/31);  total (22/31)

CACTCTCCTTGCTGCATTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCC  >  NC_000913/1941088‑1941356                                                                                                                                          |                                                                                                                                    cACTCTCCTTGCTGCATTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGttgctg                                                                                                                                   >  2:459622/1‑140 (MQ=255)  aCTCTCCTTGCTGCATTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGa                                                                                                                                  >  2:165507/1‑140 (MQ=255)   ctctCCTTGCTGCATTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGAc                                                                                                                                 >  1:789667/1‑140 (MQ=255)   ctctCCTTGCTGCATTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGAc                                                                                                                                 >  1:208368/1‑140 (MQ=255)         ttGCTGCATTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTTACCGAttt                                                                                                                           <  2:846465/140‑1 (MQ=255)          tgctgcATTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGAtttt                                                                                                                          <  2:789667/140‑1 (MQ=255)                aTTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGtt                                                                                                                    <  2:208368/140‑1 (MQ=255)                         gTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTc                                                                                                           <  1:115868/140‑1 (MQ=255)                                       gTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAg                                                                                             <  1:268387/140‑1 (MQ=255)                                        ttAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGGCACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGc                                                                                            >  2:827053/1‑140 (MQ=255)                                        ttAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGGCACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGc                                                                                            <  1:827053/140‑1 (MQ=255)                                        ttAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGc                                                                                            <  2:823679/140‑1 (MQ=255)                                        ttAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGc                                                                                            >  1:95671/1‑140 (MQ=255)                                        ttAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGc                                                                                            >  2:318677/1‑140 (MQ=255)                                          aGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCa                                                                                          >  2:141998/1‑140 (MQ=255)                                                 gTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTggg                                                                                    <  1:301514/139‑1 (MQ=255)                                                 gTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTggg                                                                                    >  2:301514/1‑139 (MQ=255)                                                     ggTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGata                                                                               <  1:571556/140‑1 (MQ=255)                                                     ggTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGata                                                                               <  1:252407/140‑1 (MQ=255)                                                      gTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGatat                                                                              <  1:165507/140‑1 (MQ=255)                                                       tCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATc                                                                             >  2:167289/1‑140 (MQ=255)                                                        cATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCa                                                                            <  2:854301/140‑1 (MQ=255)                                                          tcggTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCAcc                                                                          <  2:801004/138‑1 (MQ=255)                                                          tAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCAcc                                                                          <  1:719957/140‑1 (MQ=255)                                                          tAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCAcc                                                                          <  1:39358/140‑1 (MQ=255)                                                            ggTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCC‑CGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCat                                                                       >  1:330517/1‑140 (MQ=255)                                                               tCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCatat                                                                     >  1:314623/1‑140 (MQ=255)                                                                cGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGatat                                                                              <  1:578298/130‑1 (MQ=255)                                                                cGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGatat                                                                              >  2:578298/1‑130 (MQ=255)                                                                cGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCatat                                                                     >  1:701173/1‑139 (MQ=255)                                                                cGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCatat                                                                     <  2:701173/139‑1 (MQ=255)                                                                 gggTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTggg                                                                                    >  1:645925/1‑123 (MQ=255)                                                                 gggTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTggg                                                                                    <  2:645925/123‑1 (MQ=255)                                                                 gggTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCatat                                                                     <  1:577798/138‑1 (MQ=255)                                                                 gggTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCatat                                                                     >  2:577798/1‑138 (MQ=255)                                                                 gggTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCa                                                                   >  2:75660/1‑140 (MQ=255)                                                                     ttCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATAc                                                               <  1:141998/140‑1 (MQ=255)                                                                            cgaGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTg                                                        <  2:155929/140‑1 (MQ=255)                                                                                            tGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCtt                                        <  1:167289/140‑1 (MQ=255)                                                                                              aGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGaa                                                >  1:785061/1‑130 (MQ=255)                                                                                              aGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGaa                                                <  2:785061/130‑1 (MQ=255)                                                                                              aGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAg                                      >  2:758651/1‑140 (MQ=255)                                                                                              aGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAg                                      >  1:691579/1‑140 (MQ=255)                                                                                              aGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAg                                      >  2:194619/1‑140 (MQ=255)                                                                                                gcgcTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCtt                                        >  2:242317/1‑136 (MQ=255)                                                                                                gcgcTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCtt                                        <  1:242317/136‑1 (MQ=255)                                                                                                            gCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTgg                        <  1:75660/140‑1 (MQ=255)                                                                                                            gCCGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTgg                        <  1:758651/140‑1 (MQ=255)                                                                                                             ccGTCCGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGa                       <  2:691579/140‑1 (MQ=255)                                                                                                               cgtcccgGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGaaa                     <  2:330517/135‑1 (MQ=255)                                                                                                                 ccGCGGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAaaca                   >  1:831486/1‑140 (MQ=255)                                                                                                                  cgcgGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAaacaa                  <  2:314623/140‑1 (MQ=255)                                                                                                                   gcgGTTAATGTCCAGGAGAGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAaacaac                 <  2:831486/140‑1 (MQ=255)                                                                                                                                  gagaGTTGCTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGcc  <  1:459622/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                    CACTCTCCTTGCTGCATTTCGCTGGTCATTTTAAAACCGTTAGCGGCAGTACGGTCATAGGTTCGGGTTTCACGACGAGACACTTCCCAGGTGAGGCGCTGCAGTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCC  >  NC_000913/1941088‑1941356

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