Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,984,040	C→T	A361A (GCC→GCT)	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,984,040	0	C	T	100.0%	200.0 / NA	54	A361A (GCC→GCT)	aslB	putative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/23);  total (31/23)

AGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGA  >  NC_000913/3983892‑3984185                                                                                                                                                     |                                                                                                                                                  aGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAg                                                                                                                                                 <  2:463749/151‑1 (MQ=255)  gCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGc                                                                                                                                                >  2:686499/1‑151 (MQ=255)      aCAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCa                                                                                                                                            <  2:617724/151‑1 (MQ=255)        aaTTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCaaa                                                                                                                                          <  1:207199/151‑1 (MQ=255)         aTTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAAc                                                                                                                                         >  2:530170/1‑151 (MQ=255)           tGCAGACATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCg                                                                                                                                       >  2:128442/1‑151 (MQ=255)            gCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCgg                                                                                                                                      >  2:30000/1‑151 (MQ=255)            gCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCgg                                                                                                                                      <  2:707715/151‑1 (MQ=255)                     aTCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAatt                                                                                                                             >  2:514131/1‑151 (MQ=255)                      tCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAatt                                                                                                                             >  2:465615/1‑150 (MQ=255)                      tCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAatt                                                                                                                             >  2:674953/1‑150 (MQ=255)                       cGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAattat                                                                                                                           >  1:469428/1‑151 (MQ=255)                              ccGCAACCGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAACGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTtgtgtg                                                                                                                    >  1:22008/1‑151 (MQ=255)                               cGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGc                                                                                                                   <  1:761931/151‑1 (MQ=255)                               cGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGc                                                                                                                   <  2:521417/151‑1 (MQ=255)                                 cAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCg                                                                                                                 <  1:674417/151‑1 (MQ=255)                                     agcagGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAATACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTa                                                                                                             <  2:171950/151‑1 (MQ=255)                                     agcagGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTa                                                                                                             <  1:53685/151‑1 (MQ=255)                                            gTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGt                                                                                                      >  2:295797/1‑151 (MQ=255)                                               tGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTa                                                                                                    <  1:328485/150‑1 (MQ=255)                                                ggTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTAtt                                                                                                  <  2:61704/151‑1 (MQ=255)                                                    aaGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCg                                                                                              <  1:316905/151‑1 (MQ=255)                                                         aaaTTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATc                                                                                         >  2:614184/1‑151 (MQ=255)                                                           aTTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTa                                                                                       <  2:145322/151‑1 (MQ=255)                                                               aaGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGc                                                                                   >  2:459543/1‑151 (MQ=255)                                                                    gTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCatat                                                                              >  1:611034/1‑151 (MQ=255)                                                                     ttACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATc                                                                             >  2:380709/1‑151 (MQ=255)                                                                      tACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCt                                                                            >  2:195354/1‑151 (MQ=255)                                                                                 tgtCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATgg                                                                 >  2:765544/1‑151 (MQ=255)                                                                                  gccgcgtttGTAACGTGTAAAAAGCGTGCTGGGGAGGCCGCCCAAAACCCCGCTTCATGCTCGATGCTAGCGGCAACCCGGGATTGAATTATTTGTGGGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATgg                                                                 <  1:625729/143‑1 (MQ=255)                                                                                    cGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTg                                                              >  1:96472/1‑151 (MQ=255)                                                                                         ttGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTg                                                         >  2:755222/1‑151 (MQ=255)                                                                                         ttGTAACGTGTTAAAAGCGTGCGGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTTTTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGAtt                                                           <  1:346623/149‑1 (MQ=255)                                                                                         ttGTAACGTGTAAAAAGCGTGCTGGGGAGGCTGCCCGAAACCCCGCTTCATGCTCGATGCTAGCGGCAAACCGGGATTGAATTATTTGTGTGCCGGGTATCAGCTTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGAttt                                                          <  1:335761/150‑1 (MQ=255)                                                                                          tGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGc                                                        >  2:292509/1‑151 (MQ=255)                                                                                             aaCGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTgg                                                     >  1:224274/1‑151 (MQ=255)                                                                                                       aaGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCg                                           >  1:619767/1‑151 (MQ=255)                                                                                                          ggTGCGGGGGAGGCTGCCCGAAACCCCGTTTCATGCTGGATGCTAGCGCCAACCCGGGACTGAATTATTTGTTTGCCGGGTATCAGCTTTTTTCCCCCCTTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGccc                                        <  2:97392/150‑1 (MQ=255)                                                                                                           gTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCg                                       <  1:258477/151‑1 (MQ=255)                                                                                                             gCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGc                                     >  2:346904/1‑151 (MQ=255)                                                                                                             gCTGGGGAGGCGGCCAAAAACCCCGCTTCATGCTCGATGCTAGCGGCAAACCGGGATGAAATTATTTGTGTGCCGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCg                                       <  2:633900/149‑1 (MQ=255)                                                                                                              cTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGcc                                    >  1:352954/1‑151 (MQ=255)                                                                                                                ggggAGCCTGCCGGAAACACCGCTTCATGCTCGAGGCTGGCGGCAACCCGGGACTGAATTATTTGTGTGCCGGGTATAAGCGTTATTTCCGCCATCTCCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCgg                                      <  2:449471/147‑1 (MQ=255)                                                                                                                 gggAGGCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGc                                 >  1:153422/1‑151 (MQ=255)                                                                                                                      gCTGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACAt                            >  1:702493/1‑151 (MQ=255)                                                                                                                        tGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACAtt                           >  1:726101/1‑150 (MQ=255)                                                                                                                        tGCCCGAAACACCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACAtt                           >  1:446882/1‑150 (MQ=255)                                                                                                                              aaaCACCGCTTCATGCTCGATGCTAGCGGCAACCCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCTTTTTCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTatgcat                     <  1:530007/150‑1 (MQ=255)                                                                                                                                  acCGCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATgcgc                 >  1:493965/1‑150 (MQ=255)                                                                                                                                     gCTTCATGCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCAttt             >  2:702332/1‑151 (MQ=255)                                                                                                                                            gCTCGATGCTAGCGGCAACCCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCt          <  1:20402/147‑1 (MQ=255)                                                                                                                                            gCTCGATGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCtggtg      >  2:753937/1‑151 (MQ=255)                                                                                                                                                gATGCTAGCGCAAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCTTTATTCCCCCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGa  <  2:611421/151‑1 (MQ=255)                                                                                                                                                 aTGCTAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGgtggtg   >  1:466531/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGGTGAAGATAAATTTAAGCAGTTACCGGCGCAGTGTCGCAGTTGTAACGTGTTAAAAGCGTGCTGGGGAGGCTGCCCGAAACACCGCTTCATGCTCGATGCCAGCGGCAAACCGGGACTGAATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGA  >  NC_000913/3983892‑3984185

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