Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,985,729	C→T	L93L (CTG→CTA)	aslA ←	putative sulfatase AslA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,985,729	0	C	T	100.0%	133.4 / NA	37	L93L (CTG→CTA)	aslA	putative sulfatase AslA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/21);  total (16/21)

CAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTGCGCTGGCCATGAGAGGAAGAGCGGCGGCGACA  >  NC_000913/3985483‑3985969                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 caACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGc                                                                                                                                                                                                                                               >  1:322121/1‑250 (MQ=255)      gCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGGTGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAATTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGa                                                                                                                                                                                                                                           <  1:431780/249‑1 (MQ=255)        ggCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAaac                                                                                                                                                                                                                                       <  2:63580/251‑1 (MQ=255)            gCGTGGTTAACCCTTGCAGCCCGACCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGGCCATCCAGCCCACATCGTCTAGCAAGAAAACAAcc                                                                                                                                                                                                                                   >  1:117652/1‑251 (MQ=255)                          tGCAGCCGGCCGGGTTGCCGTTACATGGGCGGCTTAAGAATGCCGGGGTGGATGGAGTATTGTCGGG‑TGAGAAGGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCGGGCTGGCACCGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTc                                                                                                                                                                                                                       <  2:91966/249‑1 (MQ=255)                           ccAGCCCGCCGGGTTGCCCGTCCATGGGGGGCATCAGAAT‑CCGGGGGGGGTGGGATAATTTCCGGGTAAGAATGTGGGCGCGGTTTGGGAAAG‑GGTTGGTTTGAAAATCCCCCAAATTTAAAATCACCCCCTGCTTGGCAACGCCGCCGATTTCTGTGGAAGGTTTCCCCACCGCCACGCCGC‑ACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTa                                                                                                                                                                                                                    <  1:240474/249‑1 (MQ=255)                                  gcccgGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGttt                                                                                                                                                                                                             >  2:29238/1‑251 (MQ=255)                                            ccGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCa                                                                                                                                                                                                     >  1:121696/1‑249 (MQ=255)                                              gTACATTGGCGGCATCGGAATGCCGTGGTGGATGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGtt                                                                                                                                                                                                 <  1:145294/251‑1 (MQ=255)                                                ccATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGT‑TGAGAATGGTGGCGCGGTGTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTCCCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTc                                                                                                                                                                                                <  1:1894/249‑1 (MQ=255)                                                               gagtggcggggtggAGGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGAAAGAGCTGGGTT‑GAGAATACGCCGAAGTTAAAACCAGCCCCGGGCTGGCAACGGCGTCGATATCGGGTGTAGGGTTCCCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTgg                                                                                                                                                                                <  2:26430/242‑1 (MQ=255)                                                               gAATGCCGGGGGGGAGGGAGTATTGGCCGG‑TGAGAATGGTGGTGCGGGTTGGGGAAGAGCTTGGTG‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCGGGCAACGGCGTCGATATCTGGTGGAGGGTTACCCCCCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGGTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTgg                                                                                                                                                                                <  1:418159/251‑1 (MQ=255)                                                                          tggatggaGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATc                                                                                                                                                                     >  1:459655/1‑251 (MQ=255)                                                                                 aGTATTGTCCGG‑TGAGAATGGTGGCGCGGGATGGGGACGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGcc                                                                                                                                                              >  1:195711/1‑251 (MQ=255)                                                                                      tGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGCGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCggctg                                                                                                                                                         >  2:61735/1‑248 (MQ=255)                                                                                         ccGG‑TGAGAATGGTGGAGCGGGTTGGGGAAGAGCTGGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTGTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGAtgct                                                                                                                                                       <  2:282317/250‑1 (MQ=255)                                                                                                             gggTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTCCCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTg                                                                                                                                   <  1:76491/251‑1 (MQ=255)                                                                                                             gggTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCACTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTg                                                                                                                                   >  1:326542/1‑251 (MQ=255)                                                                                                              ggTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTGTTTCTCAAGTTCTGCCAGCGTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGt                                                                                                                                  <  2:56644/251‑1 (MQ=255)                                                                                                               gTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTc                                                                                                                                 >  1:365930/1‑251 (MQ=255)                                                                                                                 tGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTAAAAACCAGCCCCGGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCGTCATAGTGTCg                                                                                                                                <  2:263435/250‑1 (MQ=255)                                                                                                                   gggAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGGGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCa                                                                                                                             <  1:11842/251‑1 (MQ=255)                                                                                                                    ggAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTCTCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCa                                                                                                                             >  1:220248/1‑250 (MQ=255)                                                                                                                        gagCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCACTAgg                                                                                                                        >  1:59639/1‑250 (MQ=255)                                                                                                                             tGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACAGACGTCCATCCAGCCCACATCGTCTAGCAAGAAACCAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCCTATTGTCGGCAATAGTGGTCg                                                                                                                   >  2:82743/1‑251 (MQ=255)                                                                                                                                               aaGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATAc                                                                                                  >  2:350745/1‑251 (MQ=255)                                                                                                                                                aGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGATTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGCTACt                                                                                                 >  2:114242/1‑251 (MQ=255)                                                                                                                                                      aaTCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTg                                                                                           >  2:250452/1‑251 (MQ=255)                                                                                                                                                                 ggctggcAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTa                                                                                 <  1:69445/250‑1 (MQ=255)                                                                                                                                                                          cGGCGTCGATACCGGGTGTAGGGTTCCCCACCGCCACGCCGCCTCCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCATATACTGGTTTGGGTGATCGTATCCGTCAAAg                                                                       <  2:19603/251‑1 (MQ=255)                                                                                                                                                                            gCGTCGATGTCTGGTGGAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGAGACTGGTTTGGGTGATCGTATCCGGCAAAGcc                                                                     <  2:47725/251‑1 (MQ=255)                                                                                                                                                                             cGTCGATATCTGGGGTAGGGTTCCCCCCCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCt                                                                    <  2:59639/251‑1 (MQ=255)                                                                                                                                                                                         gtgtAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTTCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTg                                                        <  1:284833/251‑1 (MQ=255)                                                                                                                                                                                           gtAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTgg                                                       <  1:393810/250‑1 (MQ=255)                                                                                                                                                                                                                 cgccACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTgcg                                >  2:229905/1‑251 (MQ=255)                                                                                                                                                                                                                                cGACGTCCATCCTGCCCACATCGTCTAGCAAGAAAGCAACCACTGTCGTTTTCTTGCCGGTTTTTTGCTCAAGTTCTACCGGCTTCTGCTGGGTTACTTTATTCTGCGCCGGATGCTGCATTACTGGCATCTTATTGTCGGCAATAGTTGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTGTGCGTGCGGTTGCAGTTGACGGGGTATCTGCTGCGCTGGCCATGAGAGGa                 <  1:285276/251‑1 (MQ=255)                                                                                                                                                                                                                                               ccACATCGTCTAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTGCGCTGGCCATGAGAGGAAGAGCGGCGGCGACa  <  1:250452/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 CAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGG‑TGAGAATGGTGGCGCGGGTTGGGGAAGAGCTTGGTT‑GAGAATACGCCGAAGTTAAAATCAGCCCCTGGCTGGCAACGGCGTCGATATCTGGTGTAGGGTTACCCACCGCCACGCCGCCACCGTTAAAACCGACGTCCATCCAGCCCACATCGTCCAGCAAGAAAACAACCACATTCGGTTTCTTACCGGTTTTTTTCTCAAGTTCTGCCAGCTTCTGCTGGGTTTCTTTATCCTGCGCCGGATGCTGCATTACTGGCATCATATTGTCGGCAATAGTGGTCGCCGGTTTAACCAGATACTGGTTTGGGTGATCGTATCCGGCAAAGCCTTTGCGTGCGGTGGCAGTTGACGGGGTATCTGCTGCGCTGGCCATGAGAGGAAGAGCGGCGGCGACA  >  NC_000913/3985483‑3985969

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