Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,672,251	G→A	V98M (GTG→ATG)	ydgD →	putative serine protease YdgD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,672,251	0	G	A	100.0%	111.6 / NA	31	V98M (GTG→ATG)	ydgD	putative serine protease YdgD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/13);  total (18/13)

TAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAA  >  NC_000913/1672025‑1672486                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             tAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAAt                                                                                                                                                                                                                      >  1:360039/1‑250 (MQ=255)       aGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTc                                                                                                                                                                                                               >  1:387910/1‑251 (MQ=255)                   aTCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCCAAGGGTAAAGCCGATAAAGCCATGGCGCTTCGTTTTGTGTC‑AAATAAAGGACTTTGGCGCTAtc                                                                                                                                                                                                   >  2:315018/1‑250 (MQ=255)                       gCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGAt                                                                                                                                                                                               <  2:181304/251‑1 (MQ=255)                                gAGGTCAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGGTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTCTGAGATCCACGACAt                                                                                                                                                                                      >  2:321620/1‑251 (MQ=255)                                gAGGTCAGCACCCTGTTTTTTGGTCATGATGATCGCGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGGCTTTGGCGCTATGAGATCCACGACAt                                                                                                                                                                                      >  2:62501/1‑251 (MQ=255)                                  ggTCAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACAt                                                                                                                                                                                      >  1:224305/1‑249 (MQ=255)                                     cAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCTCTATGAGATCCACGACATAGAAg                                                                                                                                                                                 <  2:28118/251‑1 (MQ=255)                                       gCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGGGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGc                                                                                                                                                                               >  2:255990/1‑251 (MQ=255)                                           ccTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCcgcg                                                                                                                                                                           <  2:44399/251‑1 (MQ=255)                                                tttttGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCGAGCGGCAATTTATGTACGGTGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATGGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGATTGCCACGACATAGAAGGCCGCGTTGa                                                                                                                                                                       <  2:392144/250‑1 (MQ=255)                                                      gTCATGATGATCGTGTGCCAGTGAATGCCACGACGCAATCACCGTGGGATGCGGTTGGGCAGCTGGAAACGGCCAGCGGCATTTTATGTACGGCGACGCTGCTTGCACCCAATCTGGCATTAACGGCAGGAAACTGGTTCTTGACACCTCCCAAGGGGAAAGCCGATAAAGCAATGGCGCTGCCTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTAAGAGGTCCACCACATAGAAGGCCGCGTTACTCCGac                                                                                                                                                                 >  2:269923/1‑250 (MQ=255)                                                                                      aCCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATgg                                                                                                                                >  1:87391/1‑251 (MQ=255)                                                                                             cACCGTGGGATGAGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGAAGATGGGGATGGTTGGAtt                                                                                                                         >  2:18295/1‑251 (MQ=255)                                                                                                           gTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGAAACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGGTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGc                                                                                                           >  2:299496/1‑251 (MQ=255)                                                                                                                      ggAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTggg                                                                                                <  1:272725/251‑1 (MQ=255)                                                                                                                       gAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGa                                                                                               <  2:224305/251‑1 (MQ=255)                                                                                                                                aGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGgatt                                                                                      <  1:255990/251‑1 (MQ=255)                                                                                                                                                 ggCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATc                                                                     <  1:321620/251‑1 (MQ=255)                                                                                                                                                                aCCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGGAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCAtt                                                       >  2:340838/1‑250 (MQ=255)                                                                                                                                                                        tGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACgccgttg                                              >  2:325060/1‑251 (MQ=255)                                                                                                                                                                               aaCGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGGTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTAt                                       >  1:213671/1‑251 (MQ=255)                                                                                                                                                                                   gCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAg                                                                                                                                                                                 >  2:279754/1‑109 (MQ=255)                                                                                                                                                                                   gCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAg                                                                                                                                                                                 <  1:279754/109‑1 (MQ=255)                                                                                                                                                                                      ggACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAggg                                <  1:62501/251‑1 (MQ=255)                                                                                                                                                                                                tATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGc                       <  1:340838/250‑1 (MQ=255)                                                                                                                                                                                                   tGACACCTCCCATGGGTAAAGCCGATCAAGCAATGGCGCTGCGTTTTGCGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGAGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCcgcg                   >  2:252578/1‑251 (MQ=255)                                                                                                                                                                                                          tCCAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTACTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTAccg            >  2:179254/1‑251 (MQ=255)                                                                                                                                                                                                           ccAAAGGGTAAAGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTAcc             <  1:315018/249‑1 (MQ=255)                                                                                                                                                                                                                   aaaaGCCGATAAGGCAATGGCGCTGCGGTTTG‑GTGGAAATAAAGGTCTTTGGCGGTAGGAGAGCCACGACTTAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCGGAGGGGGATGGTTGGATTGTCCCTCCCGCAGCTGGCCCGTGGTACTTCGGATTGGTTGTGCTACGTAAGCCCCCTTCTGGCATTACGGTGTTGCCGTTGGTTGATGGAGATAAAGCCGCGCTAGCCGCCTCATTaa   <  2:211074/250‑1 (MQ=255)                                                                                                                                                                                                                      aGCCGATAAAGCAATGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTaaa  >  2:329622/1‑249 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                             TAAACCAGACGTTGCCAGATCGGCAAACGATGAGGTCAGCACCCTGTTTTTTGGTCATGATGATCGTGTGCCAGTGAATGACACGACCCAATCACCGTGGGATGCGGTTGGGCAACTGGAAACGGCCAGCGGCAATTTATGTACGGCGACGCTGATTGCACCCAATCTGGCATTAACGGCAGGACACTGTTTATTGACACCTCCAAAGGGTAAAGCCGATAAAGCAGTGGCGCTGCGTTTTGTGTC‑AAATAAAGGTCTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAA  >  NC_000913/1672025‑1672486

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