Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,610	G→A	L336L (CTG→CTA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,610	0	G	A	94.6%	293.3 / 4.2	93	L336L (CTG→CTA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (35/53);  total (37/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

GGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGT  >  NC_000913/4620364‑4620859                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          ggtTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTAcgcg                                                                                                                                                                                                                                                       <  1:616848/251‑1 (MQ=255)           ggTAAAATTGCGGACATCTACGCCATCTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACAt                                                                                                                                                                                                                                             <  2:490134/251‑1 (MQ=255)              aaaaTTGCGGTCATCTCCGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGATGTTTGACGCCACCATCAAAGAGAGGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGTGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCt                                                                                                                                                                                                                                          <  2:163055/251‑1 (MQ=255)                aaTTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGGCGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGGCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGa                                                                                                                                                                                                                                        <  1:274823/251‑1 (MQ=255)                 aTTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGCTGTCTCTGCTACGCGGTTACGACATCCTgtt                                                                                                                                                                                                                                       >  2:398065/1‑249 (MQ=255)                     cGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGTCCTGGACGCGCTGTTTGTCGCCACCATCAAAGAGATGAAATAAGCGGGTGATAACACCATCGTCTTCACCACCTTCTTTGACTTCGACTCTTCCTGGGGCCCCCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGCCATCCTGATCCTc                                                                                                                                                                                                                                   <  2:132180/251‑1 (MQ=255)                          aTCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCCCCGc                                                                                                                                                                                                                              >  2:293342/1‑251 (MQ=255)                                cccAACTGCGGTATCCCCAAAAAAGGGAAAGCGACGGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCGTCACCAACTTCTTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCGGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCa                                                                                                                                                                                                                        <  2:412398/250‑1 (MQ=255)                                  cAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACg                                                                                                                                                                                                                      >  1:603480/1‑251 (MQ=255)                                         ggTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCg                                                                                                                                                                                                                >  1:739797/1‑250 (MQ=255)                                            aTCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATc                                                                                                                                                                                                             <  1:602356/250‑1 (MQ=255)                                            aTCACCAAAAAAGGGAAAGCGACGGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGCTAACCCCATCGTCTTCACCAACTTCGTTGACTTCGCCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCGCTGCTACGCGATGACGACATCCTGATCCTCCCCGCTGACCACGGTTGCGATcc                                                                                                                                                                                                            <  2:170034/251‑1 (MQ=255)                                                 cAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCt                                                                                                                                                                                                       <  1:128626/251‑1 (MQ=255)                                                    aaaaGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGCTCCTCACCGCTGCCCACGGTTGCGATCCGACCTggc                                                                                                                                                                                                    >  2:178086/1‑250 (MQ=255)                                                     aaaGTGAAAGCGACTGGGCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTGCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGAc                                                                                                                                                                                                   <  1:350987/251‑1 (MQ=255)                                                     aaaGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTAACCGCTGACCACGGTTGCGATCCGACCGGGAc                                                                                                                                                                                                   >  2:607079/1‑251 (MQ=255)                                                      aaGTGAAAGCGACTGGCCGGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGAcc                                                                                                                                                                                                  <  2:83441/251‑1 (MQ=255)                                                       aGTGAAAGCGACTGGCCTGGTCGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCg                                                                                                                                                                                                 <  1:358210/251‑1 (MQ=255)                                                          gAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGt                                                                                                                                                                                               >  2:485621/1‑250 (MQ=255)                                                                gACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCCCCAACTTCGTTGACTTCTACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGCCGACATCCTGATCCTCACCGCTGCCCACGGTTGCGTTCCGACCTGGCCCGGTACTGAcc                                                                                                                                                                                        >  2:365029/1‑251 (MQ=255)                                                                 aCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACca                                                                                                                                                                                       <  1:586635/251‑1 (MQ=255)                                                                    ggCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCacac                                                                                                                                                                                    >  1:559670/1‑251 (MQ=255)                                                                     gcgtGGACGCGCTGTTGGACGCCACCAGCAAAGAGATGAAAGAAGCGGGTGATACCACCCTCGTCTCCACCAACTTCGTTGACTTCGTCTCTTCCTGGGGCCCCCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCCCTGCTACGCGATGACGACGTCCGGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACGGACCACAcg                                                                                                                                                                                   <  1:342868/248‑1 (MQ=255)                                                                         ggACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAAATGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTg                                                                                                                                                                               >  1:458238/1‑251 (MQ=255)                                                                          gACGCGCTGTTTGACGCCACCAGCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTGCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGa                                                                                                                                                                              <  2:164263/251‑1 (MQ=255)                                                                            ggcgcTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTGCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAac                                                                                                                                                                            <  2:356686/250‑1 (MQ=255)                                                                                  gtgtGACGCCACCATCAAAGAGAGTAAAGAAGCGGGTGATAACACCATCGTCTGCACCAACTGCGTTGACTTCGACTCTTCCGGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCGGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTc                                                                                                                                                                      <  2:482643/248‑1 (MQ=255)                                                                                  gTTTGACGCCCCCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCCCCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACAtt                                                                                                                                                                       <  2:211727/250‑1 (MQ=255)                                                                                     tGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCg                                                                                                                                                                    <  2:559670/250‑1 (MQ=255)                                                                                      gACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGCCCTGGACCGGTCCTGACCACACGCGTGAACACATTCCGGt                                                                                                                                                                  >  2:575383/1‑251 (MQ=255)                                                                                      gACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCgg                                                                                                                                                                   <  1:365949/250‑1 (MQ=255)                                                                                          ccaccaTCAAAGAGATGAAAGAATCGGGTGATAACCCCATCTTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGGCGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGCCCACACGCGTGAACCCATTCCGGTACTg                                                                                                                                                              <  2:520791/251‑1 (MQ=255)                                                                                          ccaccaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTAGGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTg                                                                                                                                                              <  1:414486/251‑1 (MQ=255)                                                                                           caccaTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTgg                                                                                                                                                             <  1:269161/251‑1 (MQ=255)                                                                                               aTCAAAGAGATGAAAGAAGCGGGTGATAACACCGTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCAACGTCGCCGGTTAAGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtat                                                                                                                                                          >  1:340292/1‑250 (MQ=255)                                                                                               aTCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTTGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGTAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtata                                                                                                                                                         <  1:575383/251‑1 (MQ=255)                                                                                                acAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTGCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACGGACCACACGCGTGAACACATTCCGGTACTGGTatat                                                                                                                                                        <  2:458238/250‑1 (MQ=255)                                                                                                 cAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATg                                                                                                                                                       <  1:485621/251‑1 (MQ=255)                                                                                                  aaaGAGAGGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATgg                                                                                                                                                      <  2:579233/251‑1 (MQ=255)                                                                                                     gagaTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGccc                                                                                                                                                   >  1:635800/1‑251 (MQ=255)                                                                                                      agaggAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGTCTTCGACTCTTCCGGGGGCCCCCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCg                                                                                                                                                  <  2:739797/247‑1 (MQ=255)                                                                                                      agaTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCg                                                                                                                                                  >  1:366098/1‑251 (MQ=255)                                                                                                        aTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGaa                                                                                                                                                <  1:331149/251‑1 (MQ=255)                                                                                                        aTGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGa                                                                                                                                                 >  2:791938/1‑250 (MQ=255)                                                                                                          gAAAGAAGCGGGTGATAACACCCTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAg                                                                                                                                              <  1:65955/251‑1 (MQ=255)                                                                                                               aagCGGGTGATAACACCATCGTCTTCACCAACTTCTTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTaaaa                                                                                                                                         <  2:182441/251‑1 (MQ=255)                                                                                                                     gTGATACCCCCATCGCCTTCCCCAACTTCGTTGACTTCGACTCTTCCGGGGGCCACCGTCGCGCCGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGc                                                                                                                                   <  2:715155/251‑1 (MQ=255)                                                                                                                     gTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACCTTCCGGTACTGGTATATTGCCCGCAAGTAAAACCGGGc                                                                                                                                   >  2:310544/1‑251 (MQ=255)                                                                                                                       gATAACGCCATCGTCTTCACCATCTTCGTTGCCTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCt                                                                                                                                  <  1:507765/250‑1 (MQ=255)                                                                                                                              ccATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGt                                                                                                                          >  2:244298/1‑251 (MQ=255)                                                                                                                               cATCGTCTCCTCCAACTTCTTTTCTTTCGCCTCTTCCGGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGAGGACTACATCCTGATCCTCACCGCTGACCACGGTTGCGATCTGACCTGGACCGGTACTGTCCTCACGCGTGAACACATTCCGGTACTGGTATATGGCCTGAAAGTAAAACCGGGCTCACTGGGTc                                                                                                                         <  2:701199/251‑1 (MQ=255)                                                                                                                                aTCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACCTCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACATGGACCGGTACTGACCACAAGCTTGAACACATTCCGGTACTGGTATGTGTCCCGAAAGAAAAACAGGGCTCACTGGGtt                                                                                                                         >  2:571174/1‑249 (MQ=255)                                                                                                                                  cGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTCATCCTCACCGCTGCCCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACCCCTTCCGGTACTGGTCTATGGCCCGAAAGTAAAACCGGGCCCACTGGGTcctc                                                                                                                      >  2:807299/1‑248 (MQ=255)                                                                                                                                             cctctgtgtcatacgcctctTCCTGTGGCCACCGTCGCGCCGTCGCCGGTTAGGCCGCTGGTCTGGAACTGGTCGTCCGCCGTCGTCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGGCCACTCGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTc                                                                                                           <  2:664710/235‑1 (MQ=255)                                                                                                                                             ccgccgtggatttCGGCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGAGGACGACAGCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCTCTGGGTCATCGTGAAACCTTc                                                                                                           <  2:519907/240‑1 (MQ=255)                                                                                                                                             aCTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTc                                                                                                           <  1:417660/251‑1 (MQ=255)                                                                                                                                                 cGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTcgcg                                                                                                        <  1:411179/250‑1 (MQ=255)                                                                                                                                                     gACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGatat                                                                                                   <  1:295956/251‑1 (MQ=255)                                                                                                                                                     gACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACATTCGCGGatat                                                                                                   <  1:375130/251‑1 (MQ=255)                                                                                                                                                          cgtctctACCTGGGGCCTCCGTCGCGACGTCGCCGGGTAGGCCGCGGGTCTGGAACGGGTCGACCGCCGTCTGCCGGAGCTGATGGCTCCGCTACGCGATGACGACCTCCTGATCCTCACCGCGGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACAAATTCCGGTACTGGTATATGGCCCGCAAGTAAAACCGGGCTGACTGGGTCATCGTGAAACCTTCGCGGATATCgtc                                                                                               <  2:93692/247‑3 (MQ=255)                                                                                                                                                          cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGcc                                                                                              <  1:293342/251‑1 (MQ=255)                                                                                                                                                           gACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGcc                                                                                              >  1:156259/1‑250 (MQ=255)                                                                                                                                                            ccctttCCTGGGGCCACCGTCGCAACGTCGCCGGTTAAGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAg                                                                                            <  2:340292/247‑1 (MQ=255)                                                                                                                                                                ttCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGActc                                                                                        >  1:465603/1‑251 (MQ=255)                                                                                                                                                                    tGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGGGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCAGCCAGACTCTGGc                                                                                    >  1:755681/1‑251 (MQ=255)                                                                                                                                                                        gCCACCGTCGCGACGTCGCCGGTTATTCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGCCATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGGCCTGGACCGGTACTGACCACACGCGTGAACCCATTCCGGTACTGGTATATGGCCCGAACGGAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaca                                                                                 >  2:311399/1‑248 (MQ=255)                                                                                                                                                                         ccACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAat                                                                               <  2:465603/251‑1 (MQ=255)                                                                                                                                                                           cccGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAatat                                                                             <  1:388691/250‑1 (MQ=255)                                                                                                                                                                            ccGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAtt                                                                            >  2:396753/1‑251 (MQ=255)                                                                                                                                                                                     cGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGCCCGCCGTCTGCCGGAGCTGATGTCGCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACtt                                                                   <  1:58714/251‑1 (MQ=255)                                                                                                                                                                                        cGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCt                                                                 <  1:365029/250‑1 (MQ=255)                                                                                                                                                                                         gCCGGTTTTGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGa                                                               <  2:773713/251‑1 (MQ=255)                                                                                                                                                                                          ccGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGat                                                              >  1:577808/1‑251 (MQ=255)                                                                                                                                                                                              ttATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATgg                                                          >  1:666148/1‑251 (MQ=255)                                                                                                                                                                                                        ggTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCa                                                 >  1:314263/1‑250 (MQ=255)                                                                                                                                                                                                             ggAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGcc                                            <  1:791938/250‑1 (MQ=255)                                                                                                                                                                                                             ggAACTGTTCGACCGCCGGCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGGACTTCTGATATGGAATATGGCAAAGcc                                            <  1:310544/250‑1 (MQ=255)                                                                                                                                                                                                                       gACCGCCGTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATg                                 >  1:524413/1‑251 (MQ=255)                                                                                                                                                                                                                         ccgccgTCTGCAGGAGCTGATGTCTATGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGAGATCCGACCTGGCCCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCAATGGGTCATCGTGAAACCTTCGCGGGTATCGGCCAGACTCTGGCAAAATATTTTGGTACTTATGATATGGCATATAGCAAAACCATGTTCTGTTgg                                >  2:457088/1‑250 (MQ=255)                                                                                                                                                                                                                           gccgTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGa                                                               <  2:156254/217‑1 (MQ=255)                                                                                                                                                                                                                           gccgTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGa                                                               >  1:156254/1‑217 (MQ=255)                                                                                                                                                                                                                            ccgTCTGCCGGAGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACCCGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGCCTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGAttt                            >  2:325324/1‑251 (MQ=255)                                                                                                                                                                                                                             cgTCTGCCGGATCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTg                           <  2:635800/251‑1 (MQ=255)                                                                                                                                                                                                                                 tGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATACTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACCGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCg                       >  1:115591/1‑251 (MQ=255)                                                                                                                                                                                                                                    cggggCTGATGTCTCTGATGCGCGATGACGAAATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAg                    <  2:16369/247‑1 (MQ=255)                                                                                                                                                                                                                                    cGGAGCTGAAGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTAATTCTGCTATGGAATATGGCAAAGCCATGTTCTGGTGGATTTGGGCgggg                    >  2:515753/1‑249 (MQ=255)                                                                                                                                                                                                                                       aGCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGt                 >  1:55599/1‑251 (MQ=255)                                                                                                                                                                                                                                        gCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGaa                                                                                                                                                                             <  1:420494/94‑1 (MQ=255)                                                                                                                                                                                                                                        gCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGaa                                                                                                                                                                             >  2:420494/1‑94 (MQ=255)                                                                                                                                                                                                                                        gCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGtt                <  1:396753/251‑1 (MQ=255)                                                                                                                                                                                                                                        gCTGATGTCTCTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGtt                >  1:768545/1‑251 (MQ=255)                                                                                                                                                                                                                                           gATGTCTCTGCTACGCGATGACGACAGCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGAc             <  1:460388/251‑1 (MQ=255)                                                                                                                                                                                                                                               tctcTGCTACGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGAGAACACATTCCGGTACTGGAATATGGCCCGAAAGTAAAACAGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATAGGCCAGACTCTGGCAAAATATTTTGGTATTTCTGATATGGAATATGGCAAAGCAATGTACTGATGGATTTGGGCGGAGCGtttaatccg         >  2:203629/1‑244 (MQ=255)                                                                                                                                                                                                                                                      tgcGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAACCTATTTTGGTACTTCTGATCTGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGt  >  2:296083/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          GGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGT  >  NC_000913/4620364‑4620859

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