Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	500,604	G→A	G160G (GGG→GGA)	gsk →	inosine/guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	500,604	0	G	A	100.0%	155.9 / NA	45	G160G (GGG→GGA)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/23);  total (22/23)

GAACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGGGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCA  >  NC_000913/500359‑500846                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   gAACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACg                                                                                                                                                                                                                                               >  2:330477/1‑251 (MQ=255)                 aCCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGt                                                                                                                                                                                                                               <  1:357345/251‑1 (MQ=255)                       ttACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGCTCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGc                                                                                                                                                                                                                         >  1:159686/1‑251 (MQ=255)                         aCCCATAAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACAACTCGGTGCTCGCGGACGACCGTTCGGGGCTGCTGGGAGTCATGTGCAGAAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGcc                                                                                                                                                                                                                        >  1:191089/1‑250 (MQ=255)                          cccATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCcac                                                                                                                                                                                                                      <  1:282248/251‑1 (MQ=255)                          cccATCAGTGTGCGGTTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTGCCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCcac                                                                                                                                                                                                                      <  1:185267/251‑1 (MQ=255)                               cAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGaa                                                                                                                                                                                                                 <  2:237914/251‑1 (MQ=255)                                aGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGTTATCAGTCCAGGCCACATGAAc                                                                                                                                                                                                                >  1:234757/1‑251 (MQ=255)                                  tttGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGCTTGGTCGTTGCTTTACGATGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCCGGCCACATGAAcc                                                                                                                                                                                                               >  2:449590/1‑250 (MQ=255)                                      ggggTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGAAATATTGAAATTGGCAGGTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAgatg                                                                                                                                                                                                          <  2:159686/250‑4 (MQ=255)                                                           ccATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGATGCGGGCTGAAAGCATTCCGGa                                                                                                                                                                                      >  1:243996/1‑250 (MQ=255)                                                                 acaACTACTCGGTGCTCGAGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGAtgtg                                                                                                                                                                               >  1:362238/1‑251 (MQ=255)                                                                 acaACTACGCGGTGCTCGCGGACTACCGTCCGGGTCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGTCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGTCGTCCCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAACTGCGGGCTGAAAGCATTCCGGAAGAtgtt                                                                                                                                                                               <  2:225874/251‑2 (MQ=255)                                                                       actCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTTTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGcc                                                                                                                                                                         <  1:312298/251‑1 (MQ=255)                                                                                         acCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGtt                                                                                                                                                       >  1:11480/1‑251 (MQ=255)                                                                                            gTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTtctc                                                                                                                                                    >  1:341742/1‑251 (MQ=255)                                                                                                ggTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCt                                                                                                                                                >  2:320197/1‑251 (MQ=255)                                                                                                        tgGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTg                                                                                                                                        >  2:82218/1‑251 (MQ=255)                                                                                                          ggCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGt                                                                                                                                      <  1:103883/251‑1 (MQ=255)                                                                                                               cATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGtt                                                                                                                                 <  1:151867/251‑1 (MQ=255)                                                                                                                 tgtgCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGCTTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTAGGCACTGGTTCTCACCTCATATCTGGTGCGTTGc                                                                                                                               >  1:455204/1‑251 (MQ=255)                                                                                                                       gcaATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCg                                                                                                                         <  1:330477/251‑1 (MQ=255)                                                                                                                        caATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCgg                                                                                                                        >  1:282237/1‑251 (MQ=255)                                                                                                                            atTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGGGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTg                                                                                                                     <  1:449590/250‑1 (MQ=255)                                                                                                                                     ggCATTTATGCCTATCGTTACCTGGGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGcc                                                                                                           <  2:2027/251‑1 (MQ=255)                                                                                                                                           tATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGGGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGc                                                                                                     >  2:133853/1‑251 (MQ=255)                                                                                                                                                cTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAAcc                                                                                                 >  1:392272/1‑250 (MQ=255)                                                                                                                                                       tACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGana                                                                                           >  2:419970/1‑249 (MQ=255)                                                                                                                                                       tACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGaaa                                                                                           <  1:419970/249‑1 (MQ=255)                                                                                                                                                        aCCTGTGTAACACTTCCAGCCGTACCGATCTTAACCATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGcc                                                                                        <  2:239305/251‑1 (MQ=255)                                                                                                                                                                 acacTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTAc                                                                               <  1:82218/251‑1 (MQ=255)                                                                                                                                                                  cacTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAGCCATCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCTGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTAcg                                                                              <  1:252379/251‑1 (MQ=255)                                                                                                                                                                       ccAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGCAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCga                                                                           >  2:67783/1‑249 (MQ=255)                                                                                                                                                                              tACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAatat                                                                   >  1:261617/1‑250 (MQ=255)                                                                                                                                                                               aCCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGCAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGAGACGAAAtata                                                                  >  2:322009/1‑250 (MQ=255)                                                                                                                                                                                  gATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGt                                                              >  2:345762/1‑251 (MQ=255)                                                                                                                                                                                       tAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAAc                                                                <  2:341742/244‑1 (MQ=255)                                                                                                                                                                                                    ggCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCt                                            <  1:320197/251‑1 (MQ=255)                                                                                                                                                                                                         ggATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGc                                        <  1:359072/250‑1 (MQ=255)                                                                                                                                                                                                                                  aCGCTGATTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATcc              <  1:17958/251‑1 (MQ=255)                                                                                                                                                                                                                                        aTTGGCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAgtg        <  2:11480/251‑1 (MQ=255)                                                                                                                                                                                                                                            gCGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGGGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCGTGAAAGCCATTGAGTATGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATGGCCGAGAATCCGCAGtggtgg    <  2:167895/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGgc   >  2:80149/1‑251 (MQ=255)                                                                                                                                                                                                                                              gAGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCTAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGgca  <  1:80149/251‑1 (MQ=255)                                                                                                                                                                                                                                               aGTCCGGAGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGgc   >  1:407450/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   GAACTAAAACAGAAAAACCTGATTACCCATCAGTTTGCGGGTGGCACCATTGGTAACACCATGCACAACTACTCGGTGCTCGCGGACGACCGTTCGGTGCTGCTGGGCGTCATGTGCAGCAATATTGAAATTGGCAGTTATGCCTATCGTTACCTGTGTAACACTTCCAGCCGTACCGATCTTAACTATCTACAAGGCGTGGATGGCCCGATTGGTCGTTGCTTTACGCTGATTGGCGAGTCCGGGGAACGTACCTTTGCTATCAGTCCAGGCCACATGAACCAGCTGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCA  >  NC_000913/500359‑500846

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