BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	578,811	0	T	G	66.1%	119.6 / 122.3	120	Q28P (CAG→CCG)	borD	prophage lipoprotein BorD
Reads supporting (aligned to +/- strand): ref base T (17/23); new base G (31/47); total (50/70)
Fisher's exact test for biased strand distribution p-value = 8.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-04
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 123 total aligned reads displayed.

CAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCAC  >  NC_000913/578562‑579050

cAGAGAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTGCCAGCGGAGTAGAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGGGCTACTGctgtctg                                                                                                                                                                                                                                                 <  2:665624/251‑1 (MQ=255)

   aTAGAGCTGCCCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATACATGACTAAAGTCATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGCAACGAAGAAATGTTGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCtgttt                                                                                                                                                                                                                                              >  1:677408/1‑251 (MQ=255)

   aTAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGAAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCtgttt                                                                                                                                                                                                                                              >  2:192641/1‑251 (MQ=255)

          tGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTg                                                                                                                                                                                                                                        >  1:615178/1‑250 (MQ=255)

           gACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGaa                                                                                                                                                                                                                                      <  1:125611/251‑1 (MQ=255)

           gACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGaa                                                                                                                                                                                                                                      <  2:23521/251‑1 (MQ=255)

              cATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAAc                                                                                                                                                                                                                                     >  1:430727/1‑249 (MQ=255)

                 atCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTa                                                                                                                                                                                                                                 >  1:119401/1‑250 (MQ=255)

                        ggCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACACATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTg                                                                                                                                                                                                                         >  1:442191/1‑251 (MQ=255)

                        ggCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTg                                                                                                                                                                                                                         >  2:441376/1‑251 (MQ=255)

                               cATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAg                                                                                                                                                                                                                    <  1:588839/249‑1 (MQ=255)

                                  gCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACAt                                                                                                                                                                                                               >  2:323284/1‑251 (MQ=255)

                                      ttattTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTg                                                                                                                                                                                                           >  1:176450/1‑251 (MQ=255)

                                       tattTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGt                                                                                                                                                                                                          >  2:128057/1‑251 (MQ=255)

                                          tttGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGtaa                                                                                                                                                                                                        <  1:52937/250‑1 (MQ=255)

                                           ttGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAACCCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGtaat                                                                                                                                                                                                       <  1:144202/250‑1 (MQ=255)

                                             gAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCTTTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAg                                                                                                                                                                                                    <  1:192641/251‑1 (MQ=255)

                                              aGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGc                                                                                                                                                                                                   <  2:20685/251‑1 (MQ=255)

                                                 aaTACACCCGCGCTTCCAGCGGAGTATAAATGCCTAAATTAATAAAACCGATCAATCCATTTTCGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTTCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAgg                                                                                                                                                                                                <  2:442191/251‑1 (MQ=255)

                                                  aTACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAggg                                                                                                                                                                                               >  2:197917/1‑251 (MQ=255)

                                                   tACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGc                                                                                                                                                                                              >  1:153071/1‑251 (MQ=255)

                                                     cacacGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCa                                                                                                                                                                                            <  1:183098/251‑1 (MQ=255)

                                                        acGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAgcg                                                                                                                                                                                         >  2:444202/1‑251 (MQ=255)

                                                        acGCGCTTCCAGCGGAGTAGAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTGTGCTGGGTTTCTTTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTGTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTGTCCTTTGGTGCTACTGCTGCCGGTTTGTATTGAACAGTAAACGTCTGGTGAGCACATCCTGTAAGAAGCAGGGCCAgcg                                                                                                                                                                                         <  1:631694/251‑1 (MQ=255)

                                                                ccAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGc                                                                                                                                                                                  <  1:276634/250‑1 (MQ=255)

                                                                ccAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGc                                                                                                                                                                                  <  1:113622/250‑1 (MQ=255)

                                                                ccAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATATACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCg                                                                                                                                                                                 >  1:384059/1‑251 (MQ=255)

                                                                        gTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAt                                                                                                                                                                         <  1:247604/251‑1 (MQ=255)

                                                                        gTATAAATGCATAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAt                                                                                                                                                                         <  2:615178/251‑1 (MQ=255)

                                                                          ataAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAg                                                                                                                                                                                           <  2:69172/231‑1 (MQ=255)

                                                                          ataAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAtt                                                                                                                                                                        <  1:233454/250‑1 (MQ=255)

                                                                            aaagtccTAGAGTAATAAAACCGAGCAAGCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAATATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAAGTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTGTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTTGCAtgttt                                                                                                                                                                     <  2:226179/246‑1 (MQ=255)

                                                                            aaaTGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAttttt                                                                                                                                                                     >  1:119413/1‑251 (MQ=255)

                                                                             aaTGCCTAACGTAATAAAACCGGGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTCGCTGCATCGACAGTTTTCTTCTGCCCACTTCCAGAAACGAAGAAATGATGGGTGCTGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAtttttt                                                                                                                                                                    >  2:309512/1‑251 (MQ=255)

                                                                             aaTGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCAtttttt                                                                                                                                                                    <  1:655332/251‑1 (MQ=255)

                                                                              ctGCCTAAAGTAATAAAACCTAGCAATCCATTTACGCATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTTTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCAttttttt                                                                                                                                                                   >  2:92398/2‑251 (MQ=255)

                                                                               tGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTc                                                                                                                                                                  >  1:454950/1‑251 (MQ=255)

                                                                                  cTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATg                                                                                                                                                               >  1:437178/1‑251 (MQ=255)

                                                                                    aaaGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATgg                                                                                                                                                              <  2:156029/250‑1 (MQ=255)

                                                                                    aaaGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCGGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTGGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATgg                                                                                                                                                              <  1:337420/250‑1 (MQ=255)

                                                                                     aaGTAATCAACCCGAGCAATCCATTTCCGAATGTTTGCTGGGTTTCTGTTTTATCAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTGTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGTGCAGTAGCGAGTAGCATTTTTTTCATGgtg                                                                                                                                                            <  2:148770/251‑1 (MQ=255)

                                                                                        taatacAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGtt                                                                                                                                                          >  2:40233/1‑250 (MQ=255)

                                                                                            aaaaCCGAGAAATCCAGTTACGAATGTTTGCTGGGTTTCTGTTTTAACAGCATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAGCGAAGAAATGATGGGGGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGTAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑cc                                                                                                                                                    <  1:96736/251‑1 (MQ=255)

                                                                                                     cAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCttt                                                                                                                                           >  2:174245/1‑251 (MQ=255)

                                                                                                      aaTCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCtttt                                                                                                                                          >  1:350409/1‑251 (MQ=255)

                                                                                                         ccATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGGTATT‑CCCGATGCTTTTTg                                                                                                                                        <  1:510803/250‑1 (MQ=255)

                                                                                                              tACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTc                                                                                                                                  >  1:581670/1‑251 (MQ=255)

                                                                                                               aCGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCg                                                                                                                                 >  2:181161/1‑251 (MQ=255)

                                                                                                               aCGAATGTTTGCGGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTc                                                                                                                                  <  2:487818/250‑1 (MQ=255)

                                                                                                                 gAATGTTTGCTGGTTTTCTGTTTTAACAACATTTTCGGCGCCTCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGTGTGATGGTTGCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCTTGTAAGAATCATGGCCAGCGCAGTATCGAGTATCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCa                                                                                                                               <  1:428113/251‑1 (MQ=255)

                                                                                                                  aaTGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAg                                                                                                                              <  2:347910/251‑1 (MQ=255)

                                                                                                                  aaTGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAg                                                                                                                              <  2:523906/251‑1 (MQ=255)

                                                                                                                    tGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTC‑TTTGGTGCCACTGCTGCCCGGTTGTTTTGGACCGTAAAAGTCTGTTTGGCCCATTCTGTAATAAGCAGGGCCAGCGCCGTAGGGGGTAGCATTTTTTTCATGGGGGTATTTCCCG‑TGCTTTTTGAAGTTTGCAGAAt                                                                                                                           >  1:633492/1‑251 (MQ=255)

                                                                                                                     gTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAAt                                                                                                                           <  2:476695/251‑1 (MQ=255)

                                                                                                                       ttGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATcc                                                                                                                                                                                                             >  1:62309/1‑168 (MQ=255)

                                                                                                                         gCTGGGTTTCTGTTGTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGt                                                                                                                        <  2:119413/250‑1 (MQ=255)

                                                                                                                           tGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTAtg                                                                                                                     >  2:44487/1‑251 (MQ=255)

                                                                                                                              gttccTGTTTTAACAACTTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTTTTCTGCGCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCGGTGGTGCGACTGCTGCCGGTTTGTTTTGGACAGTAAACGTCTGTTGAGCACATTCTGTAGTAATCAGGGTCATGGCAGTTGCGAGGAGCATTTTTGTCAGGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATgtgt                                                                                                                  <  1:36309/247‑1 (MQ=255)

                                                                                                                                        tAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTa                                                                                                         <  1:309512/250‑1 (MQ=255)

                                                                                                                                        tAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTa                                                                                                         >  2:317617/1‑250 (MQ=255)

                                                                                                                                         aacaacATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATg                                                                                                                                                               <  2:167998/196‑1 (MQ=255)

                                                                                                                                         aacaacATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATg                                                                                                                                                               >  1:167998/1‑196 (MQ=255)

                                                                                                                                         aacaacATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTaaa                                                                                                       <  2:229536/251‑1 (MQ=255)

                                                                                                                                           caacaTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGaa                                                                                                                            >  1:456234/1‑228 (MQ=255)

                                                                                                                                           caacaTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTGTGAAGTTCGCAGaa                                                                                                                            <  2:456234/228‑1 (MQ=255)

                                                                                                                                              caTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGCTGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGAATGTGTAGAAAATTaaacaaac                                                                                                  >  1:137754/1‑251 (MQ=255)

                                                                                                                                              caTTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTaaacaaac                                                                                                  <  1:44487/251‑1 (MQ=255)

                                                                                                                                                  ttCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTGTTCTTCTGCCCAATTCCTGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACAGCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAATCCCt                                                                                               <  2:677408/250‑1 (MQ=255)

                                                                                                                                                     tGCGCAGCCACACCTTTTCGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAAc                                                                                           >  2:355350/1‑251 (MQ=255)

                                                                                                                                                         ccgccACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATTGGTGATGTGTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGAGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATg                                                                                       <  2:564896/251‑1 (MQ=255)

                                                                                                                                                           gccACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAg                                                                                     <  1:524452/251‑1 (MQ=255)

                                                                                                                                                            ccACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGt                                                                                    <  1:617856/251‑1 (MQ=255)

                                                                                                                                                                aaaTTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGaa                                                                                <  2:176450/251‑1 (MQ=255)

                                                                                                                                                                        gCTGCATCGCCAGTTTTCTTCTGCCCAATTCCCGGAACGAAGGAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGATTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGTCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCTCAGACTCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCat                                                                         >  2:271841/1‑250 (MQ=255)

                                                                                                                                                                        gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGTGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCata                                                                        <  1:617750/251‑1 (MQ=255)

                                                                                                                                                                        gCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCat                                                                         >  2:557758/1‑250 (MQ=255)

                                                                                                                                                                         cTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTTATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCatat                                                                       >  1:728227/1‑251 (MQ=255)

                                                                                                                                                                           gCATCGACAGTTTTCTTCTGCCCAATTCCCGAAACGACGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTg                                                                     <  2:710932/251‑1 (MQ=255)

                                                                                                                                                                            cATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGt                                                                    <  1:182767/251‑1 (MQ=255)

                                                                                                                                                                                   aGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAAtatt                                                             <  1:323284/251‑1 (MQ=255)

                                                                                                                                                                                    gTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAACAAATGAGTTGAAATTTCATATTGTTAAtattt                                                            >  1:651294/1‑251 (MQ=255)

                                                                                                                                                                                             tGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTACTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGt                                                    <  1:92398/250‑1 (MQ=255)

                                                                                                                                                                                              gCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTAt                                                  <  2:601524/251‑1 (MQ=255)

                                                                                                                                                                                              gCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACGGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTGCGCAGAAGCGTATGTGTAGAAAATTTAATAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATAGTTATTAATGTAt                                                  <  2:370796/251‑1 (MQ=255)

                                                                                                                                                                                                         gAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGACCAGCGCAGTAGCGAGTAGCATTTATTTCATGGTGTTATT‑CCCGCTGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGATGCg                                        >  2:625608/1‑250 (MQ=255)

                                                                                                                                                                                                               aagaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGa                                                                                                               <  1:686558/173‑1 (MQ=255)

                                                                                                                                                                                                               aagaaATGATGGGTGATGGTTTCCTTTGGAGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGa                                                                                                               >  2:686558/1‑173 (MQ=255)

                                                                                                                                                                                                                 gaaATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTTTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGt                               <  1:730968/251‑1 (MQ=255)

                                                                                                                                                                                                                       atgGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGt                         >  2:194117/1‑251 (MQ=255)

                                                                                                                                                                                                                          ggTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCATCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTACGCCAGGTGCGATGAATCGTCATTGTAtt                      <  1:681086/251‑1 (MQ=255)

                                                                                                                                                                                                                              aTGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGGTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCg                  <  2:251011/251‑1 (MQ=255)

                                                                                                                                                                                                                               ggggttCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCgg                 <  2:241196/247‑1 (MQ=255)

                                                                                                                                                                                                                                ggTTTC‑TTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGa                <  2:633492/250‑1 (MQ=255)

                                                                                                                                                                                                                                  tttCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGACTCGTCATTGTATTCCCGGAtt              >  1:99434/1‑251 (MQ=255)

                                                                                                                                                                                                                                    tCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTa             >  1:510701/1‑250 (MQ=255)

                                                                                                                                                                                                                                       tttGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACt          >  2:577292/1‑250 (MQ=255)

                                                                                                                                                                                                                                          ggTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGt      <  1:187023/251‑1 (MQ=255)

                                                                                                                                                                                                                                          ggTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGt      >  1:109685/1‑251 (MQ=255)

                                                                                                                                                                                                                                              cTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAACCCAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCcac  >  2:356725/1‑251 (MQ=255)

                                                                                                                                                                                                                                              cTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCcac  >  1:335569/1‑251 (MQ=255)

CAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAGCATTTTTTTCATGGTGTTATT‑CCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGAAAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCAC  >  NC_000913/578562‑579050

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