BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,073	G→A	P91P (CCG→CCA)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,073	0	G	A	98.9%	286.7 / NA	91	P91P (CCG→CCA)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand): ref base G (0/1); major base A (40/49); minor base C (1/0); total (41/50)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCG‑ATGACCCTGGAAAAGAA‑TATGCC‑GATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGC  >  NC_000913/2786‑3366

aGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCAGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCCGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGc                                                                                                                                                                                                                                                                                            >  2:120717/1‑300 (MQ=255)

aGTTAGGAGTCCGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGATGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATCTCGTTTATCAGTGCTGGGAGCGTTTTTGCAAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑CTGACTCTGGAAAAGAC‑TACGCC‑AATCAGTTCGGGc                                                                                                                                                                                                                                                                                            >  2:119739/1‑300 (MQ=255)

 gTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCt                                                                                                                                                                                                                                                                                           <  1:524215/300‑1 (MQ=255)

 gTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCAAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCt                                                                                                                                                                                                                                                                                           <  1:43139/300‑1 (MQ=255)

      gAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGNTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGc                                                                                                                                                                                                                                                                                      <  1:506980/300‑1 (MQ=255)

      gAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGc                                                                                                                                                                                                                                                                                      <  2:422344/300‑1 (MQ=255)

       aGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAggat                                                                                                                                                                                                                                                                                     <  1:464652/300‑3 (MQ=255)

       aGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCt                                                                                                                                                                                                                                                                                     >  1:95795/1‑300 (MQ=255)

        gTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGTAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGCA‑TATGCC‑AATCGGTTCGGGCTTAGGCTc                                                                                                                                                                                                                                                                                    >  2:290561/1‑300 (MQ=255)

                tGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAACTTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCt                                                                                                                                                                                                                                                                            >  2:110172/1‑300 (MQ=255)

                 ggTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCt                                                                                                                                                                                                                                                                            >  1:399667/1‑299 (MQ=255)

                             tGCCCTGGCTTCCAGTGCCAATATGAGCTTCGGGTTTGATGTGCTCGTGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATTCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTcgc                                                                                                                                                                                                                                                               <  2:368428/300‑1 (MQ=255)

                                    gCTTCCAGTGCCAATATGATCGTCGGGTTTGATGGGCCCTGGGTGGCGTTGACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGTTTGAGGCGGCATAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAACTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCGGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTcgc                                                                                                                                                                                                                                                               <  2:368379/293‑1 (MQ=255)

                                       tCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATg                                                                                                                                                                                                                                                     <  2:410003/300‑1 (MQ=255)

                                          aGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCTGAGATGTAGTCACGGTTGAGGCGGCATAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACTGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGc                                                                                                                                                                                                                                                   <  1:312203/299‑1 (MQ=255)

                                          aGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGc                                                                                                                                                                                                                                                   <  1:21992/299‑1 (MQ=255)

                                                 gcgagatgtgtataagagacaggCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGatgaat                                                                                                                                                                                                                                            >  1:473509/23‑299 (MQ=255)

                                                      agcttCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTAGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAacac                                                                                                                                                                                                                                      <  1:10909/296‑1 (MQ=255)

                                                      aGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAacac                                                                                                                                                                                                                                      <  1:113721/300‑1 (MQ=255)

                                                      aGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCATAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAacac                                                                                                                                                                                                                                      <  1:112822/300‑1 (MQ=255)

                                                      aGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAacac                                                                                                                                                                                                                                      <  1:70337/300‑1 (MQ=255)

                                                          tCGGGTTTGATATGCTCGGGGCGGCGGTTACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGGTTGAGGCGACAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTTGCG‑ATGACTCTGGAAAAGAA‑GATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCACTTGTTCGGTGGTcgc                                                                                                                                                                                                                                                               <  2:368250/271‑1 (MQ=255)

                                                             ggTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCa                                                                                                                                                                                                                               >  2:630940/1‑300 (MQ=255)

                                                              gTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATAGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCaa                                                                                                                                                                                                                              <  1:258703/300‑1 (MQ=255)

                                                                       gCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTaa                                                                                                                                                                                                                     >  1:460798/1‑300 (MQ=255)

                                                                       gCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTAGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTaa                                                                                                                                                                                                                     >  2:87981/1‑300 (MQ=255)

                                                                         tCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGGTTGAGGCGGCATAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATg                                                                                                                                                                                                                   <  1:94666/300‑1 (MQ=255)

                                                                         tCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATg                                                                                                                                                                                                                   <  1:94161/300‑1 (MQ=255)

                                                                           ggggCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGNGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGa                                                                                                                                                                                                                  <  2:508181/299‑1 (MQ=255)

                                                                            gggCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGaca                                                                                                                                                                                                                >  2:591952/1‑300 (MQ=255)

                                                                              gcggcgGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAacacctgtctcttatacacatctcgcag                                                                                                                                                                                                              <  1:190451/300‑25 (MQ=255)

                                                                                ggcgtTGACACCTGTTGATGGTGCATTGCTCGTAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGTGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGacac                                                                                                                                                                                                               <  2:384329/297‑1 (MQ=255)

                                                                                             gTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCttt                                                                                                                                                                                                >  1:696663/1‑299 (MQ=255)

                                                                                             gTTGATGGTGCATTGCTAGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCA‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCttt                                                                                                                                                                                                >  2:627299/1‑299 (MQ=255)

                                                                                                 aTGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCATGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGAt                                                                                                                                                                                             >  1:633622/1‑298 (MQ=255)

                                                                                                    gTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGc                                                                                                                                                                                         >  1:395229/1‑299 (MQ=255)

                                                                                                           gCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGa                                                                                                                                                                                 >  2:403898/1‑300 (MQ=255)

                                                                                                           gCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGa                                                                                                                                                                                 >  1:175722/1‑300 (MQ=255)

                                                                                                               gtagaTGTAGTCACGGTTGAGGCGGCATAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACTGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAgga                                                                                                                                                                             <  2:398959/298‑2 (MQ=255)

                                                                                                                tagaTGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAggac                                                                                                                                                                            >  1:561126/2‑298 (MQ=255)

                                                                                                                gagaTGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAggac                                                                                                                                                                            >  1:105834/1‑298 (MQ=255)

                                                                                                                  gaTGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAggacga                                                                                                                                                                          >  1:452935/1‑296 (MQ=255)

                                                                                                                              gTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGTTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCAGgcagc                                                                                                                                                              >  2:136214/1‑300 (MQ=255)

                                                                                                                              gTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGgccgc                                                                                                                                                              >  2:136267/1‑297 (MQ=255)

                                                                                                                              gTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGgcagc                                                                                                                                                              >  1:454414/1‑300 (MQ=255)

                                                                                                                                      ggcAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGTGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTa                                                                                                                                                      <  1:216236/300‑1 (MQ=255)

                                                                                                                                          gagaCATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCATGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGAc                                                                                                                                                  <  2:395229/300‑1 (MQ=255)

                                                                                                                                           agaCATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGAc                                                                                                                                                  <  2:410962/299‑1 (MQ=255)

                                                                                                                                              cATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAAc                                                                                                                                               >  2:603531/1‑299 (MQ=255)

                                                                                                                                                    ttctcAACAACCTCTTACGCTTTGCCGATAAGCTTCCTTCAGAACCACGGTAAAATATCGTTTATCAGTGCTGGGATCGTTTTTGCCAGGAGCTTGGCAGGCAAATGTCGGTGGCG‑ATGACCCTGGAAAAAAA‑TATGCC‑AATCGGTTCTGGCTTAGGCTCCAGCGCCTGTGCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTTCGGCAGCATTCATGACGACAACGTGGCAc                                                                                                                                        <  2:390753/299‑1 (MQ=255)

                                                                                                                                                    gTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGGCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGGAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCAc                                                                                                                                        <  2:648081/300‑1 (MQ=255)

                                                                                                                                                        ccacaacCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGNCACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCgt                                                                                                                                     >  1:135923/3‑299 (MQ=255)

                                                                                                                                                        caacaaCCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCgt                                                                                                                                     >  1:53487/1‑299 (MQ=255)

                                                                                                                                                        caacaaCCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCgt                                                                                                                                     >  2:184423/1‑299 (MQ=255)

                                                                                                                                                               cTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTtctt                                                                                                                             >  1:489019/1‑299 (MQ=255)

                                                                                                                                                                   gACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTggtg                                                                                                                         <  1:348333/300‑1 (MQ=255)

                                                                                                                                                                         ttGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGc                                                                                                                   <  1:913/300‑1 (MQ=255)

                                                                                                                                                                           gCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCa                                                                                                                  >  1:235846/1‑299 (MQ=255)

                                                                                                                                                                              gATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGGGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTtg                                                                                                              <  1:426808/300‑1 (MQ=255)

                                                                                                                                                                              gATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTtg                                                                                                              <  1:426903/300‑1 (MQ=255)

                                                                                                                                                                              gATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATAAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGt                                                                                                                <  1:695291/298‑1 (MQ=255)

                                                                                                                                                                                        cGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCga                                                                                                      >  2:69869/1‑298 (MQ=255)

                                                                                                                                                                                                     ggAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGAcatcat                                                                                       <  1:436179/300‑1 (MQ=255)

                                                                                                                                                                                                          atatCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATcagc                                                                                   >  2:488849/1‑299 (MQ=255)

                                                                                                                                                                                                                       ttgctagaagcgcccttGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑CATCGGTTCGGGCTTCGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGACTGCACACTGTGGCAAGCCGCTTAATGACACTCGTTTGTTGGCTTTGATGGGCGAGCTGGAAGGACGCATCTCCGGCCGCACTCATTACGCCAACGTGGCACCGTGTTTTCTCGGTGGTATGCCGTTGATGATCGAAGAAAACGACATCATCGGCCAGCAAGTGCCa                                                                       >  2:679497/16‑298 (MQ=255)

                                                                                                                                                                                                                       gTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATTGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCAAGCAAGTGc                                                                         >  2:468591/1‑296 (MQ=255)

                                                                                                                                                                                                                       gTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGCCAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGc                                                                         >  2:198176/1‑296 (MQ=255)

                                                                                                                                                                                                                       gTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGCTCGAAGAAAACGACATCATCAGCCAGCAAGTGCCa                                                                       >  2:29265/1‑298 (MQ=255)

                                                                                                                                                                                                                       gTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCa                                                                       >  1:126819/1‑298 (MQ=255)

                                                                                                                                                                                                                       gTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCATGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAAAGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGCTCGAAGAGAACGACATCATCATCCGGCAAGTGCCa                                                                       >  2:30784/1‑298 (MQ=255)

                                                                                                                                                                                                                        tacTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAggg                                                                    <  1:519027/298‑1 (MQ=255)

                                                                                                                                                                                                                        tGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAggg                                                                    <  1:518965/300‑1 (MQ=255)

                                                                                                                                                                                                                                       tGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCtg                                                    <  1:571289/300‑1 (MQ=255)

                                                                                                                                                                                                                                             gagcttggcaaGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGt                                               >  2:438674/10‑299 (MQ=255)

                                                                                                                                                                                                                                               gcttggcaaGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AAGCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCt                                            <  2:558439/293‑1 (MQ=255)

                                                                                                                                                                                                                                                   ggcaaGCAAATTCCAGTTGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTTGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCTGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCg                                        <  2:489019/297‑1 (MQ=255)

                                                                                                                                                                                                                                                   ggcaaGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGNTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCg                                        <  1:149383/297‑1 (MQ=255)

                                                                                                                                                                                                                                                   ggcaaGCAAATTCCAGTGGCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCg                                        <  1:599649/297‑1 (MQ=255)

                                                                                                                                                                                                                                                      agacgaaTTCCAGTGGCNCATGACTCTGGAAAAGAANTATGCCAGATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTNGGCTGTGGGTGCTGGc                                         <  1:581687/295‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ggCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGTGGATTAAAGt                       <  1:120717/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ggCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTCAAGt                       >  2:614327/1‑300 (MQ=255)

                                                                                                                                                                                                                                                                    ggCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGACAGCAAGTGCAAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGt                       >  2:458130/1‑300 (MQ=255)

                                                                                                                                                                                                                                                                    ggCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGAGGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTTGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGt                       <  1:119723/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                     gCG‑ATGACTCTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGt                       <  1:119739/299‑1 (MQ=255)

                                                                                                                                                                                                                                                                             ctcTGGAAAAGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGAc                 >  1:189942/3‑298 (MQ=255)

                                                                                                                                                                                                                                                                                     aaGAG‑TATGCC‑AATCGGTTCGGGCGTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGAGGGCGATGAATGAACACTTCGGCAAGCCGCTTAATGACACTCGTTGGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACTTGGCACCGTGTTTTCTTGTGGGTATGCAGGGGATGATCGAAGAAAACGACACCAGCAGCCAGCAAGGGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGGATCCTGGGATTAAAGTCTCGACGGCAGAAGcc       <  1:331504/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                                     aaGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGTGGATTAAAGTCTCGACGGCAGAAGcc       <  1:391586/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                                     aaGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGcc       <  1:392749/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                                     aaGAA‑TATACC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGcc       <  1:108031/300‑1 (MQ=255)

                                                                                                                                                                                                                                                                                      aGAA‑TATGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTTGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCGGGAAGCACGAATCTCCGGCAGTATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGcc       <  2:561126/299‑1 (MQ=255)

                                                                                                                                                                                                                                                                                           tatGCC‑AATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCAAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGTGGATTAAAGTCTCGACGGCAGAAGCCAGGGc  <  1:87981/300‑1 (MQ=255)

AGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCG‑ATGACCCTGGAAAAGAA‑TATGCC‑GATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGT‑GGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGC  >  NC_000913/2786‑3366

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