BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,970	T→C	A456A (GCA→GCG)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,970	0	T	C	97.2%	439.3 / ‑2.0	149	A456A (GCA→GCG)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (50/91); total (56/93)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

Only 100 of 149 total aligned reads displayed.

CAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTC  >  NC_000913/66723‑67209

cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCa                                                                                                                                                                                                                                                 <  1:929328/251‑1 (MQ=37)

          aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCAT‑GGCTGa                                                                                                                                                                                                                                         >  1:461821/1‑248 (MQ=37)

            gCCGGATGTGCTTTGCTTTTCCGGACTACAAAATCGCAGCGTGTAGGCTTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCGGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACg                                                                                                                                                                                                                                    <  2:809136/251‑1 (MQ=255)

               ggATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTa                                                                                                                                                                                                                                 <  1:173395/251‑1 (MQ=38)

                 attagcTTTGGTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTCCTTAGCGACGAAATTCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATg                                                                                                                                                                                                                               <  2:327825/247‑1 (MQ=37)

                   gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATgg                                                                                                                                                                                                                              <  1:137351/250‑1 (MQ=38)

                             tATCCGGCCTACAAAATCTCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGcca                                                                                                                                                                                                                     <  1:77961/249‑1 (MQ=255)

                              aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGGGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGccacca                                                                                                                                                                                                                  <  2:810264/251‑1 (MQ=255)

                                  ggCCTAAAAAATCGTAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTGACACACCTGATGCGACTACTCAGCGACGAATTCCGTAATACACTTCGTTCCA‑GCGCATCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGa                                                                                                                                                                                                              <  2:285197/251‑1 (MQ=255)

                                   gCCTACAAAATCGCAGAGTGTAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGCATCGTTTTCAATCACCGTGATTTCAAATTCGTGCATCTCGGCGAACTGGGGCGTAACGTTGAGGTTCCACGCAC‑GGCTTGAGACGGTATGGTGCGCGccccccgaggg                                                                                                                                                                                                             >  2:595543/1‑242 (MQ=255)

                                    ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCTCATATCGTTGAGGGTCCGAGCAT‑GGCTGCAGACGGTATGGTGCGCGCCCCCCGCGAgg                                                                                                                                                                                                            >  2:347618/1‑251 (MQ=255)

                                       aCAAAAGCGGAGCGTGTAGGCCTGATAAGCCGCGTCAGCAGAGCAGCAGGCGTTACACAACGGAGGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTGCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTGATCAATCACCGTGATTTCAATGTCGTGCATCGCGGCGAACTGGCGAATATCGTTGAGGTTCAGCGCAG‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑Tc                                                                                                                                                                                                        <  2:319158/251‑1 (MQ=255)

                                         aaaaACGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCa                                                                                                                                                                                                      <  2:135556/251‑1 (MQ=255)

                                          aaaTCGCCGCGTGTAGGCCTGATAAGACGCGACAGCGTCGCAACAGGCGTTACGCACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCA‑TCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGGTGTCAATCACCGTGATTGCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGATGA‑TCCAc                                                                                                                                                                                                     <  2:205031/251‑1 (MQ=255)

                                              cGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACCCTTCGTTCCA‑GCGCATCGCGTCTATAAACGCTGGCAGGCTGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGATGTTCAGCGCAT‑GGCTGAAGACGGTTTGGTGCGCGCCACCAGCGAGGA‑TCCACGCtt                                                                                                                                                                                                 <  1:584342/251‑1 (MQ=255)

                                               gCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTc                                                                                                                                                                                                <  1:496030/251‑1 (MQ=255)

                                                cAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCg                                                                                                                                                                                               <  1:6759/251‑1 (MQ=255)

                                                 aGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGACGGGTCTCGTTATCAATCACCGTGATTTACATGTCGTGCAACTCGGCGAACTGGAGCATGACGTTGATGGTCAG‑GCCTAGGCTGAAGAATGTATGGGTCGAGCCACACGCGAGGA‑TCCACGccttgg                                                                                                                                                                                              >  2:213804/1‑246 (MQ=255)

                                                        aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGt                                                                                                                                                                                      <  2:632408/251‑1 (MQ=255)

                                                           cctgATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑AGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTgg                                                                                                                                                                                   <  1:241907/248‑1 (MQ=255)

                                                               aTAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGc                                                                                                                                                                                  >  1:702692/1‑248 (MQ=255)

                                                                           aGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCt                                                                                                                                                                   <  1:269609/251‑1 (MQ=255)

                                                                              gTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCCCCGTGCTTTCAATGTCGTGCATCTCGGGGAACTGGCGCATATCGTTGAGGTTCCGCGCAT‑GGCTGAAGACGGTATGGTTCGCGACACACGCGAGGA‑TCCACGCTTCGGAA‑GAAGTTGGCAGATCCGGTTGCGCtttt                                                                                                                                                                >  1:645861/1‑250 (MQ=255)

                                                                                 gCATTAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCcac                                                                                                                                                             <  1:717794/251‑1 (MQ=255)

                                                                                 gCATCAGGCGTTACACACCTGAGGCGATTAATTAGCGACGAAATCCGTAATACCCTTCGTTCCC‑GCGCAGCGAGGCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCTCCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGTCAGATCCGGTTGCGCTTTCcac                                                                                                                                                             <  2:108475/251‑1 (MQ=255)

                                                                                        gCGTTACACCCCTGATGCGTCCACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCAt                                                                                                                                                      <  2:519981/251‑1 (MQ=255)

                                                                                        gCGTTACACACCTGATGCGACTTCTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGCCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCCATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCAT‑GGCTTAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTTGCAGATCCGGTTGCGCTTTCCACAGCGCAt                                                                                                                                                      <  2:104429/251‑1 (MQ=255)

                                                                                         cGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCAtt                                                                                                                                                     <  2:167406/251‑1 (MQ=255)

                                                                                          gTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAAATGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGGATTGGGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTc                                                                                                                                                    >  1:868891/1‑251 (MQ=255)

                                                                                           ttACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTTTCGTTATCAATCACCGTGATTTCAATGTCTTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCCACGCCT‑GGCTGAAGACGGTATGGTGCGCGCCCCCCGCGGGGA‑TCCACGCTTCGGAA‑GCCGTTTGCAGATTCGGTTGCGCCTTCCACAGCGCATTct                                                                                                                                                   >  1:349276/1‑250 (MQ=255)

                                                                                             acacacctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGcc                                                                                                                                                 <  1:921060/243‑1 (MQ=255)

                                                                                                  cctgATGCGACTACTTCGCGACGAAAGCCGTAATCCACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCgg                                                                                                                                            <  2:632705/248‑1 (MQ=255)

                                                                                                     gATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCACCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGTA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                                                                                                                                         <  2:406851/251‑1 (MQ=255)

                                                                                                       ggcgacttactAGCGACGAAATCCCTAATACACTTTGGTCCAAGC‑CAGGGCGGCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTGCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGtt                                                                                                                                       <  2:954197/241‑1 (MQ=255)

                                                                                                         cgacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGTGGTTCAGCGCAT‑GGCTGGAGGCGGGATGGTGCGCGCCCCCAGCGAGGA‑TCCACGCCTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCCTTTGCCACCGGCAGtttt                                                                                                                                     >  2:255382/4‑250 (MQ=255)

                                                                                                         cgacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGTAGACGGTATGGGGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGGAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAggttt                                                                                                                                     >  1:704961/4‑247 (MQ=255)

                                                                                                         cgacTACTGAGCGACTAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTAGCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc                                                                                                                                     <  2:611681/248‑1 (MQ=255)

                                                                                                          gacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTCTCGTTATCCATCACCGGGATTTCAATGTCGAGCATCCCGGCGAACTGGCGCATCTCGGTGAGGGTAAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCCCCAACGGGGG‑TCCACCCTTCGGAACGGAGT‑GGCAGGTCCGGGTGGGCTTTCCACCGCGCCATCGCCACCAGCaacttcg                                                                                                                                    >  2:486716/3‑245 (MQ=255)

                                                                                                           acTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTTTCAATCACCGGGGTTTCCATGTCGTGCCTCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGTGCAT‑GGCTGCAGACGGTATGGTGTCCGCCACCGGCGGAGT‑GCCCCGCCTCGGGA‑GGAGTTGGCAGAGCCCGTTGCGCTTTCCACCGCCCATTCTCCCCCGGccgattccc                                                                                                                                   >  2:274558/2‑243 (MQ=255)

                                                                                                               cTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGGTTTCAATGTCGTGCCTCTCGGCGAACTGGCGAATATCGTTGAGGTTCAGAGCGT‑GGCTGGAGACGGTAGGGGGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑ACAGTTGGCAGATCCGGTTGCGCATTCCACAGCGCCATCGCCACAGGCCGTTTCGGcgcg                                                                                                                               >  1:295358/1‑248 (MQ=255)

                                                                                                                ttAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGATGGCAGGCGGGGATCGCTAACAATCACCGCGACTTCAATGACCAGCATCTCGGAGAACCGGCGAATCAAGTTGAGGTTCAGCACAT‑GGATGAAGACGGTAAGGTGCGGGCCACCAGCGAGGA‑TCCACGATTCGGAA‑GGCGTTGGCAGAACCGGTTCCGGTTTCCACAGCGCATTCGCCCACGGGAGTTTCGGCaagg                                                                                                                              >  1:467742/1‑248 (MQ=255)

                                                                                                                  aGCGACGAAATCCTGAATACACGTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGa                                                                                                                             <  2:296435/250‑1 (MQ=255)

                                                                                                                        gaAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGgt                                                                                                                      <  1:420463/251‑1 (MQ=255)

                                                                                                                        gaAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGgt                                                                                                                      <  1:36424/251‑1 (MQ=255)

                                                                                                                         aaatccGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCgggg                                                                                                                     >  1:172407/5‑249 (MQ=255)

                                                                                                                          aatccGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGTCTGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCGCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGTTTTCCACAGCGTATTCGCCACCGGCAGTTTCGGCATGTAGTGCGgtgt                                                                                                                    <  1:204740/248‑1 (MQ=255)

                                                                                                                          aatccGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCCGCGAGGA‑TCCACGCTTCGGAA‑GCCGTTGGCAGATCCGGTTGCGCTTTCCCCAGCGCCTTCGCCACCGGCAGTTTCGGCAGGGAGTGCggggt                                                                                                                    >  1:391553/4‑248 (MQ=255)

                                                                                                                                   tACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGGA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCCGTTTCGGCCGGGGGTGCGGTGTTTTCACCgg                                                                                                           >  1:332450/1‑250 (MQ=255)

                                                                                                                                     cacTTCTTTCCA‑GCGCAGCGCGATTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGGATGGTGCGCTCCACCAGCGAGGA‑TCCTCGCTTCGGAA‑GCAGTTGGCATATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCgtgt                                                                                                         <  2:559808/251‑1 (MQ=255)

                                                                                                                                      acTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTc                                                                                                        <  1:412330/251‑1 (MQ=255)

                                                                                                                                        ttCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGGATGGGGCGCGCCCACAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTTCGCTTTCCCCAGCGCATTCGCCACCGGCCGTTTTGGCAGGGGGTGGGGGGTTTTCAcccgggcgg                                                                                                      >  2:401323/1‑244 (MQ=255)

                                                                                                                                           gTTCCA‑GCGCAGCGCGTCTTTAAACGCGGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGAtac                                                                                                   <  1:662150/251‑3 (MQ=255)

                                                                                                                                                a‑gcgcagcgcGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGGGCTTTCAATGTCGTGCCTCTCTGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑TGCTGAAGGAGGTTTGGATCGCGCCACCAGCCAGGGGTCC‑CGCTTTGGGA‑GCAGTTGGCCGATCCGGTTGCGCTTTTCACAGCGCATTTCCCACCGGCAGGTTCGGCAGGGAGGGCGGGGTTTTTTCCGTGTCGGTGCactgt                                                                                              >  2:963081/1‑247 (MQ=255)

                                                                                                                                                a‑gcgcagcgcGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCATCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTa                                                                                              <  1:551869/251‑1 (MQ=255)

                                                                                                                                                  gcgcagcgcGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCCCCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGAGCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTaa                                                                                             <  2:504026/251‑1 (MQ=255)

                                                                                                                                                     cagcgcGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCa                                                                                          <  2:454029/251‑1 (MQ=255)

                                                                                                                                                         gcgTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAgcag                                                                                      <  1:652328/251‑4 (MQ=255)

                                                                                                                                                               tgaaaCGCTGGCAGGCGGGTATCGTGATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACgg                                                                                  <  1:2236/247‑1 (MQ=255)

                                                                                                                                                                 aaactcgggCAGGCGGGTATCTTTATCAATCACCTGGGTTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGAGCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTaa                                                                               <  2:513033/243‑1 (MQ=255)

                                                                                                                                                                 aaaCGCTGGCAGGCGGGGATCGGTATCAATCACCGAGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCAAATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGGGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGGTGCGCTTTCCCCAGCGCATTCGCCACCGGCAGGTTCGGCAGGGGGGGGGGGGGTTTTACCGTGTCGATAAAGGTAACCAGCAGAAGGTaa                                                                               >  1:389155/1‑250 (MQ=255)

                                                                                                                                                                    cGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGAt                                                                           <  1:878655/251‑1 (MQ=255)

                                                                                                                                                                     gCTGGCAGGTGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATc                                                                          <  2:233832/251‑1 (MQ=255)

                                                                                                                                                                     gCTGGCAGGCGGGTATCTGTAGCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTGGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATc                                                                          <  2:511647/251‑1 (MQ=255)

                                                                                                                                                                        ggcaggctggtatCGTTATCAATCTCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGGTGGCAGATCCGGGTGCGCTGTCCACAGCGCATTCGCCACCGGCAGTTGCGGCATGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGcc                                                                       <  1:403323/239‑1 (MQ=255)

                                                                                                                                                                            ggcgggtatCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGAc                                                                                    <  1:611677/226‑1 (MQ=255)

                                                                                                                                                                             gcgggtatCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCg                                                                      >  1:258063/8‑247 (MQ=255)

                                                                                                                                                                              cgggtatCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatc                                                                 <  1:404430/245‑1 (MQ=255)

                                                                                                                                                                                   atCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatc                                                             >  1:829869/2‑250 (MQ=255)

                                                                                                                                                                                     cGTTATCAATCACCGTTATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatcagg                                                          <  2:387439/251‑3 (MQ=255)

                                                                                                                                                                                     cGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGCGTGCGGTGTTTTAACCGTGTCGGACCAGTTAACCAGCAGACGGAACCGGTCGCCATGatcaatcagg                                                          >  2:829669/1‑249 (MQ=255)

                                                                                                                                                                                     cGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatcagg                                                          <  1:489726/251‑3 (MQ=255)

                                                                                                                                                                                     cGTTATCAATCACCGTGATTTCAATGTCGTGCATCGCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGagcaatcag                                                           <  2:391553/250‑2 (MQ=255)

                                                                                                                                                                                       ttATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAggg                                                                                                                              <  2:690638/181‑1 (MQ=255)

                                                                                                                                                                                       ttATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAggg                                                                                                                              >  1:690638/1‑181 (MQ=255)

                                                                                                                                                                                        tATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCt                                                        <  1:365829/250‑1 (MQ=255)

                                                                                                                                                                                           caatcaCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCg                                                    >  1:977750/1‑251 (MQ=255)

                                                                                                                                                                                            aatcaCCGTGATT‑CAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGGTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGGTAACCAGCAGAAGGTAACGATCGCCGAGAGCAATCAGGCTGGCGAc                                                  >  1:655792/1‑251 (MQ=255)

                                                                                                                                                                                                     gATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTg                                          <  1:908744/251‑1 (MQ=255)

                                                                                                                                                                                                         tcattgtctttcATCCCGGCGCACTGGCGCATATCGTTGATGTCCAGCGCAT‑GGCTGAAGATGGTATGGGTCGCGCCCCCATCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGAAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCCCCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGCCGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggac                                       <  2:751383/239‑3 (MQ=255)

                                                                                                                                                                                                         tCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCCGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggacc                                      >  1:525423/1‑248 (MQ=255)

                                                                                                                                                                                                         tCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggacc                                      >  1:709733/1‑248 (MQ=255)

                                                                                                                                                                                                          cAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggacc                                      <  1:145590/250‑4 (MQ=255)

                                                                                                                                                                                                          cAATGGCGTGCATCGCGGCGAACGGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggacc                                      <  2:342012/250‑4 (MQ=255)

                                                                                                                                                                                                           aaTGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCgg                                    >  1:955330/1‑251 (MQ=255)

                                                                                                                                                                                                           aaTGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGGTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGGTACAGTTAACCAGCAGACGGGAACGGTCGCCGAGATCAATCAGGCTGGCGACAAACCCTGGACCgg                                    >  2:310691/1‑251 (MQ=255)

                                                                                                                                                                                                              ggcGTGCATCGCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGtt                                  <  2:709733/248‑1 (MQ=255)

                                                                                                                                                                                                                cGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCGT‑GGCTGAAGACGGTATGGTGAGCGCCACCAGCGAGGA‑TACACGCTTAGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGCTACAGGTAACCAGCAGACGGTAACGATCGCCGCGATAAATCAGGATGGCGACAATCGCTGGACCGGTTTgg                               >  1:586180/1‑251 (MQ=255)

                                                                                                                                                                                                                 gggCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGAAT‑GTCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTggg                              <  2:612003/249‑1 (MQ=255)

                                                                                                                                                                                                                 gTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTggg                              >  1:931356/1‑251 (MQ=255)

                                                                                                                                                                                                                  ggCATCGCGGCGAAGTGGCGCATATCGTTGATGGTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑GCCTCGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTGCGGCAGGGAGGGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGAGCGCCGAGAGCAATCAGGCTGGCGACAATCGCTGGATCGGTGTGGGt                             <  2:168354/250‑1 (MQ=255)

                                                                                                                                                                                                                       ctcGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGCGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCCCCGGCAGTTTCGGCAGGGAGTGCGGGGCTTTCACCGTGTCGCTCCCGTTAACCAGCAGCCGGTTCCGAGCGCCGTGCTCCATCCGGCTGGCGACACTCGCTGGCCCGGtttccgagttgg                        >  2:185765/1‑242 (MQ=255)

                                                                                                                                                                                                                         cGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCGTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCATGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAg                      <  1:36119/251‑1 (MQ=255)

                                                                                                                                                                                                                         cGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGaa                       <  2:369577/250‑1 (MQ=255)

                                                                                                                                                                                                                             gaactGGCGCATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCa                  <  1:13493/247‑1 (MQ=255)

                                                                                                                                                                                                                               actGGCGCATATCGTTGAGGGTCAGCGGAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCATCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTGCGGCAGGGAGTGCGTTGTTTTCACCGTGTCGATACAGTTACCCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAgt                <  1:892071/249‑2 (MQ=255)

                                                                                                                                                                                                                                 tGGCACATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGGGCGGTGTTTTCACCGTGGCGATACAGTTAACCAGCAGACGGTGACGATCGCCGAGAACAATCCGGCTTGCGACAATACCCGGACCGGATTGGGAATTGAAAATCAgacg              >  2:893266/1‑248 (MQ=255)

                                                                                                                                                                                                                                    cgcATATCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGc           <  1:595543/251‑1 (MQ=255)

                                                                                                                                                                                                                                        tatCGTTGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGGGGGCGGTGGTTTCACCGGGTCGATACCGTTAACCAGCAGACGGTAACGATCGCCGGGGTCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGGATTGAAGATccatcgggcaggg       >  2:695448/1‑239 (MQ=255)

                                                                                                                                                                                                                                             ttGAGGTTCAGCGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGTGCGCGGTGTTTTCACCGTGTCGATCCCGTTAACCAGCAGACGGTAACGATCGCCGCGATCAATCAGGCCGGCGAACATCGCTGGACCGGTTTGGGTAATGAAGATCAAGCGGGGAggcacggc  >  2:362259/1‑245 (MQ=255)

                                                                                                                                                                                                                                                aGGTTCAGCGCAT‑GGCTGAAGAGGGTATGGTGCGCGCCACCATCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATCCAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTc  <  1:211415/248‑1 (MQ=255)

CAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCAT‑GGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGA‑TCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTC  >  NC_000913/66723‑67209

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