Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,201,820	A→G	K125K (AAA→AAG)	yehP →	VWA domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,201,820	0	A	G	100.0%	91.7 / NA	26	K125K (AAA→AAG)	yehP	VWA domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/10);  total (16/10)

GCAGGTGATTGAGCGGCTGGAAAGCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAACATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGACTGGCAAAGGAAGTTCGTCAGGCTTTTTCTGGTGTCCGC  >  NZ_CP009273/2201691‑2201934                                                                                                                                  |                                                                                                                   gcagGTGATTGAGCGGCTGGAAAGCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCtgat                                                                                                            >  1:132725/1‑138 (MQ=255)      tGATTGAGCGGCTGGAAAGCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAAccc                                                                                                      <  2:335338/139‑1 (MQ=255)        aTTGAGCGGCTGGAAAGCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCg                                                                                                    <  1:332474/139‑1 (MQ=255)           gAGCGGCTGGAAAGCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGaa                                                                                                  <  1:332201/138‑1 (MQ=255)               ggCTGGAAAGCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACt                                                                                             >  1:115816/1‑139 (MQ=255)               ggCTGGAAAGCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCGTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACt                                                                                             >  1:57637/1‑139 (MQ=255)                        gCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTgccgc                                                                                    <  2:449194/139‑1 (MQ=255)                           aTGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTgccgc                                                                                    <  1:175310/136‑1 (MQ=255)                           aTGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTgccgc                                                                                    >  2:175310/1‑136 (MQ=255)                           aTGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGccgcccg                                                                                 <  1:349243/139‑1 (MQ=255)                                           cGGAATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGGTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGAACTGGCTGCCGCCCGTCGGATAGTGTGCCAg                                                                 >  1:129661/1‑139 (MQ=255)                                                        gtggtgACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGaaa                                                    >  2:535590/1‑139 (MQ=255)                                                        gtggtgACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGaaa                                                    >  1:460905/1‑139 (MQ=255)                                                          ggtgACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAAt                                                   >  2:579274/1‑138 (MQ=255)                                                           gtgACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTa                                                 >  1:183622/1‑139 (MQ=255)                                                             gACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATg                                               >  1:406925/1‑139 (MQ=255)                                                                   tctcGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGa                                         >  1:321096/1‑139 (MQ=255)                                                                       gACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGa                                         <  2:183622/135‑1 (MQ=255)                                                                       gACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGACt                                       <  2:115816/137‑1 (MQ=255)                                                                                     tATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGACTGGCAAAGGAAGTTCGt                       >  1:241031/1‑139 (MQ=255)                                                                                      aTGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGa                                         >  1:554171/1‑120 (MQ=255)                                                                                      aTGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGa                                         <  2:554171/120‑1 (MQ=255)                                                                                               tCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGAt                                                                           <  1:452343/77‑1 (MQ=255)                                                                                               tCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGAt                                                                           >  2:452343/1‑77 (MQ=255)                                                                                               tCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGACTGGCAAAGGAAGTTCGTCAGGCttttt             >  2:568488/1‑139 (MQ=255)                                                                                                          gCTACGCGCTGTTTTGCACACCAAGCATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGACTGGCAAAGGAAGGTCGTCAGGCTTTTTCTGGTGTCcgc  >  2:175436/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                   GCAGGTGATTGAGCGGCTGGAAAGCGATGCCGTGCTGCGCTACGGCATTGAAGATGTGGTGACGAATCTCGACGTGCTGGAACGTATGCAGCCTTCTGAAAGCCTGCTACGCGCTGTTTTGCACACCAAACATCTGATGAACCCCGAAGTACTGGCTGCCGCCCGCCGGATAGTGTGCCAGGTTGTTGAAGAAATTATGGCTCGACTGGCAAAGGAAGTTCGTCAGGCTTTTTCTGGTGTCCGC  >  NZ_CP009273/2201691‑2201934

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