Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,941,326	C→T	D104N (GAT→AAT)	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,326	0	C	T	95.7%	136.3 / NA	46	D104N (GAT→AAT)	mepM	peptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (18/26);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGCGAACGAATT  >  NC_000913/1941192‑1941462                                                                                                                                       |                                                                                                                                         ttCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTgg                                                                                                                                  >  1:73123/1‑143 (MQ=255)   cGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGaa                                                                                                                                <  2:796105/143‑1 (MQ=255)    gCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTGTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATACAACCCATATCAATACCCTACTGAGTGAGAACACTGCCTAGCCTATCGCCCCTGGacc                                                                                                                               >  2:1027779/1‑141 (MQ=255)       gTCCGCGGTTAATGTCGAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAaacaa                                                                                                                            <  1:252254/143‑1 (MQ=255)       gTCCGCGGTTAATGTCCATGAGAGTTTTTTACCTATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAaacaa                                                                                                                            <  1:746959/143‑1 (MQ=255)       gTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAaacaa                                                                                                                            <  1:640702/143‑1 (MQ=255)        tCCGGGGTTGATGTCCAGGAGGGGCGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCGTAGCGTATTGCCAGTGGAaacaac                                                                                                                           <  2:995398/143‑1 (MQ=255)              tttAATGTCCGGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATCTCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTc                                                                                                                     <  2:973728/142‑1 (MQ=255)              gTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTc                                                                                                                     <  2:336825/143‑1 (MQ=255)                tAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCAt                                                                                                                   <  2:980086/143‑1 (MQ=255)                     tCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTACCGTATTGCCAGTGGAAACAACATATTCGTGCACg                                                                                                              >  1:64323/1‑143 (MQ=255)                             agTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTc                                                                                                      >  1:13248/1‑143 (MQ=255)                                 gttgACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCg                                                                                                  >  1:585806/1‑143 (MQ=255)                                  ttgACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCCATACCATACTGATTGAGAATGCTGCCTCGCGTATCGCCAGTGGAAACACCATACTCATGCACGCCCGCCTCCCCgc                                                                                                 >  2:252708/1‑142 (MQ=255)                                       cgtttttCATGGTACCCGATTCTTTGTCGGCCGCAGCCAGTTGTGTTATATCACCCATATCAATTCCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGAtt                                                                                            <  2:590405/140‑1 (MQ=255)                                       cGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATGCCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATAGTCAGGCACGGCCGGGTCGCCGGCGagg                                                                                            >  1:902194/1‑141 (MQ=255)                                         aTTTTCAGGTTCCGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTCACATCACCCATATCAATCCCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGAtttt                                                                                          <  1:521341/143‑1 (MQ=255)                                           tttCAGGTTACGCAATTCTTTGCCGGCCGCAGCCAGTTGGGGGATATCACCCATATCAATACCATACTGATTCAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCCTGCACGCCCGCTTCACCGGCGATTTTGt                                                                                        <  1:73248/143‑1 (MQ=255)                                                ggTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATcc                                                                                   >  1:469243/1‑143 (MQ=255)                                                ggTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATcc                                                                                   >  1:825710/1‑143 (MQ=255)                                                   tACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCCTACTTAGTGAGACTGTGGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGCATTTGTCATCCAGt                                                                                >  2:124288/1‑143 (MQ=255)                                                    aCGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGtt                                                                               <  1:318387/143‑1 (MQ=255)                                                      gCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGCTCg                                                                             >  1:451923/1‑143 (MQ=255)                                                            cTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTgg                                                                       <  1:469239/143‑1 (MQ=255)                                                                gTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATCTGGTCATCGAGTTCGTCCTggggca                                                                   >  2:51549/1‑141 (MQ=255)                                                                  cGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATg                                                                 >  1:89197/1‑143 (MQ=255)                                                                    gCCGCAGCCGGTTGGGTGATCCCACCCATATCAATCCCCTACTGATTGAGAATGCTGCGTAGCGTATTGCCAGTGGAAACAACCTATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGc                                                               <  1:588674/143‑1 (MQ=255)                                                                              gTTGGGTGTTCTCTCCCTTTGCAGTCCCATACTGATTGGGATTGCTGCTTGGCGTCTTGCCAGTGGAAATAATATATTCATGCTCGCCCGCTTCACCGGCGTTTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCAtcttct                                                     <  2:515647/143‑1 (MQ=255)                                                                               ttGGGTGATATCACCCATATCAATACCATACTGATTTACAATGCTGTTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATcttctt                                                    <  1:702463/143‑1 (MQ=255)                                                                                   gTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGCTTTTGTCATCCAGTTCGTCCTAGGGAATGGCTTCATCTTCTTGTGc                                                >  1:635857/1‑143 (MQ=255)                                                                                   gTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGAACGCCCGCTTCACCGGCGAGTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGc                                                >  1:808362/1‑143 (MQ=255)                                                                                     gATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTCGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTCTGTCATCCAGTTCGTCCCGGGGAATGTGTTCATCTTCTTCTGccg                                              >  1:902641/1‑141 (MQ=255)                                                                                         tCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGTAATGGCTTCATCTTCTTGTGCAGCTTg                                          <  2:570551/143‑1 (MQ=255)                                                                                           aCCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGAt                                        <  1:94206/143‑1 (MQ=255)                                                                                             ccATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCa                                      >  1:964250/1‑143 (MQ=255)                                                                                                  tCAATACCTGACTGATCGAGAACGCTGCATCGCGTATTGCCCGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCAGCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCgg                                 <  2:409396/143‑1 (MQ=255)                                                                                                  tCAAGACCATACTGATCGAGAATGCTGCTCAGCGTATTGCCAGTGGAAACAACGTATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCCGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCgg                                 <  2:839944/143‑1 (MQ=255)                                                                                                    aaTACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCt                               <  1:73865/143‑1 (MQ=255)                                                                                                      tACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTTATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCa                             >  2:794881/1‑143 (MQ=255)                                                                                                                aTTGCGAATGCTGCGTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAg                   <  2:15319/143‑1 (MQ=255)                                                                                                                    agaATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGGTCAATCGGCTCACTGGCTTCAGCTac               >  1:126278/1‑142 (MQ=255)                                                                                                                      aaTGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCGCCGGCGATTTTGTCATCCAGTTCGTCCTGTGGAATGACTTCATCTTCTTGTGCCGCTTGATCAATCGTCTCACTGGCTTCAGGTAcga             >  2:284422/1‑143 (MQ=255)                                                                                                                           tgctTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGcgaa        <  2:825710/143‑1 (MQ=255)                                                                                                                           tgctTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGcgaa        <  2:469243/143‑1 (MQ=255)                                                                                                                              ttAGCGTATTGCCCGGGGAAACAACAGATTCCCGGCCGCCCGCGTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCAGCTTCTGGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGcgaacga     <  1:253479/143‑1 (MQ=255)                                                                                                                                 gCGTATTGCCAGCGGAAACAACATCTTCATGCACGCCCGCTTCACCGGGGTTTTGGCCAGCCACTTCGTCCGGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGCGAACGAAtt  <  1:1029039/143‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         TTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGCGAACGAATT  >  NC_000913/1941192‑1941462

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