BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,761,230	G→A	S7F (TCC→TTC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,761,230	0	G	A	99.3%	428.8 / NA	137	S7F (TCC→TTC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base G (1/0); major base A (59/76); minor base C (0/1); total (60/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

GGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGG  >  NC_000913/3760983‑3761474

   tCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAACGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAATTGACTATAGAGGAAAcact                                                                                                                                                                                                                                                 >  2:315832/1‑248 (MQ=255)

   tCCCGTCAGGGTCATCACCGTACGCTGCGCGCTCTGCTCTGCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACTCAGGCAGCGTCTGGCTGCCACTAATCACTTCTCGCGCTTCGTCGATCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAATGAAGGAGCTATCGCGCAATAAGCGATCAATATCCGGAAGTTGACTATAGAGGAAAcgc                                                                                                                                                                                                                                                  <  2:115427/250‑1 (MQ=255)

     ccGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGtt                                                                                                                                                                                                                                               >  1:72751/1‑251 (MQ=255)

      cGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCTCTTCTTGCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACCCCACGCAGTCAGCGTCTGTCTGCCACGAAGCACTTCGCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGATATCGCGCAATAAGCGATCAATAGCCGGAAGTTTACTATAGAGGAAAcgcg                                                                                                                                                                                                                                                 <  2:654394/248‑1 (MQ=255)

      cGTCAGATTGATCACGGGACGCAGCGCGCTTTGCGCTTCTTTCGTCAGCCGGGCCTCGACTTCTTGCGCCAAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAACAACTTCTTGCGCTTCGTCGAGCATCTGACGCAACACTTCCCCCACGCGGCTGTGACCATAAGTTTCACGACAAGAAAGGAAGGAGCTATGGCGCAATAAGCGATGAATGGCCGGACGTTGACTAGAGAAGAAACCCGttt                                                                                                                                                                                                                                              >  2:132044/1‑251 (MQ=255)

       gTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAACGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTc                                                                                                                                                                                                                                             >  1:5081/1‑251 (MQ=255)

       gTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTc                                                                                                                                                                                                                                             <  1:79120/251‑1 (MQ=255)

                         cGCAGCGCGCTCTGCGCTGCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCATCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCt                                                                                                                                                                                                                           <  2:325488/251‑1 (MQ=255)

                           cagtgcgcTTTGCGCTGTTTTCTTCATCCGGGCATCGCCTTCTTGCGCCCAGTTTCCACTCCACGCTGGAGGCGTCTGGCTGCCACGAATCCCTTTTCGCTCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAATGAGCTATCGCGCAATTAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTgg                                                                                                                                                                                                                         <  2:525622/247‑1 (MQ=255)

                              cgcgcTCGTCGCTTCTTTCGTCAACCGGGCATCGACTTCTTTCCCCCAGTTTTCACCCCACGCAGGCAGCGTCTGGCGGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTGGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTg                                                                                                                                                                                                                      <  2:1987/251‑1 (MQ=255)

                                     tGCGCTTCTTTCGTCAACCGGGCATAGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCGGGCTGATCaa                                                                                                                                                                                                                 <  2:363863/249‑1 (MQ=255)

                                     cgccctgcttTCGCCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCTGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAgg                                                                                                                                                                                                               <  1:535569/244‑1 (MQ=255)

                                       cgcTGCTTTCGGCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGGCGGGCTGCCACGAATCACTTCGCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTAGCACGCAAAGAAAGGAAGGAGCTAGTGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATATAGGAAACGCGTTTCGGTTGTCAGAGGAAACCTGGCTGATCAAGGcc                                                                                                                                                                                                             <  2:347822/251‑1 (MQ=255)

                                                  tCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACgg                                                                                                                                                                                                  >  1:117756/1‑251 (MQ=255)

                                                    aaCCGGGCATCGACTTCTTGCGCCCAGTTTTCACCCCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGgag                                                                                                                                                                                                <  2:258798/251‑1 (MQ=255)

                                                        gggCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAgg                                                                                                                                                                                            <  1:19524/251‑1 (MQ=255)

                                                        gggCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAgg                                                                                                                                                                                            <  1:137989/251‑1 (MQ=255)

                                                              cGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGCATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGttt                                                                                                                                                                                      >  1:432396/1‑251 (MQ=255)

                                                               gACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGTCTGATCAAGGCCCTCTCACACGGAGAAGGGCGttt                                                                                                                                                                                      >  1:309828/1‑250 (MQ=255)

                                                                   tcttGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACa                                                                                                                                                                                  <  1:397150/250‑1 (MQ=255)

                                                                   tcttGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACAt                                                                                                                                                                                 >  1:218016/1‑251 (MQ=255)

                                                                        cgcCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCa                                                                                                                                                                            <  2:290909/251‑1 (MQ=255)

                                                                         tcccAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGACGGGCGTTTAACATAACCAc                                                                                                                                                                           >  2:566231/2‑251 (MQ=255)

                                                                         gcCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCAc                                                                                                                                                                           >  1:845079/1‑251 (MQ=255)

                                                                                   cacaCCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTaa                                                                                                                                                                  <  1:37682/250‑1 (MQ=255)

                                                                                      accacGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTg                                                                                                                                                              >  1:439496/1‑251 (MQ=255)

                                                                                      accacGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTg                                                                                                                                                              <  1:123059/251‑1 (MQ=255)

                                                                                       tcacGCAGGCAGCGTCTGGCTACCACGAATCAATTATCGCGCTTCGTCGAGCATCTGAAGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGAAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGAAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGCTCAAGGCCCTCTCACACGGCGAAGGGCGTTTAACATAACCACGGCTTTTAACg                                                                                                                                                                >  2:554710/2‑248 (MQ=255)

                                                                                          cGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTg                                                                                                                                                              >  1:582159/1‑247 (MQ=255)

                                                                                                  gCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCagg                                                                                                                                                  <  1:864754/251‑1 (MQ=255)

                                                                                                     tCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCaggagg                                                                                                                                               >  1:360103/1‑251 (MQ=255)

                                                                                                         gCTGCCACGAATCACCTCTCGCGCCTCGGCGAGCAGCTGACGCAACAATGCCACCACGCGGGTGTGACCATAAGGATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGAACAATAGCCGTAAGTGGACTATAGAGGAAACGCGTGTTGGTTGTCATAGGAACCCTGGCTGATCAAGGCCCGCTCACACGGAGAAGGGCGTGTAACATAACCACGGATTTTAACGTGAGATGGGTCATGAGGACag                                                                                                                                           <  2:465/251‑2 (MQ=255)

                                                                                                           tGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCAATCAATAGCCGGAAGTTGACTATAGAGAAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATAAAGGCCCTCTCACACGGAAAAGGGAGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCCGGAGGACatat                                                                                                                                         >  2:451498/1‑251 (MQ=255)

                                                                                                              cACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGCTGGGTCGGGAGGACATATcgc                                                                                                                                      >  2:498711/1‑251 (MQ=255)

                                                                                                               aCGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATcgcg                                                                                                                                     >  1:824748/1‑251 (MQ=255)

                                                                                                                  aaTCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCAt                                                                                                                                  <  1:451190/251‑1 (MQ=255)

                                                                                                                      ccTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAg                                                                                                                              <  1:275427/250‑1 (MQ=255)

                                                                                                                          ctcGCGCTTCGTCGCGCCTCTGACGCAACGATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGACAGGAAGGCGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGAGAGATGGGTCAGGAGGACATATCGCACATCAAGCCtt                                                                                                                          >  1:701609/1‑251 (MQ=255)

                                                                                                                          ctcGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTa                                                                                                                                                                                     >  1:56/1‑192 (MQ=255)

                                                                                                                          ctcGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGAAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGAGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGAAATATCGCGAATCAAGCCtt                                                                                                                          >  2:762480/1‑251 (MQ=255)

                                                                                                                           ccgcgCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCttt                                                                                                                         <  1:235479/250‑1 (MQ=255)

                                                                                                                             gcgcTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTgg                                                                                                                       >  1:374759/1‑251 (MQ=255)

                                                                                                                             gcgcTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGCCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGCGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTgg                                                                                                                       >  1:95405/1‑251 (MQ=255)

                                                                                                                               gcTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCg                                                                                                                     <  1:850774/251‑1 (MQ=255)

                                                                                                                                  tcgtcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCCCGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCCCACGGAGGAGGGCGTTTACCATAACCACGGCTTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGtt                                                                                                                  >  2:824190/1‑251 (MQ=255)

                                                                                                                                  tcgtcgAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCgg                                                                                                                    <  1:369159/249‑1 (MQ=255)

                                                                                                                                    gtcgAGAATATGAAGCAACAATTCAACCACGCGGGTGTGACCATAAGTATCACGCAAAGCAAGGCAGGAGCTATCGCGCAATCAGCGATAAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGACACCTGGCTGATCAAGGACCTCCAACACGGAGAAGGGCCTTTAACATAACCACGGATTTGAACATGAGATGGGTCAGGAGGACATATCGAGCATCAAACATTTggagggtcg                                                                                                                >  2:818904/1‑244 (MQ=255)

                                                                                                                                      cgatcATCTGACGCAACCATTCCACCCCGCGTGTGCGGCCATAAGTACCACGAATAGAAAGGCAGTAGCTTGCTCGCTATAATCGGTCAATAGCCCGAAGTTGACTATAGAGGAGACGCGTTTCGGTTGTCATAGGAAACCTGGCTGAGCAAGGCCCTATCACACGGAGAAGTGCGTTTATCTTACCCCCGGATTGTACCGTGAGATGGGTCAGGAGGCCATATCGCGCATCAATCCATTGTCGGttttgt                                                                                                              <  2:657821/247‑5 (MQ=255)

                                                                                                                                       gAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATAGTAACGTGAGATGGGTCAGGAGGAAATATAGCGCATCAAGCATTTGGCGGTTGGGTg                                                                                                             >  2:620537/1‑251 (MQ=255)

                                                                                                                                        aGCATCGGACGCAACAATGCCACCACGCGGGGGTGACCAGAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTGCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGc                                                                                                            <  2:5081/251‑1 (MQ=255)

                                                                                                                                          cATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGGGCgg                                                                                                          >  2:319764/1‑250 (MQ=255)

                                                                                                                                           aTCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGACCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCg                                                                                                           >  1:186396/1‑249 (MQ=255)

                                                                                                                                             cggACGCAACAATGCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGAGCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCa                                                                                                       <  2:72751/249‑1 (MQ=255)

                                                                                                                                               gACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGTCTGATCAAGGCCCTCGCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAg                                                                                                         <  2:309828/247‑1 (MQ=255)

                                                                                                                                               gACGCAACAATTCCACCACGCGGGTGTGACCATAAGTAGCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGc                                                                                                        <  1:138029/248‑1 (MQ=255)

                                                                                                                                                 cGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGAGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGGGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGCGCaaaa                                                                                                    >  1:91597/1‑250 (MQ=255)

                                                                                                                                                 cGCAACAATTCCACCACGCGGGTGTGACCATAAGGATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCGAGCCTTTGGCGGTTCGGTGCGAGCAAAAc                                                                                                   <  1:267372/251‑1 (MQ=255)

                                                                                                                                                 cGCAACAAGTCCCCCACTCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAAc                                                                                                   <  2:470632/251‑1 (MQ=255)

                                                                                                                                                  gCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACt                                                                                                  <  1:412983/251‑1 (MQ=255)

                                                                                                                                                   aaacaaTCCCCCCACGCGGGTGTTACCAACAGGATCACTCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTATCGGTTGTCATAGGAAACCTGGCTGATCAAGGCTCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCAGCAAGCCTTTGGCTGTTCGGTGCGAGTAACACtt                                                                                                 <  2:759629/250‑1 (MQ=255)

                                                                                                                                                          tCCACCTCGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCgg                                                                                           <  2:872060/250‑1 (MQ=255)

                                                                                                                                                          tCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCgg                                                                                           <  2:653799/250‑1 (MQ=255)

                                                                                                                                                              ccacGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTaaa                                                                                      <  2:321474/251‑1 (MQ=255)

                                                                                                                                                               cacGCGGGTGTTACCATAAGTATCAAGCAAAGAAAGGAAGGAGGTATCGCGCAATAAGCGATCAAGAGCCGGAAGTTGACTATAGAGGAAACGCGTGTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCCCGGATTGTAACGTGAGATGGGTCAGGAGGACATATCTCGCATCAAGCCATTGGCGGTTCGGTGCGAGCAAAACTTTCGCGTCTAAAc                                                                                     <  2:764162/251‑1 (MQ=255)

                                                                                                                                                                 cgcgTGTGTGACCAGAAGGAGCACGCAAAGAAAGGAAGGAGCTATAGCGCAATAAGGGATCAATAGCTGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAg                                                                                   <  2:720236/251‑1 (MQ=255)

                                                                                                                                                                   cgGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGt                                                                                 <  1:571601/251‑1 (MQ=255)

                                                                                                                                                                      gtgtGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGtttt                                                                              <  1:498711/251‑1 (MQ=255)

                                                                                                                                                                        gtgcccATAAGTAGCACGCAAAGAAAGGAAGGAGATATCGCGCAATATGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTAGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATACCCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCg                                                                            <  1:451498/247‑1 (MQ=255)

                                                                                                                                                                             cATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAg                                                                       <  1:330343/251‑1 (MQ=255)

                                                                                                                                                                              aTAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCGCACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTGGGGGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTGCGACCAGt                                                                      <  1:39570/251‑1 (MQ=255)

                                                                                                                                                                               tAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATCGGAAACCTGGCTGATCCAGGCCCTCTCACACGGAGAAGGCCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTATGGCGGTTCTGTTGGAGCCAAACCTTCGACGCCAAACAGGTTTTTGAccgagt                                                                     >  2:297387/1‑247 (MQ=255)

                                                                                                                                                                                aaGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGttt                                                                    <  2:360103/251‑1 (MQ=255)

                                                                                                                                                                                 aGTATCACGCCAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCACTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCACATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACATTCGCGGCAAAACCGGTTTTCGACCAGtttt                                                                   >  2:176318/1‑250 (MQ=255)

                                                                                                                                                                                       aCGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGCCCAGTTTGACTa                                                               >  1:704094/1‑249 (MQ=255)

                                                                                                                                                                                       aCGCAAAGAAAGGAAGGAGATATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGGTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAg                                                             <  2:218016/251‑1 (MQ=255)

                                                                                                                                                                                        cGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAgt                                                            <  1:512135/251‑1 (MQ=255)

                                                                                                                                                                                         gCAAAGAAAGGAAGGAGCTATCGAGCAATAAGCGATCAATAGCCGGACGTTGACTATAGATGCAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGTGAAGAGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAATCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACTAGTTTGACTAAgtg                                                           <  2:80224/251‑1 (MQ=255)

                                                                                                                                                                                            aagaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTgg                                                        <  1:315832/251‑1 (MQ=255)

                                                                                                                                                                                              gaaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGCAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTAACACGGAGAAGGGCGTTTAACATAACCACGGCTTGTAACGTGAGATGGGTCCGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCTAGCAAAACTTTCGCGGCTATACAGGTTTTCGACCAGGTTGACTTAGTGggggg                                                      >  2:417673/1‑250 (MQ=255)

                                                                                                                                                                                                 aggaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTGCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGTGCGGt                                                   <  1:718814/251‑1 (MQ=255)

                                                                                                                                                                                                   gaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCg                                                 <  1:290482/251‑1 (MQ=255)

                                                                                                                                                                                                    aaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAAAAGGGTTTCGAACAGGTTGGCCAAGGGTGGGCGGTCgg                                                >  2:621591/1‑250 (MQ=255)

                                                                                                                                                                                                     aggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCATGTTGACTAAGTGTGGGCGGGCGac                                               >  1:433661/1‑251 (MQ=255)

                                                                                                                                                                                                        agcaaTCGCGCAATAAGCGATCCATAGCCGGAAGTTGACTATAGAGGAAACGCTTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGTGATGGGTCAGGAGGACATATCGCGCATCAATCCTTTGGCGGTTCGGTGCTAGCAAAACTTTCGCGGCTAAACAGGTTTTCAACCAGTTTGACTAAGTGTGGGCGGTCTacaca                                            <  2:477517/247‑1 (MQ=255)

                                                                                                                                                                                                          cTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGTCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGTCGGTTCGGGGCGAGCAAAACTTTAGCGGATAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACAcc                                          >  1:174917/1‑251 (MQ=255)

                                                                                                                                                                                                          cTATCGCGCAATAAGCGATCAATAGACGGAAGTTGACTATAGCGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACAcc                                          >  1:704112/1‑251 (MQ=255)

                                                                                                                                                                                                            aTCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACAcc                                          <  1:101392/249‑1 (MQ=255)

                                                                                                                                                                                                            aTCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑aa                                       <  1:779234/251‑1 (MQ=255)

                                                                                                                                                                                                            aTCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑a                                        <  2:15661/250‑1 (MQ=255)

                                                                                                                                                                                                              cgcgCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AAc                                      <  1:405318/250‑1 (MQ=255)

                                                                                                                                                                                                                cgcAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AACCTg                                   <  1:637191/251‑1 (MQ=255)

                                                                                                                                                                                                                  cAATAAGCGATCAATAGCCGGAAGTGGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AACCTGgc                                 <  2:408919/251‑1 (MQ=255)

                                                                                                                                                                                                                       aGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGACACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCCGTTTGACCACGTGTGGGCGGTCGAC‑CACCCAACCTGGCGCAAc                            >  2:401314/1‑251 (MQ=255)

                                                                                                                                                                                                                              aaTAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGGGTGGGCGGTCGACACACC‑AACCTGGCGCAACGcggcgg                     >  1:41334/1‑251 (MQ=255)

                                                                                                                                                                                                                                 aGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AACCTGGCGCAACGCGGCGGaaa                  <  1:781126/251‑1 (MQ=255)

                                                                                                                                                                                                                                    cGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTg                                                                                                             >  2:802112/1‑158 (MQ=255)

                                                                                                                                                                                                                                    cGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AACCTGGCGCAACGCGGCGGAAATTg               <  1:303074/251‑1 (MQ=255)

                                                                                                                                                                                                                                             aCTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGAGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCACGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTAGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AACCTGGCGCAACGCGGCGGAACTTGAGATATCCg      >  1:514688/1‑251 (MQ=255)

                                                                                                                                                                                                                                               tataGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AAACTCGCGCAACGCGGCGGAAATTGAGATGTCCgcc    >  2:264097/1‑249 (MQ=255)

                                                                                                                                                                                                                                                 taGAGGAAACGCGTTTCGGTTGTCATAGGACACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAAGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACgg  >  1:249796/1‑251 (MQ=255)

GGTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACC‑AACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGG  >  NC_000913/3760983‑3761474

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