BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,383	0	A	C	31.5%	143.2 / 33.9	112	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base A (54/21); new base C (2/35); total (56/56)
Fisher's exact test for biased strand distribution p-value = 7.33e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCT  >  NC_000913/3760146‑3760611

gCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCAc                                                                                                                                                                                                                                     <  1:1243963/239‑1 (MQ=255)

     gATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAAATCACCAg                                                                                                                                                                                                                                  >  1:975243/1‑237 (MQ=255)

     gATTGGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGCCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGa                                                                                                                                                                                                                                 <  2:51363/238‑1 (MQ=255)

                ccTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGAcgcccgc                                                                                                                                                                                                                      >  1:748027/1‑238 (MQ=255)

                ccTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCCAACAACTTGTCGCCGGAGAAAATCACCAGACTGAcgcccgc                                                                                                                                                                                                                      >  1:348646/1‑238 (MQ=255)

                ccTCGGCGCTGGGGGAAAGCAGGAGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAAAAAAGAC‑GCAACGTG‑GCTCCCAGCGCCGCGAGGGTCTTTTTATCCGCGCG‑TAATGCACGCTTCAGGGGGTGGCTTTGCG‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCCCCAGACTGAcgcccgc                                                                                                                                                                                                                      <  2:1282224/238‑1 (MQ=255)

                     gcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCTCCGGAGAAACTCACCAGACTGCCGCCcgccgcca                                                                                                                                                                                                                 >  1:1221443/1‑236 (MQ=255)

                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GTTCCTCTCTTTTTTTCC‑ACCAATAATTCCTGCCTGCGGCCCGCCTAAAAACTCTTGCCCGTTGAAACCCACCAGACCGCCGCCCCCCCCCATCCAcc                                                                                                                                                                                                           >  2:782285/1‑237 (MQ=255)

                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GGTCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAAAAACTTGTCGCCGGAGAAACTCAACCAACTGACGCCCCCCCCCATCCACt                                                                                                                                                                                                           >  2:1201719/1‑238 (MQ=255)

                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGGGAAACTCACCCGCCTGACGCCCGCCGCAATCAACt                                                                                                                                                                                                           >  2:378501/1‑238 (MQ=255)

                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACt                                                                                                                                                                                                           >  1:942655/1‑238 (MQ=255)

                           cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCCTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACCACTTGTCTCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCCAcc                                                                                                                                                                                                           >  1:916855/1‑237 (MQ=255)

                                agacAGGCGCGGGGGCGGTATTTTTCCCCTCAGCGTTTCAGGGTTAAAATAAAGAC‑GCACGGTG‑GCTCCCAGCGCCGCGAGGGTCTTTTTACCCGCGCG‑TAATGCACGCTTCAGCGGTTGGTTTTGCG‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAACCTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGc                                                                                                                                                                                                      <  2:555900/235‑1 (MQ=255)

                                       cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTg                                                                                                                                                                                               >  1:1389973/1‑238 (MQ=255)

                                        gcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GGTCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGGGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTg                                                                                                                                                                                               >  1:294276/1‑237 (MQ=255)

                                                 gTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTTGCCGGAGAAACTCACCAGACCGCAGCCCGCCGCAATCAACCCCTGCGGCCTTGGCCCTTTCgg                                                                                                                                                                                     >  2:1106405/1‑238 (MQ=255)

                                                        ttCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCCTCTCTTTTTTACC‑AACAATAATTCCCGGCTGCGGCCCGCCTAACCACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCCACCCCTGGCGCAATGGCTCTTTCTGCAAAccc                                                                                                                                                                              >  2:90491/1‑237 (MQ=255)

                                                        ttCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCCTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCCCAATCCACTCCCGCCGCATTTGCTCTTTCGGCAAAcccg                                                                                                                                                                             >  1:459866/1‑237 (MQ=255)

                                                         tgtttataaggaacgggttgaaaaaaaaGCC‑GCAAGGTG‑GCTTCCAGCGCCGCGGGGGTCATTTTACCCGCGCG‑TAATGCACGCTTCAGCGGGGGGCTTTGCG‑GGCGGGC‑GATCATCTCTTTTTTACC‑AAAAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCGGGAGAAACCCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  2:1393047/218‑1 (MQ=255)

                                                         tccctaagagCTTCAGGGTGAAAATAAAGAC‑GCACGGTG‑CTTCCCAGCGCCGCGAGGGTCATTTTATCCGCGCT‑AATTGCACGCTTCAGGGGGGGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCGGCCTGCGGCCCGCCTAACACCTTGTCGCCGGAAAACCTCCCCAGACTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  2:377831/233‑1 (MQ=255)

                                                         gcccTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGGGGGTCATTTTATCCGCGCG‑TAATCCACCCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTCCC‑AACAATAATTCCGGCCGGCGGCCCGCCTAACAACTTGCCGCCGGAGAAACTCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  1:643646/236‑1 (MQ=255)

                                                         ccccTCAGAGCTTCAGGGTGAAAAAAAGGCC‑GCACGGTG‑GCTTCCGGCGCCGGGGGGGTCTTTTTCCCCGCGCG‑TAATGCCCGCTTCAGGGGGTGGCTTTGCG‑GGCGGGC‑AATCATCTCTTTTTTACA‑AAAAATAATCCCTGCCGGCGGCCCCCCTAACAACTTGTCGCGGGAGAAACCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  2:975243/236‑1 (MQ=255)

                                                         ccACTCAGAGCTGCGGGGTTAAAAAAAAGAC‑ACAAGGGG‑GCCTCCAGCGCCGCAGGGGTCATTTTATCCGCGCG‑TAATCCCCGCTCCAGGGGGTGGGTTTGCG‑GGGGGGC‑GATCATCTTTTTTTTACA‑AAAAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  2:1164621/238‑1 (MQ=255)

                                                         ccACCCAGAGCTTAGGGTTGAAAAAAAAGAC‑GCAACGTG‑GCTTCCGGCGCCGGGAGGGTCATTTTATCCGCGCG‑TAATGCACCCTTCGGCGGGTGGCTTGGCA‑GGCGGGC‑GATCATCTCTTTTTCACA‑AACAATAATCCCTTCCTGCGGCCCGCCTAACAACTTGTCGCCGGAAAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  2:670485/238‑1 (MQ=255)

                                                          cACCCAGAGCTTCAAGGTTGAAAAAAAGGCAGCACCGTG‑GC‑TCCGGCGGCGCGGGGGGCAGTTTTACCCCGCGTTAATGCCCGCTCAACGGGGTGGCTTTGC‑‑GGCGGGCGGATCATCT‑TTTTTTACCAAAAAATAATCCCTTCCGGCGGCCCGCCTAACAACTTGTCGCGGGAGAACCCCCCCAGACTGAGGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  2:14403/239‑1 (MQ=255)

                                                             tCAGGGCCTGAGGGTGTAAAAAAAGAC‑GCACCGTG‑TCTCCCGGCGCGGGGGGGGTCTTTTTTCCCGCCGG‑TAATGCAGGCTCCAGGGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATCCCTGCCGGCGGCCCGCCTAACAACTTGTCCCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTgg                                                                                                                                                                        <  2:1386473/239‑1 (MQ=255)

                                                                  ggttCGGGGTGTAAAAAAAGAC‑GCCACGGGCGCGTCGGGCCCCGGGGGTGTATTTTTCGCCGCGGG‑TAATGCACGTTCCAGCGGGGGGCTTGGCA‑GGGGGGC‑GACCATCTCTTTTTTACC‑AACAATAATCCCTGCCTGGGCCCCCCCTAACAACTTGTCGCCGAAAAAACTCACCAGACTACCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTa                                                                                                                                                                     <  2:782924/236‑1 (MQ=255)

                                                                           tgtAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGGCAACCGTACTTGCTAAGATCCAccc                                                                                                                                                           >  2:401362/1‑237 (MQ=255)

                                                                           tgtAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGAccc                                                                                                                                                           >  1:563564/1‑237 (MQ=255)

                                                                             taaataaaGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCCCCAGACTGACGCCCGCCGCAATCCACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGc                                                                                                                                                         >  1:532573/1‑238 (MQ=255)

                                                                                   aaaacccataCGTG‑CCTTCCAGGCCCGGCTGGTTCTTTTTATCCGCGCG‑TTAGCCACCTTTCAGCGGTCGGTTTGCCG‑GGCGGGC‑AATCTCCTCTTTTTCACC‑ATAAATATTCCCGCCCTGCGGCCCCCTTACCACCTGGCCGCCGGATAACCCCCCCAGCCTACCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                   <  2:1335494/230‑1 (MQ=255)

                                                                                   aaGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCCGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCa                                                                                                                                                   >  1:1296038/1‑238 (MQ=255)

                                                                                    aGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAgc                                                                                                                                                                  >  1:1393047/1‑221 (MQ=255)

                                                                                     gAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GGTCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCATCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGGCCTTTCGGCCAACCCGACCGGCTAAAATCGACCAGCCAGCCACt                                                                                                                                                 >  1:479138/1‑238 (MQ=255)

                                                                                     gAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCCACCAGCGAGCCACt                                                                                                                                                 >  1:1382938/1‑238 (MQ=255)

                                                                                     gAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCCCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACt                                                                                                                                                 >  1:299815/1‑238 (MQ=255)

                                                                                                  gCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAAATCACCAGACTGCCGCCCGCCGCCATCACCTCCTGCGGCATTGGCCCTTTCGGCGAACCGTACAGGCTAAGCTCGCCCCGGGAGCCACGTCCTAAATCAg                                                                                                                                      >  2:801573/1‑238 (MQ=255)

                                                                                                   cTTCAAGCGCCGCGAGGTTCATTTTCCCCGGGCG‑TAATGCCCGCTTCGGGGGTTGGCTTGGCA‑GGCGGGC‑GATCATTTCTTTTATACC‑AAAAATAATTCCTGCCTGCGGCCCCCCTACCATCTTCTCGCGGGAGAAACCCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGt                                                                                                                                     <  2:449540/238‑1 (MQ=255)

                                                                                                    ttCCAGCCCCGCGGGGGTCATTTTACCCGCGCG‑TAATGCCCGCGTCAGCGGGTGGCTTTGCG‑GGCGGGC‑GATCATCTCTTTTTCACC‑AAAAATAATCCCTGCCGGCGCCCCCCCTAACAACTTGCCGCCGGAAAACCTCACCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTc                                                                                                                                    <  2:1389973/238‑1 (MQ=255)

                                                                                                     tCCAGCGCCGGGTGGTTCATTTTATCCGCCCG‑TAATGCCCGTTCCGGGGGGGGTCTTTGCG‑GGCGGGG‑TATCATCTCTTTTTTACC‑AACAATAATCCCTGCCTGCGGCCCGCCTAACACCTTGTCGCCGGAGAAACCCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCAc                                                                                                                                  <  2:942655/239‑1 (MQ=255)

                                                                                                      ccAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            >  2:1222271/1‑196 (MQ=255)

                                                                                                      ccAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTa                                                                                                                                                                            <  1:1222271/196‑1 (MQ=255)

                                                                                                             cgcgGGGTTCATTTTACCCGCGCG‑TAATGCGCGCTTCAG‑GGGTTGTCTTGCAGGGCGGGC‑GACCACCTCTTTTTTCCC‑AACAATATCTCCTCCCTGCGGCCCGCCAAACACCTCGCCGCCGGAAAAACTCCCCAGACGGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGCCACTACGGGa                                                                                                                           <  2:539177/238‑1 (MQ=255)

                                                                                                             cgcgGGGGACATTTTGCCCGGGCC‑TAAGGCCCTCGTCGGGGGGGGGTTTTGCG‑GGCGGGC‑GATCATCTCTTGTTTACA‑AACAATAATCCCCGCCCGCGCCCCGCCAAACAACTTGTCGCCGGAGAACCCCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGa                                                                                                                           <  2:294276/238‑1 (MQ=255)

                                                                                                              gcgAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCCACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAAc                                                                                                                         >  1:109222/1‑239 (MQ=255)

                                                                                                               cgAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GGTCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACa                                                                                                                        >  2:159665/1‑239 (MQ=255)

                                                                                                                       aTTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCCTACTgg                                                                                                                                                                        >  1:185648/1‑183 (MQ=255)

                                                                                                                       aTTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCCGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCCTTGGCTCTTTCGGCAAACCGTACTGGCCAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACCCAGGGAACATCCAGc                                                                                                                  >  2:340676/1‑237 (MQ=255)

                                                                                                                             tCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCCCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCCAACCGGTCTGGCTCAGAACGACCAGCGAGCCCCTTCCTAAATTCATCACTACCGGGACATcccccccttttc                                                                                                           >  2:1332683/1‑228 (MQ=255)

                                                                                                                              ccGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAAAAACCTGTCGCCCGGGAAACTCACCAGACTGACCCCCGCCGCAATCCACTCCTGCGGCATTGGACCCTTCGGCCAACCGTACCGGCTAACGTCGAACAGGGAGACCTCGCCTTAAACAATCACACAGGGGACATCCAGCGCCTTGCCg                                                                                                         >  2:211277/1‑239 (MQ=255)

                                                                                                                                 cgcg‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GGTAATCTCTTTTTTTCC‑AACAATAATTCCTGCCTGCGGCCCGCCAAACAACTTGTCGCCGGTGAAAATCACCAGGCTGGCCCCCGCCGCAATCAACTCCTGCGGCATTTGCCCCTTCGGCAAACCGGACTGGCTAAAAACGACCAGCGGGCCACCGCCCAAATTAGACACTACGGGAACAGccgacccttttccatac                                                                                                      >  2:511187/1‑224 (MQ=255)

                                                                                                                                                       ggtggtTTGCCA‑GGGGGGC‑GTTCATCCCTTTTTCACC‑AAAAAAATTTCCTGCCGGCGGCCCGCCTAACAACTTGCCGCCGGAAAACCTCCCCAGACTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATc                                                                                  <  2:1382938/238‑1 (MQ=255)

                                                                                                                                                        gTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCCGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTa                                                                                >  1:455911/1‑239 (MQ=255)

                                                                                                                                                          ggCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑CACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTTTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGGACTGGCTAAGATCGACCAGCGAGCCACTGCCCAAATCAGTCCCCACGGGGAAATCCAGCTCTTTTCCGAGCGCCACCAGTTCCCCTTCCCCTAt                                                                               >  2:613453/1‑238 (MQ=255)

                                                                                                                                                          ggCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCCGCGAGCCACTGCCCAAATCAGTCACTACGGGAACATCCAGCCCTTTGCCGCGCGCCCCCAGGTCCGccacaactaa                                                                               >  2:422004/1‑229 (MQ=255)

                                                                                                                                                          ggCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTAt                                                                               >  1:318392/1‑238 (MQ=255)

                                                                                                                                                                tggcgggggggcactcatcTTTTTTGTACC‑AACAATAGTTCCTGCCGGCGGCCCCCCTACCACTGTGGCGCGGGAGAACCCCCCCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCtt                                                                          <  2:300329/225‑1 (MQ=255)

                                                                                                                                                                     gcgggc‑gATCATCTTTTTTTTCCC‑AACAATAATCCCGGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                     <  1:243229/238‑1 (MQ=255)

                                                                                                                                                                      cgggc‑gATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGa                                                                    >  1:1198136/1‑238 (MQ=255)

                                                                                                                                                                        ggtcgtttatttctTTTTTACC‑AACAATAGTTCCTGCCGGCCGCCCGCCTAACACCTTGCCCCCGGAAAAACGCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                   <  1:90491/227‑1 (MQ=255)

                                                                                                                                                                        ggggattcctcttTTTTTCACC‑AACAAAAATCCCTGCCCGCGCCCCCCCTAACATCTTGCCGCGGGAAACACTCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                   <  2:1256886/229‑1 (MQ=255)

                                                                                                                                                                            gaccatcTCTTTTTACCC‑AAAAAAAATCCCTCCCGGCGCCCCGCCTAACACCTTTCCGCCGAAGAACCCCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                              <  2:563564/236‑1 (MQ=255)

                                                                                                                                                                            gATCATCTCTTTTTCACC‑AAAAATATTCCCGCCCTGCGCCCCGCCAAACAACTTGTCGCGGGAAAACCCCCCCAACCGCACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                              <  2:455911/239‑1 (MQ=255)

                                                                                                                                                                              caatcTTTTTTTTACC‑AACAATAATCCCTGCCTGCGGCCCGCCTAACAACTTGCCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGa                                                            <  1:401362/237‑1 (MQ=255)

                                                                                                                                                                                 tcttTTTTTAACC‑AACAAAATTCCCTGCCGGCGCCCCGCCTAACAACTTGCCGCGGGAGAAACCCCCCAAACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATg                                                         <  1:1106405/239‑1 (MQ=255)

                                                                                                                                                                                     tttttaACC‑AACAATAATCCCCGCCTGCGCCCCGCCAAACACTGTGCCCCCGGAGACCCTCCCCAGCCTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGt                                                     <  1:511187/239‑1 (MQ=255)

                                                                                                                                                                                         ttcccaaaaaaaactccctgccgcCGCCCCGCCTAACACCGTGCCGCCGGAAAACCCCCCCAGCCGCACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGt                                                  <  2:318392/223‑1 (MQ=255)

                                                                                                                                                                                          tccc‑AACAATAATTCCTGCCCGCGCCCCTCCTAACACCTTGCCGCCGGAGAAACTCACCACACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                <  1:1201719/237‑1 (MQ=255)

                                                                                                                                                                                          tccc‑AAAAAAAATCCCCGCCTGCGCCCCGCCTAACACCTTGTCGCCGGAGAAACCCCCCAGACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                <  1:378501/237‑1 (MQ=255)

                                                                                                                                                                                          caCC‑AACAATAATCCCTGCCGGCGGCCCCACTAACACCTTGTCCGCGGAAAACCTCCCCAGACCCACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                <  2:532573/238‑1 (MQ=255)

                                                                                                                                                                                           tccaaaaaatattCCCTGCCGGCCGCCCGCCCACAAACTGGCCGCCGAAAAAACCACCCAAACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                <  2:348646/232‑1 (MQ=255)

                                                                                                                                                                                           caccaaaaatatccccgccccgcGGCCCGCCAAACAACTTGTCCCCGGAGAAACCCCCCAGACGGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTCTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                <  2:144225/218‑1 (MQ=255)

                                                                                                                                                                                            gccataaataagTCCTGCCTGCGGCCCGCCTAACAACTTCTCCCCGGAAAAACTCCCCATGCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTAc                                               <  1:1020794/233‑1 (MQ=255)

                                                                                                                                                                                                acaataccccctgcctgCGCCCCGCCGAACACCTTGTCGCCGGAAAACCTCACCAGACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                           <  1:497996/230‑1 (MQ=255)

                                                                                                                                                                                                acaaTACTTCCACCCTGCGGCCCGCCGAACAACTAGTCGCCGGAGAAACGCACCAGACTGCGGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                           <  1:771967/239‑1 (MQ=255)

                                                                                                                                                                                                acaaTAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGt                                           >  1:393818/1‑239 (MQ=255)

                                                                                                                                                                                                     actccccgcctgcGGCCCGCCTAACATCTTGCCGCGGGAAAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  2:1198136/232‑1 (MQ=255)

                                                                                                                                                                                                     aatTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAACCTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  1:1256449/239‑1 (MQ=255)

                                                                                                                                                                                                     aatTCCGGCCTGCGGCCCGCATAACAACTCGACGCCGGAAAAACTACCCAGACTCACGCCCGCCGCAATCAACTCCTTCGGCATTGGCTCGTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  1:176962/239‑1 (MQ=255)

                                                                                                                                                                                                     aatCCCTCCCGGCGGCCCGCCTAACAATTTGCCGCCGGAGAACCGCACCAGACTGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  1:392242/239‑1 (MQ=255)

                                                                                                                                                                                                     aactCCTGCACGCGCCCCCCCTCACACCTTGCCCCCGGAAAATCCCCCCAGCCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  2:1411672/236‑1 (MQ=255)

                                                                                                                                                                                                     aacgcctgccccggccccGACTAAAACCTTCGCGGCGGAAAACCTCCCGACACTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                      <  2:916855/224‑1 (MQ=255)

                                                                                                                                                                                                       ttCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACt                                    >  1:454057/1‑239 (MQ=255)

                                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGc                                                                                                           >  1:593556/1‑163 (MQ=255)

                                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                <  1:918362/222‑1 (MQ=255)

                                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                >  2:918362/1‑222 (MQ=255)

                                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                >  1:913674/1‑238 (MQ=255)

                                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                >  1:849760/1‑238 (MQ=255)

                                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                >  1:153331/1‑238 (MQ=255)

                                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                >  1:674040/1‑238 (MQ=255)

                                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAACGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGAATGTACTTtca                                >  1:1318715/1‑238 (MQ=255)

                                                                                                                                                                                                                      ccccTAACACCTGTCCGCCGGAGAAACTCCCCAGACTAACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCg                     <  1:136345/237‑1 (MQ=255)

                                                                                                                                                                                                                            acaacTTGTCGCCGAAGAACCCCCCCAGGCCGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCATTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTAttt               <  1:175338/239‑1 (MQ=255)

                                                                                                                                                                                                                                   gTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc         >  1:447086/1‑238 (MQ=255)

                                                                                                                                                                                                                                   gTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc         >  2:317807/1‑238 (MQ=255)

                                                                                                                                                                                                                                      gCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCAcc      >  1:1343094/1‑238 (MQ=255)

                                                                                                                                                                                                                                      gCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGAAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattcaac      >  1:1320679/1‑236 (MQ=255)

                                                                                                                                                                                                                                        cGGAAAAACCCCCCAGCCTCCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGc    <  2:459866/238‑1 (MQ=255)

                                                                                                                                                                                                                                         ggAGAAACCCCCCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCt  <  1:991146/239‑1 (MQ=255)

                                                                                                                                                                                                                                         ggAAAAACTCCCAAGTCCGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCt  <  1:782285/239‑1 (MQ=255)

GCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGAC‑GCAACGTG‑GCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCG‑TAATGCACGCTTCAGCGGGTGGCTTTGCA‑GGCGGGC‑GATCATCTCTTTTTTACC‑AACAATAATTCCTGCCTGCGGCCCGCCTAACAACTTGTCGCCGGAGAAACTCACCAGACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTTCGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCT  >  NC_000913/3760146‑3760611

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