Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,985,771	C→T	E79E (GAG→GAA)	aslA ←	putative sulfatase AslA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,985,771	0	C	T	100.0%	36.1 / NA	13	E79E (GAG→GAA)	aslA	putative sulfatase AslA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/6);  total (7/6)

GGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTGCGCTGGCCATGAGAGGAAGAGCGGCGGCGACAGCAACAACAAGACGTTTG  >  NC_000913/3985554‑3985987                                                                                                                                                                                                                          |                                                                                                                                                                                                                          ggtggATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTTT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCttt                                                                                                                                                                                          >  1:380084/1‑250 (MQ=255)                gTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTGTTTTTT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGAtgc                                                                                                                                                                          <  1:473270/251‑1 (MQ=255)                                  cgcgGGTGGGGGAAGAGCTTGGTTGAGAAATCGCCGAAGTTAAAATCAGCCCCGGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTTTTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTgcgn                                                                                                                                                                                 <  1:306532/227‑2 (MQ=255)                                                                              aTCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGATACCCACCGCCACGCAGCAACCGTTAAAACCGACCTCCATCCAGCCCACATCGTCCAGCAAGAAACCAACAACCTTCGGCTTCTTACCGCTTTTTTTTT‑CACGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGcca                                                                                                                                                                               >  1:250629/1‑183 (MQ=255)                                                                                                                                          ncACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTTT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCCTTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATc                                                                                                      >  1:52099/2‑197 (MQ=255)                                                                                                                                           cACCGTTAAAACCGACGTCTATCCACCCCACCTCGTCCAGCAAGAAAACAACCACATTCGGATTCTTACCGGTTTTTTTTTTCAAGTTCTGCAAGCTTCCGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCCTCATATTGTCGGCAATAGTGGTCGCGGGTTTAACCAGATACTGGTTTGGGTGAGCGTATCCGGCACAGCCTTTGCGTGCGGGGGCAGTTGAGGGGGTATCTGCTgcgcg                                                >  1:412870/1‑250 (MQ=255)                                                                                                                                                             ccatccaGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTGTTTTT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTa                                                                                                                                                                                         <  1:477249/95‑1 (MQ=255)                                                                                                                                                                     cccACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTTT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATGCTGGTTGTGTTGATCGTTGCCGGCAAAGCCTGTTCGTGCTGTGGCAGTTGACGGGGTATCTGCTGCGCGGGCCATGAGAGGAAGAGCGGCGGCGAc                     <  1:473454/251‑1 (MQ=255)                                                                                                                                                                         caTCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTTTTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGNCGCNNNNNTannnnnn                                                                                                                       >  2:305054/1‑144 (MQ=255)                                                                                                                                                                             gTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTTT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTGCGCTGGCCATGAGAGGAAGAGCGGCGGCGACAGcaaca              >  1:41225/1‑250 (MQ=255)                                                                                                                                                                             gTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTTT‑CAAGTTATGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCCTTCCTGGCATCATATTGTCGGCAATAGCGGTCGCCGGTTTAACCAGATa                                                                                                                     >  1:476713/1‑147 (MQ=255)                                                                                                                                                                                        aaaCAACCACATTCGGTTTCTTACCGGTTTTTTTTT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTGCGCTGGCCATGAGAGGAAGAGCGGCGGCGACAGCAACAACAAGACGTTTg  <  1:482880/251‑1 (MQ=255)                                                                                                                                                                                                                 ggTTTTTTTTT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTGCGCtn                                               <  1:97666/181‑2 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                          GGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTTGAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCT‑CAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTGCGCTGGCCATGAGAGGAAGAGCGGCGGCGACAGCAACAACAAGACGTTTG  >  NC_000913/3985554‑3985987

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