Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,940,584	C→A	R351L (CGT→CTT)	mepM ←	peptidoglycan DD‑endopeptidase MepM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,940,584	0	C	A	96.5%	65.4 / NA	30	R351L (CGT→CTT)	mepM	peptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (19/9);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGAT  >  NC_000913/1940454‑1940722                                                                                                                                   |                                                                                                                                           cAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCCTATAACGCGTGGTGTAGATGCGACCATGAAGAAGAGCCACat                                                                                                                                >  2:50621/1‑143 (MQ=255)                    cGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATCGCCACATACTACCCTGCGGCGCCGCt                                                                                                             >  2:158620/1‑143 (MQ=255)                       ccGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATCTAACGCGTGGTGTAGCTGCGACCGTGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCg                                                                                                          >  2:592912/1‑143 (MQ=255)                        cGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTACCGGTATCACCAGACTCTTGCGCAAGTGCATATAACGCGGGGCGCAGCTGCGACCATGCAGAATAGCCACATACAAACCAGCTGAGCCCGTGCGt                                                                                                         >  2:168912/1‑143 (MQ=255)                               aCCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCaa                                                                                                  <  2:396948/143‑1 (MQ=255)                               aCCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACAGAATAACCTGCTGCGCCACTGCGTTTGGcga                                                                                                  >  1:466683/1‑141 (MQ=255)                                  ggAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGCTTTTACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAGTCGCCACATAATAACCTCCTGCTCCACTGCGGTCAGCAacca                                                                                               >  2:571737/1‑143 (MQ=255)                                   gAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAGTCTTGCGCAAGTGCCTATAACGCGTGGTGTAGCTGCGACCCTGAAGAATAGCCACATAATCACCTGCTGCGCCACGGCGGTCGGCAaccac                                                                                              >  2:525164/1‑143 (MQ=255)                                      agcggtaTACGGCCGTGACGTTTCGCCGTCTGTCCCGGTTTGCCCAGAATGTTGCGTAAGTGCATATAACGCGTTGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACcaccac                                                                                           <  2:89455/137‑1 (MQ=255)                                                  tCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCAcc                                                                               >  1:245943/1‑143 (MQ=255)                                                     ccACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGCGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTTGGCCACTGCGATTCGCAACCACCACTTACCCGTCCGCCAc                                                                            >  2:337704/1‑143 (MQ=255)                                                       aCGTCTCACCTTCTGGCTCGGGTTCACGAGAAGCCTGCGCAAGTGCATATCACGCGTGGTGTACCTGCGCCCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTg                                                                          <  2:290027/143‑1 (MQ=255)                                                          tttCACCTTCTGTCCCGGTTTCACCAGCATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGaaa                                                                       <  2:295376/143‑1 (MQ=255)                                                              aCCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATACCGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGCAAGCAc                                                                   >  2:276977/1‑143 (MQ=255)                                                               ccTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTCGCAACCACCACTTCACCGTCACCCACTGAAAGcccc                                                                  >  2:45456/1‑140 (MQ=255)                                                                      tCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTCGCAGCCACCACTTCACCGTCACCCACTGCAAGCACTGGCGTAc                                                           >  2:596838/1‑143 (MQ=255)                                                                                  cAGATTCTTGCGCAAGTTCATATCACGCGTTGTGTAGCTGCGCCCATGAAGAATAGCCACATAATATCCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATgg                                               <  2:198048/143‑1 (MQ=255)                                                                                     cttCTTGCGGTAGCGCGCATAGCGCGGGGTGTAGCTGCGACCATGAAGATTAGCCACATAATAACCTGCTGCGCCCCTGCGTTTGGCAACCACCACGTCTCCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGa                                            <  2:191054/141‑1 (MQ=255)                                                                                      aTCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCTACCGTGTAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGaa                                           <  1:168863/143‑1 (MQ=255)                                                                                      aTCTGGCGCAAGGGCATATAACGCGTGGTGTAGCTAAGTGGTTGATGGATAGCCACATAATACCCTGCTGCGCCACTGCGTTTGGCAACCCCCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGaa                                           <  2:129773/143‑1 (MQ=255)                                                                                            cgcAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAac                                     >  1:455301/1‑143 (MQ=255)                                                                                                     atatAACGCGTGGTGTAGCTGCGACCATGAAGCATAGCCACATAATAACCTGCTGCGCAACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTCCGGCATGGCGAAATCCACACCTCtgcg                            >  2:142373/1‑141 (MQ=255)                                                                                                      cacaaCGCGTTGCGCAGCTGCGACCATGAAGAACACCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCtgtgt                           <  1:50621/140‑1 (MQ=255)                                                                                                        taACGCGTGGTGTAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTCGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGGCGCCTGGCCAAATCAACCCCTCCGTgtgg                         >  2:471144/1‑143 (MQ=255)                                                                                                           cgcgTGGTGTAGCTGCGACCATGAAGAATAGCTACATAATAACCGGCTGCGCAACTGCGGTTGGCGACCACCACTTCACCGTCACGCACTGAAAGCACTGGCGTACCCGGCGGCATGGCGAGCTCAACACCGCAGTGCGGGGc                      >  1:242919/1‑143 (MQ=255)                                                                                                                   gtAGCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGAc              >  1:561309/1‑143 (MQ=255)                                                                                                                      gCTGCGACCATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACAGTCACCCACTCAAAGCACTGGCGTACCATGCGGCATGGCGAAATCAACACCTCTGGGGGGCGCAACGCAACCgg           >  2:241560/1‑143 (MQ=255)                                                                                                                       cTGCGACCATGCAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTACACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGCCGAAATCAACACCTCTGCGTCGCGCAGcgcgccagac          >  2:595154/1‑137 (MQ=255)                                                                                                                              cATGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCCCTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCCTGCCGCGATCAACACCTCTCTCTGCTGCACCGCGACCCGTCACCCGa   >  1:557541/1‑143 (MQ=255)                                                                                                                               aTGAAGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGAt  <  1:307703/143‑1 (MQ=255)                                                                                                                                   |                                                                                                                                           CAGATGCGGCCCGGTTGAACGTCCGGTATTACCGGAAAGCGCGATACGGTCGCCACGTTTCACCTTCTGTCCCGGTTTCACCAGAATCTTGCGCAAGTGCATATAACGCGTGGTGTAGCTGCGACCATGACGAATAGCCACATAATAACCTGCTGCGCCACTGCGTTTGGCAACCACCACTTCACCGTCACCCACTGAAAGCACTGGCGTACCCTGCGGCATGGCGAAATCAACACCTCTGTGTGGTGCAACGCGACCGGTCACCGGAT  >  NC_000913/1940454‑1940722

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