Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,984,193	G→A	*412* (TAG→TAA)	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,193	0	G	A	100.0%	134.8 / NA	42	*412* (TAG→TAA)	aslB	putative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/20);  total (22/20)

GTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGG  >  NC_000913/3984069‑3984322                                                                                                                             |                                                                                                                                  gtgCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgc                                                                                                                     >  2:377576/1‑139 (MQ=255) gtgCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgc                                                                                                                     >  1:413924/1‑139 (MQ=255)           tCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGa                                                                                                           <  1:662312/139‑1 (MQ=255)                   aTTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGa                                                                                                   >  1:337985/1‑139 (MQ=255)                      tCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgc                                                                                                <  2:751695/139‑1 (MQ=255)                       ccGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcg                                                                                               <  2:188554/139‑1 (MQ=255)                             tCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTaa                                                                                         <  2:648409/139‑1 (MQ=255)                                 ccGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGc                                                                                     <  2:276123/139‑1 (MQ=255)                                    ccATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCtt                                                                                  >  2:664677/1‑139 (MQ=255)                                     cATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTa                                                                                 <  1:573246/139‑1 (MQ=255)                                       tatCTTAATGTATTGGCGGATTTGCTGGCGCACGGTCGCCCGGCCAGCCACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATCCGGCGTAAACGCCTTATc                                                                               <  2:275776/139‑1 (MQ=255)                                       tatCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                               >  1:594360/1‑139 (MQ=255)                                       tatCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                               <  1:664677/139‑1 (MQ=255)                                         tCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             <  2:399378/139‑1 (MQ=255)                                           ttAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                           <  2:378621/139‑1 (MQ=255)                                              aaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                        >  1:492088/1‑139 (MQ=255)                                              aaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                        >  2:191575/1‑139 (MQ=255)                                               aGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             >  2:428001/1‑133 (MQ=255)                                               aGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             <  1:428001/133‑1 (MQ=255)                                                   aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                   >  2:679951/1‑139 (MQ=255)                                                    tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTa                                                                                 >  2:78992/1‑124 (MQ=255)                                                    tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTa                                                                                 <  1:78992/124‑1 (MQ=255)                                                    tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGAt                                                                  >  1:298329/1‑139 (MQ=255)                                                           tttGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             <  1:319470/121‑1 (MQ=255)                                                           tttGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             >  2:319470/1‑121 (MQ=255)                                                                        ggTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTg                                              <  2:1163/139‑1 (MQ=255)                                                                             cccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCg                                         <  1:191575/139‑1 (MQ=255)                                                                               cGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGtt                                       >  1:272493/1‑139 (MQ=255)                                                                                    aGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             <  1:468644/96‑1 (MQ=255)                                                                                    aGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             >  2:468644/1‑96 (MQ=255)                                                                                                catgcGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgc                                                                                                                     <  1:680958/44‑1 (MQ=255)                                                                                                catgcGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgc                                                                                                                     >  2:680958/1‑44 (MQ=255)                                                                                                 atgcGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             >  2:238632/1‑83 (MQ=255)                                                                                                 atgcGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             <  1:238632/83‑1 (MQ=255)                                                                                                   gcgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             <  1:199185/81‑1 (MQ=255)                                                                                                   gcgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                             >  2:199185/1‑81 (MQ=255)                                                                                                    cgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCGGGATGCGCCGTAAACGCTTTATCGGCCCTACTTGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCg                  >  2:218050/1‑139 (MQ=38)                                                                                                    cgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCg                  >  1:561525/1‑139 (MQ=255)                                                                                                    cgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCg                  >  1:340236/1‑139 (MQ=255)                                                                                                        tttGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAg              >  2:389125/1‑139 (MQ=255)                                                                                                         ttGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgc             >  1:269699/1‑139 (MQ=255)                                                                                                                    gAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCgg  <  2:561525/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                  GTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGG  >  NC_000913/3984069‑3984322

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