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	95.5%	117.5 / ‑2.2	45	E266E (GAA→GAG)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (19/23);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

CTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCT  >  NC_000913/4618783‑4619263                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ctctATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGCCATGAATCAGGTACTGGCCTCCAGTGCCGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGt                                                                                                                                                                                                                                         >  1:35575/1‑250 (MQ=255)              aagaaACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGtt                                                                                                                                                                                                                           >  1:337035/1‑251 (MQ=255)               agaaACTGGCGGAAGGTCTGGTCGCGCTGGTGATGGACGTGAAAGTGGGTATCGGCGCGTTTATGCCGACCGACGAACTTTCTGAAGCCTGTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTTGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCTTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGttt                                                                                                                                                                                                                          <  1:353364/251‑1 (MQ=255)                   aCTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGTCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAATATCGTAACCCGCGTCTGTTTGATg                                                                                                                                                                                                                      >  1:196819/1‑251 (MQ=255)                        cGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACg                                                                                                                                                                                                                 <  1:223886/251‑1 (MQ=255)                         ggAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑CCGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACg                                                                                                                                                                                                                 >  1:302222/1‑250 (MQ=255)                                   gacctggttgTGAGGGAAGTGTAAGTTGGTGGCGGTGCGTTTATGCCGACCTACGAACTCTTTGTAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTATCGTCGCTGTATTGGGCACCACCGCTCTGCTAATCGATATGAATCAGGTACTGGTTTTCAGTGCAGGTAACGCGGTTGAAGTGTGTTTAGCGGTGCA‑TTCCTTTGCTGGGGAGTATCGTATCCTGCGTCTGTTTGATGTCACGATGGCGCtgtg                                                                                                                                                                                                      <  2:232834/242‑1 (MQ=255)                                            gTGATGGACGTGAAAGTGGGTGGCGGCGCGTTTATGCCGTCCTATGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCGGTGCGTGGAGAt                                                                                                                                                                                             <  1:82203/251‑1 (MQ=255)                                            gTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGAt                                                                                                                                                                                             <  1:231589/251‑1 (MQ=255)                                              gATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGc                                                                                                                                                                                           <  1:113301/251‑1 (MQ=255)                                                  gACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCCGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCCCGATGGCGCTGTGCGGGGAGATggttat                                                                                                                                                                                       >  2:75343/1‑246 (MQ=255)                                                      tGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTcc                                                                                                                                                                                   >  1:158906/1‑251 (MQ=255)                                                       gAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCg                                                                                                                                                                                  >  1:142576/1‑251 (MQ=255)                                                              ggTAGCGGCGCGTTTATGCCGCCCTACGAACTCTCTGAAACCCTTGCCGACGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGAGCACACCCGCGCTGCTCACCGACATGAATCAGGCACTAGCCTCCAATGCTGGTAACGCGGTTGAAGGTCCTGCAGCTGTGCAGTTCCTG‑ACGGGTGAGAATCGTAACCCGCGTCTGATTAAAGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAAcc                                                                                                                                                                           >  1:62790/1‑250 (MQ=255)                                                              ggTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACt                                                                                                                                                                           >  1:4183/1‑251 (MQ=255)                                                                 tgctgctcgTTGATGCCGACCTACGAACACTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCGGGCGTGCGCGCCCCCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTgg                                                                                                                                                                         <  2:62722/243‑1 (MQ=255)                                                                    ggcgCGTTTATGCCAACCTACGAACTCGCTGAAGCCCTTGCCGAAGCGATGGTTGGCGTGGCTAACGGCGCTGGCGGGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGaa                                                                                                                                                                     <  2:197078/251‑1 (MQ=255)                                                                             aTGCCGACCTACGAACTCTCTGAAGTCCTGGCCGACGCGATTGTTGGCGTGGCTAACGGCGCTGGCGGGCGCACCACCGCTCTTCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAATGTGGGTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACg                                                                                                                                                             <  1:152774/250‑1 (MQ=255)                                                                                  gACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAg                                                                                                                                                       <  1:62176/251‑1 (MQ=255)                                                                                   aCCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACAGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGATCGTGAAGCGGTACAGTTCCTG‑ACGGGAGAGTAGCGTAACCGGCGTCTGTTTGATGTCACGATGGGGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGCAgc                                                                                                                                                      >  1:215494/1‑251 (MQ=255)                                                                                          aaCTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGcgcgcgcg                                                                                                                                               >  1:265970/1‑251 (MQ=255)                                                                                                   acagcctggCCGAAGCGATTGTTGGCTTGGCTAACGGCGCGGGCGTGCGCCCCACCGCGCTGCTCAACGACATGAATCAGGTACTGGCCGACAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCaa                                                                                                                                      <  2:116673/243‑2 (MQ=255)                                                                                                     gCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCaagc                                                                                                                                    >  1:303845/1‑248 (MQ=255)                                                                                                      cccTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCaagc                                                                                                                                    >  1:209973/1‑247 (MQ=255)                                                                                                       ccTTGCCAAAGAGATTGTTGGAGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGGACGGGCCGCCAGGGCAGGGAACGCGGTTGAAGTTCGTGAAGGGGTGAAGGTCCTG‑ACGGGAGAGTATCGTAACCCGCGTCTGTTTGCTGGCACGATGGCGCTGGGCGTGGAGATGCTTATGTCCGGGAAACTGGCGAAAGATGAAGCCGAAGAGCGTGCGAAATTGCaagg                                                                                                                                    >  2:382663/1‑246 (MQ=255)                                                                                                             cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGGTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTAACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAAATGACGCCGAAGCGCGCGCGGAAtttcaagcagagctgg                                                                                                                            >  2:355008/1‑237 (MQ=255)                                                                                                               aaGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGAc                                                                                                                          <  1:81720/251‑1 (MQ=255)                                                                                                                aGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCg                                                                                                                                                                                  <  1:232084/194‑1 (MQ=255)                                                                                                                  cGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTTCAAGCAGTGCTGGACAAc                                                                                                                       <  1:102146/251‑1 (MQ=255)                                                                                                                               ggCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCa                                                                                                           <  1:225853/250‑1 (MQ=255)                                                                                                                               ggCTAACGGCGCGGTCGTGCGCACCACCGCTCTGCTCACAGACATGAATAAGGTACGGGCCACCAGTGCAGGTAACGCGGTTGAATTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTTAGTAGCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGATGGACACCGGTAAAGCGGCAg                                                                                                          <  2:26565/251‑1 (MQ=255)                                                                                                                                   aacggggctggctggctccccaccGTGTTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTTCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCTGTAGAAGt                                                                                                      <  2:265970/232‑1 (MQ=255)                                                                                                                                      ggCGCTGGCGGTCGCACCGCCGCGCTGCTCACCGACAGTAATCAGTGACGGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGTAGCGGGGCAGTTCCTG‑ACGGGTGAGGATCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGGCAACTGTAAAGCGGCTGAAGTCtt                                                                                                   <  2:131375/251‑1 (MQ=255)                                                                                                                                             gCGTGCGCACCGCCGCGCTGCCTACCGACAGGAATCACGTCGTAGCCCCCAGTGCGGGTAACGCGGTTGAAGTTCTTGAAGCGGTGCATGTCCTG‑ACTGGTTATTTTCGCAACCCGCGTCTGTTTTATGGCACGATGGTGCTTTGCGTGGAGATGGTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAATCGCTCGCGAAATTGCAAGGAGTGCTTGACAACGGTAAAGCGGccgttttccttggtcgt                                                                                            <  2:195598/251‑17 (MQ=255)                                                                                                                                                 gcgcACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATg                                                                                         <  1:75343/250‑1 (MQ=255)                                                                                                                                                 gcgcACCACCGCGCTGCTCACCGAAATTAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATgg                                                                                        <  1:8588/251‑1 (MQ=255)                                                                                                                                                   gcACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGt                                                                                       >  1:154844/1‑250 (MQ=255)                                                                                                                                                                           agAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGACGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCcaaccgac                                                              <  1:323068/250‑1 (MQ=255)                                                                                                                                                                                              cAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTAGACGTATGGTAGCGGCACAAAAAGGCCCAACAGACTTCGTTGAGAACTAcgcg                                            >  1:10913/1‑250 (MQ=255)                                                                                                                                                                                                     ggTAACGCGTTTGAAGTTCGTGAATCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCt                                    <  1:305478/251‑1 (MQ=255)                                                                                                                                                                                                      gTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTg                                   >  1:244688/1‑251 (MQ=255)                                                                                                                                                                                                                            aGCGGTGCAGTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAg             <  1:294395/251‑1 (MQ=255)                                                                                                                                                                                                                                  gCAGTTCCTG‑GCGGGTGTGTATCGTAACCCGCGTCTGTTTTATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGCCAACGGTAAAGCGGCAGAAGTCTTTGTTCGTATGGTAGCGTCACAAAAATGCCCAACCGTCTTCCTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCt       <  1:190678/251‑1 (MQ=255)                                                                                                                                                                                                                                     gTTCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATg    >  1:227065/1‑251 (MQ=255)                                                                                                                                                                                                                                       tCCTG‑ACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCt  >  1:81265/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              CTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTG‑ACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCT  >  NC_000913/4618783‑4619263

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