BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,610,860	G→T	14.1%	A53S (GCT→TCT)	hyfI →	hydrogenase 4, Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,610,860	0	G	T	14.1%	218.9 / 5.4	103	A53S (GCT→TCT)	hyfI	hydrogenase 4, Fe‑S subunit
Reads supporting (aligned to +/- strand): ref base G (25/61); new base T (14/2); total (39/64)
Fisher's exact test for biased strand distribution p-value = 1.64e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

AGCGTCTGCCCGGAATGCAAACAACGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGT  >  NC_000913/2610624‑2611098

agcgTCTGCCCGGAATGCAAACAACGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTtg‑c                                                                                                                                                                                                                                                     <  1:943545/238‑1 (MQ=255)

 gcgTCTGCCCGGAATGCAAACAACGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTtg‑ct                                                                                                                                                                                                                                                    <  1:819580/238‑1 (MQ=255)

       gcgcGGAATGCAAACAACGCGCGGCGCTTATCAACGACGATACAGATGTCATGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTa                                                                                                                                                                                                                                            <  1:1126907/236‑1 (MQ=255)

        cccGGAATGCAAACAACGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTtg‑ct                                                                                                                                                                                                                                                    <  1:545797/231‑1 (MQ=255)

          cGGAATGCAAACAACGCGCGACGCTTATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTacac                                                                                                                                                                                                                                         <  1:768371/239‑1 (MQ=255)

          cGGAATGCAAACAACGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTacac                                                                                                                                                                                                                                         <  1:133879/239‑1 (MQ=255)

          cGGAATGCAAACAACGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTacac                                                                                                                                                                                                                                         <  1:980052/239‑1 (MQ=255)

           ggatgtcAAACTACGCGCGACGGTGTCAAACGCAGATTCAGATGTTCTGCTGGGGGCTAAGGCGGAGCTATGAGTGCA‑GAGCCTACACAACATGTCAGCCAGCCCATCACGCTGGAAGAGCAAACGCAAAAGATGACGCGGAATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑TAAATTTTTG‑CTGCCA‑TTACAcc                                                                                                                                                                                                                                        <  2:810679/233‑1 (MQ=255)

             aaTGCAAACACCGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCa                                                                                                                                                                                                                                       <  1:26249/238‑1 (MQ=255)

                  aaaCAACGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCG‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCTGCGATTGCGGCGGCTGCAACGCCGCTGAAAGC‑GAAAATTTTT‑TTGCCC‑tttacaacctctct                                                                                                                                                                                                                                 >  2:545797/1‑227 (MQ=255)

                       acGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTTCAACGCCTGTGAAATC‑GGAATTTTTT‑CTGCCA‑TTACACCAGATTTC‑gcccc                                                                                                                                                                                                                           >  1:568709/1‑235 (MQ=255)

                       acGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCACAGCCTGGTGAATC‑GAAATTTTTT‑CTTCCC‑TTAAACCaagatttgcccc                                                                                                                                                                                                                            >  1:536992/1‑228 (MQ=255)

                         gcgcgACGGTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGTTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAg                                                                                                                                                                                                                        <  1:248528/239‑1 (MQ=255)

                                  tGATCAACGACGATACAGATGGACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGtttt                                                                                                                                                                                                               <  1:494544/238‑1 (MQ=255)

                                    aTCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCACGTCGATTGCGGCGGCGGCAACGCCCTGGAAATC‑GGAATTTTTG‑CTGCCC‑TTACACCAGCATTC‑CGC‑GCAGAAACGtttgga                                                                                                                                                                                                             >  1:1136271/1‑236 (MQ=255)

                                           ctacttcAGTTGTACTGTTGGTTGCTAAGGTGCCGCTATGAGTCCC‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATTTCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAg                                                                                                                                                                                                      <  2:489867/232‑1 (MQ=255)

                                            cgaTACAGGGTTTCTTCTGGTGGCTAAGGAGCAGCTTGGGGTCCA‑GTGCTTCCAAAAAATGGCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATACCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGgt                                                                                                                                                                                                    <  2:94358/239‑1 (MQ=255)

                                            cgaTACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGGGAAATC‑GAAATTTTTT‑CTGCCC‑TTACACCAGTATTC‑GAC‑ACAG‑AAAGTTTTGGGATTAAggg                                                                                                                                                                                                    >  1:507255/1‑238 (MQ=255)

                                            cgaTACAGATGCACGTCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGTAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGgt                                                                                                                                                                                                    <  1:612912/239‑1 (MQ=255)

                                            ccaTACAGATGTTCTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GGTCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGgt                                                                                                                                                                                                    <  1:524857/237‑1 (MQ=255)

                                            agaTACAGATGTACGGCTGGTGGGTAAGGAACAACTATGAGGCCA‑GTGCTTACACAACATGGCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGgt                                                                                                                                                                                                    <  1:600931/238‑1 (MQ=255)

                                             aattcacaTGTACTGCGGGTGGGTAAGGGGCAGCTTGGGATCCA‑TGGTTGAAACAACATGGCAGCCCGGCCATCCCGCTTGACAGGCAAACGCAAAAGATGACGCGGCCTTTGG‑TACAGGATATCCGTCGCTTGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGgt                                                                                                                                                                                                    <  2:82829/231‑1 (MQ=255)

                                              aTACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGTCGGCTGCAACGCCTGTGAAAAC‑GAAATTTTTG‑TCGCCC‑TTACACCCGCATTT‑GAC‑GCGG‑AAAGTTTTGGCATTAAGgttg                                                                                                                                                                                                  >  1:1373139/1‑239 (MQ=255)

                                              aTACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGGGAAAAC‑GAAAATTTTT‑GTGCCC‑TTACACCAGGAATT‑GAA‑GCAG‑aaaagtttgggagttgtcttt                                                                                                                                                                                                  >  2:806541/1‑220 (MQ=255)

                                              aTACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCCTGGGAATC‑GAAATTTTTG‑CTGCCC‑TTTCACCAGAATTC‑GAA‑CCCG‑AAAGGTTTGGCATTAAGgttt                                                                                                                                                                                                  >  1:528412/1‑238 (MQ=255)

                                              aTACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAAAGCCGTTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACCCCAAGATTCTCAC‑G‑AG‑AACCGTTTGGCCTTAAgggtt                                                                                                                                                                                                  >  1:822263/1‑236 (MQ=255)

                                                    aTGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTtg‑ct                                                                                                                                                                                                                                                    <  2:436718/187‑1 (MQ=255)

                                                    aTGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTtg‑ct                                                                                                                                                                                                                                                    >  1:436718/1‑187 (MQ=255)

                                                    aTGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGGCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCa                                                                                                                                                                                                                                       <  1:971703/199‑1 (MQ=255)

                                                       tACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGCCGCTTGCGGCGGGCGCAACGCCCGGGCAATCAGAATTTTTTT‑GTGCCC‑TT‑CACCAGTCTTT‑GGC‑GCGA‑AACGGTTTTTGATTATGGGTGTTtcaaca                                                                                                                                                                                          >  2:1092290/1‑238 (MQ=255)

                                                         ctgctgGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGGCACCCCCGTTAAATC‑GCAATTTTTT‑CTGCCA‑TTAAACCAGACTTC‑GAC‑CCAC‑AACGTTTTGTAGTTAAGGTTGTTtcaaaacc                                                                                                                                                                                        >  2:1173771/1‑233 (MQ=255)

                                                         ctgctgGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGAATTC‑GAC‑GCAG‑AACGCTTTGGCATTAAGGTTGTTtcaccctgt                                                                                                                                                                                       >  1:489867/1‑233 (MQ=255)

                                                         ctgctgGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGGGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GCA‑CCAG‑AACGGTTTGGCATTAAGGTTGTTTCATCAcc                                                                                                                                                                                        >  1:1167843/1‑238 (MQ=255)

                                                         ctgctgGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCCGTGCAATC‑CAAATTTTTG‑TTGCCC‑TTAAACCAGTATTT‑GCC‑GCCG‑GACGTTTTTGGATTAAGGTTGTTTCATCACcg                                                                                                                                                                                       >  1:359652/1‑239 (MQ=255)

                                                         ctgctgGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAAAGCCTGTGAAATC‑GGAATTTTTG‑CTTCCA‑TTACACCAGTATTT‑GAC‑GCCG‑AACGTTTTGGCATTAAGGGTGTTTCATCAcc                                                                                                                                                                                        >  1:36107/1‑238 (MQ=255)

                                                         ctgctgGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCCCGAACGCCCGCGAAATC‑GCAATTTTTG‑CTGCCC‑TTTCACCAGTATTC‑GAC‑GCCA‑AACACTTTGGGATTAAGGTTGTTTCATAAcct                                                                                                                                                                                       >  2:600931/1‑238 (MQ=255)

                                                              ggtggCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑TGAAATTTTG‑CTGCCA‑TTAAACCAGTATTC‑GCC‑GCAG‑AAAGTTTTGGCATTAAGGTTGTTTCTTAACCCCGTc                                                                                                                                                                                   >  1:1331072/1‑238 (MQ=255)

                                                              ggtggCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGCCGGCTGGAACGCCTGTGGAATT‑CAAATTTTTG‑CTGCCC‑TTACACCCCTATTC‑GAC‑GACG‑AACGTTTTTGCGTTAAGGTTCTTTTACCCCCGCGTc                                                                                                                                                                                   >  1:1268160/1‑238 (MQ=255)

                                                               gtggCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCAGCAAACCCTCTGAGAAC‑GCCACTTTTG‑CTGCCA‑TTTCACCCGAATTC‑CCC‑GTCG‑AAGGTTTTGGCATTTAGTTGGTTtctctcccccgtcc                                                                                                                                                                                  >  2:935326/1‑226 (MQ=255)

                                                                 ggCTCAGGAGCCTCTTTGGGTCCC‑GGGCTTACACAACATGGCAGCCAGCCCATCACGCTGGACGGGCAAACGCAAAAGGTGAAGCGGCATTTGC‑TACAGGATTTCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGc                                                                                                                                                                               <  2:576165/239‑1 (MQ=255)

                                                                 ggCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACg                                                                                                                                                                                <  1:935326/238‑1 (MQ=255)

                                                                 ggCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAAAGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCa                                                                                                                                                                                                                                       >  1:1006620/1‑186 (MQ=255)

                                                                  gCTAAGGGGCAGCTATG‑GTCCATGTGCTTAAACAAACTGCCAGCCAGCCCAGCACCCTGGACGAGCAAACCAAAAAGAGGAAGCGGCATTTGC‑TACAGGATACCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGc                                                                                                                                                                               <  2:822263/238‑1 (MQ=255)

                                                                   taagagccaggcTAGGAGTCCA‑ATGCTTACACAACATGCCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCg                                                                                                                                                                             <  1:325456/229‑1 (MQ=255)

                                                                      aGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTTCAAAGCCTGTGAAATC‑CAAATTTTTG‑CTGCCA‑TTACCACAGCATTT‑GGC‑GCGG‑AACGATTTGGCATTAAGGTTGTTTTATCACACCGTCCCGCCGata                                                                                                                                                                          >  1:735365/1‑239 (MQ=255)

                                                                         agcagcTAAGAGTCCA‑GTTCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGgt                                                                                                                                                                                                    <  1:436329/210‑1 (MQ=255)

                                                                             gcTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCCACACCCGTGAACAC‑CAAATTTTTT‑CTGCCC‑TTACACCAGTAGTA‑CAC‑CCAA‑CCAGTTTTTGCATTAAGGTTGTTTCATCTCCCCGCCACGCCGGTATTTTGtt                                                                                                                                                                   >  2:1034207/1‑239 (MQ=255)

                                                                              cTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTtg‑ct                                                                                                                                                                                                                                                    <  1:42981/161‑1 (MQ=255)

                                                                                    atCCG‑GTTCTTACAAAACAGTTCACCCGGCCCATCCCCTGGGACGGGCAAATCCAAAAGGAGAAGGGGCTTTTGG‑TGCAGGATTTCCGTGGCTCGGCTTAGGTTTATGTCGTCGATTTCGGGGGGTGCAACGCCTTTGAAACC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTAc                                                                                                                                                             <  2:507255/237‑1 (MQ=255)

                                                                                      ccA‑GTGCTTACACAACATGTCCGCCAGCCCATCACGGTGGACCAGCAAACGCAAAAGGTTAAGGGGCATTGGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTg                                                                                                                                                           <  1:1198884/238‑1 (MQ=255)

                                                                                      ccA‑GTGCTTACACAACATGGCCGCCAGCCCATCACGCCGGACGAGCAAACGCAAAAGGTGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTgg                                                                                                                                                          <  1:740526/239‑1 (MQ=255)

                                                                                      ccA‑GTGCTTACAAAACATGTCAGCCAGCCCATCCCGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTg                                                                                                                                                           <  1:265572/238‑1 (MQ=255)

                                                                                        a‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCC‑TTACACCAGTATTC‑GAC‑GCAG‑AACGGTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCTGATATTTTGTTTTTTACTGgc                                                                                                                                                         >  1:326151/1‑238 (MQ=255)

                                                                                              ttaaccaacaGGTCCGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATAGCCGGCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGt                                                                                                                                                    <  1:412096/233‑1 (MQ=255)

                                                                                                 accacaatgtcccccctcccTTCCCGCTGGAAGAGAAACAGCAAAAGCGGAAGCGGCTTTTTC‑TACAGTATACCCGGCGCTCGGCTTAAGTTTATCGCGGCGATTGGGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCAcc                                                                                                                                                <  2:578210/222‑1 (MQ=255)

                                                                                                      caTGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTAtt                                                                                                                                                                                                                                  <  1:1105328/155‑1 (MQ=255)

                                                                                                       ctggcagccagcCCTTCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCg                                                                                                                                          <  1:613131/235‑1 (MQ=255)

                                                                                                       aTGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGGTTTGGCATTCAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTTGCGCCGTCACCCGGGCg                                                                                                                                          >  1:46574/1‑239 (MQ=255)

                                                                                                       aTGGCAGCCAGCCCAACCCGCGGGACGAACAAACGAAAAAGATGAAGCGGCATTTGC‑TTCAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCg                                                                                                                                          <  1:247538/239‑1 (MQ=255)

                                                                                                         gTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTTCAACGCCCGTGAAATC‑GGAATTTTTG‑CTGCCA‑TTATCCCCATATTC‑GAC‑GCAG‑AAAGTTTTTGCATTAAGGTTGTTTCATCCCCCCGCCACCCCCGAATTTTTTTTTTTATTTGCGCCGGCACCCggggca                                                                                                                                         >  2:461484/1‑235 (MQ=255)

                                                                                                              ccagccCCTCACCGTGGAAGAGAAAACGCAAAAGATGAAGGGGCATTTGC‑TACAGGATTTCCGCCGCCCGGCTTACGTTTAGCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑Ta                                                                                                                                  <  2:36107/239‑1 (MQ=255)

                                                                                                               cagccCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCg                                                                                                                                     <  1:181463/236‑1 (MQ=255)

                                                                                                                   ccctccccatGGGCGATCAAACGAAAAAGATGAAGCGGCATTTGC‑TCAAGGATTTCCTCTGCCGGGCGTACGTTTATGGGGTGGATTGGGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCt                                                                                                                             <  2:476441/230‑1 (MQ=255)

                                                                                                                    cATCCCGCTGGACGGGCAAACGCAAAAGGTGGAGCGGCATTTGG‑TACAGGATTTCCGGCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTg                                                                                                                            <  1:1307754/239‑1 (MQ=255)

                                                                                                                          ggtGGACGGGAAAACGAAAAAGGTGAGGGGGCATGTGC‑AACAGGATCTCCGTCGCTGGGCTTAGCTTTTCCGGGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTc                                                                                                                      <  2:57673/237‑1 (MQ=255)

                                                                                                                            tttcccgcctacctcAAAATGTGCAGCGGCATTTGC‑TCCAGGGACCCCGTCCCCTGGCTTACGTTTTGCGAGTACATTGCGGTGGCGGCAACCGCTGGAAAACC‑GAAATTTTTC‑CTGCCA‑TAACCCCAGTAGTC‑GTC‑GCAG‑ACAGGTTTGGCGTTAAGGTTGTTCCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGAGGCG‑TATGCCTTCTCTTCgg                                                                                                                    <  2:716429/225‑1 (MQ=255)

                                                                                                                            tGGACGAGCAAACGCAAAAGATGAAGGGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCgg                                                                                                                    <  1:876493/239‑1 (MQ=255)

                                                                                                                            tGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTacac                                                                                                                                                                                                                                         >  1:160084/1‑126 (MQ=255)

                                                                                                                            tGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTacac                                                                                                                                                                                                                                         <  2:160084/126‑1 (MQ=255)

                                                                                                                            tGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCCG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTTGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCgg                                                                                                                    >  1:778389/1‑239 (MQ=255)

                                                                                                                            tGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCCCGTCACGCCGATATTTTGTTATTTACTGGGGCAGTCACCCGGGGGATGCG‑GATGCCTGCACTTCgg                                                                                                                    >  1:896808/1‑239 (MQ=255)

                                                                                                                             ggACGAGCAAACGAAAAAGATGAAGCGGAATTTGC‑TACAGGAAACCCGTCGCCCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCggg                                                                                                                   <  1:1034207/239‑1 (MQ=255)

                                                                                                                                  agaaaccGCAAAATATTAAGCGGCATTTGC‑TACTGGATTCCCTTACCTCGGCTTGCGTTTATCGCTTCGATTTCGGCGGCTTCACCGCCTTTGAAATC‑AAAATTTTTG‑CTGCCA‑TTACCCCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCAGCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTa                                                                                                               <  2:885648/232‑1 (MQ=255)

                                                                                                                                   gCAAACGCAAAAGTTGAACCGG‑ACTTGCGTACCGTATATCCGTCCCGCGTCTTCCGTTTATCGGTTCGATTGCGGCCGCTGAAACGCCTGTGAAATC‑AAAATTTGTG‑CTGCCA‑TTACCCCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCGCCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑t                                                                                                             <  2:1033976/238‑1 (MQ=255)

                                                                                                                                             aagttgagcgggCATTTTT‑TACAGGTTTCCCGTCGGTCGGCTTACGTTTTCCGCTCCGTTTGCGCGGGCGACAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCAA‑TTACACCAGTATGC‑GGC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCCCCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGcccc                                                                                                  <  2:1270105/229‑1 (MQ=255)

                                                                                                                                             aaggttaaGCGGCATTTGC‑TACAGGTGATCCGTCCCCCGGCTTACGTTTGTCGCGTCGATTGCGGCGGCTGAAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGcccc                                                                                                  <  2:1331072/233‑1 (MQ=255)

                                                                                                                                             aaGATGAAGGGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGcccc                                                                                                  <  1:12778/239‑1 (MQ=255)

                                                                                                                                                 ttaaGCGGCTTTTGC‑TACGGGATTTCCGTTGTTCGGTTTACTTTTTTCGCGTCGATTGCGCGGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTCGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGAt                                                                                              <  2:528412/237‑1 (MQ=255)

                                                                                                                                                    aGCGGCATTTTC‑TAAGGGATATCCGTCGCTGGTCTTACGTTTATCGCGTCGTTTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCAt                                                                                           <  2:1269712/239‑1 (MQ=255)

                                                                                                                                                         cATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCCTCAACGCGTCCCGCCGATATTTTTTTATTTACTGGGGCAGTCACCCCGGCGATGGG‑TATCCCCGCACTTTGGGGGTAGTGAGT‑TGCCCCCGATCATaaaaa                                                                                      >  2:152498/1‑238 (MQ=255)

                                                                                                                                                              tggctcAGGATTCCCGGTGCCCGGCTTACGTTTATCGGGTCGATTGCGGCGGCTGCAACGCCTGGGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTtgt                                                                                  <  1:461484/234‑1 (MQ=255)

                                                                                                                                                                                gCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTT‑CTGCCA‑TTACCCCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATCTTTTTTTTTTTTCTGGCGCCGTCACCCGGCCGAGTCC‑TATGCCTGCACTTCGGGCGGA‑TGAGGCTTCCCCCGCCCAAAAAActttctcttataaccactctca                                                                 >  2:1226887/1‑219 (MQ=255)

                                                                                                                                                                                gCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCCGTGGAATC‑GAAATTTTTT‑CTGCCA‑TTACACCAAGATTC‑CCC‑CCAC‑AACGTTTTGGCCTTAAGGTTGTTTCATCCCCGCGCTCCGCCCATATTTTGTTTTTTACTTGCGCAGGCCCCCGGGCGCTGCG‑TATGCCCTCACCTTGGGCGCA‑TGAGTTTGCCCCCGATCCTAAAATTT‑TGtttttcttcgcgcggggt                                                               >  2:386550/1‑224 (MQ=255)

                                                                                                                                                                                      cTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGt                                                                                                                                                                                    <  1:1022867/119‑1 (MQ=255)

                                                                                                                                                                                      cTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGt                                                                                                                                                                                    >  2:1022867/1‑119 (MQ=255)

                                                                                                                                                                                           ttttctccggTCGTTTGGGGCGGTTGCAAGGGCTTTGAAATG‑GAAATTTTTG‑CTGCCA‑TTACCCCAGAATTC‑GAA‑CCAG‑AAAGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCCCTCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTcggcgg                                                     <  2:1268160/230‑1 (MQ=255)

                                                                                                                                                                                              tATCGTGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCCCTTCGGGCGTA‑TGAGGCTGCCCCCGCTCATAAAATTTTTGTTTCCTTCGGCGCGTGCGGTGTcggcggcgg                                                  >  1:426297/1‑239 (MQ=255)

                                                                                                                                                                                                tCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTT‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGGTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTTTTTACTTGAGCAGGCCCCCCGGCGATGCGTTATCCCT‑CACTTCGGGGGTA‑TGAGTCTTCCCCCGCTCATAAAATTTGTGTTTCCTACGgccggggcggggtggggcgggggg                                                >  2:1013847/1‑217 (MQ=255)

                                                                                                                                                                                                tCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCAcc                                                                                                                                                <  2:620466/145‑1 (MQ=255)

                                                                                                                                                                                                tCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCAcc                                                                                                                                                >  1:620466/1‑145 (MQ=255)

                                                                                                                                                                                                   cgTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCCGGCGATGCG‑TTTGCCTGCACTTCGGGCGTA‑TCAGTCTTCCCCCGATCATAAAATTTGTGTTTCCTACGGCGGGTTCGGTGGcagcggcggggatt                                             >  1:293169/1‑234 (MQ=255)

                                                                                                                                                                                                      cgtttgcggcgcttggacgcccTGTGAAATC‑GAAATGTTGG‑CGGCCA‑TTCCACCAGTATTC‑GAC‑GCAG‑AAAGTTTTGGCATTAAGGTTGTTTCAGCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTcc                                          <  2:231603/219‑1 (MQ=255)

                                                                                                                                                                                                        tttGGGGCGGCTGAAAAGCCTTTGAAATG‑AAAATTTTTG‑CTGCCATTTACACC‑GTATTC‑GAC‑GCAG‑AACTTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCTCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCAc                                        <  2:426297/238‑1 (MQ=255)

                                                                                                                                                                                                        cttgcgcgtcttccagcGCCTGTGAAATC‑GAATTTTTTG‑CTGCAA‑TTACCCCAGTCTTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTTTTTCATCCCCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCAc                                        <  2:1373139/223‑1 (MQ=255)

                                                                                                                                                                                                         tggggggtgctcgaaagccgttgAAATC‑GAAATTTTGG‑CGGCCA‑TTACCCCCGTATTC‑GAC‑GCAG‑AACGGTTTGGCATTAAGGTTGTTTCATCACGGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCAc                                        <  2:1135214/217‑1 (MQ=255)

                                                                                                                                                                                                         tggcggcgGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCAc                                        <  1:948034/236‑1 (MQ=255)

                                                                                                                                                                                                                                       taaTTTTTGTCTGCCC‑TTACCCAAGTATTC‑GGC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCCTCACCGCGTCACTCCGATATTTTGTTATTTACTGGCGCAGTCCCCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGa           <  2:896808/238‑1 (MQ=255)

                                                                                                                                                                                                                                        aaTTTTTT‑CTTCCA‑TTTACCCCGTTTTC‑GACGGCAG‑‑ATGTTTTGGCATTAAGGTTGTTTCATCCCCGCCGCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGaca         <  2:568709/239‑1 (MQ=255)

                                                                                                                                                                                                                                        aaTTTTTT‑CTGCCA‑TTGCACCCATTGTC‑AAC‑GCAG‑AACCTTTGGGCCTTAAGGTTTTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCTGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGaca         <  1:244404/239‑1 (MQ=255)

                                                                                                                                                                                                                                        aaTTTTGG‑CTGCCA‑TTACCCCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGaca         <  1:1092290/239‑1 (MQ=255)

AGCGTCTGCCCGGAATGCAAACAACGCGCGACGCTGATCAACGACGATACAGATGTACTGCTGGTGGCTAAGGAGCAGCTATGAGTCCA‑GTGCTTACACAACATGTCAGCCAGCCCATCACGCTGGACGAGCAAACGCAAAAGATGAAGCGGCATTTGC‑TACAGGATATCCGTCGCTCGGCTTACGTTTATCGCGTCGATTGCGGCGGCTGCAACGCCTGTGAAATC‑GAAATTTTTG‑CTGCCA‑TTACACCAGTATTC‑GAC‑GCAG‑AACGTTTTGGCATTAAGGTTGTTTCATCACCGCGTCACGCCGATATTTTGTTATTTACTGGCGCAGTCACCCGGGCGATGCG‑TATGCCTGCACTTCGGGCGTA‑TGAGTCTGCCCCCGATCATAAAATTTGTGTTTCCTACGGCGCGTGCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGT  >  NC_000913/2610624‑2611098

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