Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,721,897	G→A	D649N (GAT→AAT)	patZ →	peptidyl‑lysine N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,721,897	0	G	A	90.0%	20.8 / ‑6.1	10	D649N (GAT→AAT)	patZ	peptidyl‑lysine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (5/4);  minor base T (0/1);  total (5/5)
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.78e-01

GCAAAGTATGGCTAACCGTGCTGGCGCTCAGGAGTTGCGGGTTGTGGTTGAGCACGATCCGGTTTTCGGGCCGTTGATCATGCTGGGTGAAGGCGGTGTGGAGTGGCGTCCTGAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGT  >  NC_000913/2721671‑2722121                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 ncAAAGTATGGCTAACCGTGCTGGCGCTCAGGAGTTGCGGGTTGTGGTTGAGCACGATCCGGTTTTCGGGCCGTTGATCATACTGGGTGAAGGCGGTGTGGAGTGGCGTCCTGAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGAATGTTGCAGGCTTGAGCCAGCTTc                                                                                                                                                                                                          >  1:203839/2‑251 (MQ=255)                                                       cGATCCGGTTGTCGGGCCGTTGATCATGCTGGGTGAAGGCGGGGTGGAGTGGCGTCCTTAAGATCAAGCCGTCGTCGCACTGCCGCCGGTGAACATCAACCTGGTCCGCTATCTGGTGATTCAGGGGATCAAAAGTAAAAAGATTCGGGCTCGCAGTGCGCTACGGCCAGTGTATGTTGCAGGCTTGAGCCAGTGTCTGGTGCAGGGTGCCAACTTGATTGTCGATTGCCCg                                                                                                                                                                       <  1:217094/232‑1 (MQ=255)                                                                                     gggTGAAGGCGGTGTGGAGTGGCGTCCTGAAGATCAAGCCGCCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGAATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAAt                                                                                                                      >  1:219854/1‑251 (MQ=255)                                                                                              cGGTGTGGAGTGGCGTCCTGAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGACCATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTGCGCCCATTGAATGTTGCAGGCTTGAGCCAGCTTCTTGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGc                                                                                                                                    <  1:130471/228‑1 (MQ=255)                                                                                                                         aGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGAATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTTATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTg                                                                                                                               <  1:42692/206‑1 (MQ=255)                                                                                                                                 gcgttacgttctcctgaagtaTGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGAATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGc                                                                          <  1:20554/231‑1 (MQ=255)                                                                                                                                                          ncTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCACATTGAATGTTGCAGGCtt                                                                                                                                                                                                                     >  1:28589/2‑87 (MQ=255)                                                                                                                                                                                     naTCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGAATGTTGCAGGATTGAGCCCGCTTCCGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATCTTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAAt                      >  1:28335/2‑251 (MQ=255)                                                                                                                                                                                          aaaGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGAATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGTCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTa                 <  1:3215/251‑1 (MQ=255)                                                                                                                                                                                                         gTGCGCGCAGTGCGCTACGCCCATTGAATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGt  >  1:222480/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 GCAAAGTATGGCTAACCGTGCTGGCGCTCAGGAGTTGCGGGTTGTGGTTGAGCACGATCCGGTTTTCGGGCCGTTGATCATGCTGGGTGAAGGCGGTGTGGAGTGGCGTCCTGAAGATCAAGCCGTCGTCGCACTGCCGCCGCTGAACATGAACCTGGCCCGCTATCTGGTTATTCAGGGGATCAAAAGTAAAAAGATTCGTGCGCGCAGTGCGCTACGCCCATTGGATGTTGCAGGCTTGAGCCAGCTTCTGGTGCAGGTTTCCAACTTGATTGTCGATTGCCCGGAAATTCAGCGTCTGGATATTCATCCTTTGCTGGCTTCTGGCAGTGAATTTACCGCGCTGGATGTCACGCTGGATATCTCGCCGTTTGAAGGCGATAACGAGAGTCGGCTGGCAGTGCGCCCTTATCCGCATCAGCTGGAAGAATGGGTAGAATTGAAAAACGGT  >  NC_000913/2721671‑2722121

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