BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	64,692	C→T	V363V (GTG→GTA)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	64,692	0	C	T	98.7%	226.2 / ‑4.6	75	V363V (GTG→GTA)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (35/39); total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

AGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAG‑ATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGAT  >  NC_000913/64446‑64935

aGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCCCT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCg                                                                                                                                                                                                                                                   >  1:12678/1‑251 (MQ=255)

aGGCTGCGCCATACCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCg                                                                                                                                                                                                                                                   >  1:704933/1‑251 (MQ=255)

       gcCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGAGCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGcc                                                                                                                                                                                                                                            <  1:134592/251‑1 (MQ=255)

        ccATGCCTTCCACCAGCCCGACGGGATCACTAAGAAAGGTGCGGATGATCGACGGGTCCAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGACGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTCGGCGCGACATAACCAGCGCTGTGCATACGAGGAAAATAGAG‑ATGACCACACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCacggcg                                                                                                                                                                                                                                           >  1:554159/1‑246 (MQ=255)

           tGCCTTCCACCAGCCCGACGGGATCAATCCGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCCCT‑TCGCCGGGATTAGGCGCGACATAACCAGCGCGGGGCATACGCGGAAAATAGGG‑ATGACCAAACGCCGCCACCGATCCCCCGGGGCGGTCTACCGGCGGGCCggtc                                                                                                                                                                                                                                        >  2:983521/1‑248 (MQ=255)

            gCCTTCCACCAGCCCGACGGGAGCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCa                                                                                                                                                                                                                                       <  1:945218/251‑1 (MQ=255)

             ccTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGAGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGGGGCGGTCTACAGGCAGGCCGTTCAc                                                                                                                                                                                                                                      >  1:950443/1‑251 (MQ=255)

                  caccaGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGtcg                                                                                                                                                                                                                                 >  1:954131/1‑249 (MQ=255)

                          cGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACCAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGGGGCGGTCTACCGGACGGCCGTTAACCgcagccacgtcg                                                                                                                                                                                                                         >  2:351493/1‑240 (MQ=255)

                               ggagtaatcaGAAAGGTGCGGATGCTCGACGGGTACATGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACTTAGCCGCCAGGGCTGGCGTGTGGCGGCACG‑TCGCCGAGATTAGGTGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTaa                                                                                                                                                                                                                    <  2:584245/246‑1 (MQ=255)

                               ggATCAATCAGAAAGGTGCGGCTGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACAGGCAAGCCgatccccggcccccgttcgggga                                                                                                                                                                                                                     >  2:935608/1‑228 (MQ=255)

                               ggATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTaa                                                                                                                                                                                                                    <  1:983521/251‑1 (MQ=255)

                                  acaatcaGAAAGGGGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCAAAAAGCCCTGTCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTaaaaa                                                                                                                                                                                                                 <  2:950443/250‑1 (MQ=255)

                                   caatcaGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAAt                                                                                                                                                                                                                <  1:987014/251‑1 (MQ=255)

                                         gAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATg                                                                                                                                                                                                          >  1:394736/1‑251 (MQ=255)

                                              gTGCGGATGAGCGACGGGTACAGGCTTTTATAGTCCAGCTCCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCAGCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTc                                                                                                                                                                                                     <  2:171959/251‑1 (MQ=255)

                                                gCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGCGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAg                                                                                                                                                                                                   <  1:323638/251‑1 (MQ=255)

                                                  ggATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGtt                                                                                                                                                                                                 >  1:800480/1‑251 (MQ=255)

                                                         cGACGGGTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACCTAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGGTCGAGTACAAATGACATGATTTCAGTTTTATggg                                                                                                                                                                                          >  1:936446/1‑250 (MQ=255)

                                                             gggTACAGGCTTTTATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGaaaat                                                                                                                                                                                      <  1:667280/251‑4 (MQ=255)

                                                              ggTACAGGCTTTTATAGGTCAGCACCAGTACCGAATCAGAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATc                                                                                                                                                                                     <  1:533195/251‑1 (MQ=255)

                                                                       ttttATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCCCCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCAcc                                                                                                                                                                            >  1:419660/1‑251 (MQ=255)

                                                                        tttATAGTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTAACGGCAGGCCGTTCACAGACGCTCGTTCGAGTAAAAATGGCATGATTTCTGTTTTCTGGaaaaactggggaccc                                                                                                                                                                            >  2:520436/1‑237 (MQ=255)

                                                                              gTCCAGCACCAGTACCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCCGCCGGGCGTTCACCCTCCCCCGGTCGAGTAAAAATTGCATGGTTTTAGTTTTTTGGAAAATCttctttacccgctagc                                                                                                                                                                     >  2:684579/1‑236 (MQ=255)

                                                                                   gcaccagtacCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGGTCCCCGTCGCTCGTTCGAGGAAAACTGGCATGATTTCAGGTTTATGGAAAATCtcgctccccagcttcgaatct                                                                                                                                                                >  2:13467/1‑231 (MQ=255)

                                                                                   gcaccagtacCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGAACT‑TCGCCGCGATTAGGCGCGACATAACCAGCGCGGGGCATACGGGGAAAATAGAG‑ATGACCAAACACCGCCACCGATCCCCCGTGGCGGTATACCGGCAGGCCGTGCACCGTCGCTCGTTCGAGTAAAAAAGGCTTGGTTTTAGTTTTATGGAAAATCTggctaccccggccgcaaact                                                                                                                                                                >  2:825712/1‑232 (MQ=255)

                                                                                     cacagtaccgAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCttt                                                                                                                                                             <  2:534394/242‑1 (MQ=255)

                                                                                     accagtacCGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCttt                                                                                                                                                             <  1:840602/251‑1 (MQ=255)

                                                                                        agtaccgAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAg                                                                                                                                                          <  1:903952/245‑1 (MQ=255)

                                                                                        agtaacgAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGTG‑ATGACCAAACGCCGCCACCGATCCCCCGTAGCGGTCTACCCGCCGGCCGTTCCCCGTCGCTCGTTCGCGTAAAAAAGCCCAGATTTCAGATTTATGGAATATCTG‑GTCACCCAGCCCGCaaccgttc                                                                                                                                                            >  1:667554/7‑243 (MQ=255)

                                                                                         gtaccgAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCCGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGGGCATACGAGGAAAATAGGG‑ATGCCCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCCCCGGCGCTCGTTCGAGGAAAAATGGGATGATATCAGTTTTATGGAAAATCTTGGTCCCC‑AGTCCTCAATCTTTCAgg                                                                                                                                                         >  2:74351/6‑251 (MQ=255)

                                                                                          taccgAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGGTCACCGTCGCTCGTTCGAGTAAAAATAGCATGTTTTCAGTTTTATGGAAAAACTGCGTCACC‑AGATCGCAAACTTTCAGGt                                                                                                                                                        >  2:930517/5‑251 (MQ=255)

                                                                                           accgAATCATAAAGCCCTGGCCGGGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGt                                                                                                                                                        <  1:52860/247‑1 (MQ=255)

                                                                                              gAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGGTCACCCTCGCGTGGTAGTGTAAAAATGGCATGGTTTCAGTTTTATGGAAAATCTGAGGTACC‑AGCTCGCAATCTTTTAGGttca                                                                                                                                                     >  2:676429/1‑248 (MQ=255)

                                                                                                aTCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAg                                                                                                                                                  <  1:811190/251‑1 (MQ=255)

                                                                                                aTCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAg                                                                                                                                                  <  1:454149/251‑1 (MQ=255)

                                                                                                  cATAAAGCCCTGGCCGTGAATCCATCACGGAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGt                                                                                                                                                 <  2:534656/250‑1 (MQ=255)

                                                                                                      aaGCCCTGGCCGGGAATCCATCACGGAGCCGCCAGGGCGGGCGTGTGCCGGCCCA‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGACACCGATTCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTTGCAATCTTTCAGGTTATAAGtcgtca                                                                                                                                            <  2:712413/251‑6 (MQ=255)

                                                                                                         cccTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTAGCTCGTTCGAGTAAAAAAGGCCTGATTTCCGTTTTATGGGAAATCTGCGTCACC‑AGCTCGCAATCTTTTAGGGTATaaatcgacagcg                                                                                                                                         >  2:728905/1‑241 (MQ=255)

                                                                                                            tGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGCGCTTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCCTCGCTCGTTCGTGGAAAAATGGCATTATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCCATCTTTCATGTTATAATTCGCCAgagccg                                                                                                                                      >  2:887563/1‑246 (MQ=255)

                                                                                                               ccGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGCGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCCGGTTATAAGTCGCCAGCGCAGGtt                                                                                                                                   >  1:228135/1‑251 (MQ=255)

                                                                                                                        cATCACGTAGCCGCCAGGGCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCgg                                                                                                                          <  1:826186/251‑1 (MQ=255)

                                                                                                                                aGCCGCCATGGCTGTGGTGTTGCGGCACT‑ATGCAGAGATTAGGCGAGACATAACCAGCGCTGTGCATTCGCGGAAAAAAGAG‑ATGACCAAACGCCGCCACCGATTCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGGCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTGATGGAAAATCTGCGTCACC‑TACTCGCAAGCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAAc                                                                                                                    <  2:779621/249‑1 (MQ=255)

                                                                                                                                         ggCTGTCGTGTGGCGGC‑CTGTCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGGGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAAtt                                                                                                        <  2:817406/251‑1 (MQ=255)

                                                                                                                                          gCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAAttt                                                                                                       <  1:92365/251‑1 (MQ=255)

                                                                                                                                          gCTGGCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAAttt                                                                                                       <  1:492444/251‑1 (MQ=255)

                                                                                                                                             ggCGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCCACGATCCCCCGTGGCGGGCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTc                                                                                                      <  1:520436/249‑1 (MQ=255)

                                                                                                                                              gggtgcgggcggaccgacGCCGAGAGTAGGCGCGACATAACCAGGGCGATGAATTCGAGGAAAAT‑GAGGATGACCAAATGCCGCCGCCGATCCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGTCATGATTTCAGTTTTGTGGAAGATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGt                                                                                                    <  2:353533/234‑1 (MQ=255)

                                                                                                                                               cGTGTGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCggcc                                                                                                  >  1:824414/1‑248 (MQ=255)

                                                                                                                                                     gcggcACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGa                                                                                            <  1:559907/251‑1 (MQ=255)

                                                                                                                                                      gggcACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGAt                                                                                           <  1:755995/250‑1 (MQ=255)

                                                                                                                                                               cGCCGAGATTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGT‑TCAATTTCGTCCATTCGATCCCACgg                                                                                    <  1:369568/250‑1 (MQ=255)

                                                                                                                                                                    agatctggcgagataTAACCAGCGCGCGGCATGCTAGGAAAATAGAG‑ATGACCAATCGGCGCCACCGATCCCACGTGGCGTTCTACCGGCAGGCCGTTCGCCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCTTCACC‑AGCTCGCAATCCTTCAGGTTATAAGTCGCCAGCGCCGGTTGAACTTCGTCGATACGGCGG‑TCATTTTCGTCGATTCGATCCCACGGTTTATc                                                                              <  2:482222/237‑1 (MQ=255)

                                                                                                                                                                    agaTTAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAAATTTGTCCATTTGATCCCACggggtttc                                                                              >  1:171737/1‑246 (MQ=255)

                                                                                                                                                                       gtAGGCGCGCTATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATTGCCAATCGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGGCGCTCGTTCGAGTAAAAATGGCATGAGTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGAt                                                                           <  2:608198/250‑1 (MQ=255)

                                                                                                                                                                        tAGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGAt                                                                           <  1:247377/250‑1 (MQ=255)

                                                                                                                                                                         agggccaaCATAACCAGCGCGGTGAATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGATa                                                                          <  2:695411/243‑1 (MQ=255)

                                                                                                                                                                         aGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCtt                                                                                                                                                              >  1:881385/1‑167 (MQ=255)

                                                                                                                                                                         aGGCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCtt                                                                                                                                                              <  2:881385/167‑1 (MQ=255)

                                                                                                                                                                          ggCGCGACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCCTTCGATCCCACGGGGTATCGAt                                                                           >  1:885295/1‑248 (MQ=255)

                                                                                                                                                                              cgACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGATTATCTTCGGCGCAACGGCGG‑TCAAATTCGTCCATTCGATCCCCCGGGTTATCGAttagtttt                                                                    >  1:399849/1‑244 (MQ=255)

                                                                                                                                                                               gACATAACCAGCGCGGTGCATACGAGGAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTc                                                                   <  1:261697/251‑1 (MQ=255)

                                                                                                                                                                                                    acgaggAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTATACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGTAAAATCTGCGTCAAC‑AGCTCGCAATCTTTAAGGTTATAAGACGCCAGCGCAGGGTTATCTTCGGCGAAACGGCGG‑TCCATTTCTTCCCTTCTATCCCACGGGTGCTTGATAGATTTTCTTTCGCCTCATaaatcccg                                                >  2:775717/5‑242 (MQ=255)

                                                                                                                                                                                                     cgaggAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGGTTTCGGTAGATTTTTCTTCGCCTAATAACTCCTGAGc                                             >  1:361865/4‑251 (MQ=255)

                                                                                                                                                                                                      gaggAAAATAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑GCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGAATTTCCTTCGCCTAATAACTTCTGAGCg                                            >  1:313509/3‑251 (MQ=255)

                                                                                                                                                                                                          aaaaTAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGAcgg                                        <  1:770950/251‑3 (MQ=255)

                                                                                                                                                                                                           aaaTAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGAcggt                                       <  1:684579/251‑4 (MQ=255)

                                                                                                                                                                                                           aaaTAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGGCTACCGGCAGGCAGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAAGCTTTCAGGTTATAATTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCAGTCGAGCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGAcggt                                       <  1:956647/251‑4 (MQ=255)

                                                                                                                                                                                                             aTAGAG‑ATGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTTTGGAAAATCTGCGTCACA‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAACAGGCGG‑TCAATTTTGCCAATTCGGTCCACCGGGTTTTCGGTAGATTTTCCTTCGCCTCATaaatcctgtgtgaaggggt                                     >  1:416747/1‑233 (MQ=255)

                                                                                                                                                                                                                gag‑aTGACCAATCGTCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCTGTTCACCGTCGCTCTTTCGAGTAAAAATGGCATGATTTCAGTTTTATGTAAAAACTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTGCGTCCATTCGATCCCACGGGTTAGCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCa                                  <  1:618466/251‑1 (MQ=255)

                                                                                                                                                                                                                  g‑aTGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGCGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGGCCAATTT‑GTTCATTCGATCCCCCGGGTGTTCGATAGATTTTCCTTCCCCTAATAACTCCTggaccacggtttccacc                                >  2:805450/1‑235 (MQ=255)

                                                                                                                                                                                                                     tGACCAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCAGGCCGTTCCCCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCCGTTTTATGGAAAATCTGCGTCACC‑AGCTAGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTTGATCCCACGGGTTATCGATAGATTTTCCTTTCCCTAATAACTcccgaacgccgggttccagcgg                              >  1:97478/1‑231 (MQ=255)

                                                                                                                                                                                                                         cAAACGCCGCCACCGATCCCCCGTGGCGGTCTACCGGCCGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGACACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGGTTATCTTCGGCGCAACGGCGG‑ACAATTTCGTCCTTTTGATCCCCCGGGCTATCGATAGATTTTCCTTCCCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGaa                           >  1:898768/1‑250 (MQ=255)

                                                                                                                                                                                                                                           ccccgTGGCGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATTATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCGTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAAc        <  1:513782/249‑1 (MQ=255)

                                                                                                                                                                                                                                                 ggcGGTCTACCGGCAGGCCGTTCACCGTCGCTCGTTCGAGTAAAAATGGCATGATTTCAGTTTTATGGAAAATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTCGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGAt  <  2:800480/249‑1 (MQ=255)

AGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTGAATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACT‑TCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAG‑ATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACC‑AGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGG‑TCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGAT  >  NC_000913/64446‑64935

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