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	96.4%	72.9 / NA	28	D104N (GAT→AAT)	mepM	peptidoglycan DD‑endopeptidase MepM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/13);  minor base G (1/0);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGCGAACGAATTTCG  >  NC_000913/1941191‑1941465                                                                                                                                        |                                                                                                                                            tttCGCCGTCCGCGGTTAATGTCCAGGAGTGTTGTTGACCGATTTTTAGGTTATGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTg                                                                                                                                      <  2:247684/142‑1 (MQ=255)  ttCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTgg                                                                                                                                     >  1:297488/1‑143 (MQ=255)     gCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGaaa                                                                                                                                  >  2:241199/1‑143 (MQ=255)                ttAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGCGAATGATGCTTAGCGTATTGCCAATGTAAACAACATATTCa                                                                                                                       >  2:103463/1‑143 (MQ=255)                   ctGTCCAGCAGAGTTGTTGACCGATTTACAGGTTACGCAATTCCTTGTCGGCCGCCGCCAGATCAAACATATCACCCATATCAATGCCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGc                                                                                                                    <  2:188620/142‑1 (MQ=255)                        cAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCTAGTGGAAACAACATATTCATGCACGcc                                                                                                               <  2:76941/143‑1 (MQ=255)                                 tgttgACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCAcc                                                                                                      <  1:355368/143‑1 (MQ=255)                                      aCCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATACCACTACCCTACTGATTGAGAATGCTGCTTAGCGTATTGCCAGGGGAAACCCCATATGTTTGCAGCCCCGCTTCGCCGGGGa                                                                                                 >  2:127165/1‑143 (MQ=255)                                        cGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGACACAACATATTCATGCACGCCCGCTTCACCggagctt                                                                                               >  1:578208/1‑138 (MQ=255)                                         gATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGAttt                                                                                              >  1:89661/1‑143 (MQ=255)                                             ttCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTc                                                                                          >  2:617197/1‑143 (MQ=255)                                                aGGTTATGGTATTCTTTGTCGTCCGCATCCAGTTGGTTCATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATc                                                                                       <  2:556657/143‑1 (MQ=255)                                                      cGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTc                                                                                 >  2:262112/1‑143 (MQ=255)                                                         aaTTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATGGCCAGTGGAAACAACATATTCATGCACGGCCGCTTCACCGGCGAGCTTGTCATCGACTTCgac                                                                              >  2:295126/1‑141 (MQ=255)                                                                           gCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCAtc                                                            <  1:105596/143‑1 (MQ=255)                                                                              aGTTGGGTGACATCACCCATATCGATACCATACTGATTGAGAATGCTGCTTAGCGCATTGCGAGTGGAGCCAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTAATCCAGTTCGTCCTGGCGAATGGCTTCAtcttc                                                         <  2:509405/143‑1 (MQ=255)                                                                                ttGGGTGGTTTCGCCTCTATCAATACCATACTGTTTTAGACTGCTGCTTTGCGTATTGCCAGTTGTATCAACATATTCTTGCTCGCTCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGTGATTGGCTTCATcttctt                                                       <  2:404704/143‑1 (MQ=255)                                                                                    gTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGc                                                   >  2:230576/1‑143 (MQ=255)                                                                                          tCAGCCATATCAATACCATACTGAGTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTg                                             <  2:194897/143‑1 (MQ=255)                                                                                                atatCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAAt                                       >  1:294336/1‑143 (MQ=255)                                                                                                atatCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTCGTCAAt                                       >  1:109389/1‑143 (MQ=255)                                                                                                 tagcAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATc                                      <  1:241199/140‑1 (MQ=255)                                                                                                  atCAATAGCATACTGATTGAGAATGCTGCTTAGCGTAGTGCCAGTGGAAACAACATATTCAGGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCg                                     <  1:117629/143‑1 (MQ=255)                                                                                                   tCAATACCATACTGATTGAGAATGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCgg                                    <  1:105001/143‑1 (MQ=255)                                                                                                    cAACACCATACTGATTGATACTGCTGCTTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCCGCGATTTTGTCATTCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGc                                   <  1:94826/143‑1 (MQ=255)                                                                                                                            tgctTAGCGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGcgaa           >  1:223158/1‑143 (MQ=255)                                                                                                                                   cGTATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCCGGGGAATGGCTTCATCTTCAAGCGCAGCTCGATCAATCGGCGCACTGGCTCCCGGCACCAGcgaacgaaaac    >  2:409570/1‑140 (MQ=255)                                                                                                                                     tATTGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGCGAACGAATTTCg  >  1:612303/1‑143 (MQ=255)                                                                                                                                        |                                                                                                                                            GTTCGCCGTCCGCGGTTAATGTCCAGGAGAGTTGTTGACCGATTTTCAGGTTACGCAATTCTTTGTCGGCCGCAGCCAGTTGGGTGATATCACCCATATCAATACCATACTGATTGAGAATGCTGCTTAGCGTATCGCCAGTGGAAACAACATATTCATGCACGCCCGCTTCACCGGCGATTTTGTCATCCAGTTCGTCCTGGGGAATGGCTTCATCTTCTTGTGCAGCTTGATCAATCGGCTCACTGGCTTCAGGTAAGAGCGAACGAATTTCG  >  NC_000913/1941191‑1941465

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: ‑