BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,171	G→A	G389G (GGC→GGT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,171	0	G	A	95.6%	408.6 / 1.2	137	G389G (GGC→GGT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base G (3/3); new base A (64/65); total (68/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCG  >  NC_000913/66928‑67417

gTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGttt                                                                                                                                                                                                                                                 <  1:490953/251‑1 (MQ=255)

gTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGtt                                                                                                                                                                                                                                                  <  1:169781/250‑1 (MQ=255)

 gcGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTg                                                                                                                                                                                                                                                <  1:625289/250‑1 (MQ=255)

           cGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTAtt                                                                                                                                                                                                                                          <  1:55003/247‑1 (MQ=255)

            ggCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGa                                                                                                                                                                                                                                        >  1:263108/1‑248 (MQ=255)

             gCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGAt                                                                                                                                                                                                                                    <  1:512616/251‑1 (MQ=255)

              cGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGa                                                                                                                                                                                                                                        >  1:568041/1‑246 (MQ=255)

                aactGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATc                                                                                                                                                                                                                                   >  1:204926/4‑249 (MQ=255)

                 actGGCGCATATAGTTGCGGTTCAGCGCATGGATGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCAAAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGAGGAGTTTTCAACGTGTCGCTACAGTTAACACGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACACTCGCTGGACCGGTTTGGGTATTGAAGATCAgt                                                                                                                                                                                                                                >  1:77013/3‑250 (MQ=255)

                    ggCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCgg                                                                                                                                                                                                                             <  1:432537/251‑1 (MQ=255)

                        catagcggtGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGTGCTCCACCAGCGAGGATCCACGCTGAGGAAGCAGTTGGCAGATCCGGTTGCGCTTGCCGCAGCGCATTCGCTACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTACCGATCGGCGATATCAATCAGGTTGGCGACAAGCGCTGGACCGGTTTGGGTAGTGAATATCAGTCGGGCAg                                                                                                                                                                                                                         <  2:555911/243‑1 (MQ=255)

                          tatCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGGGTTTTCACCGTGTCGGTACAGTTAACCAGCAGACGGTAACGATCGCCGAGGTCAATCAGGCTGGCGACAATCGCTGGCCCGGTTTGGGGATTGAAGATCAGGCGGGCAggg                                                                                                                                                                                                                       >  2:108214/1‑250 (MQ=255)

                            gcgtggAGGTGCAGCGCAGGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTGGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTAGTGAAGATCAGTCGGGCAGGAtg                                                                                                                                                                                                                     <  2:66449/246‑2 (MQ=255)

                                   ggTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCCGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGATGATCAGACGGGCAGGATCGTCCtt                                                                                                                                                                                                               >  1:162335/1‑250 (MQ=255)

                                       cagcgCATGGCTGAAGCCGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTGTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTACCCAGACGGCGGTACCACTCCCGGGGATCGATAAGGCTGGCGACAAACGCAGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTACttcacac                                                                                                                                                                                                          >  2:249466/5‑246 (MQ=255)

                                       cagcgCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGGAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCCGGATCGGCCttcaccc                                                                                                                                                                                                          >  2:504431/5‑246 (MQ=255)

                                                   gAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGAt                                                                                                                                                                                              <  1:204136/251‑1 (MQ=255)

                                                   gAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGAt                                                                                                                                                                                              <  1:160544/251‑1 (MQ=255)

                                                       gcggtagtgTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTg                                                                                                                                                                                          <  1:572/243‑1 (MQ=255)

                                                       aCGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTg                                                                                                                                                                                          >  1:614736/1‑251 (MQ=255)

                                                        cGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCCATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTg                                                                                                                                                                                          >  2:374255/1‑250 (MQ=255)

                                                        cGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGCATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTg                                                                                                                                                                                          >  1:331106/1‑250 (MQ=255)

                                                         ggTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGaa                                                                                                                                                                                        <  1:379952/251‑1 (MQ=255)

                                                            aTGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGGGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGCTCAATCAGGCTGGCGACAATCGCTGGACCGGGTTTGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGTGATGCTGACCGt                                                                                                                                                                                     >  2:190454/1‑251 (MQ=255)

                                                                 gcgcgcCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAg                                                                                                                                                                                 <  1:504431/250‑1 (MQ=255)

                                                                         ccaGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCCCCGGCAGATTCAGCAGGGAGTGAGGAGGTTTCACCGTGTCGCTACAGTTAACAAGCAGGCGGTAACGATCGCCGAGATAAATCAAGCTGGCGACAATCAATGGAACGGTTTGGGTAATGAACATCAGGCGGGCAGGATAGTCCTTAACACCCATACCCAGAAGCTGAACATCGGGGAACGGttt                                                                                                                                                                        >  2:407891/1‑251 (MQ=255)

                                                                           cgCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCCGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGGTTGGGGATTGAAGATCAGACGGGCAGGATCGTCCTTACCACCAATCCCGGAATGCTGAACGTCGAGGATCGGTTtct                                                                                                                                                                      >  2:116081/2‑251 (MQ=255)

                                                                            gCGAGGATCCACGCTTCGTAAGCAGTTGGCAGACCCGGTTGCGCTTTCCACAGCGCATTCGCCCCCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTtctc                                                                                                                                                                     <  1:249466/251‑1 (MQ=255)

                                                                                   cccACGCTTCGGAATCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAAGACCGAGATGCTGAACGTCGAGGATCGGTTTCtcttctacg                                                                                                                                                              <  1:9541/250‑4 (MQ=255)

                                                                                       cGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGGGCGGGGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCtcttatacggca                                                                                                                                                           >  1:22993/1‑242 (MQ=255)

                                                                                         cTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGAAAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATc                                                                                                                                                        >  1:48996/1‑251 (MQ=255)

                                                                                             ggAAGCAGTTGGCAGATACGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATAAAGTTAACCAGCAGACGGTAAAGGTCGCAGAGATCAATCAGTCTGGCGACAATCGGTGGACCGGTTTGGGGATTGAAGATCAGTCGGGCAGGATAGTACTTACCACCAATACCGAGATGATGAAAGTCGGGGAGAGGTTTCtcttcatcggcaaacgacg                                                                                                                                                    >  2:360152/1‑237 (MQ=255)

                                                                                               aaGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTgg                                                                                                                                                                                                                                               <  1:184594/158‑1 (MQ=255)

                                                                                               aaGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACggg                                                                                                                                                  >  1:565057/1‑251 (MQ=255)

                                                                                                 gCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGGCGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCCATACCGAGATGCTGAACGTCGCGGAACGGTTTCTCTTCTACGGCAATCGACGGGc                                                                                                                                                 >  2:655671/1‑250 (MQ=255)

                                                                                                    gTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCTCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTCACCAGCAGACTGTAACGATCGCCGAGATCAATCCGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGCAGATCAGTCGGGCAGGATCCACCTTCCCACCAATACCGAAATGCTGAACGTCGAGGATCGGTTTCTCTTCTTCGGCAATTGCCGGGCAGAc                                                                                                                                             >  2:569582/1‑251 (MQ=255)

                                                                                                       ggCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACtt                                                                                                                                           <  1:454739/250‑1 (MQ=255)

                                                                                                        gCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTCACGATCGCCGCGATCAATCAGGCTGGCGCCAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCTATACCGCGATGCTGAACGCAGAGGATCGGTTTCTCTTCTACGGCACTCGACGGGCAGGCTTcc                                                                                                                                         >  2:157323/1‑251 (MQ=255)

                                                                                                              ccGGTTGCGCTTTCCACAGCGCCTTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGGCAATCGCTGGACCGGTTTGGGTATTGACGATCAGTCGGGCCGGATCGTCCTTACCACCAATACCGAGGTGATGAACGTCGAGGGTCGGGTTCTCTTCTACGGCAATcggcggaccggcttcccgcatt                                                                                                                                   >  2:531907/1‑235 (MQ=255)

                                                                                                               cGGTTGCGCTTTCCACAGCGCATGCGTCACCGGCAGTTTCGGCCGGGAGTGCGGGGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCatat                                                                                                                                  <  2:527212/251‑1 (MQ=255)

                                                                                                                         ttCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCg                                                                                                                          >  1:472589/1‑249 (MQ=255)

                                                                                                                          cccACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGc                                                                                                                       <  1:185829/250‑1 (MQ=255)

                                                                                                                            cacaGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCAc                                                                                                                     <  1:655671/251‑1 (MQ=255)

                                                                                                                             acaGCGCAGTCGCCACCGGCCGTTGCGGCAGGGAGTGCTGTGTTTTCGCCGTGTGGATACAGTTAACCCGCCGACGGTAACGATCGCCGAGAGCAAGCAGGCGGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTGTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGTATATGGGATCCGAGcgcc                                                                                                                    <  2:127941/251‑4 (MQ=255)

                                                                                                                                    aTTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCCATCGACGGGCAGACTTCCAGCATATGGGAGCCCAGCACCAGGTCAt                                                                                                             >  2:163333/1‑251 (MQ=255)

                                                                                                                                      tCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCAtt                                                                                                            <  1:167797/250‑1 (MQ=255)

                                                                                                                                       cGCCACCGGCATTTTCGGCAGGGATTGCGGTGGTTTCACCGTGTCGATACAGTTAACCAGCAGGCGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGGATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTa                                                                                                           <  1:286857/250‑1 (MQ=255)

                                                                                                                                       cGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTTGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTa                                                                                                           <  2:484964/250‑1 (MQ=255)

                                                                                                                                        gCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGCGCACCAGGTCATTAcc                                                                                                         >  1:710408/1‑251 (MQ=255)

                                                                                                                                         ccACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATg                                                                                                                                 >  1:698547/1‑226 (MQ=255)

                                                                                                                                         ccACCGGCAGTTACGGCAGGGAGTGCGGTGGTTGCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCGGGCGACAATCGCAGGACCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCGACGTCAATCGACGGGCAGACTTCCAGCATATg                                                                                                                                 <  2:698547/226‑1 (MQ=255)

                                                                                                                                          cACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGGCGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCtt                                                                                                       <  1:677786/251‑1 (MQ=255)

                                                                                                                                             cGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGGCATTACCTTtct                                                                                                    <  1:603160/251‑1 (MQ=255)

                                                                                                                                             cGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTAGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGAATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCCGCATATGGGAGCCGAGCACCAGGTCATTACCTTtct                                                                                                    >  1:126826/1‑251 (MQ=255)

                                                                                                                                                cAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGCCGGGCAGACTTCCAGCATATGGGCGCCGAGCACCAGGTCATTACCTTTCTCg                                                                                                  >  1:125102/1‑250 (MQ=255)

                                                                                                                                                  gTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTCCCTTTCTCGa                                                                                                 >  1:353336/1‑249 (MQ=255)

                                                                                                                                                    ttCGGCAGGGAGGGCGGTGTTTTCCCCGTGTCGATGCAGGTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATTGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGt                                                                                              <  2:630660/250‑1 (MQ=255)

                                                                                                                                                      cGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGa                                                                                            <  1:163333/250‑1 (MQ=255)

                                                                                                                                                       ggCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATa                                                                                          <  1:116528/251‑1 (MQ=255)

                                                                                                                                                          aGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAggg                                                                                       >  1:588929/1‑250 (MQ=255)

                                                                                                                                                          aGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGAAGACGGTAACGATCGCCGAGATCAATCAGGCTTGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACTGTCAGACTTCCAGCATATGGTCTCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAgtt                                                                                       >  1:257188/1‑249 (MQ=255)

                                                                                                                                                                 gCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGtc                                                                                 <  2:565057/250‑1 (MQ=255)

                                                                                                                                                                  cgggtttttCACCGTGGCGATGCAGGTAACATGGAGACGGTAACGATCGCCGAGATCAAGCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGGTCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCTAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGGCATTACCTTTCTCGAAGTGATAGGTGTAGtccg                                                                               <  2:381160/246‑2 (MQ=255)

                                                                                                                                                                      gtTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGCGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTCCTCGaatggttaggtgatgtcctccct                                                                           >  2:604908/1‑230 (MQ=255)

                                                                                                                                                                      gtTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGATTCTCTTCTACGGCAATCGACGGGCAGCCTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCAt                                                                           >  1:656407/1‑251 (MQ=255)

                                                                                                                                                                      ggtctCACCGTGTCGATACTGTGAACCAGCAGACGGTAACGATCGCCAAGATCAATCAGGCGGGCGACAATCGCTGGATCTGTTTGGGTATTGAAGAGCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCGCTTCTACTGCAAACTACGGGCACACGTCCAGCATATGGGAGCCTAGCACCAGGTCATTACCTTTCTCGAAGTGATAgtcgttgacctccat                                                                           <  2:401792/247‑15 (MQ=255)

                                                                                                                                                                           cacaggggcGATACAGTTAACCAGCAGACGGTAACGATCGCGGAGATCAATCAGGCGGGCGACAATCGCTGGTCCGGTTTGGGTATTGAAGATCAGTCTGGCAGGATCGTCCTTACCGCCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCTCCAGGTCATTTCCTTTCTCGAAGTGAGAGGTGTAGTCCTCCATAAagg                                                                      <  2:657410/243‑1 (MQ=255)

                                                                                                                                                                                ggtCGATCCAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCCGTCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGAAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTg                                                                 <  2:428156/250‑1 (MQ=255)

                                                                                                                                                                                ggtCGATACAGTAACCCAGCAGCCGGTACCGACCGCCGAGATCAATCAGGCGGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGGCGGGCAGTATCGTCCTTACCTCCAATACCGAGATGCTGAACGTCGAGGATCGGTTGCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTg                                                                 <  2:366582/250‑1 (MQ=255)

                                                                                                                                                                                     atacAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGAAGGCTCGTCCTTACCACCAATACCGAGATGCTGAACGTAGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGAATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTaccgcc                                                            >  1:103998/4‑251 (MQ=255)

                                                                                                                                                                                       acAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTTCTCCATAAAGGAGGTGCCGcccg                                                          <  2:296606/250‑2 (MQ=255)

                                                                                                                                                                                        cAGTTAACCAGCAGACGGTAACGATCGCCGAGATCACTCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATCCCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCATCCATCAAGGAGGTGCCGCCCTg                                                         >  1:288936/1‑251 (MQ=255)

                                                                                                                                                                                         aGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGc                                                        >  1:179960/1‑251 (MQ=255)

                                                                                                                                                                                             aaCCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGAc                                                    <  1:407059/251‑1 (MQ=255)

                                                                                                                                                                                              aCCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGTGGATCGGTTTCTCTTCTACGGCAATCGACGGACAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGCTAGGTGTAGTCCTCCATCAAGGAGGGGCAGCCCTGCAGAcc                                                   >  2:37222/1‑251 (MQ=255)

                                                                                                                                                                                               ccAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCg                                                  <  1:108214/251‑1 (MQ=255)

                                                                                                                                                                                                 agcaGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGCGGTGCCGCCCTGCAGACCGGt                                                >  1:238008/4‑251 (MQ=255)

                                                                                                                                                                                                   caGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCggctg                                              >  2:514711/2‑248 (MQ=255)

                                                                                                                                                                                                    aGACGGTAACGATCTCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTg                                              >  1:142729/1‑250 (MQ=255)

                                                                                                                                                                                                    aGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTg                                              >  1:667920/1‑250 (MQ=255)

                                                                                                                                                                                                         gTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCGGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGGCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCa                                        <  2:689093/251‑1 (MQ=255)

                                                                                                                                                                                                           aaCGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGc                                                                                                                                      <  2:228849/155‑1 (MQ=255)

                                                                                                                                                                                                           aaCGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGc                                                                                                                                      >  1:228849/1‑155 (MQ=255)

                                                                                                                                                                                                             cGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCCATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTCCCTTTCTCGAAGTGATAGGGGGAGTCCTCCCTAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCtt                                    >  2:432504/1‑251 (MQ=255)

                                                                                                                                                                                                                 cGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGCAGAGATGGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCatg                                >  1:479236/1‑251 (MQ=255)

                                                                                                                                                                                                                    cGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCatgatg                             <  1:705563/251‑1 (MQ=255)

                                                                                                                                                                                                                       gaTCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGa                          <  2:204926/251‑1 (MQ=255)

                                                                                                                                                                                                                        atcaatcaGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCa                                                                                                                                     >  2:535039/1‑143 (MQ=255)

                                                                                                                                                                                                                        atcaatcaGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCa                                                                                                                                     <  1:535039/143‑1 (MQ=255)

                                                                                                                                                                                                                          caatcagGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAg                        <  1:531907/250‑1 (MQ=255)

                                                                                                                                                                                                                                ggCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGAGCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGACCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCACAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGggcg                 <  2:162335/250‑1 (MQ=255)

                                                                                                                                                                                                                                aGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGAGCGTCCTTACCGCCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGACGTGATAGGTGTAGTCCTCCATAAAGGGGGTGCCGCCCTGCCGACCGGTTGACATCACCTTCATGATTCGAAGCAGggcg                 >  1:541017/1‑251 (MQ=255)

                                                                                                                                                                                                                                     gCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGggcggcggt            >  1:682960/1‑248 (MQ=255)

                                                                                                                                                                                                                                      cGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGtt           >  2:498733/1‑251 (MQ=255)

                                                                                                                                                                                                                                          aaTCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTgccg                                                              >  1:329338/1‑196 (MQ=255)

                                                                                                                                                                                                                                          aaTCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCGCCATAAAGGAGGTgccg                                                              <  2:329338/196‑1 (MQ=255)

                                                                                                                                                                                                                                           aTCGCTGGACCGGTTTGGGAATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCCTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGGGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTcc       >  2:281478/1‑250 (MQ=255)

                                                                                                                                                                                                                                             cGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGGTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGGGCCGCCCTGCATACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGAGGTTTTCCAGt    >  2:18783/1‑251 (MQ=255)

                                                                                                                                                                                                                                              gcttggCCGGTTTGGGTATTTAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGGTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGCGCACCTGGACATTACCTTTCTCGAAGTGATAGGTGTAGGCCTCCATAAAGGAGTTGCCGCCCTGCTGAACGGTTGAAATCACCTTCATGCTTCGAGGCGGGGCGGAGGTTTTACAGTc   >  2:22423/5‑251 (MQ=255)

                                                                                                                                                                                                                                               cTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATACGAAGCAGGGCGGCGGGTTTCCAgacg  >  2:482941/1‑248 (MQ=255)

GTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCG  >  NC_000913/66928‑67417

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