Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,581,323	G→A	T401I (ACA→ATA)	ydeN ←	putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,581,323	0	G	A	100.0%	124.4 / NA	36	T401I (ACA→ATA)	ydeN	putative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/20);  total (16/20)

TATCCAGGTCAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCC  >  NC_000913/1581205‑1581452                                                                                                                       |                                                                                                                                  tATCCAGGTCAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCg                                                                                                               >  2:496245/1‑139 (MQ=255) tATCCAGGTCAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCg                                                                                                               >  2:417655/1‑139 (MQ=255)   tCCAGGTCAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGa                                                                                                             >  1:605104/1‑139 (MQ=255)          cAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGc                                                                                                      <  1:425787/139‑1 (MQ=255)           aGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCt                                                                                                     <  2:703886/139‑1 (MQ=255)                  tATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATa                                                                                              >  2:312491/1‑139 (MQ=255)                       ggCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTAc                                                                                         >  2:243080/1‑139 (MQ=255)                                            tGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTtttcctttcc                                                                    >  1:226684/1‑139 (MQ=255)                                            tGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTtttcctttcc                                                                    <  1:312491/139‑1 (MQ=255)                                            tGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTtttcctttcc                                                                    <  1:417655/139‑1 (MQ=255)                                            tGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCACCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTtttcctttcc                                                                    <  1:590313/139‑1 (MQ=255)                                             gCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTcca                                                                   <  1:496245/139‑1 (MQ=255)                                              cAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTAc                                                                                         >  1:205488/1‑116 (MQ=255)                                              cAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTAc                                                                                         <  2:205488/116‑1 (MQ=255)                                              cAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTccac                                                                  <  1:243080/139‑1 (MQ=255)                                                   aGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACAt                                                             <  2:172916/139‑1 (MQ=255)                                                     ggCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACAt                                                             >  1:444204/1‑137 (MQ=255)                                                     ggCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACAt                                                             <  2:444204/137‑1 (MQ=255)                                                     ggCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATaa                                                           <  2:226684/139‑1 (MQ=255)                                                           aaGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTg                                                     <  2:164133/139‑1 (MQ=255)                                                            aGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTgg                                                    <  1:351832/139‑1 (MQ=255)                                                               aaaCGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGTAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGgt                                                 >  2:480232/1‑139 (MQ=255)                                                               aaaCGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGgt                                                 >  1:126705/1‑139 (MQ=255)                                                                                      tttGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTccac                                                                  >  2:425836/1‑99 (MQ=255)                                                                                      tttGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTccac                                                                  <  1:425836/99‑1 (MQ=255)                                                                                      tttGGAATGCTGATATCGGCTGCATCAAGAGCTATC‑GGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTg                         >  2:762999/1‑139 (MQ=255)                                                                                       ttGGAATGCTGATATCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTg                         <  1:408173/139‑1 (MQ=255)                                                                                                  atatCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGc                                      >  1:682220/1‑115 (MQ=255)                                                                                                  atatCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGc                                      <  2:682220/115‑1 (MQ=255)                                                                                                  atatCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGc              >  2:275456/1‑139 (MQ=255)                                                                                                   tatCGGCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGcc             <  2:33618/139‑1 (MQ=255)                                                                                                       ggCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGt                                           <  1:695142/105‑1 (MQ=255)                                                                                                       ggCTGCATCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGt                                           >  2:695142/1‑105 (MQ=255)                                                                                                         cTGCATCAACAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACAt                                                             >  1:525285/1‑85 (MQ=255)                                                                                                         cTGCATCAAAAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTTCTTTCCACCACAt                                                             <  2:525285/85‑1 (MQ=255)                                                                                                              tCAAGAGCTATCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGccc  <  1:594623/139‑1 (MQ=255)                                                                                                                       |                                                                                                                                  TATCCAGGTCAGATTTTTATGTGGCTCGCCTTGTTTCTTATCTTGCAACCAGGGCAGCAAGGAAACGCCATCCAGCTTAAGGTCTTTTGGAATGCTGATATCGGCTGCATCAAGAGCTGTCGGGTAGAAATCCATTGCGGAAATCAGCTTGTCATAATTACCGGGTTGAAGTTTTCCTTTCCACCACATAAACATTGGGGTGTGAGTACCGCCAGGATAGGTCTGACTCTTATAGCCTTTTTGCGCCC  >  NC_000913/1581205‑1581452

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