Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,026	A→G	E266E (GAA→GAG)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,026	0	A	G	100.0%	151.5 / NA	53	E266E (GAA→GAG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/31);  total (22/31)

ACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCG  >  NC_000913/4618801‑4619270                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     actggCGGAAGGTCTCGACGAGCTGGGGATGGACGTGACAACGGGTAGCGGAGCTTTTAGGACGACCTACGAACTCTCTAGAGGACTTGCCGGAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGCGCGCACCACCGCGCTGCTCACCGACATGAATAAGGTACTTGCCACCAGTGCCGGTAACGCGGTTGAAGTTCGTGAAGCGGGGCCGGTCCTGAAGGGGGAGGATCGTAAcacgggtcgggttgagg                                                                                                                                                                                                                             >  2:141156/5‑235 (MQ=255)   tggCGGAAGGTCTGGACTCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCCGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTc                                                                                                                                                                                                                           >  2:251714/3‑251 (MQ=255)     gCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCAc                                                                                                                                                                                                                         >  1:261985/1‑251 (MQ=255)           ggTCTGGACGCGCGGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGTCCTACGAACTCTCTGAAGCCCTTGCCGAAGCTATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCGCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGgc                                                                                                                                                                                                                   <  1:120889/251‑1 (MQ=255)                ggACGCGCGGGTGATGGACGTGAAAGGGGGTAGCGGCGCGTTTATGCCCACCGACGACCTCTCTGATGCCCTTGCCGAAGCGGTTGGTGGCGTGGCTAACGGCGCTGGCGTGCGCACCGCCGCGCTGCTCACCGACATGAAACAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGGCTTGTTGAGGTCACGATGgcgatg                                                                                                                                                                                                               <  2:149722/250‑4 (MQ=255)                  aCGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCtgtgg                                                                                                                                                                                                            >  1:13513/1‑250 (MQ=255)                  aCGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGACGACCTACGAACTCTCTGAAGCCCTCGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGAACCACCGCGCTGCTCACCGACATGAATAAGGTACTGGCCCCCAGTGCCGGTAACGCGGTTGAAGTTCGTGAAGCGGGGCAGTTCCTGACGGGTGAGTATCGTAACCCGAGTCTGTTTGATGTCAAGAAGGCGCtgcgc                                                                                                                                                                                                            >  1:121193/1‑248 (MQ=255)                   cgcgCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGGTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGCTGGCGCTGTGCg                                                                                                                                                                                                           >  2:138219/1‑251 (MQ=255)                            gATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGc                                                                                                                                                                                                  <  1:238092/251‑1 (MQ=255)                             aTGGACGTGAAAGTGGGTAGCGGCGCGTTGAGGCCGACCTACGAAATCTCTGAAGGCCGTGCCGAAGCGATTGTTGGGGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTAATGGGATCCAGTGCGGGTAACGCTGTTGAAGTTCGTGAGGCGGTGCAGTTCCTGACGGGTGAGGATCGTAATCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCt                                                                                                                                                                                                 <  1:155378/251‑1 (MQ=255)                                     gAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCg                                                                                                                                                                                         >  2:125709/1‑251 (MQ=255)                                             gTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTg                                                                                                                                                                                 <  1:76317/251‑1 (MQ=255)                                                    cgcgTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAACCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGAGGAGATGCTTATCTCCGGCAAACTGGCGAAAg                                                                                                                                                                          <  2:31577/251‑1 (MQ=255)                                                         ttATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGAGTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGa                                                                                                                                                                      <  1:58979/250‑1 (MQ=255)                                                               cGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCCCCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCg                                                                                                                                                                 >  2:35008/1‑249 (MQ=255)                                                                   ctacgaactccctgccgtaattGCCGAAGCGCTTGTTGGCGTGGCTAACGGCGCTGGCGAGCGCACCCCCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGGGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATTTCACGATGGCGCTGTGCGTGGAGGTGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAacgc                                                                                                                                                           >  2:38462/21‑251 (MQ=255)                                                                            ctctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAt                                                                                                                                                  >  2:214890/1‑251 (MQ=255)                                                                                gAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGc                                                                                                                                               <  2:265862/250‑1 (MQ=255)                                                                                 aaGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTc                                                                                                                                                                                                                           >  2:191740/1‑173 (MQ=255)                                                                                 aaGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGGAACCCGCGTCTGTTTGATGTc                                                                                                                                                                                                                           <  1:191740/173‑1 (MQ=255)                                                                                  aGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCa                                                                                                                                              >  2:44294/1‑249 (MQ=255)                                                                                    cccTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCaagc                                                                                                                                           >  2:40009/1‑247 (MQ=255)                                                                                     ccTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCaagcag                                                                                                                                         <  1:25640/251‑6 (MQ=255)                                                                                            aaaGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTTCAGGTAACGCGGTTTAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCTTAACCCGCGTCTGTTTGATGTCACGATTGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGa                                                                                                                                  <  1:176887/250‑1 (MQ=255)                                                                                                    tgttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGt                                                                                                                           <  1:214890/250‑1 (MQ=255)                                                                                                        ggCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGc                                                                                                                      >  1:3961/1‑251 (MQ=255)                                                                                                                               gcgcACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATgg                                                                                               <  1:35008/251‑1 (MQ=255)                                                                                                                                         gcgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGcaca                                                                                     <  1:125709/251‑1 (MQ=255)                                                                                                                                                      aCATGAATCAGGTACTGGCCTCCAGTGCCGGTAGCGCGGTAGAAGTTCGTGAAGCGGTGCAGTTACTGACGGGTGAGGTTCGTAGCCCGCGTCTTTTTGATTACACGATGGCGCTGTGCGTGGAGATGCGTATCTCCGGCAAACTGGCTAAATATGACGCCGAAGCGCGCGCGAAATTGCAAGCATTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGGCGTATGGTAGCGTCACAAAAAGGCccaacc                                                                        <  1:54896/251‑1 (MQ=255)                                                                                                                                                         ggaaccAGGTACTGGCCTCCAGTGCAGTTAACGGGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGAATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTG‑AGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACt                                                                    <  1:29872/246‑1 (MQ=255)                                                                                                                                                           aaTCAGGTACGGGTCTCCATTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGCCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACTCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGTCAGAAGTCTTTGGTCGTATTGTAGCGGCACAAAAAGGACCAACCGACtt                                                                   <  2:163296/251‑1 (MQ=255)                                                                                                                                                            aTCAGGTACTGGCCTCCATTGCAGGTAACGCGGTTGAAGTGCGTGAAGCGGTGCAGTGCCTGACGGGGTAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTc                                                                  <  2:167958/251‑1 (MQ=255)                                                                                                                                                                   aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTgaga                                                           <  2:202468/251‑1 (MQ=255)                                                                                                                                                                   aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTgag                                                            <  1:40009/250‑1 (MQ=255)                                                                                                                                                                    cTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGAACAAAAAGGCACAACCGACTTCGTTgaga                                                           >  1:52639/1‑250 (MQ=255)                                                                                                                                                                              gTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGa                                                  >  1:92664/1‑249 (MQ=255)                                                                                                                                                                               ggCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGt                                               <  2:77050/250‑1 (MQ=255)                                                                                                                                                                                  aGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTGCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATc                                            <  1:104550/251‑1 (MQ=255)                                                                                                                                                                                      agcgcgGTTGAAGTGCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGAGGGCGCTGTGCGTGGAGATGCGTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACTGTAAAGCGGCAGAAGGCTTTGGTCGTATGGTATCTTCTCAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                        <  1:250517/249‑1 (MQ=255)                                                                                                                                                                                      aaCGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTGGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTGGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                        <  2:18693/251‑1 (MQ=255)                                                                                                                                                                                            gTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAACGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGCGAACTACGCGAAGTATCTGCCGACAGc                                  >  1:9/1‑251 (MQ=255)                                                                                                                                                                                               gAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAGTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGa                                <  1:33209/250‑1 (MQ=255)                                                                                                                                                                                                    tCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTg                           <  1:251714/250‑1 (MQ=255)                                                                                                                                                                                                      gTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACg                        >  1:68623/1‑251 (MQ=255)                                                                                                                                                                                                             ggTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGc                   >  1:5692/1‑249 (MQ=255)                                                                                                                                                                                                               tGCAGATCCTGACGGGTGAGTATCGAAACAAGAGTCTGTTTGCTGTTACGATGGAGCTGTGCGTGGAGATGCTTATCCCCGGCAAACTGACGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGAAGAGCTGGACAACGGTAAAGCGGCAGCAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGACCAACCGACTTCGTTGAGAACTCAGCGAAGCATCTGCCGACAGCGATGCGGACGAAAGCCGTc               >  2:129810/1‑251 (MQ=255)                                                                                                                                                                                                                 caggtccggACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCTCGAAATTGCAAGCAGTGCTGGATAACGGTAAAGCGGCAGAAGGCTTTGGTCTTATGGTAGCGTCACAAAAAGGCCCAACCGACTTCGTTGAGGACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTa             <  1:172523/243‑1 (MQ=255)                                                                                                                                                                                                                  aGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCTAAATTGCAAGCAGTGCTGGACCACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAt            <  2:68695/251‑1 (MQ=255)                                                                                                                                                                                                                   gTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTGGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATg           <  1:138219/251‑1 (MQ=255)                                                                                                                                                                                                                     tCCTGACGGGTGCGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCt         >  2:133100/1‑251 (MQ=255)                                                                                                                                                                                                                      ccTGACGGGTGAGTAGCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTTTCTCCGGCAAACTGGCGAAAGATTACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCTCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTg        <  1:141649/251‑1 (MQ=255)                                                                                                                                                                                                                         gACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGa       <  2:261985/249‑1 (MQ=255)                                                                                                                                                                                                                            gggTGAGTATCGTAACCCGCGTCTGTTTGATGTCCCGATGGCGCTGTGCGTGGAGATGCTTATCTCAGGCAAACTGGCGAAAGATGACGCCGCAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGAAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCg  >  2:182724/1‑251 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     ACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCG  >  NC_000913/4618801‑4619270

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