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%	103.7 / NA	33	*412* (TAG→TAA)	aslB	putative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/15);  total (18/15)

ATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAG  >  NC_000913/3984060‑3984327                                                                                                                                      |                                                                                                                                       attattTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATc                                                                                                                                   >  1:769092/1‑139 (MQ=255)          gtgCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgc                                                                                                                          <  2:769092/139‑1 (MQ=255)               gggTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATcggcgg                                                                                                                              <  1:461925/130‑1 (MQ=255)               gggTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATcggcgg                                                                                                                              >  2:461925/1‑130 (MQ=255)                           tATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCgg                                                                                                         <  2:289446/139‑1 (MQ=255)                                  gCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcggc                                                                                                  >  1:327015/1‑139 (MQ=255)                                        tACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAAc                                                                                            >  1:104506/1‑139 (MQ=255)                                                tatCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                    >  1:655148/1‑139 (MQ=255)                                                  tCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                  <  2:670602/139‑1 (MQ=255)                                                            aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGGCATTAGGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                        >  2:537957/1‑139 (MQ=255)                                                            aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTa                                                                                               >  1:581037/1‑116 (MQ=255)                                                            aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTa                                                                                               <  2:581037/116‑1 (MQ=255)                                                            aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                        >  1:296967/1‑139 (MQ=255)                                                            aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                        >  2:644362/1‑139 (MQ=255)                                                             tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                  >  1:793429/1‑128 (MQ=255)                                                             tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                  <  2:793429/128‑1 (MQ=255)                                                             tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGAt                                                                       >  1:296177/1‑139 (MQ=255)                                                              ggCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                  >  2:661245/1‑127 (MQ=255)                                                              ggCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                  <  1:661245/127‑1 (MQ=255)                                                                                   tCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCa                                                 <  1:644362/139‑1 (MQ=255)                                                                                     gCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGc                                               >  1:714115/1‑139 (MQ=255)                                                                                       ccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATcggcgg                                                                                                                              >  2:742008/1‑58 (MQ=255)                                                                                       ccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATcggcgg                                                                                                                              <  1:742008/58‑1 (MQ=255)                                                                                           ccAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTg                                         <  2:104506/139‑1 (MQ=255)                                                                                           ccAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTg                                         >  1:710922/1‑139 (MQ=255)                                                                                                         catgcGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGa                           <  2:710922/139‑1 (MQ=255)                                                                                                              gcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGt                      <  2:296177/139‑1 (MQ=255)                                                                                                                 tttGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGt                                                                    <  2:340834/90‑1 (MQ=255)                                                                                                                 tttGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGt                                                                    >  1:340834/1‑90 (MQ=255)                                                                                                                              agtaagtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTg                                                   >  2:246272/1‑94 (MQ=255)                                                                                                                              agtaagtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTg                                                   <  1:246272/94‑1 (MQ=255)                                                                                                                              agtaagtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGc      >  1:782982/1‑139 (MQ=255)                                                                                                                                  agtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAg  <  2:296967/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       ATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAG  >  NC_000913/3984060‑3984327

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