BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,173,361	Δ2 bp	100%	pseudogene (1‑2/442 nt)	gatC ←	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,361	0	C	.	84.0%	330.7 / 44.5	100	pseudogene (2/442 nt)	gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand): ref base C (0/16); new base . (54/30); total (54/46)
Fisher's exact test for biased strand distribution p-value = 7.37e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*	NC_000913	2,173,362	0	C	.	83.3%	312.4 / 30.9	100	pseudogene (1/442 nt)	gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand): ref base C (0/11); new base . (54/29); total (54/46)
Fisher's exact test for biased strand distribution p-value = 3.35e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 105 total aligned reads displayed.

GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCAC  >  NC_000913/2173123‑2173591

||

ccATTAGCCGCCACTTGGGGGGGCAAACTAATAGTTGGCGTGGCGATCCAAAGGTTGATGATAAAAATGATGACCCCCGAGATTAAGTTGCGAAACAGATTTCCACGATGCCCGGCCA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCTAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGT‑TTGCCGATaccac                                                                                                                                                                                                                                        <  2:256131/238‑1 (MQ=255)

  aTTACCCGCAAGTTGGGTGTGCAAACAAATAGTTTGCGTCGCGACCCACAGGGGGAGGCCCATAAGGATGACACCCGAGAGTAAGGTGCGAAACAGATTCCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGGGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATaccaccg                                                                                                                                                                                                                                      <  2:818681/239‑2 (MQ=255)

      gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATaccaccgccgt                                                                                                                                                                                                                                  >  1:1084875/1‑234 (MQ=255)

      gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATaccaccgccgt                                                                                                                                                                                                                                  >  2:805208/1‑234 (MQ=255)

      gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATaccaccgccgt                                                                                                                                                                                                                                  >  1:297263/1‑234 (MQ=255)

           cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGCGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                                                                                                                            >  2:448059/1‑238 (MQ=255)

                    gtgCAAACCAATAGTTTGGGGCGCGATCCACAGGGTGATGCTCAAAATGAGGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                                                                                                                   <  2:1320606/238‑1 (MQ=255)

                       cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                                                                                                                >  1:813869/1‑238 (MQ=255)

                       cAAACCAATAGTTGGCGTCGGGATCAACAGGGTGATTCTAATAAGGATGACCCCCGAGATTAAGGTGCGGAACATATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACg                                                                                                                                                                                                               <  2:1008437/239‑1 (MQ=255)

                        aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                              <  1:423779/239‑1 (MQ=255)

                                   ttGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                   <  1:1457269/239‑1 (MQ=255)

                                   ttGCGTCGCGATCCACAGGGTGAGGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGTTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                    <  1:19811/238‑1 (MQ=255)

                                     gCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCa                                                                                                                                                                                                 <  1:1234337/239‑1 (MQ=255)

                                      cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                                                                >  2:9436/1‑239 (MQ=255)

                                      cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                                                                >  1:271243/1‑239 (MQ=255)

                                      cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                                                                >  1:1527071/1‑239 (MQ=255)

                                           cgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                           >  1:748246/1‑239 (MQ=255)

                                                        gATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                                                               >  1:1129566/1‑238 (MQ=255)

                                                        gATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTCAccg                                                                                                                                                                               >  2:696176/1‑238 (MQ=255)

                                                        gATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCAAACAAACAGCAATTAAAAAGGGGAGTGGGATAAAAATCA‑GGC‑TTGCCGAAACCA‑‑CCGCCGTATGTCCCAGCGGCAACGCGGGATCAAGGCCAATCAGGAGCTCCTGAccg                                                                                                                                                                               >  2:1065441/1‑238 (MQ=255)

                                                         aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGAccgcc                                                                                                                                                                             >  2:783296/1‑239 (MQ=255)

                                                          tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                                                                                                                         >  1:551327/1‑226 (MQ=255)

                                                               aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                                                                                                        >  2:170022/1‑238 (MQ=255)

                                                               aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                                                                                                        >  1:472804/1‑238 (MQ=255)

                                                               aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                                                                                                        >  1:1416995/1‑238 (MQ=255)

                                                               aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACCAGGAACTCCTGACCGCCGAAcc                                                                                                                                                                        >  2:1341274/1‑237 (MQ=255)

                                                                  acgatgaCCCCCGATATTGAGGGGCGGAACAGATTCCCACGAGGCACGGCGC‑CCGCCACCCCCACGAAGAAGCCGATGGTGGCAAGACCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                                                                                                    <  2:1084875/237‑1 (MQ=255)

                                                                      tgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                                >  2:1012076/1‑239 (MQ=255)

                                                                       gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                               >  1:909313/1‑239 (MQ=255)

                                                                       gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                               >  2:1321822/1‑239 (MQ=255)

                                                                         cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                              >  1:249239/1‑238 (MQ=255)

                                                                          ccccGAAATAAAGTGGCGAAAGAGCTTTCCACGATGCCCGCGCC‑CCCCCCTCGCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                            <  2:249239/238‑1 (MQ=255)

                                                                           cccGAGATTAAGGTGCGGGACAGATTTCCACGATGCACGGCGC‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGGACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaaac                                                                                                                                                           >  2:1276687/1‑237 (MQ=255)

                                                                           cccGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                            >  2:1099339/1‑238 (MQ=255)

                                                                               agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                                                                                                                       >  1:689880/1‑239 (MQ=255)

                                                                                    aggggggaaaaCAGATTGCCGCGGGGCACGGCGA‑CCCCCATCCCCAGAAAGAAGCCAATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGCCCCACAAACAACAATTAAAATGGTGGGTGGAATAAAAATCA‑GGC‑TTGCCGTTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  2:566580/230‑1 (MQ=255)

                                                                                    agggggCGAAACAGTTTCCCGCGAGGCCCGGCGC‑CCGCCATCCCCACGAAGAAGCCAAGGGTGGCAAGATCGCCAAACGGCAGCCCCTAATTCCCCGCCCCACAAACAGAAATAAAAATGGGGAGGGGGAAAAAAACCA‑GGT‑TTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  2:297263/234‑1 (MQ=255)

                                                                                    agagtccGGAGCAGATTCCCACGAGGCACGCCAC‑CCGCCACCGCCACGAAGAAGCCGAGGGTGGCAAAATCGCAAAACGGCAGCACGTGATTACCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  2:533077/233‑1 (MQ=255)

                                                                                    aaggggggaaccAGTTTCCCACGAGGC‑CGGCGAGCCCCCATCGCCAGGAAAAGGCCGATGGGGGCAAGAGCCCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TGGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  2:580748/228‑1 (MQ=255)

                                                                                    aaGGTGGGGAACAGATTTCCACGATGCACGGCAC‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  1:1139676/239‑1 (MQ=255)

                                                                                    aaGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                                                                                                   >  1:1239869/1‑238 (MQ=255)

                                                                                    aaGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                                                                                                   >  2:310586/1‑238 (MQ=255)

                                                                                    aaGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  1:170022/239‑1 (MQ=255)

                                                                                    aaGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  1:9436/239‑1 (MQ=255)

                                                                                    aaGGTGCGAACAAGATTCCCCGGATCCCCGACGA‑CCGCCATCGCCCCAAAGAACCGAAGGGTGTAAGGACCGCAAACCGACAGCACCAGATTCCCGGGCACAAAAACGAAAATAAAACGGGAGGGAGGGATAAAAACCA‑GCC‑TT‑CAG‑GACCCCACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                                  <  2:689880/239‑1 (MQ=255)

                                                                                      ggTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                                                                                                                   >  1:226406/1‑236 (MQ=255)

                                                                                           ggAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTCCGACTGTAACCGACTACGAGCCTGCTTa                                                                                                                                           >  2:1479623/1‑239 (MQ=255)

                                                                                                 gATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGAGGGAGCCAAGATCGCCAAACGGCAGCACCTGAGTACCCGGCGCACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCA‑GGT‑TTGCCGAT‑CC‑CTCCGCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                              <  1:117986/166‑1 (MQ=255)

                                                                                                     tCCACGAAGCACGGCGA‑CCGCCACCCCCACGAAGAGGCGGAGGGTGGAAAGATCGCCAACCGGCAGCACCTGATTCCCCGGCACACAAAAAGAAATTAAAAGGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                                                                                                                 <  2:1532994/239‑1 (MQ=255)

                                                                                                         ggaggCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                                                                                                             <  1:1517544/235‑1 (MQ=255)

                                                                                                           aTGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                    >  2:117986/1‑166 (MQ=255)

                                                                                                            tGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                          >  1:365671/1‑239 (MQ=255)

                                                                                                            tGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                                                                                                                           >  2:1299874/1‑238 (MQ=255)

                                                                                                                     gcccGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACCCAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                                  <  1:1473511/237‑1 (MQ=255)

                                                                                                                     a‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCTGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                 >  1:913404/1‑239 (MQ=255)

                                                                                                                     a‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                 <  1:1321822/239‑1 (MQ=255)

                                                                                                                      tctccaaccGCAACGAAGAGGCCGATGGTGGCAAAAACGCCAAACGGCAGCACTAGATTCCCCGGCACAAAAACCGCAATAAAAAGGGTGAGTGGAATAAAAATCA‑GGC‑TGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                 <  2:568045/231‑1 (MQ=255)

                                                                                                                       ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                >  1:1295278/1‑239 (MQ=255)

                                                                                                                       ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                >  1:1249124/1‑239 (MQ=255)

                                                                                                                              cGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                                                                         >  1:395630/1‑239 (MQ=255)

                                                                                                                               gCCGCGAAGAGGCAGATGGTGGCAAGTTCGCAAACCGGCAGCCCCGGATTACCCGGCACACAAACACCAATTAAAATGGTGAGGGGGATAAAAATCA‑GCG‑TGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                        <  2:1475834/239‑1 (MQ=255)

                                                                                                                                    gaagagGCCGGTGGTGGCAAGATCGCAAAGCGCCAGCACCTGATTCCCAGGCCCACAAACAGCAATAAAAAGGGTGAGGGGAAAAAAAATGA‑GGC‑TTG‑CG‑TACCCCACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  2:485011/239‑1 (MQ=255)

                                                                                                                                    aaagaagCCAATGGGGGCAAGATCGCAAACCGGCAGCACCTGATTACCCGGCCCAAAAACAGCAATTAAAAGGGTGAGTGGGAAAAAAATCA‑GCT‑TTCCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  2:226406/238‑1 (MQ=255)

                                                                                                                                     aagaggCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGGACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  1:604229/238‑1 (MQ=255)

                                                                                                                                     aagaagGGGAGGGGGGAAAGATCCCCAACGGGCGACACATGATTACCCGGCACACAAACAGCAATAAAAATGGGGAGTGGGAAAAAAACAA‑GGC‑TGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  2:204155/238‑1 (MQ=255)

                                                                                                                                     aagaagCCAAGGGTGGAAAGACCGCCAAACGGCAGCACCTGATTACCCGGCCCAAAAACAGCAATTAAAAGGGTGAGTGGAATAAAAATCA‑GCC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  1:1479623/238‑1 (MQ=255)

                                                                                                                                            gcagggtggCAAGATCCCCAACCGCCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGTTGAGTGGGATAAAAACCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                           <  2:371845/235‑1 (MQ=255)

                                                                                                                                               tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                        >  2:817992/1‑239 (MQ=255)

                                                                                                                                               tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                        >  1:45176/1‑239 (MQ=255)

                                                                                                                                               tggggccAAAATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGAAATAAAAATGGGGAGGGGGAAAAAAATCA‑GGC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag                                                                                         <  2:442399/232‑1 (MQ=255)

                                                                                                                                                         aTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGACCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGTCATTAGCAGCATAACt                                                                               >  2:114209/1‑238 (MQ=255)

                                                                                                                                                           cGCCAAAGGGCAGCCCCTAATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                            <  1:454277/239‑1 (MQ=255)

                                                                                                                                                           cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                                  >  2:875510/1‑201 (MQ=255)

                                                                                                                                                           cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                                  <  1:875510/201‑1 (MQ=255)

                                                                                                                                                                          ccTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGAAAAAAAATCA‑GGC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACAGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                             <  1:785690/239‑1 (MQ=255)

                                                                                                                                                                          ccTGATTACCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                             <  1:1099339/239‑1 (MQ=255)

                                                                                                                                                                            gaaTTCCCCGCCCCACAAACACAAATAAAAAGGGGGAG‑GGGAAAAAAAACAGGGC‑TCGCGGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAg                                                           <  2:278860/237‑1 (MQ=255)

                                                                                                                                                                             gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCa                                                            >  2:501657/1‑237 (MQ=255)

                                                                                                                                                                                   ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                    <  1:805208/239‑1 (MQ=255)

                                                                                                                                                                                    ggtcacacaAACAGCAAGAAAAGGTGGGAG‑GGGAAAAAAAACA‑GGCTTTGCCTCTCCCACGCC‑CGATA‑GCCCCAGCAGCACCGCCGGAAAAAGCCCAATCGGGACCTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATa                                                   <  2:45176/236‑1 (MQ=255)

                                                                                                                                                                                     ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                  >  1:1314224/1‑239 (MQ=255)

                                                                                                                                                                                     ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATa                                                   >  2:1067879/1‑238 (MQ=255)

                                                                                                                                                                                      gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                  >  1:278860/1‑238 (MQ=255)

                                                                                                                                                                                        acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATaca                                                 >  2:32221/1‑237 (MQ=255)

                                                                                                                                                                                                 gCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                                  <  1:42223/163‑1 (MQ=255)

                                                                                                                                                                                                 gCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                                  >  2:42223/1‑163 (MQ=255)

                                                                                                                                                                                                 gCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAACCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                      >  1:1176739/1‑239 (MQ=255)

                                                                                                                                                                                                    aTTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                   >  1:847569/1‑239 (MQ=255)

                                                                                                                                                                                                            ggggAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACCAACCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCACCAGCTGCAATACACCTTTGACATCGAAACCcgcaa                           >  1:1092112/4‑237 (MQ=255)

                                                                                                                                                                                                            ggggAGTGGAAAAAAAACGG‑GGC‑TTGCCAAT‑CC‑CCCCGCCGAATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                            <  2:1370858/235‑1 (MQ=255)

                                                                                                                                                                                                               gAGGGGAATAAAAACCA‑GGT‑TTGCCG‑‑ACCCCACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                        <  1:1055682/239‑1 (MQ=255)

                                                                                                                                                                                                               aggTGGAAAAAAAACCA‑GGT‑TTGCCAAT‑CCA‑CCCGCCGTTGGTCCCAGCAGCACCGCCGGACCAAGCCCAATCAGTACCGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                        <  2:1295278/237‑1 (MQ=255)

                                                                                                                                                                                                                  tGGGATAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                         >  1:1207620/1‑219 (MQ=255)

                                                                                                                                                                                                                     taTAAAAACCG‑GCC‑TTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                   <  1:448059/237‑1 (MQ=255)

                                                                                                                                                                                                                       tAAAAATCA‑GGC‑TTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATCCACCTTTGACATCGTAACCCGCGAGGATGCCGATa                >  2:472710/1‑239 (MQ=255)

                                                                                                                                                                                                                            aTCG‑GGC‑TTGCAGATACC‑CCCC‑CCGTATGTCCCAGCAGCAACGCGGGATCAAGGCCAACCAGGATCCCCGTACCGCCGAACTTCGCCTGTAAACGACTCCGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAg           <  1:198355/239‑1 (MQ=255)

                                                                                                                                                                                                                                   c‑TTGCCAATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcccc       <  1:332774/237‑1 (MQ=255)

                                                                                                                                                                                                                                       gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac  >  2:1237971/1‑239 (MQ=255)

                                                                                                                                                                                                                                       gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac  >  1:80652/1‑239 (MQ=255)

                                                                                                                                                                                                                                       gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac  >  2:1285097/1‑239 (MQ=255)

||

GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGA‑CCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCA‑GGC‑TTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCAC  >  NC_000913/2173123‑2173591

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