Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,983,841	A→G	Y295C (TAC→TGC)	aslB →	putative AslA‑specific sulfatase‑maturating enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,983,841	0	A	G	100.0%	184.2 / NA	51	Y295C (TAC→TGC)	aslB	putative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/20);  total (31/20)

GACCACTGGATCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAA  >  NC_000913/3983693‑3983985                                                                                                                                                     |                                                                                                                                                 gACCACTGGATCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCg                                                                                                                                                >  1:521441/1‑151 (MQ=255)        ggATCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGa                                                                                                                                         <  2:740353/151‑1 (MQ=255)              aaCGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACt                                                                                                                                    >  1:479200/1‑150 (MQ=255)                   gACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGttt                                                                                                                              >  1:266507/1‑151 (MQ=255)                       tCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATcc                                                                                                                          >  1:209034/1‑151 (MQ=255)                            aaGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAat                                                                                                                     >  1:608045/1‑151 (MQ=255)                             aGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAat                                                                                                                     >  2:664554/1‑150 (MQ=255)                                 tttCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCg                                                                                                                 <  2:640609/150‑1 (MQ=255)                                   tCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCgg                                                                                                                <  2:622115/149‑1 (MQ=255)                                   tCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGc                                                                                                               <  2:591351/150‑1 (MQ=255)                                    cGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTg                                                                                                             >  1:633937/1‑151 (MQ=255)                                               tCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGc                                                                                                   >  2:133532/1‑150 (MQ=255)                                                  aaTGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCa                                                                                               <  1:213809/151‑1 (MQ=255)                                                       gCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCaa                                                                                           >  2:322718/1‑150 (MQ=255)                                                       gCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCa                                                                                            >  1:449492/1‑149 (MQ=255)                                                            tGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAAt                                                                                      <  2:160098/150‑1 (MQ=255)                                                             gCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTg                                                                                    >  2:335747/1‑151 (MQ=255)                                                              ccAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGc                                                                                   <  2:98066/151‑1 (MQ=255)                                                               cAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCa                                                                                  <  1:524053/151‑1 (MQ=255)                                                                aaTTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGc                                                                                   >  1:107577/1‑149 (MQ=255)                                                                 aTTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGa                                                                                <  2:534715/151‑1 (MQ=255)                                                                   tttGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGaa                                                                               <  2:72101/150‑1 (MQ=255)                                                                     tGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATg                                                                            >  2:760828/1‑151 (MQ=255)                                                                      gTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATg                                                                            >  2:761885/1‑150 (MQ=255)                                                                       tCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGTGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATg                                                                            >  2:665572/1‑149 (MQ=255)                                                                          gTGCGCCGGGTGCGGTTTTCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCGGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGa                                                                           >  2:767198/1‑147 (MQ=255)                                                                           tGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGAt                                                                      >  1:394782/1‑151 (MQ=255)                                                                                         tttGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAAcagc                                                         >  1:703459/1‑150 (MQ=255)                                                                                              catcatCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGt                                                   <  1:567865/151‑1 (MQ=255)                                                                                                 catcaGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGAACGATTCCCCGCAACAGCAGGCGttt                                                 <  1:632481/150‑1 (MQ=255)                                                                                                   tcaGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGt                                              >  2:523706/1‑151 (MQ=255)                                                                                                    caGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGt                                              <  2:132262/150‑1 (MQ=255)                                                                                                    caGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTg                                             >  1:38361/1‑151 (MQ=255)                                                                                                     aGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTg                                             >  2:674376/1‑150 (MQ=255)                                                                                                     aGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGa                                            <  2:685874/151‑1 (MQ=255)                                                                                                         aaCCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGa                                            >  1:181496/1‑147 (MQ=255)                                                                                                            cTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGGTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACATCAGGCGTTTGGAGAAGATaa                                      >  2:201618/1‑150 (MQ=255)                                                                                                              gtgGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATaaa                                     <  2:214124/149‑1 (MQ=255)                                                                                                                gggCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTa                                 <  2:617215/151‑1 (MQ=255)                                                                                                                    gcTCGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGc                              >  1:569506/1‑150 (MQ=255)                                                                                                                       cGGTGATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGtt                          >  1:669280/1‑151 (MQ=255)                                                                                                                           gATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTAcc                       <  1:373408/150‑1 (MQ=255)                                                                                                                           gATTGTTGAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCg                      >  1:617761/1‑151 (MQ=255)                                                                                                                             ttgttgAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGgc                    >  2:421862/1‑151 (MQ=255)                                                                                                                               gttgAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGgcgc                  <  1:76305/151‑1 (MQ=255)                                                                                                                                  gAGCACAACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAgt               >  1:463770/1‑151 (MQ=255)                                                                                                                                      acaACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCg           >  1:82412/1‑151 (MQ=255)                                                                                                                                       caACGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAACATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGc          <  1:756251/151‑1 (MQ=255)                                                                                                                                          cGGCGACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGt       >  2:20402/1‑151 (MQ=255)                                                                                                                                              gACGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGGGTCGCAGTTGt    >  2:352357/1‑150 (MQ=255)                                                                                                                                               aCGTTTGCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTaa  <  1:155830/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GACCACTGGATCAAACGCGACGTCGGCAAGATTTTCGTGATGAATATCGAATGGGCGTTTGCCAATTTTGTCGGTGCGCCGGGTGCGGTTTGCCATCATCAGCCAACCTGTGGGCGCTCGGTGATTGTTGAGCACAACGGCGACGTTTACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAA  >  NC_000913/3983693‑3983985

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