BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,878,536	G→A	Y528Y (TAC→TAT)	gyrB ←	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,878,536	0	G	A	98.8%	524.8 / ‑3.4	170	Y528Y (TAC→TAT)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (66/99); total (69/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 170 total aligned reads displayed.

AGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTT  >  NC_000913/3878289‑3878781

agataAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAg                                                                                                                                                                                                                                                     <  1:187268/251‑1 (MQ=255)

     aGCTCTTGTAGCATTGCTTGCGGATAACCACGCTCCATACGATTTATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCGTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTg                                                                                                                                                                                                                                                <  2:124440/251‑1 (MQ=255)

           ttCAGCTTTGCTTTCGGATAACGACGCTCCATACGATTTATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGAcc                                                                                                                                                                                                                                             <  1:920580/248‑1 (MQ=255)

                cATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCGGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCaa                                                                                                                                                                                                                                     <  2:156459/251‑1 (MQ=255)

                cATTGCTTTCGGATAACGACGCTCCATACGATTGAACATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCAAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCaa                                                                                                                                                                                                                                     <  1:948932/251‑1 (MQ=255)

                    gCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGttacacgat                                                                                                                                                                                                                                 >  2:65986/1‑244 (MQ=255)

                        tCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCTCCAGCCCATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGGCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGCTTTAGACGTTACCGCGTTCAAACATTTcc                                                                                                                                                                                                                             >  2:174366/1‑251 (MQ=255)

                          ggATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCgg                                                                                                                                                                                                                           <  1:127383/251‑1 (MQ=255)

                              aaCGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGCGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATc                                                                                                                                                                                                                       >  1:242090/1‑251 (MQ=255)

                              aaCGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCAt                                                                                                                                                                                                                        <  1:991487/250‑1 (MQ=255)

                              aaCGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTCCCTGCTTGCCTTCCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCAt                                                                                                                                                                                                                        <  1:757377/250‑1 (MQ=255)

                               acgacgCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCt                                                                                                                                                                                                                      <  1:549476/251‑1 (MQ=255)

                                        cATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGCTATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGCCGAt                                                                                                                                                                                                                >  2:296776/1‑248 (MQ=255)

                                            cGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATagaaga                                                                                                                                                                                                          >  1:571599/1‑250 (MQ=255)

                                                     aTTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAAcag                                                                                                                                                                                                <  1:278221/251‑1 (MQ=255)

                                                     aTTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGGGACCGCGTTCAACGGTTTCCGGCCTCTGACGATAGAAGAAGGTCAAc                                                                                                                                                                                                  >  2:645893/1‑249 (MQ=255)

                                                      ttttCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAAcagc                                                                                                                                                                                               >  1:243103/1‑251 (MQ=255)

                                                         tCTGCGTCGCGTTGTACTCAGATCCCAGTTTCCCTAACGCTCCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCGGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCGTCACTTTGTACAGCGGCGGCTGAGCGATATAGTCGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGGCAAcagcag                                                                                                                                                                                             <  2:418826/250‑1 (MQ=255)

                                                         tCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATTCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACagcagc                                                                                                                                                                                            >  1:613031/1‑251 (MQ=255)

                                                              gTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACg                                                                                                                                                                                       >  1:334757/1‑251 (MQ=255)

                                                                             gATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTc                                                                                                                                                                        >  1:238001/1‑251 (MQ=255)

                                                                              aTACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGGCAACAGCAGCGTACGAATGTGCGAGCCGTCg                                                                                                                                                                       <  1:826804/251‑1 (MQ=255)

                                                                                ccccgtttcccgaaCGCTTCGCCAGCCAAGGCCTGTGCACTGGCGTTGGTGTGCCGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTTATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGAc                                                                                                                                                                     <  2:816909/239‑1 (MQ=255)

                                                                                cccAGTTTCTCTAACCCTCCGCCATCCAATGCCGGTGCACTGGTGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGAc                                                                                                                                                                     <  2:705994/250‑1 (MQ=255)

                                                                                  cAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGt                                                                                                                                                                   >  1:39704/1‑251 (MQ=255)

                                                                                      ttCTCTAACTCTTCCCCAGCCAATGCCGGTGCACGGGCGTTGGTGTGCAGCGGTGCGCCGTCCAGCGCGATAGAGATCTGGTACGGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTTCTTGCCTTTCTTCACTTTGTACAGCGGCTGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCg                                                                                                                                                                <  2:543879/250‑1 (MQ=255)

                                                                                            aaCGCTTCTCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTGCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGt                                                                                                                                                         <  2:980794/251‑1 (MQ=255)

                                                                                             aCGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTc                                                                                                                                                        >  1:560359/1‑251 (MQ=255)

                                                                                              cGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGCCGATAGAAGAAGGTCAACAGCAGCGTCCGAATGTGCGAGCCGTCGACGTCCGCATCGGTc                                                                                                                                                        >  1:162067/1‑250 (MQ=255)

                                                                                                cTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAGTGTGCGAGCCGTCGACGTCCGCATCGGTCatg                                                                                                                                                     >  1:815859/1‑251 (MQ=255)

                                                                                                   cGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTTCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGTTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCatgat                                                                                                                                                   <  1:816929/250‑1 (MQ=255)

                                                                                                          caATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATTGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGATCAACAACAGCGGACCAATGTGCGAGCCGTCGACTTCCGCCTCGGTCatgatgatgttgg                                                                                                                                           >  2:517947/1‑247 (MQ=255)

                                                                                                              gCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTTATCCATCGCTTCGGCGTCTTGAATGTACTGTTCCTGCTTGCCTGTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCGACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCtgtg                                                                                                                                       <  2:577657/251‑1 (MQ=255)

                                                                                                               ccGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCtgtg                                                                                                                                       >  1:485666/1‑250 (MQ=255)

                                                                                                                  ttgccctgtcGTTGGTCTGCAGCGTTCCGCCTTTCAGCGCGATAGATATTTGGTACTTATCCACTGCCTCGTCGTCTTGAATGTACTGTACCTGCTTTCCTTTCTTCACTTTGAACAGCGGCGGAAGAGCGATATATACGTGCCCGCGTTCAACTATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCTACGTCCGCATCGTTCATGATGATGATGCTGTGAt                                                                                                                                     <  2:100453/240‑1 (MQ=255)

                                                                                                                      aCTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATTTACTGTTCCTGCTTGCCCTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCa                                                                                                                               <  1:313463/251‑1 (MQ=255)

                                                                                                                         ggCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGt                                                                                                                                                                                                      >  1:524132/1‑177 (MQ=255)

                                                                                                                         ggCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGt                                                                                                                                                                                                      <  2:524132/177‑1 (MQ=255)

                                                                                                                         ggCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAAAGCAGtt                                                                                                                            >  1:303603/1‑251 (MQ=255)

                                                                                                                               tggggtcAGCGTTGCGTCGGCCACCTCGATAGGGATCAGTGACTGATCCAACGCTTCGTCTTCTTTAAGGGTCGGTTCCTTCTTGCTTGTCTTCTCTGTGTACAGCGTCGGCTGAGCGATATAGACTTGATCGCGTTCAACGATTTCCGGCATCTGACGATAGATGAAGGGCAACATCAGCGTACGAATGTGCGAGCTGTTGACGTTCGCATCGGTCAAGATGATGATGCTGTGATAATGCAGTTTGTcc                                                                                                                       <  2:973832/244‑1 (MQ=255)

                                                                                                                                     cAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGt                                                                                                                >  1:971253/1‑251 (MQ=255)

                                                                                                                                       gCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGGTATAACGCAGTTTGTCCGGGTTGTa                                                                                                               >  2:770918/1‑250 (MQ=255)

                                                                                                                                        cGTTGCGCCGTCCAGCGAGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATCGACGTGACCGCGTTCAACGATTTCCGGCATCTGCCGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCCGTTTGTCCGGGTTGtcct                                                                                                             >  2:41686/1‑248 (MQ=255)

                                                                                                                                         ggtgctccGTCCAGCGCGATAGAGATCTGGTCCTGATCCATCGCTTCGTCGTCTTTAATTTACTGTTCCGGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGTTATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTCCTc                                                                                                            <  2:542352/245‑1 (MQ=255)

                                                                                                                                           tGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGt                                                                                                          >  1:990567/1‑251 (MQ=255)

                                                                                                                                              gccgccCAGCGCGATAGAGATCTGGTCCTGATCCATCGCTTCGTCGTCTTTAATGTACTGTCCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGGGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTGTGTCCGGGTTGTACTCGTCAc                                                                                                       <  2:495870/251‑1 (MQ=255)

                                                                                                                                              gcCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTCCGAATGTGCTAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCAc                                                                                                       >  2:313353/1‑251 (MQ=255)

                                                                                                                                              gcCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCAc                                                                                                       <  1:663485/251‑1 (MQ=255)

                                                                                                                                                      gcgcGATAGAGAGCTGGTACTGAGCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATa                                                                                               <  1:340467/251‑1 (MQ=255)

                                                                                                                                                         cgATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGCACAGCGGCGGCTGAGCGATATAGACGTGACAGCGTTCAACGATTTCCGGCATCTAACGATAGAAGAAGGGCAACAGCAGCGTACGACTGTGCGAGCCGTCGACGTCCGCATCGGTCATGTTGCTGATGATGTGATAACGCAGTTTGTCCGCGTTGTACTCGTCCCGCCCGATAcc                                                                                             >  1:620941/1‑250 (MQ=255)

                                                                                                                                                          gATAGAGATTTGGTTCTGAGCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTGCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACcac                                                                                           <  2:200242/251‑1 (MQ=255)

                                                                                                                                                               agaTCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACAACAGcc                                                                                       >  1:747079/1‑250 (MQ=255)

                                                                                                                                                                         cTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTg                                                                             >  2:177888/1‑250 (MQ=255)

                                                                                                                                                                            aTCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATaa                                                                         <  1:704730/251‑1 (MQ=255)

                                                                                                                                                                              ccATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCa                                                                                                        >  1:377056/1‑218 (MQ=255)

                                                                                                                                                                              ccATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCa                                                                                                        <  2:377056/218‑1 (MQ=255)

                                                                                                                                                                                  cGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGAc                                                                                                    <  1:514477/218‑1 (MQ=255)

                                                                                                                                                                                  cGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTc                                                                    <  1:939396/250‑1 (MQ=255)

                                                                                                                                                                                  cGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCg                                                                   >  1:875533/1‑251 (MQ=255)

                                                                                                                                                                                   gCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGc                                                                  >  1:265700/1‑251 (MQ=255)

                                                                                                                                                                                     ttagtcgtcTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGAGGTGATAAGCGTCGccc                                                                >  1:859421/4‑250 (MQ=255)

                                                                                                                                                                                     ttCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCCCTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGcc                                                                 >  1:859402/1‑250 (MQ=255)

                                                                                                                                                                                      tcgtctTCTTTAATGTCCTGTTCCTGCTGGCCTTTCGTCACTTGGTACAGCGGTGGCTGAGCGAGATAAGCGGGACCGCGTTCAACGATTTCCGGAATCTGACGATAGAAGAAGTTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGcaa                                                                <  2:586351/250‑3 (MQ=255)

                                                                                                                                                                                      tcgtctTCCTTAATGTACTGTTCCTTCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGcc                                                                 <  2:84719/249‑1 (MQ=255)

                                                                                                                                                                                      tcgtcgTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACGTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCa                                                                <  1:147375/250‑1 (MQ=255)

                                                                                                                                                                                        gtcgtcTTTAATGTTCTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCAGCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACtt                                                             <  2:243103/251‑1 (MQ=255)

                                                                                                                                                                                          cgtcTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTc                                                            <  1:928804/250‑1 (MQ=255)

                                                                                                                                                                                            tcTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGCTAGAAGAAGGTCAACCACAGCGTACGAATGTGCGAGCCGTCGCCGTCCGCATCGGTCATGATGCTGATACTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGAAACCACCGCCAAGCGCGGTGCTAAGCGCCGCCACttaccg                                                         >  2:650812/1‑247 (MQ=255)

                                                                                                                                                                                            gcTTTAATTTACTGTGCCGGAGTGCCTTTCTTCACGTTGTACAGGGGCGGCTGAGCAATTTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGGCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAATCGTCGCCACTTCCtt                                                         <  2:125032/250‑2 (MQ=255)

                                                                                                                                                                                                tAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTgaga                                                      <  1:529023/250‑1 (MQ=255)

                                                                                                                                                                                                tAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTgaga                                                      <  1:200620/250‑1 (MQ=255)

                                                                                                                                                                                                tAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTgaga                                                      <  1:443516/250‑1 (MQ=255)

                                                                                                                                                                                                tAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTGGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTgaga                                                      <  1:41686/250‑1 (MQ=255)

                                                                                                                                                                                                      aCTGTTCCGGCTTGCCTTTCTTCCCTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAagagag                                                <  2:268470/250‑1 (MQ=255)

                                                                                                                                                                                                         gtttcTGCTTGCCTTTCTCCACTTTGTACAGCTGCGGCTGAGCGAGATAGACGTGACCGCGTTCAACGAGTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGGACTCGTCACTACCGATACCACAGCCAAGCGCTGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CAt                                            <  2:329201/246‑1 (MQ=255)

                                                                                                                                                                                                         gtcccTGCTTGCTGTGCTTCACTGTGTACAGGGGCGGCGGAGCGATATAGAAGTGACCGCGTTCAACGATTGCCGGCATCTGACGAAAGAAGAAGGGCAACAGCAGCGAACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCGGAGAAGAGAG‑CAtt                                           <  2:502430/248‑2 (MQ=255)

                                                                                                                                                                                                         gTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGCTATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCCCGACCGATACCACCGCCACGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGat‑catc                                           >  2:452658/1‑251 (MQ=255)

                                                                                                                                                                                                         gTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CAt                                            <  1:453248/250‑1 (MQ=255)

                                                                                                                                                                                                            tcTGCTTGCCTTTCTTCACCTTGGACAGCGGCGGCTGAGCGATAGAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCATCGTACTAAGTTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCtt                                         <  2:892084/249‑1 (MQ=255)

                                                                                                                                                                                                              tGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATc                                      >  1:899196/1‑251 (MQ=255)

                                                                                                                                                                                                               gcgtGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCg                                     <  1:131097/248‑1 (MQ=255)

                                                                                                                                                                                                               gattgacttccttcactgttttgaGCGCCGGCTGAGCGAGATAGATGTGACCGCGCTCAACGATGTCCGGTATCCGTAGAAAGTAGAAGGTCAACCGAAGAGTATGAATGGGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGTTGTGTTAACGAAGTTTGTCATGGTAGTACTCGTCACCGCAGATACTACAGACAATCGCGGTGATAAGCGTCTCCACTTCCTGAGAAGAGCG‑CATCTTATCg                                     <  2:273424/228‑1 (MQ=255)

                                                                                                                                                                                                                  tGCCTTTCTTCCCTTTGCACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCTGCAGCTGACGATAGAAGAAGGTAAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTTATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACCGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATTTTATCGaa                                   <  2:594939/250‑1 (MQ=255)

                                                                                                                                                                                                                   gCCTTTCTTCACTGTGTACAGCGGCGGCTGAGCGATATATACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGGCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAgc                                 <  2:39704/251‑1 (MQ=255)

                                                                                                                                                                                                                    cctttcctcccgtttgaccgtggcggcTGAGCGATATAGACGTGACCGCGATCAACGATTTCCGGTATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGAGAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAgcg                                <  2:869034/230‑1 (MQ=255)

                                                                                                                                                                                                                       ttcttcACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAgcgcg                              <  1:913571/250‑1 (MQ=255)

                                                                                                                                                                                                                          ttcACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCCCGTTCAAAGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGCCCGATCCCCCAGCCAAGCGCGGGGATAAGCGTCTCCACTTCCTGTGAAGAGAGACA‑CCTATCGAAGCGCGCttt                          >  2:710388/1‑251 (MQ=255)

                                                                                                                                                                                                                           tcACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGTAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCt                            <  1:22596/248‑1 (MQ=255)

                                                                                                                                                                                                                            cccTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTtct                        <  2:207333/249‑1 (MQ=255)

                                                                                                                                                                                                                               tttGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCg                      >  1:776215/1‑250 (MQ=255)

                                                                                                                                                                                                                                ttGTACAGCGGCGGCTTGGCGATATAGACGTGCCCGCTTTCAACGATTTCCGGCATCTGACTATAGAAGAAGGTCAACAGCAGCGTACGAATTTGCGGGCGGGCGACGCCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACGCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGGCGCCTCTTCCTGATAATAGAG‑CATCTTATCGAATCGCGCTTTCTCGa                     <  2:717934/250‑1 (MQ=255)

                                                                                                                                                                                                                                 tGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGTAGCGCGCTTTCTCGACg                   >  2:84046/1‑251 (MQ=255)

                                                                                                                                                                                                                                    aCAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTg                <  1:602543/251‑1 (MQ=255)

                                                                                                                                                                                                                                      aGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTGAg              <  1:65986/251‑1 (MQ=255)

                                                                                                                                                                                                                                           cggcTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGAttt         <  1:84046/251‑1 (MQ=255)

                                                                                                                                                                                                                                            ggcTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATACAGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGCCGTTGAGGCtttt        >  2:148130/1‑251 (MQ=255)

                                                                                                                                                                                                                                              cTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTGGGGATTTTAc      >  2:210488/1‑251 (MQ=255)

                                                                                                                                                                                                                                                gAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCCGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTAcc     <  2:613031/250‑1 (MQ=255)

                                                                                                                                                                                                                                                gAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTAcc     <  1:171226/250‑1 (MQ=255)

                                                                                                                                                                                                                                                gAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTAc      <  1:177888/249‑1 (MQ=255)

AGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTCCATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAG‑CATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTT  >  NC_000913/3878289‑3878781

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