BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,381	T→C	15.4%	S290G (AGT→GGT)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,381	0	T	C	15.4%	162.4 / 4.3	87	S290G (AGT→GGT)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base T (42/28); new base C (3/12); total (46/41)
Fisher's exact test for biased strand distribution p-value = 8.69e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGA  >  NC_000913/3760158‑3760613

aTGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACCCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCGGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑Acg‑cc                                                                                                                                                                                                                            <  1:66783/239‑1 (MQ=255)

    ccTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑Acg‑cccgc                                                                                                                                                                                                                         >  1:1027718/1‑238 (MQ=255)

        tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTCGCTTTGCAGGCGGT‑CGGTCCTCTCTTTTTTACCAACCATTATTACTGCCTGCCGCCCGCCTTACCACTT‑TTCGCCGGCAAAAA‑CCACCAG‑ACTG‑ACG‑CCCCCCCCa                                                                                                                                                                                                                     >  2:1084389/1‑238 (MQ=255)

             tgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAAAAATAATTCCCGCCTGGGGCCCGCCTAACAACTT‑TTCGCCCGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCac                                                                                                                                                                                                                >  2:941188/1‑237 (MQ=255)

               cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTTCTGCCTGCGGCCCGCCTCACAACTT‑GTCGCCGGTGAAAC‑TCACCAG‑ACTG‑ACG‑CCCCCCGCaacctcat                                                                                                                                                                                                              >  1:234477/1‑232 (MQ=255)

               cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑AGATCATCTCTTTTTTACCAACAAAAATTTCTGCCTGCCGCCCCCCTAACAACTT‑GTCTCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCCATCAACt                                                                                                                                                                                                              >  1:70131/1‑238 (MQ=255)

                   tAAGCGGGCGCAGGGTCGTTAATTGTTCACTCAGGGCTTCAGGGTGT‑AAATAAAGACGAAACGTGCCTCCCGGCGCGGCGAGGGTCTTTTTACCCGCGCGTAATGCTCGCTCCA‑GCGGGTGGCTTGGCGGGCGGG‑CGATCATCTCTTTTTTCCCAACAATATTCCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCGGGAGAACC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTg                                                                                                                                                                                                          <  2:817594/238‑1 (MQ=255)

                           cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGGCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCCCCGCAATCAACTCCTGCGGCAttt                                                                                                                                                                                                  >  1:960411/1‑237 (MQ=255)

                           cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCCATCAACTCCTGCGGCATTg                                                                                                                                                                                                  >  1:444207/1‑238 (MQ=255)

                                ggTCGGTAATTTTTCACTCAGAGCTTCAGGGCGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGTG‑CGATTATCTCTTTTTTACTAACAATAATTCCTGCCTCCGGCCCCCCTTACATCTT‑TTCGCCCGGGAAAC‑TCCACAG‑CATC‑ACC‑CCCCGCGCAACACACTCCTCCGGCATTGGgtcct                                                                                                                                                                                            >  2:900750/1‑237 (MQ=255)

                                 gTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGGTCATCTCTTTTTTTCCACAAATAATTCCTGCCTGCGGCCCGCCTTACAACTT‑CTTGCCCGCGCAAC‑TCCCCAC‑GCTT‑GAG‑CCCCCCGCCATCCACTCCTCCGGCCTTGGCTCttt                                                                                                                                                                                           >  2:1174002/1‑239 (MQ=255)

                                    ggTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCCACAATAATTCCTTCCTGCGGCCCGCCCACCAACTT‑GTGCCCGGCGCAAC‑ACACCCG‑TCTG‑AGG‑CCCGCCGCAATCAACTCCTTCGGTAGTGGctttctct                                                                                                                                                                                         >  2:1057581/1‑232 (MQ=255)

                                      tacgtttcccgccagggctcgaGGTTGA‑AAATAAAGCCGAGAGGTGGGTCCCAGGGCCGGGAGGGTCATTTTATCCCGGCTAACTGCCCGCTACG‑GGGGGGGGCTTTGCGGCGGGG‑CATTCATCTCTTTTTTACCAAAAATAATTCCGGCCTGCGCCCCGCCTAACAACTT‑TTCGCGGGAGAACC‑CCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGc                                                                                                                                                                                      <  2:1179664/217‑1 (MQ=255)

                                        tttttcCACCAAGGGCTTGGGGGTTT‑AAATAAGGCCGCAAGGCGGCGTCCAGCGCTGGGGGGTTCTTTTTTCCCGCGCGTAATCCACGCTTCA‑GCGGGTGGCTTT‑CCGGCGGGCCGCTCATCTTTTTTTTCCCAACAATAATTCCTGCCTGCGCCCCGCCTAACAACTT‑GCCGCCGGAGAACC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCaaa                                                                                                                                                                                   <  2:305989/238‑1 (MQ=255)

                                        aTTTTTCCCCCGGGGCTTCGGGGTTT‑AAAAAAAGCCGCAACGTGGCTTCCGGCGCCGCGGGGTTCATTTTCCCCGCGCTTAATGCCCGCTTCA‑GCGGGTGGCTTTGCGGGCGGG‑CGATCATCTTTTTTTTCCAAAAAATACTCCCTGCCTGCGGCCCGCCTAACACCTT‑GCCGCCGGAAAACC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCaaa                                                                                                                                                                                   <  2:502082/239‑1 (MQ=255)

                                            ttCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGGTCATCTCTTTTTTTACAACAAAAATTCCTGTCTGCGGGCCGCCTAACACATT‑TTCTCCGGGCAAAC‑TCCCCAG‑ACTG‑CCG‑CCCGCCGCCAACAAATCCCGCGGCATTGTCTCTTT‑CGGCAAAACg                                                                                                                                                                                >  2:905426/1‑238 (MQ=255)

                                            ttCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGCTCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGGCCGCCTAACCACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCGTTGGCTCTTT‑CGGCAAACCGt                                                                                                                                                                               >  1:21135/1‑239 (MQ=255)

                                             tCACTCAGAGCTTCAGGGTAT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCTTTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCGGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAACC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTa                                                                                                                                                                              <  1:764069/239‑1 (MQ=255)

                                             acccccAGAGTTTAGGGGTGA‑AAAAAAAGACGCAACGGGGCTCCCAGCGCGGGGAGGGTCATTTTATCCGCGGGTAATGCACGCTTCA‑GCGGGGGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAAAAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTa                                                                                                                                                                              <  2:1029858/234‑1 (MQ=255)

                                              gcccgcagagtgtcggtgttaaaaaaaaGGCGCATCGTGGTTCCCAGCGGCGCGGGGGTCATTTTAGCCGCGCGTAATGCCCGCTCCA‑GCGGGTGGCTTTGGAGGCGGG‑CGATCATCTCTTTTTCACCAACAATAATCCCGGCCTGGGGCCCGCCAAACAACTT‑GCCGCCGGAGAACC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTa                                                                                                                                                                              <  2:960411/218‑1 (MQ=255)

                                              ccacgcggagtatggggTGTAAAATAAGGCCGCAGCGGGCCTCCCGGCGCGGCGGGGGTTTTTTTCTCCGCGGTTAATTCTCGCTCGA‑GGGGGGGGCTTTGCGGGGGGG‑CGATCACCTCTTTTTCCCAAAAAAAAATCCCTCCCGGCGGCCCGCCTAACATCTT‑GTCGCGGGAGAATC‑CCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTa                                                                                                                                                                              <  2:489372/225‑1 (MQ=255)

                                              cACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGGTTATCTCTTTTTTACCACAAATAATTTCTGCCTGGCGCACGGCTAACCACTT‑TTTCCCGGAGAAAC‑CCACCCG‑GCTG‑AAG‑CCCGCCGCAATCAACTTCAGCGGCATCGGCTCTTTCCGCCAAAccaa                                                                                                                                                                               >  2:327966/1‑236 (MQ=255)

                                                          cAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCCTCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCCCTTT‑CGGCAAACCGTACTGGCTAAGATc                                                                                                                                                                  >  1:364016/1‑238 (MQ=255)

                                                             ggTAT‑AAAAAAACACGCAACGGTGCCTCCGGCGCGGGGGGGGATATTTTACCCGCGCGTAATGCACGCCGCA‑GCGGGTGGCTTTGCGGGGGGG‑CGATCATCTCTTTTTCACCAAAAATAATCCCTGCCTGCGCCCCGCCTAACAACTT‑GTCGCGGGAGAACC‑CCCCCAG‑CCTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGAc                                                                                                                                                               <  2:278515/238‑1 (MQ=255)

                                                               tgt‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAA‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCCTT‑CGGCAAACCGTACTGGCTAAGATCGcccc                                                                                                                                                             >  1:272990/1‑237 (MQ=255)

                                                                gt‑AAATAAGGCCGCATCGTGCCTTCCTGGGCCGCGGGGTTATTTTTACCCGCGCGTAATGCGCGCTGCA‑GCGGGTGGTTTTGCGGGCGGG‑CTTTCATCTTTTTTTTACAAACAATACTCCCCCCCTGCGCCCCCCCAAACAACTT‑GTCGCCGGAGAACC‑CCCCCAG‑CCGG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAg                                                                                                                                                            <  2:1027718/238‑1 (MQ=255)

                                                                   aaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGGTCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCCGAGAAAC‑TCACCAG‑ACTG‑ACC‑CCCGCCCCCATCCACTCCCGCGGCATTGGCTCTTT‑TCGCAAACCGCACTGGCTAAGATccaccaacca                                                                                                                                                         >  2:1098570/1‑235 (MQ=255)

                                                                            cGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCa                                                                                                                                                                                                                     >  1:425278/1‑171 (MQ=255)

                                                                            cGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGAACGCTTCA‑GCGGGTGGCTTTGCCGGCGGG‑CGATCATCTCTTTTTTTTCAACAACAATTCCTCCCTGTGGCCCGCCTAACAACTT‑TTCGCCAGAGACAC‑CCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAATTCCTGCGCCAATGGCTCCTT‑CGGCCAACAGTCATGGCTagaatccccccacgctccaccgg                                                                                                                                                 >  2:8894/1‑216 (MQ=255)

                                                                              cAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCAGCGAGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑AGATCATCTCTTTTTTACCCACAATAATTCCTGCCTGCGGACCGCCTAACAACTT‑GTCTCCGGAGAAAC‑TCACAAG‑ACTG‑ACG‑CCCGCCGCCATCAACTCCTGCGGCCTTTGCGCTTT‑AGGCATACCACTATGGCTAACAACGACCAGCGAGCCCATGCCtt                                                                                                                                              >  2:160737/1‑238 (MQ=255)

                                                                                      cgcccggcggcgtgggtgTCATTTTACCCCGGCGTAATGCACGCTGCA‑GGGGGTGCTTGGGGAGCGGGG‑CATTCATCTTTTTTTCCCAAAAAAAAATCCCGGCCGGCGGCCCGCCTAACAACTT‑GCCGCCGGAGAACC‑TCCCCAG‑CCTG‑ACC‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                                       <  2:909890/221‑1 (MQ=255)

                                                                                        tctaggcgcgccggggTCATTTTTGCCACGAGTAATGCACCACTCA‑GCGGGTGGCTTTGCAGGCGGG‑CACTCCTCTTTTTTTAACCAAAAATAACCCCGCCCGGCGGCCCGCCAAACAACGT‑GGCGCGGGAGAACC‑TCCCCAG‑AC‑G‑ACGCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCAc                                                                                                                                    <  2:652008/226‑1 (MQ=255)

                                                                                               ccgggggggtcttttctcccGCGCTTAAT‑CCCGCCTCACGCGGGTGCTTTTGCGGCGGGG‑CGTTCTCTTTTTTTTTCCCAACAATATTCCCGCCCTGCGCCCCGCCTAACAACTT‑GTCGCGGGAGAACC‑TCACCAG‑CCGG‑CCC‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                             <  2:977269/221‑1 (MQ=255)

                                                                                                  cgAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGCTCCTCCCTTTTTTACCCACAATAATTCCTCCCTGCGGCCCGCCTAAAAACTT‑TTCGCCGGAGAAAC‑TCACCAA‑ACTG‑CCG‑CCCGCCCCAATCAATCCCTGCCGCTTATACCCTTT‑TGGCGCACCCGACCGGCTACAATCTGCCAGCCTGCCACTTCCTCAAttcgtcccaccaggaaca                                                                                                                          >  2:694710/1‑222 (MQ=255)

                                                                                                               ctCCGCGTGTAACTCTCGCTTCG‑GGGGGTGCCTTGCCGGGCGGC‑CACTCACCTCTTTTGTACCAACAAGAATCCCCGCCCGCGCCCCCCACAACACCTT‑CTCGCCGGAAAACC‑CCCCCAG‑ACTG‑CCG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGc                                                                                                             <  2:1150263/238‑1 (MQ=255)

                                                                                                               aTCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCCCTACGGGAACATCCAGCTCTTTg                                                                                                              >  1:333412/1‑238 (MQ=255)

                                                                                                                    cgcgTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑Cccgccg                                                                                                                                                                                                                       >  2:883705/1‑129 (MQ=255)

                                                                                                                    cgcgTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCCGAGAAAC‑TTACCAC‑ACTG‑ACG‑CCCGCCGCAAACAACTCCTTCGGCATTTGCTCTTT‑CGGCAAACCCTACTGGCTCAAATCCACCAGCGAGCCACCGCCTAAATCCATCACTAAGGGCACATCCAGCTCTTTGCCgcg                                                                                                         >  2:154336/1‑236 (MQ=255)

                                                                                                                    cgcgTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GCCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑Cccgccg                                                                                                                                                                                                                       <  1:883705/129‑1 (MQ=255)

                                                                                                                    cgcgTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCCCTTTTTTTCCAACAATAATTCCTGCCTTCGGCCCGCCTAACCACTT‑GTCCCCGGAGAAAT‑TCCCCAG‑ACTG‑ACG‑CCCCCCCCCATCAACCCCTCCGGCAACGGCCCTTT‑CGGCAAAACCCACTTGCTAAGATTGACCCGCGTGCGCCTCCCTACAACAGGCCCCACGGGGACATCCCGCCCTTTGCCCAg                                                                                                         >  2:644311/1‑238 (MQ=255)

                                                                                                                      cgTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGGTCATCTCTTTTTTACCAACAATAATTCCTTCCTGCGGCCCGCCTAACAACTT‑TTTGCCGGAGAAAA‑TCACCAG‑CCTG‑ACG‑CCCGCCCCCATCAACTCCTTCGGCCATGGCCCTTT‑CGGCAAACCGGACTGGCTAACATCGACCAGCGAGCCACTTCCCAAATCACTCCACACGGGGACATCCACCTCCTTGccccccc                                                                                                       >  2:410863/1‑233 (MQ=255)

                                                                                                                         tacgccacgcgacctctggTGTCTTTGCGGGGGGG‑CTCTCCTCTGTTGTTTCCAAACAATAA‑TCCCGCCCTCCGCCCCCCCAACAAACT‑GTCCCCGG‑GAAACACCCCCAG‑GCCGCACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGcca                                                                                                    <  2:726647/222‑1 (MQ=255)

                                                                                                                                  ttCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAAt                                                                                                                                           >  1:233750/1‑190 (MQ=255)

                                                                                                                                   acA‑GCGGGTGGTTGTGCAGGCGGG‑CGATCATCTCTTTTTAACAAAAAAAAATCCCTGCCGGCGGCCCGCCTAACAACTT‑GCCGCCGGAGAAAC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                                         <  1:405956/238‑1 (MQ=255)

                                                                                                                                        cGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑TTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCCCCGCAATCAACTCCTGCGGCCTTGGCGCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCCTAATAAGTCCCTACGGGGACATCCCGCTCTTTTCCGACCGCCCCCAGTTCCGCtttc                                                                                      >  2:1174111/1‑236 (MQ=255)

                                                                                                                                           gTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAA‑TCACCAG‑ACTG‑ACG‑CCCGCCGAAATCAACTCCTGCGGCATTGGCTCTTT‑CCGCAAACCGTACTGGCTAAAATCCACCAGCGAGCCACTGCCCTAATCAGTCACTCCGGGAACATCCAGCTCCTTGCCCGCCGCCCCCCGTTCCCCTTTATCTa                                                                                  >  2:788547/1‑239 (MQ=255)

                                                                                                                                             ggCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCCGCGAGCCACTGCCTAAATCAGTTACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCCGTTCCCCTTCATTTAt                                                                                 >  1:681652/1‑238 (MQ=255)

                                                                                                                                             ggCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCCCCAGTTCCGCTTCATCTAt                                                                                 >  1:667153/1‑238 (MQ=255)

                                                                                                                                             ggCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATc                                                                                >  1:974549/1‑239 (MQ=255)

                                                                                                                                                    cAGGCGGG‑CGATCATCTCTTTTTTCCCAACAATAATTCCTGCCCGCGCCCCGCCTAACAACTT‑GCCGCCGGAAAAAC‑TCCCCAG‑ACGG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAg                                                                                                                                                                     <  2:131538/147‑1 (MQ=255)

                                                                                                                                                       gcggg‑cgATCATCTTTTTTTACCCAACAATAATTCCTGCCTGCGGCCCGCCTACCAACTT‑GTCGCGGGAGAACC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                       <  1:1174111/238‑1 (MQ=255)

                                                                                                                                                       ccggg‑cgATCTTCTTTTTTTTACCACCAATAATTCCTCCCTGCGCCCCCCCTACCACCTT‑GTCGCCGGAGAACC‑TCCCCAG‑CCTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                       <  1:674227/237‑1 (MQ=255)

                                                                                                                                                         ggg‑cgATCATCTCTTTTTTACAAACAATAATCCCTGCCTGCGGCCCGCCTAACAACTT‑GCCGGCGGAGAAAC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                     <  1:941188/238‑1 (MQ=255)

                                                                                                                                                             cgATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGGACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAAAATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAAccc                                                                  >  1:202393/1‑238 (MQ=255)

                                                                                                                                                               ctcatcTCTTTTTTACAAACAATAATCCCTGCCCGCGCCCCGCCTAACAACTT‑GTCGCCGGAGAACC‑TCACCAG‑ACGG‑ACC‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                                <  1:644311/237‑1 (MQ=255)

                                                                                                                                                               atcatcTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGc                                                                                                                                                           >  1:131538/1‑147 (MQ=255)

                                                                                                                                                                   tctcTTTTTTACCAACAAAAATCCCTGCCTGCGGCCCGTCTAACAATTT‑GCCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATg                                                           <  1:410863/239‑1 (MQ=255)

                                                                                                                                                                     tgttttttACCAACAATAGTTCCTGCGTGCGCCCCGCTAAACAACTT‑GCCGCCGGAGAACC‑CCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                         <  1:1084389/237‑1 (MQ=255)

                                                                                                                                                                     tctTTTTTACCAACAATAATTCCCGCCTGCGCCCCGCCTAACACCTG‑GTCGCCGGAAAAAC‑TCCCCAG‑ACGG‑CCG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                         <  1:152371/239‑1 (MQ=255)

                                                                                                                                                                     ctttttttACAAACATTACTCCCTCCCCGCGCCCCGCCTAACAACTT‑GCCGCCGGAAAACC‑CCCCCAG‑ACGG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                         <  1:444212/237‑1 (MQ=255)

                                                                                                                                                                     cctgtagtACAAACAATACTCCCACCCTGCGGCCCGCCT‑ACAACCT‑CGGCCCGGGAAAAC‑CCCCCAGTACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGCCCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                         <  2:272990/232‑1 (MQ=255)

                                                                                                                                                                      ctTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTCAATCAGCCACTACGGTAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAAt                                                            >  2:1083922/1‑235 (MQ=255)

                                                                                                                                                                      ctTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAAt                                                            >  2:1009469/1‑235 (MQ=255)

                                                                                                                                                                          tttACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAACGCTGTAGt                                                    >  1:378630/1‑239 (MQ=255)

                                                                                                                                                                            tcccAACAATATTCCCCCCCTGCGGCCCGCCAAACAACTT‑GGCGCGGAAGAATC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                  <  1:469717/237‑1 (MQ=255)

                                                                                                                                                                                 acaaTAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAACATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCACCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGTTGAACCCCTGACTGCTGTAGTTACTGGt                                             >  2:637232/1‑239 (MQ=255)

                                                                                                                                                                                   aaaaatCCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑CCC‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTAt                                           <  1:245273/236‑1 (MQ=255)

                                                                                                                                                                                      actccctgcctcCGCCCCCCCAAACACCTG‑TTCCCCGGAGAAAC‑CCCCCAG‑CTGG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                        <  1:618753/235‑1 (MQ=255)

                                                                                                                                                                                      acccccgccctgc‑GCCCCCCCAACCACTTGGTCCCCGGAAAATC‑TCCCCAC‑GCTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                        <  2:97368/231‑1 (MQ=255)

                                                                                                                                                                                      acatgagtaaaagagccacccTACCACCTG‑GCCGCGGGAAAACC‑CCCCCAG‑ACGG‑ACC‑CCCGCCGCAATCAACCCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                        <  2:605714/220‑1 (MQ=255)

                                                                                                                                                                                      aatccctgctgtcgccccGCCAAACAACGT‑GCCGCCGGAGAAAC‑TCCCCAC‑ACGT‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                        <  1:694710/224‑1 (MQ=255)

                                                                                                                                                                                      aatCCCTGCCGGCGCCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACGG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                        <  1:1174002/239‑1 (MQ=255)

                                                                                                                                                                                      aatCCCCGCCTGCGGCCCGCCTAACAACTT‑GCCGCGGGAAAACC‑CCCCCAG‑ACTA‑CCG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                        <  1:154336/239‑1 (MQ=255)

                                                                                                                                                                                             gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                     >  2:987507/1‑203 (MQ=255)

                                                                                                                                                                                             gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                     <  1:987507/203‑1 (MQ=255)

                                                                                                                                                                                             gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                  >  2:290151/1‑238 (MQ=255)

                                                                                                                                                                                             gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                  >  1:952470/1‑238 (MQ=255)

                                                                                                                                                                                             gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                  >  1:1016477/1‑238 (MQ=255)

                                                                                                                                                                                                tgcGCCCCCCCAAACAACTT‑GCCGCGGGACAACC‑TCCCCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                              <  1:30437/239‑1 (MQ=255)

                                                                                                                                                                                                      ccGCCTAAAAACTT‑GTCCCGGGAAAAAC‑GCACAGC‑CCTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTC‑TT‑CGGCAAACCGTACTGGCTAAGACCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCg                       <  1:584192/239‑1 (MQ=255)

                                                                                                                                                                                                                     gTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCGAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGAATGTACTTTCATCaa                              >  1:605714/1‑219 (MQ=255)

                                                                                                                                                                                                                     gTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc           >  2:60411/1‑238 (MQ=255)

                                                                                                                                                                                                                     gTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc           >  2:1075598/1‑238 (MQ=255)

                                                                                                                                                                                                                     gTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc           >  1:88549/1‑238 (MQ=255)

                                                                                                                                                                                                                     gTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc           >  1:838445/1‑238 (MQ=255)

                                                                                                                                                                                                                          cGGAGAACC‑TCACCAG‑CCTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGcc     <  1:1098570/239‑1 (MQ=255)

                                                                                                                                                                                                                              gaAAC‑TCCCCAG‑ACTG‑CCG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGa  <  2:944497/238‑1 (MQ=255)

ATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGT‑AAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCA‑GCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAC‑TCACCAG‑ACTG‑ACG‑CCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGA  >  NC_000913/3760158‑3760613

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