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	98.0%	156.0 / ‑4.8	51	coding (798/1323 nt)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/30);  minor base . (0/1);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

AACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAC  >  NC_000913/4618800‑4619268                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    aaCTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGAt                                                                                                                                                                                                                            <  1:178787/251‑1 (MQ=255)           aGGTCTGGACGCGCTGGTGATGGACGTGGAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGGAGCGTTTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTCCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGGTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGAA‑TCGTAACCAGCGTCTGTTTTATGTcccgttgg                                                                                                                                                                                                                  >  2:420623/1‑244 (MQ=255)                     gcgcTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGt                                                                                                                                                                                                        <  1:129927/251‑1 (MQ=255)                      cgcTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTg                                                                                                                                                                                                       >  1:82640/1‑251 (MQ=255)                       gcTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTgg                                                                                                                                                                                                      <  2:216438/251‑1 (MQ=255)                         tgggtaTGGACGGGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGGGCGTgg                                                                                                                                                                                                      <  2:425267/244‑1 (MQ=255)                             gATGGACGGGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCGCTGAAGCCCTGGCCGAAGCGATTGGTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATg                                                                                                                                                                                                 <  2:147039/250‑1 (MQ=255)                                ggACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCtt                                                                                                                                                                                              <  1:216350/250‑2 (MQ=255)                                       aaaGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACCCTCTTAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCTCCACCGCGCTGCTCACCGACATTAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCgg                                                                                                                                                                                      <  1:468943/251‑1 (MQ=255)                                                aGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCGCCGGCAAACTGGc                                                                                                                                                                             <  1:85360/251‑1 (MQ=255)                                                        gTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATg                                                                                                                                                                     >  1:148549/1‑251 (MQ=255)                                                          ttATGCCGACCTACGATCCCTCTGAAGCCCTCGCGGAAGCGACTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTTCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACTCGGTTGAAGTTCGTGAAGCGGTGCATTTCCTGACGGGTGAGTA‑TCGTATCCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGa                                                                                                                                                                    <  1:266305/250‑1 (MQ=255)                                                             tGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCAACACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGcc                                                                                                                                                                >  1:53276/1‑251 (MQ=255)                                                                 gACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAg                                                                                                                                                            <  2:290559/251‑1 (MQ=255)                                                                      aCGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAgcgcgc                                                                                                                                                       <  2:185996/251‑1 (MQ=255)                                                                               ctGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAt                                                                                                                                                <  1:451813/249‑1 (MQ=255)                                                                                 gAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCTCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAtt                                                                                                                                               <  1:188586/248‑1 (MQ=255)                                                                                    gCCCTTGCCGCAACGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGGACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGATTATCTCCGGCAAACTGGCGAAAAATGACGCCGAAGCGCGCGCGACATTGCaagg                                                                                                                                         >  2:199447/1‑248 (MQ=255)                                                                                           ccGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCCCCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAGGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAGTGCAAGCAGTGCTg                                                                                                                                  <  2:68503/251‑1 (MQ=255)                                                                                              aaGCGATTGTTGGCGTGGCTAACGGCGCGGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGAc                                                                                                                               <  2:340312/251‑1 (MQ=255)                                                                                              aaGCGATTGGTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGAc                                                                                                                               <  1:1193/251‑1 (MQ=255)                                                                                                  gATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGAAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCGGGACCACg                                                                                                                           >  2:88953/1‑251 (MQ=255)                                                                                                        tgGCGTGGCTAACGGCGCTGGCGTGCGCAACACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGGGGAGATGCATATCTCAGGCAAACTGGAGGAAGATGAAGCCGAAGCGGGCGCGAAATTGCaagcaggtcttgcacacggtag                                                                                                                       >  2:14512/1‑228 (MQ=255)                                                                                                                    cGGCGCTGGCGTGCGCGCCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCt                                                                                                         <  1:423003/251‑1 (MQ=255)                                                                                                                            gCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTACGTGAAGAGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCATGCAGTGCTGGACAACGGTAAGGCGGAAGAAGTCTTTGGTCGt                                                                                                 >  1:86449/1‑251 (MQ=255)                                                                                                                                gcgcACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATgg                                                                                             <  1:88953/251‑1 (MQ=255)                                                                                                                                 cgcACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTTCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGt                                                                                            <  1:113511/251‑1 (MQ=255)                                                                                                                                  gcACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGAGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAAATGGCGAAAGATGACGCCGACGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCATTGGTCGTATgg                                                                                             >  2:495334/1‑249 (MQ=255)                                                                                                                                            gctgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGCAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAACGGCACaa                                                                                  >  1:366105/1‑250 (MQ=255)                                                                                                                                            gctgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACaaa                                                                                 <  1:187896/251‑1 (MQ=255)                                                                                                                                            gctgctCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGAATGGTAGCGGCACaa                                                                                  >  2:242599/1‑250 (MQ=255)                                                                                                                                                     cGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGccca                                                                         >  1:214338/1‑249 (MQ=255)                                                                                                                                                      gACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGcccaac                                                                       >  1:372800/1‑248 (MQ=255)                                                                                                                                                       acaggAATCAGGTTCTGGCCCTCAGTGCGGGTAGCGTGGTGTAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGGA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGTGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTAGGGTAGCGGCACAAAAAGGCccaacc                                                                      <  1:360786/247‑1 (MQ=255)                                                                                                                                                       aaaTGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGGCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCccaacc                                                                      <  2:148549/249‑1 (MQ=255)                                                                                                                                                        caggAATCAGGTACAGGCCTCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCATTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCGTTGGTCGTATGGTAGCGGCACAAAAAGGCccaaccg                                                                     <  2:367642/248‑1 (MQ=255)                                                                                                                                                        cATGAATCGGGTACGGGCCCCCAGTGCAGGTACCGCGGTTGAAGTTCGTGAATCGGTGCGGTTCCTGAGGGGTGA‑TACTCTTAGCCCGCGGCTGTTTGATGTCACGATGGCGCTGTGCGTGGGGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGTAGCGCGCGCGAAATGGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCACAAAAAGGCccaaccg                                                                     <  2:63378/251‑1 (MQ=255)                                                                                                                                                         aTGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGACGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGGAGCGGCACACAAAGGCccaaccg                                                                     >  2:483374/1‑250 (MQ=255)                                                                                                                                                                    aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGAGAGTA‑TCGTAACACGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGACGCGCGCGAGACATTGCCAGCAGTGCTGGCCAACGGTAAAGCGGAAGAAGTCTATGGTCGCATGGTAACGGCACAAAAAGGCCCAACCGACTTCGTTgcg                                                          >  2:19236/1‑248 (MQ=255)                                                                                                                                                                       ggCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTACCCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACt                                                      <  1:483374/251‑1 (MQ=255)                                                                                                                                                                        gCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTGGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTa                                                     <  1:78930/251‑1 (MQ=255)                                                                                                                                                                           tCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTAcgc                                                  >  1:81775/1‑251 (MQ=255)                                                                                                                                                                                  cAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTAt                                           <  1:38844/251‑1 (MQ=255)                                                                                                                                                                                  cAGGAAACGCGGTTGAAGTTCGTGAAGCGGGGCAGTTCCTGAAGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCAGTGAGTGGAGATGCTGAACTCCGGCACACTGGCGAAAGAAGACGCCGAAGCGCGAGCGAAATTGCAAGAAGTGCTGGACAAAGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCAAAAAAAGGACCAACCGTATTCGATGAGGACTACGCGaggt                                             >  2:353377/1‑246 (MQ=255)                                                                                                                                                                                   cgggaaCGCGGTTGAGGTGCGTGAAGCGGTGCAGTTCCTGAGGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTGCGTTGAGAACTACGCGAAGTATc                                          <  2:275808/247‑1 (MQ=255)                                                                                                                                                                                                 aaGTTCGTGACGCGGTGCAGTTCCGGACGGGTGAGTC‑TCGAAACCCGCGTCTGTTTGATGTCTCGATGGCGCTGTGCGTGGAGATGCTTATCGCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACATCGATg                            <  2:406497/251‑1 (MQ=255)                                                                                                                                                                                                   gTTCGTGAGGCGGTGCAGTTCCGGCCGGGGGAGTA‑TCGTAACCCGCGGCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGAGGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCt                          <  2:134121/251‑1 (MQ=255)                                                                                                                                                                                                   gTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCCCAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                      >  2:35095/1‑239 (MQ=255)                                                                                                                                                                                                   gTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                      <  1:35095/239‑1 (MQ=255)                                                                                                                                                                                                                  cAGTTCCTGACGGGAGAGTA‑TCGTAACCCGCGTATGTTTGATGACACGATGGAGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCACAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCt            >  1:167956/1‑250 (MQ=255)                                                                                                                                                                                                                           aCGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGAGGCACAAAAAGGCCCAACCGACTTAGTTGGAAACTACGCAAAGTCTCTGCAGACAGAGATGCTGACGAAAGCAGGTTATGCTGATAc  >  2:274465/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    AACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAC  >  NC_000913/4618800‑4619268

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