Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,329,743	G→T	D261Y (GAT→TAT)	dacB →	peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,329,743	0	G	T	94.4%	54.8 / NA	18	D261Y (GAT→TAT)	dacB	peptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (8/9);  total (8/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

GCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGC  >  NC_000913/3329503‑3329983                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 gCCCAAAAGCCTGGTGAGATGGTTTTGATACGCGTGGCATTTTATTCCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTCGTTATGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTAAGCTGACGGGATGCCTGCCACAACGTTTTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCAAGCTATGCCGGTGCAATTCTGAAATATGAGTTaaa                                                                                                                                                                                                                                        <  1:106788/251‑1 (MQ=255)  ccccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaaa                                                                                                                                                                                                                                        >  1:133296/1‑250 (MQ=255)                              aCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAacac                                                                                                                                                                                                           <  1:226459/251‑1 (MQ=255)                                  gTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgct                                                                                                                                                                                                       <  1:47111/251‑1 (MQ=255)                                              taCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATTACCCGGGGCGGAAAACTGCTGCGCCAGACTc                                                                                                                                                                                            >  1:150801/1‑250 (MQ=255)                                               aCCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGACTc                                                                                                                                                                                            >  1:10225/1‑249 (MQ=255)                                                                agtcatgtccgacccctccctCGTGGTTCTGCCGAAGCGCAATACTGCGAACTTGATGTGGTGCCAGTCGACCTGAACCGCTTTACGCTTACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCTGGTATCTCCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAc                                                                                                                                                                         <  1:167450/238‑1 (MQ=255)                                                                aGCCAGGTACGCCCCCCCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGTCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAc                                                                                                                                                                         <  1:227505/251‑1 (MQ=255)                                                                     ggTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGCGTTAAAACAGGCGGGGATCAACTGGAGGGGGACACTGCTGCACCCGACTccagtagacgcaacttgaaaggt                                                                                                                                                                      >  2:136313/1‑227 (MQ=255)                                                                                       tgtttCTGCCGAATCTAAATCCTGCGAACTGGAGGGGGTGCCAGGCGACCTGAACCGCTTTCCTCGGAGGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTGTTGCCTTGCAGGATGGAGCCAGCTATGCCGGTGCAAGTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTTCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTc                                                                                                                                                    <  2:150801/246‑1 (MQ=255)                                                                                                cGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGGGGAAAAATGCTGCGCCAGACCCCGGGTAACGAAACTGGGACGGTAGTTGCCAGTAAACAGTCGGCCCCGcggc                                                                                                                                         >  1:159210/1‑248 (MQ=255)                                                                                                                                cAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCg                                                                                                         <  1:41974/251‑1 (MQ=255)                                                                                                                                  ggCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGa                                                                                                       <  1:253006/251‑1 (MQ=255)                                                                                                                                                             ggagtcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCTGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCg                                                                            <  2:278174/246‑1 (MQ=255)                                                                                                                                                                    tgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGGACACTGCTGCGCCCGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCAGACCCGCTGCACGATCTGCTTAAGATTATGCTGTAAAAATAGGACAAAATGATCGCCGATACGGGTTTCCGCAttca                                                                      >  1:63337/1‑247 (MQ=255)                                                                                                                                                                                          cccGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCAACTGGAGCGGAACAATGCTGGGCCAAACTCAGGTTAAAGAAACTGGAACGGGAGTTGCCAGGAAAAAGTCGGCCCCGCGGCCCGATCTGCTTAAGATTATGATGAAAGAGGCGGACAAAATGATCGACGGTACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATTTGc                                               >  1:67629/1‑251 (MQ=255)                                                                                                                                                                                                          gttcAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACGGCTGCGCCAGAATCAGGGTAACGAAACTGGAAAGGGAGTTGACCGTAAACAGTCGGACCCGCCGCCCGATCTGCTTAAGATTAATCTGAAAAAGGCGGACAACATGCTCGCCGATACGGGTTTCCGCATTAAAGGCCATGCGCGCTTAAATGTGGCTGGAACAAGGCgggg                               >  1:245968/4‑250 (MQ=255)                                                                                                                                                                                                                                       aaTTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACCCTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTgc  <  1:50029/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 GCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGC  >  NC_000913/3329503‑3329983

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