BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,348,705	0	A	G	31.7%	121.2 / 90.7	104	intergenic (‑70/+404)	PP_3820/galU	Group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand): ref base A (37/34); new base G (17/16); total (54/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 104 total aligned reads displayed.

GCCTGAGGTTGGCCGGTGCAAGCACCCGCTCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCA  >  AE015451/4348437‑4349004

tggagttcagacgtgtgctcttccgatcttCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                                                                                                                                                                               <  4:676083/272‑1 (MQ=255)

gcCTGAGGTTGGCCGGTGCAAGCACCCGCTCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                                                                                                                                                                               <  2:652450/301‑1 (MQ=11)

 ctttccctacacgacgctcttccgatcttCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCACCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                                                                                                                                                                                                                              <  1:192696/273‑1 (MQ=11)

   ggAGGTTGGCCGGTGCAAGCACCCGCTCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGc                                                                                                                                                                                                                                                                            <  3:7704/300‑1 (MQ=14)

        ttGGCCGGTGCAAGCCCCCGCTCCATCAGCTTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCg                                                                                                                                                                                                                                                                       <  1:319995/301‑1 (MQ=17)

            ccGGTGCAAGCACACGCTCCATCAGCGTGACCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                                                                                                                                                                                   <  2:331048/301‑1 (MQ=17)

                        ccgatctcATCAGCGTGTCCGGCCCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAggg                                                                                                                                                                                                                                                       <  1:436356/294‑1 (MQ=17)

                             tCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAtt                                                                                                                                                                                                                                                  <  4:349779/301‑1 (MQ=17)

                             tCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGtagatcggaagagcgtcggtgagggacagg                                                                                                                                                                                                                                                  >  2:675861/1‑272 (MQ=255)

                             tCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCACCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTgagatcggaagagcacacgtctgaactcc                                                                                                                                                                                                                                                  >  3:192694/1‑273 (MQ=11)

                               cATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgggagagcgg                                                                                                                                                                                                                                                >  3:436358/1‑294 (MQ=17)

                                  cAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                                                                                                                                                                             >  3:235559/1‑301 (MQ=32)

                                        gtCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACTCGCGGTCGCACCGTGGCAAGCGTCCAACTCCCCCAGGGTAATTCGCCCGAcact                                                                                                                                                                                                                                       >  2:381786/1‑298 (MQ=11)

                                          ccGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                                                                                                                                                                                     >  4:177421/1‑301 (MQ=32)

                                            ggCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCCGTAGCACCGTTGCAAGCGTCCCACTCCCCCAGGGCAATTCCGGCGACCTTCCtt                                                                                                                                                                                                                                   >  2:102021/1‑300 (MQ=11)

                                              cTCCGCGTTCGTCAACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCTGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTgtt                                                                                                                                                                                                                                 <  2:681820/301‑3 (MQ=17)

                                                 tgtgctcttccgatctgCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                                                                                                                                                                                              <  3:315421/285‑1 (MQ=12)

                                                         ggagttcagacgtgtgctcttccgatcttGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCg                                                                                                                                                                                                                      <  4:381547/274‑1 (MQ=11)

                                                                gttgtgctcttccgatctcTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc                                                                                                                                                                                                               <  4:185512/283‑1 (MQ=12)

                                                                gtgtgctcttccgatctccTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc                                                                                                                                                                                                               <  4:542272/284‑1 (MQ=12)

                                                                 gCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAtagatcggaagagcgtc                                                                                                                                                                                                              >  1:315420/1‑285 (MQ=12)

                                                                   tGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCCCGGTAGCCCAGTGGCAAGCGTAAAGATACCCCAGGGTAATTCGCGCGACCTTCCTTCTTATGCCTGGCGGATCAACGTc                                                                                                                                                                                                            >  2:428033/1‑301 (MQ=11)

                                                                     ctttccctacacgacgctcttccgatctcATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTcac                                                                                                                                                                                                          <  1:231583/273‑1 (MQ=255)

                                                                      cgccgcCGGCACCTTTGCACTGTCAGGCATAACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTcaca                                                                                                                                                                                                         <  3:135392/301‑1 (MQ=17)

                                                                        ccgccgGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGc                                                                                                                                                                                                       >  1:78440/1‑301 (MQ=32)

                                                                                 ccTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCtccagatcggaagagcttcg                                                                                                                                                                                              >  2:542115/1‑282 (MQ=12)

                                                                                 ccTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGGTCCGTGc                                                                                                                                                                                              >  2:545751/1‑301 (MQ=17)

                                                                                  cTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCAACCTTCCTGCTTATCCCTGTCGGCTCTACGTCCCAGCGTTCCTTGcc                                                                                                                                                                                             >  2:97244/1‑300 (MQ=255)

                                                                                     tGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAACGATCTCCCAGGGTAATTCGCGGGGCCTTCCTGCCTATGCCTGTCgaagtcgggataggcgcgtggtgggaaaa                                                                                                                                                                                          >  2:381457/1‑273 (MQ=11)

                                                                                     tGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAgg                                                                                                                                                                                          >  3:96702/1‑300 (MQ=17)

                                                                                              ggggaTCACCCTCGGCTGCTGGGCGTGATTCCGGAGGACAGTTTCCGCTGTGGGAAAGTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCtt                                                                                                                                                                                 <  2:177379/297‑1 (MQ=255)

                                                                                                 ggagttcagacgtgtgctcttccgatctgggATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTgg                                                                                                                                                                              <  4:424714/273‑1 (MQ=255)

                                                                                                 cATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAgatcggaagagcacacgtctgaactcc                                                                                                                                                                              >  3:231581/1‑275 (MQ=255)

                                                                                                    cACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCgg                                                                                                                                                                           >  1:605308/1‑300 (MQ=17)

                                                                                                          cGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgtg                                                                                                                                                                     <  4:545908/301‑1 (MQ=32)

                                                                                                           ggCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGc                                                                                                                                                                    >  3:21415/1‑301 (MQ=32)

                                                                                                           cgctcttccgatctaTTCGGGATGACAGTTTCCACTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGc                                                                                                                                                                    <  1:177483/287‑1 (MQ=12)

                                                                                                               tcgggCCGGGAATCGGGATGACAGTTTCCGCGGTGGGGAACTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTGCTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCAcc                                                                                                                                                                <  2:148707/298‑1 (MQ=255)

                                                                                                                       ggATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCg                                                                                                                                                       <  3:78442/301‑1 (MQ=32)

                                                                                                                         aTTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACGGCGTTCCGTGCAAGTATTGGGCTTTGGCGTTTTGTCCCACCTCCACCCCGc                                                                                                                                                      >  2:590002/1‑301 (MQ=12)

                                                                                                                         aTTCGGGATGACAGTTTCCACTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGcagatcggaagagca                                                                                                                                                      >  3:177485/1‑287 (MQ=12)

                                                                                                                          cccctcaacgccgctcttccgatctgaactttCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCtg                                                                                                                                                    <  2:243845/272‑1 (MQ=255)

                                                                                                                            cGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCtgtt                                                                                                                                                  >  1:259923/1‑300 (MQ=17)

                                                                                                                             gggATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCttttgagatcggaagagctgcgtttagggaaag                                                                                                                                                  >  2:424605/1‑271 (MQ=255)

                                                                                                                                 tGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTgccgcc                                                                                                                                             >  4:630888/1‑301 (MQ=32)

                                                                                                                                    cAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTTTGCCACCT‑CACCCCCCTTTGCCGCctcc                                                                                                                                          >  2:272208/1‑300 (MQ=12)

                                                                                                                                     aGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTa                                                                                                                                         <  3:259922/301‑1 (MQ=32)

                                                                                                                                           gcctgggggAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCt                                                                                                                                   <  2:630685/295‑1 (MQ=17)

                                                                                                                                           cGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCt                                                                                                                                   <  3:605310/301‑1 (MQ=17)

                                                                                                                                            gCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCtt                                                                                                                                  >  4:384765/1‑301 (MQ=32)

                                                                                                                                            ccTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCtt                                                                                                                                  <  4:272274/300‑1 (MQ=17)

                                                                                                                                              tgtgGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTcc                                                                                                                                >  1:416076/1‑301 (MQ=32)

                                                                                                                                               gtgGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCt                                                                                                                               <  1:96699/301‑1 (MQ=32)

                                                                                                                                                 gggAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACCTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAGCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGAACTACGTCACAGCGTTCCGTGCAATGATTGGGCTTTGGCGATTTTTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGt                                                                                                                             >  4:679138/1‑301 (MQ=255)

                                                                                                                                                  ggAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGtt                                                                                                                            <  3:416078/301‑1 (MQ=32)

                                                                                                                                                   tacacgacgctcttccgatctccAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTc                                                                                                                           <  1:378562/280‑1 (MQ=255)

                                                                                                                                                   gAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTc                                                                                                                           >  4:686923/1‑301 (MQ=17)

                                                                                                                                                   gAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCtgagatcggacgagcacacgtctgaac                                                                                                                           >  4:243903/1‑276 (MQ=255)

                                                                                                                                                    aaaTTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTACAcccgatgttgcgccgccttccgcctccgtgtgt                                                                                                                           >  2:701091/1‑271 (MQ=11)

                                                                                                                                                      aTTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTc                                                                                                                        >  1:249745/1‑301 (MQ=17)

                                                                                                                                                         tCTGCATTCCGGTTACCAACGGGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAgg                                                                                                                     <  1:115221/301‑1 (MQ=17)

                                                                                                                                                         tCTGCATCACGGTTCCCAACGAGACGGCGAAGCCTACTTGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAgg                                                                                                                     <  2:384668/301‑1 (MQ=12)

                                                                                                                                                           tGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGcc                                                                                                                   <  3:609089/301‑1 (MQ=32)

                                                                                                                                                               ttACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCa                                                                                                               <  4:428147/301‑1 (MQ=32)

                                                                                                                                                                        ccAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTcagatcggaagagcacacgtct                                                                                                      >  3:378562/1‑280 (MQ=255)

                                                                                                                                                                          aaCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTTGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCgt                                                                                                    >  3:587861/1‑300 (MQ=17)

                                                                                                                                                                          aaCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCgg                                                                                                    >  4:168890/1‑301 (MQ=32)

                                                                                                                                                                           aCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCggg                                                                                                   >  3:623198/1‑301 (MQ=17)

                                                                                                                                                                            ccgtccacgacgccctcccgatctaaTAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGc                                                                                                  <  2:502404/276‑1 (MQ=255)

                                                                                                                                                                                 cGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGcttcct                                                                                             <  1:623196/301‑1 (MQ=17)

                                                                                                                                                                                 cGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTAGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGcttcct                                                                                             <  1:21415/301‑1 (MQ=17)

                                                                                                                                                                                    cgaagCCTACTGGCGGCATAACCTTTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCa                                                                                          <  2:168848/297‑1 (MQ=17)

                                                                                                                                                                                    cGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCa                                                                                          >  3:38884/1‑301 (MQ=32)

                                                                                                                                                                                                 cggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAACATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCAc                                                                             >  4:207973/1‑301 (MQ=17)

                                                                                                                                                                                                    cATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGcagatcggaagagcacacgtctgaa                                                                          >  4:502546/1‑277 (MQ=255)

                                                                                                                                                                                                     aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCcgc                                                                         >  2:418756/1‑301 (MQ=17)

                                                                                                                                                                                                     aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCcgc                                                                         >  3:379033/1‑301 (MQ=32)

                                                                                                                                                                                                        aCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGac                                                                      <  3:249744/301‑1 (MQ=17)

                                                                                                                                                                                                        aCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGac                                                                      >  1:6434/1‑301 (MQ=32)

                                                                                                                                                                                                          cTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGacac                                                                    <  4:381876/301‑1 (MQ=32)

                                                                                                                                                                                                              tCGGCCCTTGGTGGCGCCGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCtt                                                                <  4:590187/301‑1 (MQ=17)

                                                                                                                                                                                                                 gCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGcc                                                             <  1:235561/301‑1 (MQ=32)

                                                                                                                                                                                                                   ccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTc                                                           <  4:102042/301‑1 (MQ=32)

                                                                                                                                                                                                                       ggtggCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGc                                                       >  3:7203/1‑301 (MQ=32)

                                                                                                                                                                                                                              tgtgctcttccgatctactacGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACt                                               <  4:613210/286‑1 (MQ=18)

                                                                                                                                                                                                                               ggTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑ACACtt                                              <  4:418865/301‑1 (MQ=14)

                                                                                                                                                                                                                               ggTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACtt                                              <  1:38883/301‑1 (MQ=33)

                                                                                                                                                                                                                                 ttAACGGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑ACACTTcc                                            <  2:686693/301‑1 (MQ=14)

                                                                                                                                                                                                                                       cccACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGATGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACTTCCCCTTGc                                      <  4:701336/300‑1 (MQ=18)

                                                                                                                                                                                                                                        ccACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCAACCGAGACACCCTTTCCTCCGGCTACGC‑CTTTCCCTTGcc                                     >  2:379856/1‑301 (MQ=255)

                                                                                                                                                                                                                                             tactccGGCTTCGGCGGACTTCTGCACGCTCATCCTCTGGCGTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAAATCTCCCCAGGTAGATCTCGCCACCTTCCTGCTTCTGCCTGTCGGATCTACGTCACAGCGTTCCTTGCAAGTATTGGACTTTGGCGATTTGTGCCACCT‑CACCCCGCTCTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACTTCCCCTTGCCGGgtg                                <  2:678916/301‑1 (MQ=11)

                                                                                                                                                                                                                                              actacGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑CCACtagaccggaagagcttc                               >  2:613013/1‑285 (MQ=12)

                                                                                                                                                                                                                                                            tGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑ACACTTCCCCTTGCCGGGTGGGTAGAGGACTTTcc                 >  1:79067/1‑300 (MQ=14)

                                                                                                                                                                                                                                                            tGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGTCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCa                 <  1:379033/301‑1 (MQ=18)

                                                                                                                                                                                                                                                             gACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTccc                >  2:727855/1‑299 (MQ=18)

                                                                                                                                                                                                                                                                    gCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTGCGGCTC‑GCACTTCCCCTGGCCGGGGGTGTAGAGGACTTTcccaccacaa         >  2:413006/1‑292 (MQ=255)

                                                                                                                                                                                                                                                                     cACGCTCGTCCCGTCGCCTTTCGACGCGCTGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑ACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGt        >  4:64638/1‑301 (MQ=14)

                                                                                                                                                                                                                                                                      aCGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑ACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCACGTc       >  4:45242/1‑301 (MQ=14)

                                                                                                                                                                                                                                                                        gCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCAc     >  3:401107/1‑301 (MQ=33)

                                                                                                                                                                                                                                                                           cATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCAcccc  >  3:278491/1‑300 (MQ=18)

GCCTGAGGTTGGCCGGTGCAAGCACCCGCTCCATCAGCGTGTCCGGCTCCGCGTTCGTCCACGTCGCTGACGCCGCCGGCACCTTTGCACTGTCAGGCATCACCCTCGGCTTCTGGCCGGGATTCGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCT‑CACCCCGCTGTGCCGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTA‑GCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCA  >  AE015451/4348437‑4349004

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