BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,383	A→C	20.5%	V289G (GTG→GGG)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,383	0	A	C	20.5%	159.7 / 13.3	100	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base A (53/24); new base C (2/20); total (55/45)
Fisher's exact test for biased strand distribution p-value = 7.04e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCC  >  NC_000913/3760146‑3760610

gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGGGAAACTCa                                                                                                                                                                                                                                      >  2:520674/1‑238 (MQ=255)

     gATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCAcccg                                                                                                                                                                                                                                  >  2:1184312/1‑235 (MQ=255)

         tGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGGTCATCTCCTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCCCCTAACAACTTGTCGCC‑GGGGAAACTCACCAAACTgc                                                                                                                                                                                                                             >  2:956296/1‑237 (MQ=255)

           aaTGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAAAAATAATTCCTG‑CCTG‑CGGCCCGCCTAAACACTTGTCGCC‑GGAGAAACTCACCCGACTG‑aag                                                                                                                                                                                                                          >  2:1119427/1‑236 (MQ=255)

               aCCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCCGACTG‑Acgcccg                                                                                                                                                                                                                      >  2:1024758/1‑238 (MQ=255)

                ccTCTGCGCTGCGGGTAAGCGGGCGAGGGGTCGGTATTTTTTCCCTCAGAGCTTCAGGGTGTAAATAAAGACGCACCGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTACCCGCGCGTAATGCACGCTTCAGCGGGTGGTTT‑TGCGGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATCCCTG‑CCTG‑CGGCCCGCCTAACACCTTGTCGCC‑GGAGAAACTCCCCAGACTG‑Acgcccgc                                                                                                                                                                                                                     <  1:654259/238‑1 (MQ=255)

                ccTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑Acgcccgcc                                                                                                                                                                                                                    >  2:524713/1‑239 (MQ=255)

                    tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTTCTG‑CCTG‑CGGCCCGCCCAACAACTTGTCGCC‑GGAGAAAATCACCAGACTG‑ACGCCcgccgcca                                                                                                                                                                                                                >  2:181639/1‑237 (MQ=255)

                    tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑AACCACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGGGAAAATCACCAAACTG‑ACGCCCGCCCCaa                                                                                                                                                                                                                >  2:434287/1‑239 (MQ=255)

                          gggggTAAGAGGGCGCAGTGTCGGTATTTTTCCACTCAGAGCTACAGGGTGTAAAAAAAGACG‑AAAGTGGTCTCCCGGGCGCGGGGGGGTCATTTTCTCCGCGCTTAATGCACGCTTCAGGGGGTGGCTT‑TGCAGGGGGGCGATCATCTCTTTTTC‑ACCAAAAATATTTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAAc                                                                                                                                                                                                           <  1:237974/236‑1 (MQ=255)

                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGCCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATATCTTTTTT‑ACCCACAATAATTCCTG‑CCTC‑CGGCCCGCCCAACCACCTGTCGCC‑GGCAAAACTCACCACACCG‑ACGCCccccgccacccact                                                                                                                                                                                                          >  1:1034606/1‑230 (MQ=255)

                              gTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCCATCAACTCCt                                                                                                                                                                                                       >  2:419216/1‑238 (MQ=255)

                                        gcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTTCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTTACCCGACTG‑ACCCCCGCCCCAATACACTCCTCCGGCCTTgg                                                                                                                                                                                             >  2:664365/1‑238 (MQ=255)

                                        gcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTCACAACTTGTCGCC‑GGAGAAACTCCCCAGACTG‑ACGCCCGCCGCAATCAACTc                                                                                                                                                                                                         >  2:521860/1‑226 (MQ=255)

                                                    ttttttCCCCCACGGCTTGGGGGTGAAAAAAAGGAGGAACCGTGCCTCCCAGGCC‑CGCGAGGGCCTTTTTACCCGGCCGAAATCCCCGCTCGAGCGGTGGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑CCCAACAATAATCCCTG‑CCTG‑CGGCCCGCCTAACACCTTGTCGCC‑GGAGAAACTCCCCGGACGG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCaaa                                                                                                                                                                                <  1:269451/238‑1 (MQ=255)

                                                        ttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTTGCTCTTTCGGCAAACCGt                                                                                                                                                                            >  2:51402/1‑239 (MQ=255)

                                                        ttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCCG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAAATCACCAGACTG‑ACGCCCGCCGCAATCCACTCCTGCGGCATTGGCTCTTTCGGCCAAACga                                                                                                                                                                            >  2:899859/1‑238 (MQ=255)

                                                         tcccccAGATCTTAAGGGTGAAAATAAGGACGCAACGTGGCTCCCAGCCC‑CGCGAGGGTCATTTTACCCGCGCTTAATCCACGCTTCAGCGGGTGGTTT‑TGCAGGCGGCCGATCATCTCTTTTTT‑CCCAACAATAATTCCTC‑CCTG‑CGCCCCCCCTAACAACTTGTCGCC‑GGAGAAACTCCCCAGACTC‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                           <  1:744110/234‑1 (MQ=255)

                                                         tcccTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                           <  2:465687/236‑1 (MQ=255)

                                                         taccccagggcataaggggcaaaagaacGACGCAGGGTGGCAGCCGGCGC‑GGGGGGGCTCATTTTGCCGGCGCGAACTGCGCGTTTCACGGGGTGGCTT‑GCCCGCGGGGTGATCATCTCTTTTTA‑ACCAAAAAAATTCCCGG‑CTTG‑TGCCCGGCCAACCATCTTCCCGAG‑AGAGAACCACTCAAGATGG‑AGCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                           <  1:484461/218‑1 (MQ=255)

                                                         tCACTCAGAGCTTCAGGGTGTAAAGAAAGACGCAACGTGGCTTCCAGCGC‑CGCGGGGGTCATTTTACCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATCCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCG‑GGAGAAACCCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                           <  2:913145/239‑1 (MQ=255)

                                                         tCACCCAGAGCTTCGGGGTGAAAAAAAAGCCGGCACGGGGCTCCCAGCGC‑CGCGGGGGTCATTTTCTCCGCGTGTAATGCACGCTTCAGCGGGTGGTTT‑TGCAGGCGGGCGATCATCTTTTTTTT‑ACCAACAATAATCCCTC‑CCGG‑CGCCCCGCCTAACAACTTGTCGCC‑GGAAAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                           <  1:328497/239‑1 (MQ=255)

                                                         tCACCCAGACCTCGAGGGTGAAAATAAAGCCCAAACGTGGCTTCCAGCGC‑CGCGGGGGTCATTTTATCCGCGCGTAAGGCCCTCTTCAGCGGTTGGCTT‑GGCAGGCGGGCGCTCTTCTCTTTTTT‑ACCAACAATAATCCCGG‑CTGG‑CGGCCCGCCTAACAACTTGTCCCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                           <  1:245078/239‑1 (MQ=255)

                                                            cTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGGTCATCTCTTTTTT‑TCCAAACATAATTCCTG‑CCTG‑CGGCCCCCCTAACAACTTGTTGCC‑GGAGAAAATCACCCGACTG‑CCGCCCGCCGCCATCAAATCCTGCGGGCTTTGCTCTTTCCGCAAACCGtcatg                                                                                                                                                                        >  1:142854/1‑235 (MQ=255)

                                                            cTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACCCCTGCGGCAATGGCTCTTTTGGCAAACCGTa                                                                                                                                                                           >  1:195238/1‑236 (MQ=255)

                                                                           tgtAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑TCCAAAAATAATTCCTT‑CCTG‑CGGCCCGCCTAACAACTTGTCTCC‑CGAGAAACTCCCCCGACTG‑GCGCCCCCCCCAATCAACTCCTGCGGCATTTGCCCTCTCGGCACACCCTACTTGCTAAAATcccaca                                                                                                                                                          >  1:115502/1‑233 (MQ=255)

                                                                           tgtAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTTTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCa                                                                                                                                                          >  2:192171/1‑238 (MQ=255)

                                                                             aaaataaaGCCGCACCGGGGCATCAGGCGC‑GGCGGGGGTATTTTTACCCGCGCGAAAGGCCCCGTCCAGCGGGGGGTTT‑TGCGGGCGGGCGATCTTCTTTTTTCT‑CCCAAAAATATTCCCGG‑CCGG‑CGGCCCGCATAAAATCTTGGCGCG‑GGAAAAACCCCCCAGCCGC‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGc                                                                                                                                                        <  1:419216/237‑1 (MQ=255)

                                                                                 aaaaGCCCCACCGGGGCTTCCGCCGC‑CGCGGGGTTCTTTTTACCGGCGCGTAATGCACGTTCCGGGGGGTGTCTT‑TGCGGGCGGGCGTCCATCTCTTTTTC‑CCCACAAATATTCCCCC‑CCGG‑CGGCCCCCCTAACACCTTGTCGCC‑GAAAACACCCCCCAGCCTG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGcc                                                                                                                                                   <  1:990018/238‑1 (MQ=255)

                                                                                  aaggaccaaaCGGGGCTCCCTGCGC‑CGGGGGGGTCATTTTACCCCGGAGTAAGGCCCGTTTCAGCGGGTGGTGT‑TGCGGGGGGCCGACCACCTCTTTTTT‑ACAAAAAATAATCCCGG‑CCTG‑CGCCCCGCCTAACACCTTGCCGCG‑GGAAAACCCCCCCAGACGG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                  <  1:966106/231‑1 (MQ=255)

                                                                                  aaaGCCGCAACGTGGCTCCCAGCGC‑CGCGGGGTTCATTTTACCCGCGCGTAATGCCCGCTTCAGCGGGTGGTTT‑GGCAGGCGGGCGATCATCTCTTTTTT‑CCCAAAAATAATCCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                  <  2:1117461/239‑1 (MQ=255)

                                                                                   aaGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCCG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCCGACTG‑ACGCCCGCCGCAATCAACCCCTGCGGCATTGGCTCTTTCGGCAAACCGCACTGGCTAAAATCCACCCGCGAGCCa                                                                                                                                                  >  1:344121/1‑238 (MQ=255)

                                                                                     gACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCCATCAAATCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACt                                                                                                                                                >  2:744250/1‑238 (MQ=255)

                                                                                     gACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACt                                                                                                                                                >  2:966106/1‑238 (MQ=255)

                                                                                     gACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACt                                                                                                                                                >  2:671111/1‑238 (MQ=255)

                                                                                       cGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGCGAAACTCACCCGACTG‑ACGCCCCCCGCAATCAACTCCTCCGGCATTTGCTCTTTCGGCCACCCGTACTGGCTAAGATCGACTCGCGAACCCCTGc                                                                                                                                              >  1:164752/1‑238 (MQ=255)

                                                                                       cGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTg                                                                                                                                               >  2:826777/1‑237 (MQ=255)

                                                                                          aaCGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGGTCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGGCTG‑GCGCCCGCCGCAAACAAATCCTGCGGCAATGGCTCTTTCGGCAAACCGGACTGGGTAAGATCGACCAGCggacccacggctaa                                                                                                                                          >  1:286690/1‑226 (MQ=255)

                                                                                            cGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTTTCGCC‑GGAGAAAATCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGGGGCATTGGCTCTTTTGGCAAACCGTCCTGGCTAAGATCCGCCCGCGAGGCCCTGCCTaaa                                                                                                                                         >  1:1117461/1‑238 (MQ=255)

                                                                                                 tctcaggcgc‑CGCGGGGTTCTTTTGCTCCGCGTGGATTCCCCGCTCCGGGGGGTGGCTT‑TGACGGGGGCCTATCATCT‑TTTTTT‑ACCCAAAAAAATTCCCG‑CCCGCCGCCCCCCCTAACACCTTGCCCCC‑GGAAAAACCCCCCACCCTG‑CCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                                    <  1:403480/232‑1 (MQ=255)

                                                                                                  ggagggggcccgggggggTCATTTTCCCGGCGCTAAATGCACGTTCGAGGGGGTGGTTT‑CGCCGGCGGGCGATCTTCTTTTTTTT‑ACCAACAAAAATCCCTG‑CCTG‑CGGCCCGCCTAACACCTTGTCCCC‑GGAAAACCCCACCAGACTG‑ACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                                    <  1:907045/223‑1 (MQ=255)

                                                                                                      aGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTAc                                                                                                                              >  2:88345/1‑239 (MQ=255)

                                                                                                        cgc‑cgcGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGTTCCTCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCACTGGCCCCTTCGGCAAACCCTACTGGCCAAGATCGACCAGCGGCCCACCGCCTAAATCCGTCCCtcccg                                                                                                                            >  1:727132/1‑235 (MQ=255)

                                                                                                           ccggggggttCATTTTACCCGCGCTTAATGCCCGTTTCAGCGGGTGGTTT‑TGCGGGCGGGCGATCATCTTTTTTTT‑ACCAAAAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGCCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                          <  2:195238/230‑1 (MQ=255)

                                                                                                            cgcgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTTTTGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGGACTTGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAATCACTACGGGaa                                                                                                                         >  2:503824/1‑239 (MQ=255)

                                                                                                                tgggaaaaggacacagcgcaTAATGCCCGCTACGGGGGGTGGCTT‑GGCGGGCGGGGGTCAACCTCTTGTTA‑CCCAACAATAAGTCCGG‑CCGG‑AGGCCCGCCAAAAACCTTGCCGCG‑GGAAAACCTCACCAGTTAG‑AGGACCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                          <  1:428247/219‑1 (MQ=255)

                                                                                                                   gTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                          >  2:907045/1‑231 (MQ=255)

                                                                                                                           ccccgccctaaatccccgtttCAGGGGGTGCCTT‑GGCGGCGGGGCGACCATCTTTTTTTC‑CCAAACAATAATCCCGG‑CCGG‑CGGCCCCCCTAAC‑ACCTGTCCCCGGGAAAACCCCCCCAGCCGG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGc                                                                                                          <  1:671111/220‑1 (MQ=255)

                                                                                                                             ccGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACCACTTGTTGCC‑GGAGAAACTCACCCGACTG‑ACGCCCCCCGCCATCAACTTCTGCGGCATTGGCTCTTTCGGCAAAACGTACTGGCTAAGCTCGACCCGCGAGCCACCGCCCAAATCAGCCACTACGGGAACATCCAGCCATTTGCCg                                                                                                        >  1:575510/1‑239 (MQ=255)

                                                                                                                                cgcgTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCCTCTCTTTTTT‑ACCAACAATAATTACTG‑CCTG‑CGGCCCGCCCAACCACTTTTTCCC‑GGAGAAACTCACCAAACTG‑CCGCCCCCCGCAATCAACCCCTGCGGCACTTGCTCTTTCGGCAAACCGGCCTGGCTAAAAACGACCACAGAACCACTTCCTAAATACCTCACCACGGGAACCTCCCGCTCTTTGCCgcg                                                                                                      >  1:1022562/1‑236 (MQ=255)

                                                                                                                                 gcgTAATCCCCCCTCCAGCGGTTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑CCCAAAAATAATCCCTG‑CCTG‑CGGCCCGCCTAACAACTTGCCGCC‑GGAGAACCTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgc                                                                                                     <  2:344121/238‑1 (MQ=255)

                                                                                                                                  cgTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGTCGGGCGGTCATCTCTTTTTT‑TCCCACCATAATTCCCG‑GCTG‑CCGCCCGCCTAACCACTTTTCGCC‑GCAGAAACTCACCAGACTG‑ACGCCCGCCGCCATCCACTCCTCCGGCACTTGCTCTTTCGGCAAACCGTACTGGCCTAGATCGACCCCCGCGCCCCCTCCTAAATCAGTCACTACGGGCCCACCCCGCCCTTTcccccgcc                                                                                                    >  1:645574/1‑233 (MQ=255)

                                                                                                                                  cgTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCCTCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgcg                                                                                                    >  2:206075/1‑238 (MQ=255)

                                                                                                                                  cgTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACCATAATTCCTG‑CCTG‑CGGCCCGCCCAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑TCGTCCGCCTCAATAAACTCCTGCCGCATTGTCTCTTTCAGCAAACCGTACTGGCTAAGATCGAACAGCAAGCCACTGACTAAAACAACAACAACGGTACCATCCAGCTCTTTGTCGAgcg                                                                                                    >  1:919462/1‑238 (MQ=255)

                                                                                                                                  cgTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACCG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAAATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACAACCACATCTTTTCCGAgcg                                                                                                    >  1:300057/1‑238 (MQ=255)

                                                                                                                                     aaTGCACGCCTGAG‑GGGTGGGTTGTGCAGGCGGGCCATCA‑CTTTTTTTT‑ACCAAACAAAATTCCCGCCCGG‑CGCCCCGCTGAACACCTTGTCGCC‑GGAAAACCTCCCCAGCGTG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGccac                                                                                                <  1:1064150/239‑1 (MQ=255)

                                                                                                                                           cGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCCCTGCCTAAATCAGTCACTACGGGAACATCCAGct                                                                                                                >  2:428247/1‑217 (MQ=255)

                                                                                                                                             cTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATTATCTCTTTTTT‑ACCAACAATAATTTCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCCGCTCTTTGCCGAGCGCCACCCATTcc                                                                                         >  2:227662/1‑238 (MQ=255)

                                                                                                                                                 aGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAAAAACCTGTCGCC‑GGAGAAACTCACCAGACTG‑ATGCCCGCCGCAATCAACTCCTGCGGCCTTGGCTCTTTCGGCAAAACGTACTGGCTAAGATCGACCCGCGAGCCACTCCCTAAATCAGGCACTACGGGAACATCCAGCTCTTTGCCGCGCGCCACCAGTTCCGCTTc                                                                                    >  1:589954/1‑239 (MQ=255)

                                                                                                                                                   cGGGTGGTTT‑TGCAGGCGGGCGATCATCCTTTTTTT‑CCCAAAAATAATCCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAAAAACTCCCCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                  <  2:1170651/175‑1 (MQ=255)

                                                                                                                                                        tgctttggcgggcgggggAACAT‑TCTTTTTTCCCCAATAATAATCCCGC‑CCTC‑GCCCCCGCTTAACAACTTGCCCCC‑GGAAACCCCCCCAAAACGG‑GCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTa                                                                               <  1:122799/228‑1 (MQ=255)

                                                                                                                                                        ggCTT‑TGCAGGCGGGCGGTCCTCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTTCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCa                                                                                                                                  >  2:948311/1‑186 (MQ=255)

                                                                                                                                                        ggCTT‑TGCAGGCGGGCGGTCATCTCTTTTTT‑ACCAACCATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTTTCGCC‑GGAGAAACTCACCCGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTTGCTCTTTCGGCAAACCGTACTGGCCAAGATCCACCCGCGAGCCACCGCCTAAATCAGTCACTTCGGGAACATCCAGCTCTTTGCCGAGCGCCCCCAGATCCGCTTCATCTAt                                                                              >  1:524663/1‑238 (MQ=255)

                                                                                                                                                        ggCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCCAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAt                                                                              >  2:841063/1‑238 (MQ=255)

                                                                                                                                                        ggCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACCCCCGCCCCAATCAACTCCTGCGGCATTGGCTCCTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAt                                                                              >  1:1034933/1‑238 (MQ=255)

                                                                                                                                                        ggCTT‑TGCAGGCGGGCGATATTCTCTTTTTT‑ACCAAAAATAATCCCTG‑CCTG‑CGGCCCGCCTAACACCTTGTCGCC‑GGAGAAACTCACCAGCCTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAt                                                                              <  2:1034606/238‑1 (MQ=255)

                                                                                                                                                           cgtcgagggccggcgACCACCTTTTTGTT‑ACAAACAATATCCCCTC‑CCGC‑CGGCCCCCCCAACACCTTGACGCC‑GAAGAACACCCCCAGACCG‑CGCCCCGCCGCAATCAACTCCGGCGTCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCg                                                                            <  1:977515/231‑1 (MQ=255)

                                                                                                                                                                cgggcgggctACCCTCTCTTTTTT‑ACAAACAATAATCCCTC‑CCTG‑CGCCCCGCCTAAAACTTTGTCCCC‑GGAGAACCTCCCCAGACGC‑ACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCt                                                                                <  1:756849/227‑1 (MQ=255)

                                                                                                                                                                 cggcgggcaCTCATCTTTTTTTT‑ACCAACAATAATCCCTC‑CCTG‑CGGCCCGCCTAACAACTGGTCGCC‑GGAAAAACCCACCAGACTG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGt                                                                     <  2:119590/238‑1 (MQ=255)

                                                                                                                                                                 aGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGt                                                                     >  2:23191/1‑239 (MQ=255)

                                                                                                                                                                   gggggcgATCATCCCTTTTTT‑ACCAACAAAAATCCCCG‑CCTG‑CGGCCCGCCAAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                    <  2:500541/236‑1 (MQ=255)

                                                                                                                                                                   gcgggctATCATCTTTTTTTA‑CCCAACAATAATTCCTC‑CCGG‑GGGCCCCCCAAACAACTTGTCGCG‑GGAAAACCCCCCCAGACTG‑CCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                    <  2:575510/238‑1 (MQ=255)

                                                                                                                                                                   gcgggcgATCATCTCTTTTTT‑CCCAACAATAATCCCTG‑CCTG‑CGGCCCGCCTAACACCTTGTCGCC‑GGAGAAACTCCCCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                    <  2:114463/238‑1 (MQ=255)

                                                                                                                                                                     gggcgATCATCTCTTTTTT‑CCAAACAAAAATTCCGC‑CCGG‑CGCCCCGCCAAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                  <  2:300057/238‑1 (MQ=255)

                                                                                                                                                                         tatcatcTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                             <  2:1190011/238‑1 (MQ=255)

                                                                                                                                                                         gATCTTCTTTTTTTT‑ACCAACAAAATTCCCTC‑CCTG‑CGCCCCCCCTAACAACTTGTCGCC‑GGAAAAACTCACCAGACGG‑ACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                             <  1:826777/239‑1 (MQ=255)

                                                                                                                                                                          atcatcTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTgg                                                                      >  2:756849/1‑229 (MQ=255)

                                                                                                                                                                           tcatcTTTTTTTC‑ACAAAAAAAACTCCCTG‑CCGG‑CGCCCCGCTAAACACGTTGCGGGC‑GGACAACCACCCCAGCCTG‑ACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                           <  1:1024758/239‑1 (MQ=255)

                                                                                                                                                                           tcatcTCTTTTTT‑ACCAACAATAATTCCTC‑CCTG‑CGGCCCGCCTAACAACTTGGCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                           <  2:449485/239‑1 (MQ=255)

                                                                                                                                                                              cctctTTTTT‑ACAAAAAATATTCCCGG‑CCGG‑CGGCCCGCCTAACATCTTGTCGCC‑GGAGAAACCCACCAGCCTG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATg                                                        <  1:206075/238‑1 (MQ=255)

                                                                                                                                                                                tttttttt‑ACCAAAAAAAATTCCTG‑CCGG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGCCTA‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                      <  1:520674/237‑1 (MQ=255)

                                                                                                                                                                                tctTTTTT‑ACCAACAATAATTCCTC‑CCTG‑GGGCCCGCCTAACAACTTGCCGCC‑GGAAAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                      <  2:524663/239‑1 (MQ=255)

                                                                                                                                                                                tctTTTTA‑CCCAACAATAATCCCTG‑CCTG‑CGCCCCGCCTAACAACTTGTCGCC‑GGAGAACCCCACCAGACTG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                      <  2:276405/239‑1 (MQ=255)

                                                                                                                                                                                 ctTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGCAt                                                         >  2:672941/1‑235 (MQ=255)

                                                                                                                                                                                        tcccAACATTAATTCCGG‑CCTG‑GGCCCCGCCAAACACCTTGTCGCC‑GGAGAAACCCACAAGATTG‑CGGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAACCCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCTCTTTGGTGAACCCCTGAATGCTGTAGTTa                                               <  2:1147582/237‑1 (MQ=255)

                                                                                                                                                                                             cataataaactccgcccgtcgccccGCTTAACACCTGCTGCCC‑GGAGAATCTCCCCAGGCGC‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTg                                            <  1:636576/217‑1 (MQ=255)

                                                                                                                                                                                             acaaTAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                          >  1:346016/1‑239 (MQ=255)

                                                                                                                                                                                                  attccccg‑cctg‑cTCCCCGCCTACCAACTTGGCGCC‑GAACAACCTCACCAGACGG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTGCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                     <  1:524713/232‑1 (MQ=255)

                                                                                                                                                                                                  aatTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGGCGCC‑GGAGAAACCCCCCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                     <  2:1119036/239‑1 (MQ=255)

                                                                                                                                                                                                         g‑cctg‑cGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                               >  2:304789/1‑238 (MQ=255)

                                                                                                                                                                                                         g‑cctg‑cGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                               >  2:742232/1‑238 (MQ=255)

                                                                                                                                                                                                         g‑cctg‑cGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                               >  1:146529/1‑238 (MQ=255)

                                                                                                                                                                                                         g‑cctg‑cGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCAAGCTCTTTGCCGAGCGACACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATCTCCtttac                               >  1:1040943/1‑236 (MQ=255)

                                                                                                                                                                                                             tg‑cGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGAAAACCGTACTGGCTAAGATCGACCAGCAAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCAGCTTCATCTATCGCTTTGGTGAACCCATGAATGCTGTAGTTACTGGTATGAACTTTAATaaa                           >  1:731630/1‑239 (MQ=255)

                                                                                                                                                                                                             cg‑cgcCCCCCCTAACAACTTGCCCCC‑GAAG‑AACCCCCCCGACCGTACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                           <  1:744250/238‑1 (MQ=255)

                                                                                                                                                                                                                cGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACa                         >  1:70522/1‑239 (MQ=255)

                                                                                                                                                                                                                      gCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCgg                   >  1:784119/1‑239 (MQ=255)

                                                                                                                                                                                                                        cTAACACCTGGCCGCG‑GGAGAACCTCCCCAGCCTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTa                 <  2:164752/239‑1 (MQ=255)

                                                                                                                                                                                                                                  gTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc        >  1:318470/1‑238 (MQ=255)

                                                                                                                                                                                                                                       c‑GGAAAAACCCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGcc  <  2:346016/239‑1 (MQ=255)

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCC‑GGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCC  >  NC_000913/3760146‑3760610

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