Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,526,613	G→T	13.6%	C253F (TGC→TTC)	ddlB →	D‑alanine‑‑D‑alanine ligase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,526,613	0	G	T	13.6%	120.1 / 7.6	44	C253F (TGC→TTC)	ddlB	D‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base T (3/3);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

AGGTGCTGCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGA  >  NC_002947/1526466‑1526757                                                                                                                                                     |                                                                                                                                                 aGGTGCTGCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCg                                                                                                                                                >  2:465857/1‑150 (MQ=255)                   cccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTTAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTTCGATGCCAATGGCATCGCAg                                                                                                                              <  2:359869/150‑1 (MQ=255)                     cTGGTTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCCCCCAGCCTTCGATGCCAATGGCATCGCAgg                                                                                                                             >  1:146879/1‑150 (MQ=38)                             gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg                                                                                                                    >  2:403610/1‑150 (MQ=255)                              ccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATc                                                                                                                                  <  1:349513/135‑1 (MQ=255)                              ccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATc                                                                                                                                  >  2:349513/1‑135 (MQ=255)                                   cGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGGGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGa                                                                                                              <  2:12927/150‑1 (MQ=255)                                    gtgtTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGAc                                                                                                             <  2:701596/150‑1 (MQ=255)                                      gtTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGt                                                                                                           <  1:465857/150‑1 (MQ=255)                                             gACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAg                                                                                                    >  1:116881/1‑150 (MQ=255)                                              aCTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                 <  1:184252/136‑1 (MQ=255)                                              aCTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                 >  2:184252/1‑136 (MQ=255)                                                      gCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAgg                                                                                                   >  1:226426/1‑142 (MQ=255)                                                      gCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAgg                                                                                                   <  2:226426/142‑1 (MQ=255)                                                              cATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCT‑GGGGCGGTTGGAc                                                                                                             >  2:567429/1‑123 (MQ=255)                                                              cATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCT‑GGGGCGGTTGGAc                                                                                                             <  1:567429/123‑1 (MQ=255)                                                                tCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGAGCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCt                                                                                 >  2:730165/1‑150 (MQ=255)                                                                     aaCGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctg                                                                            >  2:767937/1‑150 (MQ=255)                                                                           aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACCCGTAGGGCCGGTTCTGGCTGCTCGaa                                                                      <  1:706594/150‑1 (MQ=255)                                                                             ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATg                                                                                                                                             >  1:444122/1‑77 (MQ=255)                                                                             ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATg                                                                                                                                             <  2:444122/77‑1 (MQ=255)                                                                              cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCCCCCAGCCTTCGATGCCaatg                                                                                                                                       >  2:634585/1‑80 (MQ=25)                                                                              cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCCCCCAGCCTTCGATGCCAATg                                                                                                                                       <  1:634585/83‑1 (MQ=25)                                                                               aGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCa                                                                  >  2:509184/1‑150 (MQ=255)                                                                                     gCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCg                                                            >  2:575527/1‑150 (MQ=255)                                                                                          ccGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCa                                                       <  1:767937/150‑1 (MQ=255)                                                                                           cGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAg                                                      >  1:69686/1‑150 (MQ=255)                                                                                             tGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                      >  2:674931/1‑132 (MQ=255)                                                                                             tGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                      <  1:674931/132‑1 (MQ=255)                                                                                                 ggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTTCGATGCCAATGGCATCgcag                                                                                                                              >  1:161059/1‑69 (MQ=37)                                                                                                 ggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTTCGATGCCAATGGCATCGCAg                                                                                                                              <  2:161059/72‑1 (MQ=37)                                                                                                 ggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGa                                                <  1:460623/150‑1 (MQ=255)                                                                                                    tGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAg                                             <  1:200449/150‑1 (MQ=255)                                                                                                      gATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                 <  2:224314/80‑1 (MQ=255)                                                                                                      gATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                 >  1:224314/1‑80 (MQ=255)                                                                                                       aTAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATc                                          >  2:464134/1‑150 (MQ=255)                                                                                                          gCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATa                                       <  1:623079/150‑1 (MQ=255)                                                                                                               aaGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTg                                  <  2:611235/150‑1 (MQ=255)                                                                                                                                   ccTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgct                                                                             <  2:720755/87‑1 (MQ=255)                                                                                                                                   ccTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgct                                                                             >  1:720755/1‑87 (MQ=255)                                                                                                                                        ccGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTg         <  1:509184/150‑1 (MQ=255)                                                                                                                                         cgcgagc‑gCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGt                                <  2:156690/126‑1 (MQ=255)                                                                                                                                         cgcgagc‑gCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGt                                >  1:156690/1‑126 (MQ=255)                                                                                                                                               c‑gcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTggg  <  2:353894/150‑2 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AGGTGCTGCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGA  >  NC_002947/1526466‑1526757

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