BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,388	A→C	13.1%	S287R (AGT→AGG)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,388	0	A	C	13.1%	263.1 / 3.7	126	S287R (AGT→AGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base A (56/47); new base C (3/18); total (60/66)
Fisher's exact test for biased strand distribution p-value = 7.26e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

ATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTC  >  NC_000913/3760152‑3760625

     aaTGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTT‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCCGACTGAcg                                                                                                                                                                                                                                         >  1:780484/1‑238 (MQ=255)

               gcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCaa                                                                                                                                                                                                                               >  1:227702/1‑238 (MQ=255)

                cgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGGCGTCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCCAACCA‑CTTGTCGCC‑GGAGAAACTCACCAGAATGACCCCCGCCCCAATc                                                                                                                                                                                                                             >  2:1314062/1‑239 (MQ=255)

                   ggggggTAAGCAGGCGTGGGGGTGTTAATTTTTCACGAAGAGCTGCAGGGTAT‑AAATAAAGACGCAGCGTGTCTTCCAGCGCCGCGAGGG‑TCATTTTACCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGCTCATCTCTTTTTTACCAACAATAATCCCTG‑CCTGCGCCCCCCCTAACAT‑CTTGTCGCC‑GGAGAAACCCACCAGCCTGACGCCCGCCGCAATCaa                                                                                                                                                                                                                           <  2:796582/235‑1 (MQ=255)

                     cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAAATCACCAGACTGACGCCCGCCGCCATCAACt                                                                                                                                                                                                                         >  1:498153/1‑238 (MQ=255)

                        gTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCt                                                                                                                                                                                                                      >  1:585189/1‑238 (MQ=255)

                                          gttagttttcCGCTCAGAGCTTCGGGGTGT‑AAATAAAGACGCAACGGTGCCTCCAGCGCCGGGGGGG‑TATTTTTACCGGCGCGTAATG‑CACGCTTGAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTTTTTTGTCCCAACAAAAATCCCTC‑CCTGCGCCCCGCCAAACAC‑CTGGTCGCC‑GGAGAAACGCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCg                                                                                                                                                                                                    <  2:14641/233‑1 (MQ=255)

                                           gTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGT‑CTTTGCAGGCGGGCCCTCCTCCTTTTTTTCCCAACAAAAATTCCTG‑CCTTCGGCCCGCCTAACAA‑ATTGTCGCC‑GGCGAAAATCCACAGCATGACTCCCGACCCACACAACTCCCTCGGCATCGGCCATTTctct                                                                                                                                                                                                  >  2:224037/1‑236 (MQ=255)

                                           gTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGGTCATCTCTTTTTTTCCAACAATAATTCCTG‑CCTGCGGCCCGCCTAAAAA‑CTTTTCGCC‑CGAGAAAATCACCAGACCGACGCCCGCCCCCATCAACTCCTCTGGCCTTTGCCCTTTTgg                                                                                                                                                                                                   >  2:705759/1‑238 (MQ=255)

                                                ttttCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGGTCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTTGCC‑GGAGAAACTCACCAGACTGACGCCCCCCGCAATCAACTCCTGCGGCATTGGCCCTTTCGGCCAAc                                                                                                                                                                                              >  2:239088/1‑238 (MQ=255)

                                                  ttCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCACTCAACTCCTGCGGCATTTGCCCTTTCGGCAAACCGt                                                                                                                                                                                           >  1:154885/1‑239 (MQ=255)

                                                   tCACTCAGAGCTTAAGGGTGT‑AAATAAAGCCGCACCGTGGCTCCCAGCGCCGCGAGGT‑TCATTTTACCCGCGCGTAATG‑CCCGCTTCAGCGGGTGG‑CTTGGCAGGCGGGCGATCATCTCTTTTTTCCCAAAAATAATTCCTG‑CCTGCGGCCCGCCTAAAAA‑CTTGTCGCC‑GGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                                          <  1:239088/239‑1 (MQ=255)

                                                   tCACTCAGACCTTCAGGGTTT‑AAATAAAGACGCAAGGGGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATC‑CACGCTTCAGCGGGTGG‑CTTTGCGGGCGGGCGATCATCTTTTTTTTCCCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAACCTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                                          <  1:857656/239‑1 (MQ=255)

                                                   ccccTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CCCGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTTTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                                          <  1:782402/236‑1 (MQ=255)

                                                   ccACGCAGAGCTCCAGGGTGA‑AAAAAAAGGCGCAACGGGGCTCCCCGCGGCGGGGGTG‑TCATTTTACCCGCGCGTAGTG‑CACGCTCCAGGGGGGGG‑TTTGGCAGGCGGGCGATCATTTCTTTTTTCCAAACAAAAATCCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAACCCCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                                          <  2:291162/238‑1 (MQ=255)

                                                    ccactgatggctgcgggTGTAAAATAAAGCCCCAAGGTGCTTTCCAGCCCCGGGGGGT‑TCTTTTTACCCCGGCGTACTG‑CCCCTTTCAGCGGGTGC‑TTTGGCAGGGGGGCGATCTTCTCTTTTTTCCCAACAATAATCCCTG‑CCGGCGCCCCGCCTAACAA‑CTTGTCGCC‑GGAAAAACCCACCAGCCGGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                                          <  2:241154/225‑1 (MQ=255)

                                                        cAGAGTTTCAGGGTGT‑AAAAAAAGACGCAACGGGGCTCCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGc                                                                                                                                                                                     <  1:887983/239‑1 (MQ=255)

                                                           agCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTAACAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGATCTTTCGGCAAACCGTACTGGCTaa                                                                                                                                                                                  >  1:1048005/1‑239 (MQ=255)

                                                            gCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTa                                                                                                                                                                                   <  1:1261200/237‑1 (MQ=255)

                                                                  ggggtaaaaaaaaggccGCAAGGGGGCTTCCGGCGCCGGGGGGT‑TCTTTTTACCCGCGCGTATGG‑CCCGCTTCGGGGGGTGG‑TTTGGGGGGCGGGCGATCATCTTTTTTTTCCCAACAATAATTCCTG‑CCGGCGCCCCGCCAAACAA‑CTTGTCGCC‑GGAAAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCg                                                                                                                                                                             <  2:173725/222‑1 (MQ=255)

                                                                   ggTGT‑AAATAAAGACGCAACGTGCCTTCCAGCGCCGCGAGGT‑TCTTTTTACCCGCGCTTAATG‑CACGCTTCAGGGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACAAAAAATAATCCCTG‑CCTGCGCCCCGCCTAACAA‑CTTGTCGCC‑GGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGAc                                                                                                                                                                           <  2:129161/238‑1 (MQ=255)

                                                                     tgt‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCa                                                                                                                                                                         >  1:871926/1‑238 (MQ=255)

                                                                       t‑aaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGCTCATCTCTTTTTTTCCAACAATAATTCCTG‑GCTGCGGCCCGCCCTACAA‑CTTTTCGCC‑GGAAAAAATCACCATTCTGAAGACCGCCCCAATCAACACCTGCGGCACTGGCACTTTTGGGAGACCGTActaactaaagtcaccccgag                                                                                                                                                                      >  2:882049/1‑227 (MQ=255)

                                                                         aaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGa                                                                                                                                                                     >  1:86723/1‑239 (MQ=255)

                                                                            aaagaacgaaccggggCTCCCGGCCCCGCGAGGG‑TTATTTTACCCGCGCTTAATG‑CCCCCTCCGGGGGGTGT‑CTTTGCGGGCGGGCGATCTTCTCTTTTTACCCAACAATAATCCCTG‑CCTGCGCCCCGCCTAACAA‑CTTGCCGCG‑GGAAAACCTCACCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGcc                                                                                                                                                                  <  2:739975/225‑1 (MQ=255)

                                                                             aaccaCGCAACGTGGTGTCAAGCGGCGAGGGCG‑TTATTTTCCCCGCGGAGAATG‑CCCGCGACAGCGGGGGG‑GTTGGCAGGCGGGCTAACCTCTCTTTATACCCAAAAAGAATGCCTG‑CCTGCCGCGCGCCCACCAA‑CTTGCCGCC‑GGAGAACCCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                                 <  2:544867/235‑1 (MQ=255)

                                                                             aaatccgcatgagggCTCCCGGCGCGGTGAGGG‑TATTTTCATCCGCGAG‑AATGCCCCGCTCCAGGGGGTGG‑CTTTGCGGGCGGGCGATATTCTCTTTTTTCCCAACAA‑AATCCCCGCCCTGCGCCCCGCCTAACAC‑CTTGCCGCC‑GGAAAACCCCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                                 <  2:952658/226‑1 (MQ=255)

                                                                             aaaGGCGCAGCGGGGCTTCCAGCGCCGCGAGGG‑TCTTTTTACCCGCGCGTAATC‑CACGCTTCAGCGGTTGG‑CTTGGCAGGCGGGCAATCATCTTTTTTTTACAAACAATAATTCCTG‑CCTGCGGCCCGCCAAACAA‑CTTGTCGCC‑GGAAAAACTCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                                 <  2:189804/239‑1 (MQ=255)

                                                                             aaaGCCGCACCGTGGCTCCCAGCGCCGCGAGGG‑TCATTTTACCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑TTTTGCAGGCGGGCGATCATCTCTTTTTTCCCAACAATAATCCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                                 <  1:850478/239‑1 (MQ=255)

                                                                              aaagccccaagggggttcccggcccgcGGGGGTTCTTTTTACCCGCCGTTAATG‑CCCGCTTCGGCGGGTGG‑CTTGGCAGGCGGGCGTTCATCTTTTTTTTACCAAAAAAAATCCCTG‑CCCGCGCCCCGCCTAACAA‑CTTGCCGCG‑GGAGAAACCCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                                 <  2:1162896/216‑1 (MQ=255)

                                                                              aaGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGTGTTCATCTCTTTTTTCCCCACCATAATTCCTG‑CCTCCCGCCCGACCAACAA‑CTTTTCCCC‑GGAGCACATCACCAGACTGACGCCCGCCACCAACAAATCCCGGGTAAATGGCTCTTTTGGCAAAACGGACCGGCTAAGacgcgcccacaaccca                                                                                                                                                                 >  2:1194505/1‑223 (MQ=255)

                                                                                gACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCCCACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTCAGATCCACCAGCGAGCCACt                                                                                                                                                               >  2:129739/1‑238 (MQ=255)

                                                                                gACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCCCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTCACAA‑CTTGTCGCC‑GGAGAAACTCACCAGAATGACGCCCGCCGCCATCAACTCCTGCGGCATTGGCTCTTTCGGCAAAACGTACTGGGTAAGATCGACCAGCGAGCCACt                                                                                                                                                               >  1:268460/1‑238 (MQ=255)

                                                                                        ggggtctCCCGCGCGGGGGGGG‑GCATTTTACCCGCGGTGAGTG‑CACCGGTGAGGGGGGGT‑CTTTGCGGGGGGGAGACAGTCTCTTTTTCACAAAAAAAATGCCCTC‑CCTGCGCCCCGCCAAAAAA‑CTGGTCGCG‑GGAGAACCCCCCCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTaaa                                                                                                                                                        <  2:259634/231‑1 (MQ=255)

                                                                                           gCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGAACATCTCTTTTTTACCAACAATAATTCCTG‑CCCGCGGCCCGCCCAACAA‑CGTGTCGCC‑GGCGAAAATCACCAGACTGAGGCCCCCCCAAATCCACTCCCGCGGCATTGGGCCGTTCGGCTAACACTACTTGCTAAGATCGACCAGCGAGCCGCATCCTAAAACac                                                                                                                                                    >  2:1301437/1‑237 (MQ=255)

                                                                                            cTTCCAGGGCGGGGGGGT‑TCATTTTACCCGCGCGTAATC‑CACGCTTGGGCGGGGGG‑CTTGGGAGGCGGGCGATCATCTCTTTTTACCCAACAATAATTCCTG‑CCGGCGCCCCCCATAACAA‑CTTGCCGCC‑GGAGAACCTCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                                                   <  2:928817/238‑1 (MQ=255)

                                                                                                cAGCGCCGCGAGGT‑TCATTTTATCCGCGAGTAAGG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATCCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACCCACCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTa                                                                                                                                              <  1:994292/239‑1 (MQ=255)

                                                                                                       cggggggtgctttttACCCGCGCGTAATT‑CACCCTTCGGGGGGGGG‑TTTTGCGGGCGGGCGATCCTCTTTTTTTACCCAAAAAAAATCCCTC‑CCTGCGGCCCGCCAAACAA‑CTTGCCGCG‑GGAGAAACCCCCCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                                         <  2:780484/227‑1 (MQ=255)

                                                                                                            ggttcttttttcccGCGCTTATTG‑CCCGCTTC‑GGGGGTGGTCTTTGCCGGCGGGCAACCATCTTTTTTTTACCAAAAATATTCCCTC‑CCTGCGCCCCGCCAAACAA‑CTTGCCGCG‑GGAAAACCTCACCAGACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATc                                                                                                                                    <  2:1039319/226‑1 (MQ=255)

                                                                                                              gtgaTTTTATCCGGGGTTAATG‑CCCGCTGAGGGGGGTGG‑CTTTGCGGGCGGGCGATCATCTCTTTTTTACCAAAAATAATTCCTG‑CCGGCGCCCCGCCTAACAA‑CTTGTCGCC‑GGAAAAACTCACCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCa                                                                                                                                  <  2:645528/235‑1 (MQ=255)

                                                                                                                        ccGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGGTCCTCCCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑ATTGTCGCC‑CGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCCGCAAACCGTACACGCTCAGATCCACCCTCCAGGCACTCCCACACACAATCCCTACGGCTACATCCACGTCTTTTcc                                                                                                                        >  2:1083419/1‑238 (MQ=255)

                                                                                                                        ccGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGcc                                                                                                                        >  1:173320/1‑238 (MQ=255)

                                                                                                                           cgcgTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAgc                                                                                                                    >  1:623427/1‑239 (MQ=255)

                                                                                                                            gcgTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGGgc                                                                                                                    >  1:1095781/1‑238 (MQ=255)

                                                                                                                             cgTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCAAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCcgcgcg                                                                                                                   >  1:922549/1‑238 (MQ=255)

                                                                                                                                         cTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCCTCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCCACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCCAAGATCGACCAGCGCGCCACTGCCTAAATCCGTCACTACCGGAAAATCCCGCTCTTTTCCGCGCGCCACCAGTTcc                                                                                                        >  2:575101/1‑238 (MQ=255)

                                                                                                                                           tCAGCGGGGGG‑CTTTGCGGGCGGGCCATCTTCTCTTTTTCACCAAAAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑TTTGGCGCC‑GAAAAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                                                     <  1:302413/239‑1 (MQ=255)

                                                                                                                                                ggtggg‑CTTGGAAGGCGGGCGATCATCTTTTTTTTACAAAAAATAATCCCGG‑CCGGCGGCCCGCTAAACAA‑CTTGTCGCC‑GGGAAAACCCCCAAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCAt                                                                                                 <  1:1194505/234‑1 (MQ=255)

                                                                                                                                                  gTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCCCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTATATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGAGCCACCAGTTCCCCTTCATCTa                                                                                              >  1:123471/1‑239 (MQ=255)

                                                                                                                                                    gg‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCCCCAGACTGACGCCCCCCCCAATCAACTCCTTCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCCCTTTGCCGAGCGCCACCCGTTCCGCCTCATCAAt                                                                                             >  2:1163637/1‑238 (MQ=255)

                                                                                                                                                    gg‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGCCGCCCGCCGCAATCAACCCCTGCGGCCATGGCTCTTTCGGCCAACCGTACTGGCCAAGATCGACCAGCGAGCCCCTGCCTAAATCAGTCCCCACGGGAACATCCAGCTCTTTGCCGGGCGCCCCCCGTTCCGCCTCCTCTAt                                                                                             >  2:913543/1‑238 (MQ=255)

                                                                                                                                                    gg‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCTGCCAACCGGACTGGCTAAGATCGACCACCGCGCCACTGCCTAAATCACTCACTACGGGAACATCCAGCCCTTTTCCGAGCCCCACCAGTTCCCCTCCATCTAt                                                                                             >  2:1260478/1‑238 (MQ=255)

                                                                                                                                                       cTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑ATTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCCACTCCTTCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCCGCTCTTTGCCGGGCCCCACCCGTTCCGCTTCATCTATc                                                                                            >  2:298689/1‑237 (MQ=255)

                                                                                                                                                               gggggcgAACACATTTTTTTAAACAAAAATACTCCCCG‑CCGGCGGCCCGCAAAACCA‑TGGGGCGCG‑GGAAAATCCCCCCAGACAGCGCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                                   <  2:871926/236‑1 (MQ=255)

                                                                                                                                                                 gggcgTTCATCTCTTTTTTCCCAACAATAATCCCTG‑CCGGCGGCCCCCCTACCAA‑CTTGTCGCC‑GGAAAACCCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                                 <  1:1317188/238‑1 (MQ=255)

                                                                                                                                                                 gggcgCTCTTCT‑TTTTTTCACAAAAAAAATTCCCGCCCTGCGCCCCCCCAAACAA‑TTGGTCGGG‑GGAAAACCCCCCCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                                 <  2:273304/238‑1 (MQ=255)

                                                                                                                                                                 gggcgATCCTCTCTTTTTTTCCAACAATAATTCCTG‑CCCGCGGCCCGCCTAACCA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCCATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCCGACCGGCTAAGATCCACCAGCGAGCCCCCGCCTAAATCCGTCCCCACGGGAACATCCAGCTCCTTTCCGAGCGCCACCAGATCCCCTTCATCTATCGCTTTGGCGaa                                                                                 >  2:314966/1‑238 (MQ=255)

                                                                                                                                                                 gggcgATCATCTCTTTTTTACCAACAATAATCCCTG‑CCTGCGGCCCGCATAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                                 <  1:607469/238‑1 (MQ=255)

                                                                                                                                                                 gggcaATCTCCTCTTTTTTACCAACAATATTTCCTC‑CCTGCGGCCCGCCTAACAT‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                                 <  1:1288360/238‑1 (MQ=255)

                                                                                                                                                                     gATCATCTCTTTTTTCCAAACAATATTTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                                            <  1:1314062/239‑1 (MQ=255)

                                                                                                                                                                       tcatcTTTTTTTTACAAAAAAAAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTGGCCGCC‑GGAAAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                                          <  1:23588/239‑1 (MQ=255)

                                                                                                                                                                        ctttttttttttCCCAAAAATATATCCGC‑CCGGCGCCCCGCTTAACAC‑CTTGCCCCC‑GGAAAAACTCACCAGCCGGCTCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGaa                                                                         <  2:1012718/233‑1 (MQ=255)

                                                                                                                                                                          ttttttttttACCAAAAATAATCCCT‑‑CCCGCAGGCCCCCCAACAACCTGGCCGCC‑GAAGAACCCCCCCAGCCGCCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATg                                                                       <  2:1303211/235‑1 (MQ=255)

                                                                                                                                                                            tctTTTTTACCAAAAATAATCCCCG‑CCTGCGGCCCGCCTAACAA‑CTTGTCCCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                                     <  1:1260478/239‑1 (MQ=255)

                                                                                                                                                                              ttttttACCAACAATAACTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGt                                                                   <  1:487925/239‑1 (MQ=255)

                                                                                                                                                                                   tccaaacaaTAATTCCGG‑CCTGCGGCCCGCCTAACAA‑CTGGCCGCG‑GGAGAAACCCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                              <  1:312203/235‑1 (MQ=255)

                                                                                                                                                                                   tACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGCTCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCCTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGtt                                                               >  1:629527/1‑238 (MQ=255)

                                                                                                                                                                                   aaCCAACAATAATTCCGG‑CCTGCGGCCCGCCAAACAA‑CTTGCCGCC‑GGAGAAACTCACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                              <  1:768788/238‑1 (MQ=255)

                                                                                                                                                                                    tccaaaaaaaaatcccgcccggcgcCCCGCCTAACAA‑CTTCCCGCC‑GGAAAACCCCACCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                              <  2:585189/218‑1 (MQ=255)

                                                                                                                                                                                     ccAACAATACTCCCTG‑CCTGCGGCCCTCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACt                                                            <  1:710829/239‑1 (MQ=255)

                                                                                                                                                                                        acaaTAATTCCTG‑CCTGCGGCCCGCCTAACAC‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCAAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTggc                                                         >  2:463383/1‑238 (MQ=255)

                                                                                                                                                                                        acaaTAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                                         >  1:1166611/1‑239 (MQ=255)

                                                                                                                                                                                        acaaTAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                                         >  2:684446/1‑239 (MQ=255)

                                                                                                                                                                                        acaaTAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                                         >  1:378570/1‑239 (MQ=255)

                                                                                                                                                                                        acaaTAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                                         >  2:1064392/1‑239 (MQ=255)

                                                                                                                                                                                          aataatCCCTC‑CCTGCGCCCCCCCAACCAC‑CTTCGCGCC‑GGACACACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGCCCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTAt                                                       <  2:1095781/239‑1 (MQ=255)

                                                                                                                                                                                             aatccgtg‑cttgcGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCCCCATACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                                    <  1:1236876/233‑1 (MQ=255)

                                                                                                                                                                                              atccccccccgcccgCCCGCCAAAAAA‑CTGGCCGCG‑GAAAAACCCCACCGGCCGCACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                                    <  2:748248/225‑1 (MQ=255)

                                                                                                                                                                                              aatcccgccctgcGGCCCGCCTAACAT‑CTTGTCGCC‑GGAGAAACTCCCCAGCCTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                                    <  1:1163637/231‑1 (MQ=255)

                                                                                                                                                                                               tccccg‑ccggcGCCCCGCCTAACAC‑CTTGTGCCC‑GGAAAAACCCCCCAGAGTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                                                  <  2:623427/234‑1 (MQ=255)

                                                                                                                                                                                                 cccgccccgcGGCCCCCCCAACAA‑CCTGGCCCCGGGAAACACCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                                                  <  2:227702/231‑1 (MQ=255)

                                                                                                                                                                                                  ctg‑cctgcCGCCGGCGTAACAA‑CTTCTGTAC‑GGAGAAACTCACCGGACAGACCCCAGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAAACGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCATTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtc                                               <  1:1238803/239‑1 (MQ=255)

                                                                                                                                                                                                    g‑cctgcGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                              >  1:670644/1‑238 (MQ=255)

                                                                                                                                                                                                    g‑cctgcGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                              >  2:72714/1‑238 (MQ=255)

                                                                                                                                                                                                    g‑cctgcGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                              >  1:455080/1‑238 (MQ=255)

                                                                                                                                                                                                        tgcGGCCCGCC‑AACAACCTTGTCGCC‑GGAGAACCTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                                          <  1:705759/239‑1 (MQ=255)

                                                                                                                                                                                                          cGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACa                                        >  2:964038/1‑239 (MQ=255)

                                                                                                                                                                                                           ggCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAg                                       >  2:489200/1‑239 (MQ=255)

                                                                                                                                                                                                                 cccaaccaactgcgcGCA‑GGAGAACCCCACCAAGCCGGCACCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCgg                                  <  2:450317/225‑1 (MQ=255)

                                                                                                                                                                                                                             gTCGCC‑GGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                                          <  2:606939/75‑1 (MQ=255)

                                                                                                                                                                                                                             gTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                                          >  1:606939/1‑75 (MQ=255)

                                                                                                                                                                                                                             gTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc                       >  1:454821/1‑238 (MQ=255)

                                                                                                                                                                                                                                gCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACCTCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGCAGTTACTGGTATGTCCTTTCTTCAACAGTGCGGGATTttcattccgc                    >  2:214601/1‑235 (MQ=255)

                                                                                                                                                                                                                                  cgggAGAAACCCCCCGGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGc                  <  1:472381/237‑1 (MQ=255)

                                                                                                                                                                                                                                     gagaAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCt                <  1:489200/238‑1 (MQ=255)

                                                                                                                                                                                                                                               ccAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCa     >  1:715784/1‑239 (MQ=255)

                                                                                                                                                                                                                                                  gACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                              <  1:922285/179‑1 (MQ=255)

                                                                                                                                                                                                                                                   aCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                                          >  2:1051070/1‑166 (MQ=255)

                                                                                                                                                                                                                                                   aCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                                          <  1:1051070/166‑1 (MQ=255)

                                                                                                                                                                                                                                                   aCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCAt    <  1:1095757/236‑1 (MQ=255)

ATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGG‑TCATTTTATCCGCGCGTAATG‑CACGCTTCAGCGGGTGG‑CTTTGCAGGCGGGCGATCATCTCTTTTTTACCAACAATAATTCCTG‑CCTGCGGCCCGCCTAACAA‑CTTGTCGCC‑GGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTC  >  NC_000913/3760152‑3760625

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