Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,134	G→A	A302A (GCG→GCA)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,134	0	G	A	97.5%	249.2 / ‑4.6	82	A302A (GCG→GCA)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (35/44);  total (37/45)
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.98e-01

TGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑GGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGG  >  NC_000913/4618888‑4619373                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                tGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGAAGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGc                                                                                                                                                                                                                                             >  1:229651/1‑251 (MQ=255)    cGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCAAGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCCGGTAACGCGGTTGAAGTTCGTGACGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGCTGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTgg                                                                                                                                                                                                                                          >  1:411657/1‑251 (MQ=255)     gAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGACGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTgg                                                                                                                                                                                                                                          >  1:83665/1‑250 (MQ=255)     gAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTgg                                                                                                                                                                                                                                          >  1:46069/1‑250 (MQ=255)         agATTGTTGGCGTGGCTAACGGCGCTGGCGTGAGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGACAGATGACGCAGAAGCGCGCGCGCAATTGCA‑AGCATGGCTGGACAAc                                                                                                                                                                                                                                     >  2:332586/2‑251 (MQ=255)          gATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACg                                                                                                                                                                                                                                    <  1:347686/251‑1 (MQ=255)             tgttgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGATGCTCCCCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGAGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGt                                                                                                                                                                                                                                  >  2:553892/1‑250 (MQ=255)                tgGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTaaa                                                                                                                                                                                                                               <  1:507785/250‑1 (MQ=255)                  gCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCg                                                                                                                                                                                                                            <  2:167158/251‑1 (MQ=255)                   cGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCgg                                                                                                                                                                                                                           <  1:553892/251‑1 (MQ=255)                          aaCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAg                                                                                                                                                                                                                     <  1:548592/250‑1 (MQ=255)                          aaCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGt                                                                                                                                                                                                                    <  1:356066/251‑1 (MQ=255)                           aCGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGCAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGaagaaggc                                                                                                                                                                                                                   >  2:222720/1‑249 (MQ=255)                            cGGCGCTGGCGTGCGCACCAACGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGACGCGGAGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGATGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAACGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCCGTGCTGGACAACGGGAAAGCGGCAGACGTCt                                                                                                                                                                                                                  >  2:542879/1‑251 (MQ=255)                                gcTGGCGTGCGAAACACCGAGCTGCTCACAGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGATGAAGTTCGTGAAGCGGTGCCGTTCCTGGCGGGTGCGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGGGCGTGGAGATGCTTATCTACGGCAAACTGGCGCAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGAAGAAGTATTTgg                                                                                                                                                                                                              >  2:558002/1‑251 (MQ=255)                                 cTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGCTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTTGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTCGGt                                                                                                                                                                                                             >  1:108740/1‑251 (MQ=255)                                                  gcgcTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGcaca                                                                                                                                                                                            <  2:241846/251‑1 (MQ=255)                                                           aCCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGATGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGAAGCGGGACAAAAAAGccc                                                                                                                                                                                   >  1:59242/1‑251 (MQ=255)                                                                  tGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGGTCGTGAAGCGGTGCAGTTCATGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGATGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGCCGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTAGGGGAGCGGCACaaaacagccccacggac                                                                                                                                                                            >  2:195975/1‑238 (MQ=255)                                                                         ggTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGGTCGTGAAGCGGTGCAGTTCCTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTg                                                                                                                                                                     <  2:229651/251‑1 (MQ=255)                                                                            ccTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCGGACGGGTGAGTATCGTAACCCGCGTCTGGTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTgaga                                                                                                                                                                  <  2:352220/250‑1 (MQ=255)                                                                            aCTGGCCTCCAGTGCAGGTAACGCCGTTGAAGTTCGTGACGCGGTCCAGTTCATGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCGCGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGCAAGCTGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGAATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTg                                                                                                                                                                     >  2:618561/1‑248 (MQ=255)                                                                              tGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGCTGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAAc                                                                                                                                                                >  2:42293/1‑251 (MQ=255)                                                                                                   cgGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCTTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACa                                                                                                                                           <  1:431748/251‑1 (MQ=255)                                                                                                    ggTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGAc                                                                                                                                            >  1:450275/1‑249 (MQ=255)                                                                                                    ggTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGAc                                                                                                                                            >  1:541385/1‑249 (MQ=255)                                                                                                        gAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGAt                                                                                                                                      <  1:108980/251‑1 (MQ=255)                                                                                                             tCGTGAAGCGGTGCAGTTCTTGACGTGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGGGTGGAGATGCGTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAAGCGATTTCGTTGGGAACTACGCGAAGTATCTGCCCACAGCGATGCTg                                                                                                                                  <  2:429027/250‑1 (MQ=255)                                                                                                              cGTGAAGCGGTGCAGTTCCTGACGGTTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGTTTTCGTTGAGAACTACGCGAAGTAGCTGCCGACAGCGATGCTGAc                                                                                                                                <  1:503959/251‑1 (MQ=255)                                                                                                              cGTGAAGCGGTGCAGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGCAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGAc                                                                                                                                >  1:257609/1‑251 (MQ=255)                                                                                                              agTGAAGCGGTGCAGTTCCTGAAGGGTGAGTCTCGTAACACGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAAAGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGAACATAAAGGCCCAACCGACTTCGTTGAGAACTCAGCGAAGTATAGACCGACAGCGATGCTGa                                                                                                                                 >  2:416431/2‑250 (MQ=255)                                                                                                                           aGTTCCTGACGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAt                                                                                                                   >  1:351969/1‑251 (MQ=255)                                                                                                                               ccTGACGGGGGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAATATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGTCTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTg                                                                                                               <  2:492749/251‑1 (MQ=255)                                                                                                                                   gcGGGTGAGTGTCGTAACCCGCGGCTGTGTGATGTCACGATGGCGCTGTGCGGGGAGATGCGTAGCTCCGGTAAAGTGGCGAAAGATGACGGCGAAGCGCGCGCGAAATGGCAAGGAGG‑GTTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTTGTATGGTAGCTGCACAAAAAGGCCCAAGCGTCTTCGTTGAGAACTGCGCGAAGTATCGGCCGACAGTGAGGCTGTCGAAAGCACTCTTTGCTGATAc                                                                                                           <  2:568568/250‑1 (MQ=255)                                                                                                                                    cGGGTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATAcc                                                                                                          <  2:83665/251‑1 (MQ=255)                                                                                                                                     gggTGAGTAGCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGGCTTTGGTCGTATGGTTGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATGTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCg                                                                                                         <  1:266720/251‑1 (MQ=255)                                                                                                                                      ggTGAGTATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGa                                                                                                        <  1:432227/251‑1 (MQ=255)                                                                                                                                         gattatCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAg                                                                                                      <  2:254860/247‑1 (MQ=255)                                                                                                                                                gTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTc                                                                                              >  1:482498/1‑251 (MQ=255)                                                                                                                                                     ccGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTg                                                                                          >  2:616974/1‑250 (MQ=255)                                                                                                                                                      cgcgTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGTCGAAGCTCGCGCGAAATTGCA‑TGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCATCCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATTCTGATACCGAAGGTTTTGTTAGTGaa                                                                                        <  2:123000/251‑1 (MQ=255)                                                                                                                                                       gcgTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGaaa                                                                                       <  1:190327/251‑1 (MQ=255)                                                                                                                                                       gcgTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGaaa                                                                                       <  2:153408/251‑1 (MQ=255)                                                                                                                                                           cggTTTGATGTCTCGATGGCGCGGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAGTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGTTGATACCGAAGGTTTTGTCAGTGAAATgg                                                                                    <  1:209565/248‑1 (MQ=255)                                                                                                                                                             gTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAt                                                                                  >  1:600478/1‑250 (MQ=255)                                                                                                                                                             gTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTAGGCTGATACCGAAGGTTTTGTCAGGGAAATgg                                                                                    <  1:162722/248‑1 (MQ=255)                                                                                                                                                                tGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGTCAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAccc                                                                              >  1:532582/1‑251 (MQ=255)                                                                                                                                                                 gATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAccc                                                                              <  1:540913/250‑1 (MQ=255)                                                                                                                                                                  aTGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTCAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcgc                                                                            >  2:39987/1‑251 (MQ=255)                                                                                                                                                                   tGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcgcg                                                                           >  2:503947/1‑251 (MQ=255)                                                                                                                                                                      cACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCcgcgcg                                                                         <  1:542879/250‑1 (MQ=255)                                                                                                                                                                       aCGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCt                                                                       <  1:488618/251‑1 (MQ=255)                                                                                                                                                                         gagtgcgcTGTGCGTGGAGATGCTTATCTCCGGCAACCTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGGCGTATGTTAGCGGCACAAAAAGGCCCAACCGTCTTCGTTGAGAACGACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCTGTCTATGCTGATACCTAAGGTTTTGTCAGTGAAATGGATACCCCCGCGCTgg                                                                     <  2:368932/247‑1 (MQ=255)                                                                                                                                                                            ggCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTgggg                                                                   <  1:77422/250‑1 (MQ=255)                                                                                                                                                                              cgcTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGGCGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGTTAAAGCGGCAGAAGTCTTTGGGCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATg                                                                <  1:359188/251‑1 (MQ=255)                                                                                                                                                                                 tgtgCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTg                                                                                                                                  >  2:295829/1‑182 (MQ=255)                                                                                                                                                                                 tgtgCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTg                                                                                                                                  <  1:295829/182‑1 (MQ=255)                                                                                                                                                                                   ggCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGt                                                           <  2:482498/250‑1 (MQ=255)                                                                                                                                                                                      gTGGAGATGCTTATCTACGGCAAACTGGCGAGAGATGACGGCGAAGCGCGCGAGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGACTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGACGACAGCGATGCTGCCGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATAACCGCGCGCTGGGGATGGCAGTGGt                                                        >  2:400796/1‑251 (MQ=255)                                                                                                                                                                                        ggAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGTGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAATGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg                                                      <  1:630110/251‑1 (MQ=255)                                                                                                                                                                                        ggAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg                                                      >  2:9074/1‑251 (MQ=255)                                                                                                                                                                                        ggAGATGCTTAGCTCCGGCAAACTGGCGAAAGAGGACGCCGAAGCGTGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACGACGCGAAGGATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTGGTCATTGAAATGGATACCCGCTCGCTGGGGATGGCAGTGGTTg                                                      <  2:434436/251‑1 (MQ=255)                                                                                                                                                                                                ttaTCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATgggg                                              >  2:107911/3‑250 (MQ=255)                                                                                                                                                                                                 taTCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGggcg                                             >  1:544252/2‑251 (MQ=255)                                                                                                                                                                                                  aTCTCGGGCAAACTGGCGAAAGCTGACGCCGAAGCGCGCGAGCACTTGCA‑AGCGGAGCTGGACAAAGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGATGATACCGAAGGTTTTGTCAGGGAAAAGGAAACACGCGCGCTGGGGATGGCAGTGGTTGCAGTGggctg                                            >  2:378391/1‑249 (MQ=255)                                                                                                                                                                                                  aTCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGCAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGATGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGggcgg                                            >  2:412119/1‑251 (MQ=255)                                                                                                                                                                                                      ccGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGggcggc                                           <  1:339168/248‑1 (MQ=255)                                                                                                                                                                                                      ccGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATGGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGGCGTATGGTAGCGGCGCAAAAAGGCCCAACCGACTTCGTTGAGAACTATGCGAAGTATCGGCCGGCAGCGATGCTGGCGAAAGCAGTCTATGCTGAGACCGAAGGTGTTGTGGGGGAAATGGAGGCTCGCGCTCGGTGGAtggcagtggtggcagtggttggcggt                                        <  2:411657/251‑9 (MQ=255)                                                                                                                                                                                                          cAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGcggcggtcgtc                                    <  1:618924/251‑6 (MQ=255)                                                                                                                                                                                                                ggCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCa                                <  1:537161/249‑1 (MQ=255)                                                                                                                                                                                                                 gCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGc                             <  1:321486/251‑1 (MQ=255)                                                                                                                                                                                                                   gAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGCAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCACTGGGCGGCGGTCGTCGTCAGGCAt                           >  1:215275/1‑251 (MQ=255)                                                                                                                                                                                                                        aTGACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGac                      <  1:26845/251‑1 (MQ=255)                                                                                                                                                                                                                          gACGCCGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGaca                     <  1:2064/250‑1 (MQ=255)                                                                                                                                                                                                                              ccGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGGCGTATGGTAGCGGCACAAAAAGGCCCAACCGACGTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATGCCCGCGCGCGGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCAtt                <  2:319616/251‑2 (MQ=255)                                                                                                                                                                                                                               cGAAGCGCGCGCGAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCAttg               <  2:515860/251‑3 (MQ=255)                                                                                                                                                                                                                                    cgcgcgcgAAATTGCA‑AGCAGTGCTGGACAACGTTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCGGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTAc          <  2:3506/251‑1 (MQ=255)                                                                                                                                                                                                                                      cgcgcgGAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGAAGACGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGATGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTAc          >  2:581872/1‑249 (MQ=255)                                                                                                                                                                                                                                       gcgcgAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGc       <  2:133835/251‑1 (MQ=255)                                                                                                                                                                                                                                           gAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGCGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACAGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCg   >  2:410557/1‑251 (MQ=255)                                                                                                                                                                                                                                           gAAATTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGAATTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCg   >  1:583067/1‑251 (MQ=255)                                                                                                                                                                                                                                            aaaTTGCA‑AGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCgg  >  1:428612/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                TGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACCGACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATATCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCA‑GGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGG  >  NC_000913/4618888‑4619373

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: ‑