BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	540,610	C→T	9.7%	Q16* (CAG→TAG)	ybbY →	putative uracil/xanthine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	540,610	0	C	T	9.7%	322.2 / 4.6	123	Q16* (CAG→TAG)	ybbY	putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand): ref base C (32/75); new base T (11/2); total (46/77)
Fisher's exact test for biased strand distribution p-value = 2.06e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

AGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAACCATCCTTACTATCACCCTTGGTGAAGCATCGCGCGGG  >  NC_000913/540373‑540843

aGTCACCCCTTTTTTTGGCCGTACCCTTGGCGCGCGTATCACGAAAACCATCTTACGTGGGGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCa                                                                                                                                                                                                                                             <  2:478623/239‑1 (MQ=255)

                ggCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCC‑CGTTAAGGTGGGTTTATG‑TTCAATTTTTCAGTCAGCCCCGGAA‑ACCCGTTATAAGGCTTTTAGT‑G‑GTGTTTCtttt‑tt                                                                                                                                                                                                                          >  1:729088/1‑239 (MQ=255)

                    gTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTTTG‑TTCAATTTTTCAGTCAGCCGCGCAA‑ACCCGGTATCAGGATTTTTGG‑T‑GGTTTTCtttt‑tttttt                                                                                                                                                                                                                      >  1:997606/1‑239 (MQ=255)

                      tccATTGGGGCACGTATCACGAAAATCAACTTACGTGGTGATTGGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTg                                                                                                                                                                                                                     <  1:518518/237‑1 (MQ=255)

                          ttGGTGCGGGTTTCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAaca                                                                                                                                                                                                                <  2:40282/239‑1 (MQ=255)

                           tgggcggcgTTTCACGAAAACCATCTTACGTGGTGATGTGATTTACGCCATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAacac                                                                                                                                                                                                               <  2:862854/233‑1 (MQ=255)

                           tGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCATTTCC‑AGTTCAGGGGTGTTTTTG‑TTTCAATTTGCCGTCAGCCGCCGAA‑GACCGTTTTTTGGGTTTTAGT‑G‑GGTtttttttt‑ttttttttaaaaa                                                                                                                                                                                                               >  2:220315/1‑232 (MQ=255)

                           tGGCGCCGGTTTCAAGAAAACCATTTTACGGGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAacac                                                                                                                                                                                                               <  2:729088/239‑1 (MQ=255)

                                cgcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGGTCAGGGGTGTTTATG‑TTCCATTTTGCCGTCAGCCGCGAAA‑GCCTGTTATCGAGGTTTTTGT‑G‑GTTTTTCTTTA‑TTTTTTGCCACAAGGttt                                                                                                                                                                                                          >  1:679891/1‑238 (MQ=255)

                                cgcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCAATCTCC‑CGGTAAGGGGGGGTTATG‑GTCAAATTTTCAGTCAGCCCCGAAA‑GCCTGTTTTTAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAAAACGGGtg                                                                                                                                                                                                          >  1:950964/1‑239 (MQ=255)

                                  cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCTCC‑GGTTAGAGGGGGTTTATG‑TTTAATTTTGCAGTCACCCCCGCAA‑GCCTTTTTTTCGGGTTTCTGTGG‑GTTTTTTTTTT‑TTTTTTTCAACAGGGtgtg                                                                                                                                                                                                         >  2:406543/1‑236 (MQ=255)

                                  cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTCCC‑GGGTAAGGGGGGGTTTTGTTTCATTTTTGC‑GTCAGCCCCGAAACGCCTTTTTTCGGGTTTTAGGG‑G‑GTTTTTTTTTT‑TTTTTTGAAAAAGGGTtgg                                                                                                                                                                                                         >  2:1156894/1‑238 (MQ=255)

                                  cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTTAGGGGGGTTTATG‑TTCCATTTTTCCGTCAGCCGCGCAA‑GACCGTTATCAGGGTTTTAGT‑G‑GTTTTTCTTTT‑TTTTTTTCTACACCGGtgtg                                                                                                                                                                                                        >  1:65955/1‑239 (MQ=255)

                                     aTCCCGAAAACCCTCTTGCGGGGGGATGTTATTTTCGGCATTTAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAAGCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGtt                                                                                                                                                                                                     <  2:1272716/239‑1 (MQ=255)

                                       cccGAACACCATCTTTCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTTTTGCGCGGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTtcc                                                                                                                                                                                                   <  2:674317/237‑1 (MQ=255)

                                              cccATTTTTCCTGGTGATGTGATTTACGACATTGAACGGGGCTTCCCTGGTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                                                                                                                                                            <  2:657563/238‑1 (MQ=255)

                                              aCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                                                                                                                                                            <  1:962702/239‑1 (MQ=255)

                                              aCCATCTTACGTGGGGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                                                                                                                                                            <  1:1029272/239‑1 (MQ=255)

                                                 aTCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGCAA‑GACTGTTTTCCGGGTTTCCTT‑G‑GTTTTTCTTTA‑TTTTTTGTAAAACGGTTGTGGTTCCTCCCACGcc                                                                                                                                                                                          >  1:1261000/1‑237 (MQ=255)

                                                       cGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGGGTGTTTATG‑TTCAATTTTGCAGGCAGCCCCGAAA‑GCCTGTTATCAGGGTTTTAAT‑G‑GGTTTTTTTTT‑TTTTTTTCAAAACGGTTTTGGTTCCTCCTCCGCTACttttt                                                                                                                                                                                   >  1:474000/1‑237 (MQ=255)

                                                          ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCC‑GGTTTAGGTGTGGTTTTG‑TTCAATTTTTCCGTCAGCCGCGAAA‑GCCTGTTTTCTGGGTTTCAGT‑G‑GGTTTTTTTTT‑TTTTTTTCAACCCGGTTTTGTTTTCCCCCCCGCCACTTTCTgtt                                                                                                                                                                                >  2:701133/1‑237 (MQ=255)

                                                           gtgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGGTAAGGGGTGTTTTTG‑TTCAATTTTGCAGTCAGCCCCGAAA‑GCCCGTTATCAGGATTTACGT‑G‑GGTTTTCTTTA‑TTTTTTTCAACACGGGTTTGGGTCCTCCTCCCCGACTTTCTGCtt                                                                                                                                                                               >  1:1080003/1‑239 (MQ=255)

                                                            tgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGGTAAGGGGTGTTTATG‑TTCAATTTTGCAGTTCGCCCCGCAA‑GCCTGTTATCAGGATTTTAGT‑G‑GTTTTTCTTTA‑TTTTTTTCCACACGGTTGTGGTTCCTCCCACGCTACTTTTTGCttt                                                                                                                                                                              >  1:1267629/1‑239 (MQ=255)

                                                              atgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTtcctc                                                                                                                                                                                                 <  1:934171/218‑1 (MQ=255)

                                                                      tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑CGTTAAGGTGGGTTTTTG‑TTCAATTTTTCAGTCAGCCCCGAAA‑GCCCTTTATCAGGGTTTTCGT‑G‑GGTTTTCTTTT‑TTTTTTTCAACACGGTTTGTGTTCCTCCTTACCTTCTTTCTTTTTTTCAGTTTcc                                                                                                                                                                     >  1:414693/1‑238 (MQ=255)

                                                                        ccgactttaaCCGGGCTTCCCTTTTGCGCCGAAAGGTCCATTTATCCTTAAACATCAGCAGTAATCGGGCCCCGGCAATGCCCGGCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCa                                                                                                                                                                  <  2:1317139/231‑1 (MQ=255)

                                                                        cGACATTGAACAGGGCTTCCGTGTTGCGCCGAAAGGTCAATTTATCCGTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCa                                                                                                                                                                  <  1:292796/239‑1 (MQ=255)

                                                                        cGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGCATCTTC‑GGTTTAGGTGTGTTTATA‑TTTTATTTTTGAATCAGCCCGGGAA‑GACCGGTATTCGGGTATGCGG‑G‑GTTTTTTTTTATTTTTTTGCAA‑ACGGGTGTTGTTTCTCCCACGCCCCTTTCTGGTTTTCAATTTGCGCa                                                                                                                                                                  >  2:1257318/1‑239 (MQ=255)

                                                                        cGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGt                                                                                                                                                                                                      <  1:1364628/203‑1 (MQ=255)

                                                                           cATTGAACAGGGCTTCCATGTTTCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCa                                                                                                                                                                  <  1:1005280/236‑1 (MQ=255)

                                                                               gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCGTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                                                                           <  1:406543/239‑1 (MQ=255)

                                                                               gAACAGGGCTTCCCGGTTGGGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTTCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                                                                                                                                                           <  1:527335/239‑1 (MQ=255)

                                                                                     ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGGGTGTTTATG‑TTCAATTTTGCAGTCAGCCCCGCAA‑GACCGTTATCAGGGTTTCTGG‑G‑GGTTTTTTTTT‑TTTTTTTCAAACCGGTTTTGGTTTCCCCCTCCCTACTTTTTGTTTTTTTTTTTCCGCaaaaaaagccgcgt                                                                                                                                                     >  2:662200/1‑228 (MQ=255)

                                                                                      tcTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCt                                                                                                                                                      <  1:538914/236‑1 (MQ=255)

                                                                                          cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTTCC‑GGGTTAGGGGGGTGTATT‑TTTCATTTTGCAGTCAGCCGCCGAA‑GCCTGGTTTCAGGATTTTATG‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGGTGGTGTGTCTCCCACGCTCCTTTCTTCTTTTCAGTTTCCGCCAAAGAGCCCCCTTACGc                                                                                                                                                 >  2:1197987/1‑238 (MQ=255)

                                                                                          cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCC‑TGTTAAGGTGTGTTTATG‑TTCAATTTTGCCGTCAGCCGCGGAA‑CCCTGTTATAAGGGTTTCTTG‑G‑GTTTTTTTTTT‑TTTTTTTCAACACGGTTGTGGGTCTCCCTCCGCTACTTTTTTCTTTTTTGTTTTCGCACAGGAGcactcttaccc                                                                                                                                                 >  2:495442/1‑228 (MQ=255)

                                                                                          cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGGGGGGTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTTGTATCAGGATTTCCGT‑G‑GGTTTTCTTTT‑TTTTTTTCAACACGGTTTGGGTTCTCCCTAGCACACTTATTTCTTTTTAGTTTCCGCAAAGTAGcgcgtttaacc                                                                                                                                                 >  2:679232/1‑228 (MQ=255)

                                                                                          cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑CGTTAAGGTGTGTTTTTG‑TTCCATTTTGCCGTCAGCCGCGAAA‑GCCCGTTTTTCGGATTTTAGT‑G‑GTTTTTTTTTT‑TTTTTTGCTACACGGTTTTGGGTCCTCCCCCCCCACTTTTTGCTTTTTTGTTGCCCCCAAGACGCCCGCtttccc                                                                                                                                                 >  2:842352/1‑234 (MQ=255)

                                                                                           ccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATTCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTc                                                                                                                                               <  1:305822/239‑1 (MQ=255)

                                                                                           ccTGTGGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTc                                                                                                                                               <  1:388235/239‑1 (MQ=255)

                                                                                            cTGTTGCGCCCAAAGTCTAATGTTTTCTTGAACATCAGCAGTAGTCTGGACCCTGAAATCCCCGTCGCTCCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTACG‑TGGAATTTTGCAGTCAGGCGCGAAA‑GCCTGGTATCAGGATTTCCGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTCGCCGCAAAGTAGCCTGCTTACGCTCa                                                                                                                                              <  2:1052035/239‑1 (MQ=255)

                                                                                                tGCGCCGAAGGGTCAATTTATCCTTAAACATCATCAGTAATCGGGCCCCTGCAATGCCCGTCCTTGCGGCGGTCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCGCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCa                                                                                                                                          <  1:63240/239‑1 (MQ=255)

                                                                                                 gcgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAAGGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCaa                                                                                                                                         <  1:209131/239‑1 (MQ=255)

                                                                                                            tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCtttt                                                                                                                                <  1:842352/237‑1 (MQ=255)

                                                                                                             taaTTTATCCTTAAACTTCAGCAGTAATCTGGCCCCTGCAATGCCCGTGCTTGCGGCGGGCATTCGCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCt                                                                                                                              <  1:836241/237‑1 (MQ=255)

                                                                                                                      cTTAACCCTCAGGAGTAATCTGGCCCCTGGAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCATTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                                                                                    <  1:1346637/239‑1 (MQ=255)

                                                                                                                      cTTAAACCTCAGCAGTAAGCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCTTTCTCC‑GGTTAAGGTGTGTTTATT‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                                                                                    <  1:779106/239‑1 (MQ=255)

                                                                                                                      cTTAAACAGTAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                                                                                    <  1:203772/239‑1 (MQ=255)

                                                                                                                              tCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGGGTGTTTATG‑TTCAATTTTTCAGTCAGCCGCGAAA‑GCCTGTTATCAGGGTTTCAGG‑G‑GTTTTTCTTTA‑TTTTTTTCAACACGGTTGTGGGTCCTCCTACCCTACTTTCTTCTTTTTTTTTGCCCCCAAGTAGCCTGCTTACGCTCACCCAAAATGTTTTTTTTTCTACCGCACCGGCCt                                                                                                             >  2:1098944/1‑238 (MQ=255)

                                                                                                                                   agTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCtt                                                                                                         <  1:1214556/237‑1 (MQ=255)

                                                                                                                                     tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCCCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTGCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAata                                                                                                                                       <  1:732519/205‑1 (MQ=255)

                                                                                                                                      aaTTTGGCCCCTGCATTGCCCGTCCTTGCGGCGGGCCTTCCCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCg                                                                                                       <  2:128078/236‑1 (MQ=255)

                                                                                                                                       agttGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCa                                                                                                                    <  1:677710/219‑1 (MQ=255)

                                                                                                                                        tCTGGCCCCGGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAg                                                                                                  <  1:317415/239‑1 (MQ=255)

                                                                                                                                            gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTTAGGTGTGTTTTTG‑TTCAATTTTTCAGTCAGCCGCGAAA‑ACCTGTTATCAGGATTTTAGT‑G‑GGTTTTTTTTT‑TTTTTTGCCACAAGGTTTTGGTTCCTCCTACGCTCCTTTCTTCTTTTCAGTTGCCGCAAAGTAGCCCGCTTACGCTCACCCAATATGCTTTTCTTTCTACCGCACTGGCCTGCTTCGCCTAGGCGt                                                                                              >  1:1354322/1‑239 (MQ=255)

                                                                                                                                            gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGGTAAGGTGTGTTTATG‑TTCAATTTTTCAGTCAGCCGCGGAA‑GCCTGTTTTCAGGGTTTCCGT‑G‑GGTTTTCTTTA‑TTTTTTGCAACACCGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTTCGTTGCCGCAAAGTAGACTGCTTACGCTCACCCCATATGCTTTTCTTGCTACCGCCCTGGCCTGCTTCGCTCAAGc                                                                                                >  1:162642/1‑237 (MQ=255)

                                                                                                                                                 tGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTTTG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTTCCGCAAAGGAGCCTGCTTACGCTCACGCAATATGCTTTTTTTGCTACCGCACCGGCCTGCTTCGCTCAGGCGgttttc                                                                                         >  1:826301/1‑237 (MQ=255)

                                                                                                                                                 tGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                                                                                                                                                            <  1:241919/140‑1 (MQ=255)

                                                                                                                                                      tGCCCGTCCTTGCGGCGGGCATTCTCCGGGTTAAGG‑GGGGTTATG‑TTTCATTTTGCAGTCAGCCGCGAAA‑GCCTGGTATCAGGGGTTCAGT‑G‑GTTTTTTTTTA‑TTTTTTTCCACACGGGTTTGGGTCCTCCTACGCTACTTTCTTCTTTTTAGTTGCCCCACAGTAGCCCGCTTACTCTCCCCCACTAAGCTTTTCTTGTTCCCCCACCGGCCCGCTTCGCTCCGGCGgttttctggc                                                                                     >  2:213432/1‑232 (MQ=255)

                                                                                                                                                           gTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGGGTGGTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCCCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgt                                                                                >  1:350094/1‑238 (MQ=255)

                                                                                                                                                            ttcTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                                                                              <  1:531921/237‑1 (MQ=255)

                                                                                                                                                            tCCTTGGGGCGGGCATTCTCC‑GTTTAAGGTGTGTTTATT‑TTCAATTTTGCAGTCCGCCGCGAAA‑GCCTTTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                                                                              <  2:51650/239‑1 (MQ=255)

                                                                                                                                                            tCCCTGCGGCGGGCATTCTCG‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCATCACGGTTGTGGTTCCTCCTACGCTACGTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCAtcgtcg                                                                              <  1:150251/239‑1 (MQ=255)

                                                                                                                                                             ccTTGGGGGGGGGCTTCTCC‑GGTTAAGGGTTTTTTATG‑TTCAATTTTGCAGGCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgc                                                                             <  2:997606/239‑1 (MQ=255)

                                                                                                                                                             ccTTGCGGCGGGCATTCTCC‑GGTTAAGGTTTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgc                                                                             <  1:422194/239‑1 (MQ=255)

                                                                                                                                                             ccTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCCGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgc                                                                             <  1:640780/239‑1 (MQ=255)

                                                                                                                                                              cttggcgcggGCATTCTCC‑GGTTAAGGTGTTTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                                                                            <  1:1004925/233‑1 (MQ=255)

                                                                                                                                                              cTTGCGGCGGGCCTTGTCC‑GGTTTAGGTGTGTTTATG‑TTCAATTGTGCAGTCAGCCGCGAAA‑GCCTGGTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                                                                            <  1:407384/239‑1 (MQ=255)

                                                                                                                                                              cTTGCGGCGGGCATTTTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                                                                            <  1:213432/239‑1 (MQ=255)

                                                                                                                                                              cTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                                                                            <  1:583211/239‑1 (MQ=255)

                                                                                                                                                                         cATTCCCC‑GGTTAAGGTGTTTTTATT‑TTCAATTTTGCAGGCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGa                                                                  <  1:679232/238‑1 (MQ=255)

                                                                                                                                                                                   gTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCTTCCCGCTATTATGGAAGGGCCAggg                                                        >  2:1207146/1‑238 (MQ=255)

                                                                                                                                                                                    ttAAGGTGTTTTTATT‑TTCAATTTTGCCGTCGGCCGGGAAA‑GCCTGTTATCAGTATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTGCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGGTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTg                                                       <  2:417525/239‑1 (MQ=255)

                                                                                                                                                                                      aaGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGc                                                     <  1:1257318/239‑1 (MQ=255)

                                                                                                                                                                                          tttgtTTATT‑TTAAATTTTGCATTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑T‑GTTTTTCTTTA‑TGTTTTGCAACCCGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCtg                                                  <  2:414693/236‑1 (MQ=255)

                                                                                                                                                                                          tgtttTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCtg                                                  <  1:1052581/238‑1 (MQ=255)

                                                                                                                                                                                          gttttttATG‑TTCAATTTTGCAGTCAACCCCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACTGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCtg                                                  <  1:264224/234‑1 (MQ=255)

                                                                                                                                                                                            tgtTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTgtgg                                               >  2:94667/1‑239 (MQ=255)

                                                                                                                                                                                              tttATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACg                                                                                                                                                  >  2:791385/1‑138 (MQ=255)

                                                                                                                                                                                              tttATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACg                                                                                                                                                  <  1:791385/138‑1 (MQ=255)

                                                                                                                                                                                              tttATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTgtggtg                                             >  1:719359/1‑239 (MQ=255)

                                                                                                                                                                                                 aTG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAATGGCCAGGTGGCCTGTGGTgggg                                          >  1:755244/1‑239 (MQ=255)

                                                                                                                                                                                                   g‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTTTGGGGGGGa                                         >  1:599032/1‑238 (MQ=255)

                                                                                                                                                                                                     tttaaTTTTGGAGTCAGCCGGGAAA‑GCCTGTTTGTAGGATTTCAGT‑G‑GTTTTTTTTTA‑GTTTTGGAAACCCGGTGGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCGTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAAc                                       <  2:372403/236‑1 (MQ=255)

                                                                                                                                                                                                     ttCAATTTTGCATTCATCCGCGAAA‑TCCTGTTTGTAGGATTTC‑GT‑GTGTTTTTCTTTT‑TTTTTGGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGAAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAAc                                       <  2:1080003/239‑1 (MQ=255)

                                                                                                                                                                                                       cAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAAcc                                      >  1:1008957/1‑238 (MQ=255)

                                                                                                                                                                                                         aTTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCCCAGGCGTTTTGCGGTCATCTTCGCGCTCTTATGGAAGGGCCAGGTGGCCTGGGGGGGGGAACCaa                                    >  2:183793/1‑237 (MQ=255)

                                                                                                                                                                                                               cAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGt                                                                                                                                                                         <  2:101709/99‑1 (MQ=255)

                                                                                                                                                                                                               cAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGt                                                                                                                                                                         >  1:101709/1‑99 (MQ=255)

                                                                                                                                                                                                                 gttccccccgAAA‑GCCTTTTTTTAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACCCGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAACCATCCTTACTAt                           <  2:1008957/231‑1 (MQ=255)

                                                                                                                                                                                                                     gtcgcaAAA‑GCCTTTTTTCAGGTGTGTCGT‑G‑TTTTTTCTTTA‑GTTTTGGTCACACGGGTGGGGTTCCCCCCACGATGCTTTCTGGTTTGCAGGTGCCGGTAAGTGGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACGGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAACCATCCTTACTATCAcc                       <  2:356882/237‑1 (MQ=255)

                                                                                                                                                                                                                       ggcgAAA‑GCCTTTTATCAGGATTTTGGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACCCGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGGAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAACCATCCTTACTATCACCCt                     <  2:65955/238‑1 (MQ=255)

                                                                                                                                                                                                                       cgcgAAA‑GCCTGTTATCAGGATTTCAGT‑T‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAACCATCCTTACTATCACCCt                     <  1:1098944/239‑1 (MQ=255)

                                                                                                                                                                                                                                         aGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAacac                                                                                                                                                                                                               >  2:491875/1‑36 (MQ=255)

                                                                                                                                                                                                                                         aGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAacac                                                                                                                                                                                                               <  1:491875/36‑1 (MQ=255)

                                                                                                                                                                                                                                         aGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAACCATCCTTACTTTCACCCTTGGTGAAGCAtaccccg    >  1:1288056/1‑233 (MQ=255)

                                                                                                                                                                                                                                            aTTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                                                                            <  1:684720/164‑1 (MQ=255)

                                                                                                                                                                                                                                            aTTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTcgcg                                                                            >  2:684720/1‑164 (MQ=255)

                                                                                                                                                                                                                                              ttCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCGAAGTAGCCTGCt                                                                                                                                                      <  1:1148414/88‑1 (MQ=255)

                                                                                                                                                                                                                                              ttCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCt                                                                                                                                                      <  2:1308808/88‑1 (MQ=255)

AGTCAGCCCGTATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCC‑GGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAA‑GCCTGTTATCAGGATTTCAGT‑G‑GTTTTTCTTTA‑TTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACGCAATATGCTTTTCTTGCTACCGCACTGGCCTGCTTCGCTCAGGCGTTTTGCGGTCATCGTCGCGCTATTATGGAAGGGCCAGGTGGCCTGTGGTGGGGAACCATCCTTACTATCACCCTTGGTGAAGCATCGCGCGGG  >  NC_000913/540373‑540843

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