Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,455,812	G→A	G10S (GGC→AGC)	mpl →	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,455,812	0	G	A	97.1%	112.1 / NA	34	G10S (GGC→AGC)	mpl	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (23/10);  total (23/11)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AACTCGTGCTGCTTTTCGACAATAAATTCACCTAACGTTTTCATAAAACTTTCCCTGCAATGTTTATGGAGTAAAGCGACCGCAACAATCTTAACAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTC  >  NC_000913/4455566‑4456058                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aaCTCGTGCTGCTTGTCGACAATAAATTCACCTAACGTTTTCATAAAACTTTCCCTGCAATGTTTATGGAGTAAAGCGACCGCAACAATCATAACAAATAATCTCAATGTTGCGCTCAGGGGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGcac                                                                                                                                                                                                                                                    <  2:286749/251‑1 (MQ=255)                          aTTCACCTAACGTTTTCATAAAACTTTCCCTGCAATGTTTATGGAGTAAAGCGACCGCAACAATCTTAACAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTg                                                                                                                                                                                                                           >  1:295149/1‑251 (MQ=255)                                           aTAAAACTGTCTGTGCAATGCTTATGAAGTAAAGCGCCCGCAACAATCTTCACAAATAATCTCAATGTTGCGCTCAAGCAAATCACGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGcc                                                                                                                                                                                                           >  2:220665/1‑250 (MQ=255)                                               aaCTTTCCCTGCAATGTTTATGGAGTAAAGCGACCGCAACAATCTTAACAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGa                                                                                                                                                                                                       >  1:27545/1‑250 (MQ=255)                                                    tCCCTGCAATGTTTATGGAGTAAAGCGACCGCAACAATCTTAACAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAAc                                                                                                                                                                                                 >  2:135840/1‑251 (MQ=255)                                                                   tGGAGTAAAGCGACCGCAACAATCTTAACAAATAATCTCAATGTTACGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCa                                                                                                                                                                                   >  1:329427/1‑250 (MQ=255)                                                                          aaGCGACCGCAACAATCTTAACAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTAt                                                                                                                                                                           >  2:257721/1‑251 (MQ=255)                                                                              gACCGCAACAATCTTAACAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATccgc                                                                                                                                                                       <  2:68260/251‑1 (MQ=255)                                                                                            taaaaaaTAATCTCAATGGTGCGCTCAGGTGAAACGCGCAAGCAAAATACGGAGTATCCTGAAATGAGGTTCGAACGTGCCCGAAAAAGGCGTAAAATGAGTGGCAGATTAAAAAAGGAAAGTGCCGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTAGGGGCGGGCTGGCGAGGCTGGCGCGCCAGTGAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCAc                                                                                                                                                             <  2:304474/243‑1 (MQ=255)                                                                                            tAACAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGGGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATTCTGCCCCGCCAGTTAGGCCATGAAGTAGCGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCt                                                                                                                                                           <  1:31363/249‑1 (MQ=255)                                                                                            tAACAAATAATCTAAATGTTGCGATCGGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGGTTTTCTCACTTGCCCGACATATGCGTAAAATGGGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTGCGGACGCCAATGTGTATCCGCCGATGAGCAc                                                                                                                                                             <  2:277606/247‑1 (MQ=255)                                                                                               cAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGACATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGCTGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACtt                                                                                                                                                      >  2:40816/1‑251 (MQ=255)                                                                                                                      aGGTGAATAGCGACAGCAAATTACGGATTATCCTGCAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGAt                                                                                                                               >  2:320917/1‑251 (MQ=255)                                                                                                                                     gCAAATTACGGATTATCCTGAAATGCGTTTCTAACTTGCCCGACATATCCGTAAAATGAGCGGCCGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTACCGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCCCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATgc                                                                                                                >  1:232910/1‑251 (MQ=255)                                                                                                                                        aaTTACGGATTATCCTGAAATGCGTTTCTCACTCGCCCGACATATGCGTAAAATGAGCGGCAGATTACAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATgcc                                                                                                               >  1:154305/1‑249 (MQ=255)                                                                                                                                                     cctaaaaTGCTTTTCCAACTTGCCGTACTTATGCGTAAAATGAGCGTAAGATTAAAAAAGGATAGGAACGTGTTGC‑AATTCATATTTGAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTGGGCAATGAAGTAAGGGGTTCGGACGCAAATGGGTTTCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCg                                                                                                <  2:118657/247‑1 (MQ=255)                                                                                                                                                      cTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGc                                                                                               >  1:13716/1‑251 (MQ=255)                                                                                                                                                          aaTGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGcc                                                                                           >  1:16297/1‑251 (MQ=255)                                                                                                                                                               gTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAg                                                                                                                           <  1:195696/214‑1 (MQ=255)                                                                                                                                                               gTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAg                                                                                                                           >  2:195696/1‑214 (MQ=255)                                                                                                                                                                    tcACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTg                                                                                 >  1:87646/1‑251 (MQ=255)                                                                                                                                                                    tcACTTGCCCGACATATGCGAAAACTGAGCGGCAGATTAAAAACGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTg                                                                                 >  2:355802/1‑251 (MQ=255)                                                                                                                                                                           cccGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTg                                                                          >  1:268384/1‑251 (MQ=255)                                                                                                                                                                                     gCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCAt                                                                >  1:123669/1‑251 (MQ=255)                                                                                                                                                                                         aaaaTGAGCGGCAGATAAAAAAAGGATAGTGACGT‑ATGCGCATTAATTTTTGAGGAATTGGTAGCACTTTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTGGGCCATAAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGAc                                                             <  2:207463/250‑1 (MQ=255)                                                                                                                                                                                                     aaataaaaaaaGGTTAGTGACGT‑ATGCTCATTCATTTTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTACCGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCg                                                <  2:16297/246‑1 (MQ=255)                                                                                                                                                                                                                ggATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGGGTGGAAGc                                     >  2:144205/1‑251 (MQ=255)                                                                                                                                                                                                                  aTAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCgg                                   >  1:322600/1‑251 (MQ=255)                                                                                                                                                                                                                    aGTGACGT‑AGGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTa                                 <  1:99491/251‑1 (MQ=255)                                                                                                                                                                                                                       gACGT‑ATGCGCATTCATATTTTAGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTg                              >  1:40220/1‑251 (MQ=255)                                                                                                                                                                                                                                       atatTTTAGGAATTTGTAGCACATTTATGGGCGGCCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAACAAATATCCCt               >  1:42810/1‑251 (MQ=255)                                                                                                                                                                                                                                              aGGAATTTGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATg         <  1:116114/250‑1 (MQ=255)                                                                                                                                                                                                                                                    ttGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTc  >  2:106881/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ttGTAGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTc  >  1:113781/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AACTCGTGCTGCTTTTCGACAATAAATTCACCTAACGTTTTCATAAAACTTTCCCTGCAATGTTTATGGAGTAAAGCGACCGCAACAATCTTAACAAATAATCTCAATGTTGCGCTCAGGTGAATCGCGCCAGCAAATTACGGATTATCCTGAAATGCGTTTCTCACTTGCCCGACATATGCGTAAAATGAGCGGCAGATTAAAAAAGGATAGTGACGT‑ATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTC  >  NC_000913/4455566‑4456058

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