BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,384	A→C	19.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,384	0	A	C	19.5%	155.6 / 12.2	94	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base A (52/22); new base C (2/18); total (54/40)
Fisher's exact test for biased strand distribution p-value = 1.25e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCC  >  NC_000913/3760147‑3760611

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                           tgtAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGCTCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGCCTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCa                                                                                                                                                          >  1:977515/1‑238 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                                                                        tcccAACATTAATTCCGG‑CCTG‑GGCCCCGCCAAACACCTTGTCGCCGGAGAAACCCACAAGATTG‑CGGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAACCCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCTCTTTGGTGAACCCCTGAATGCTGTAGTTa                                               <  1:1147582/237‑1 (MQ=255)

                                                                                                                                                                                             acaaTAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                          >  2:346016/1‑239 (MQ=255)

                                                                                                                                                                                                  attccccg‑cctg‑cTCCCCGCCTACCAACTTGGCGCCGAACAACCTCACCAGACGG‑CCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTGCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                     <  2:524713/232‑1 (MQ=255)

                                                                                                                                                                                                  aatTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGGCGCCGGAGAAACCCCCCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                     <  1:1119036/239‑1 (MQ=255)

                                                                                                                                                                                                         g‑cctg‑cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                               >  1:304789/1‑238 (MQ=255)

                                                                                                                                                                                                         g‑cctg‑cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                               >  2:146529/1‑238 (MQ=255)

                                                                                                                                                                                                         g‑cctg‑cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                               >  1:742232/1‑238 (MQ=255)

                                                                                                                                                                                                         g‑cctg‑cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCAAGCTCTTTGCCGAGCGACACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATCTCCtttac                               >  2:1040943/1‑236 (MQ=255)

                                                                                                                                                                                                             tg‑cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGAAAACCGTACTGGCTAAGATCGACCAGCAAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCAGCTTCATCTATCGCTTTGGTGAACCCATGAATGCTGTAGTTACTGGTATGAACTTTAATaaa                           >  2:731630/1‑239 (MQ=255)

                                                                                                                                                                                                             cg‑cgcCCCCCCTAACAACTTGCCCCCGAAG‑AACCCCCCCGACCGTACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                           <  2:744250/238‑1 (MQ=255)

                                                                                                                                                                                                                cGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACa                         >  2:70522/1‑239 (MQ=255)

                                                                                                                                                                                                                      gCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCgg                   >  2:784119/1‑239 (MQ=255)

                                                                                                                                                                                                                        cTAACACCTGGCCGCGGGAGAACCTCCCCAGCCTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTa                 <  1:164752/239‑1 (MQ=255)

                                                                                                                                                                                                                                  gTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc        >  2:318470/1‑238 (MQ=255)

                                                                                                                                                                                                                                       cGGAAAAACCCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGcc  <  1:346016/239‑1 (MQ=255)

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGC‑CGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTT‑TGCAGGCGGGCGATCATCTCTTTTTT‑ACCAACAATAATTCCTG‑CCTG‑CGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTG‑ACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCC  >  NC_000913/3760147‑3760611

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