BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,621,090	C→T	R633R (CGC→CGT)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,621,090	0	C	T	96.2%	287.2 / NA	157	R633R (CGC→CGT)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (85/66); minor base G (2/4); total (87/70)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCG  >  NC_000913/3620844‑3621334

ttATGACACTGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGAGGTATGCATACAACATCGCCGGTGACCGGCCCGCCGTCACTGCCCCGGCCGGCAGAAGAAACGGGACCCAGCACGATGCGTGGCGAAAGTCCGTCAGTACACCGCAGGGCGGGCGGACGCGGAGt                                                                                                                                                                                                                                                    >  2:205659/1‑250 (MQ=11)

 tATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCgtagta                                                                                                                                                                                                                                                   <  1:506663/250‑1 (MQ=14)

       catgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAAt                                                                                                                                                                                                                                             >  1:307254/1‑250 (MQ=14)

        atgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGT‑CCGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGaacac                                                                                                                                                                                                                                           >  2:510338/1‑248 (MQ=14)

         tgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCAGAGT‑CAGTACCGCGCATACGACAGACGTGGAAAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCCGCAGAAACGGGACACAGTACGATGCGAGGGGAAAGGCCGTCCGTAACACGCAGGGCGGGCTGACGCGTAGTATGGAAGACg                                                                                                                                                                                                                                          >  1:347204/1‑251 (MQ=11)

          gaccGCTTCGGGCAGATGCCGGCGGTGCACCGCGCGGACGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTCATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGACGTCCGGACCACGCAGGGCGGGCTGACGCGTAGTATGGAAAACg                                                                                                                                                                                                                                          >  2:260408/1‑250 (MQ=11)

            ccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCGGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACAGCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtg                                                                                                                                                                                                                                       <  1:548530/251‑1 (MQ=14)

             cGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGACGTCCGTACCACGCAGGGCGGGCTGCCGCGTAGTATGGAATACGAtgc                                                                                                                                                                                                                                      >  1:81316/1‑251 (MQ=14)

               cTTCGGGCAGATGACGGCGGTGCACCGCGAGGCAGGGCTGAGT‑CCGTACCGCGCATACGCCAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGGCTGAATACAACATCGCCGGAGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTCCCCCGCAGGGCGGGCTGACGCGGAGTATGGGATACGAtgctg                                                                                                                                                                                                                                    >  2:210104/1‑251 (MQ=11)

                      cAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAATACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCTTCATTGCCCCGGACGTCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGTCCGTCCATACCACGCCGGTCGTGCTGACGCTTATTATTGAATACGAttcttcctgac                                                                                                                                                                                                                              >  1:310970/1‑240 (MQ=11)

                      cAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAAT‑CAGTCCCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACg                                                                                                                                                                                                                             <  1:390674/251‑1 (MQ=14)

                           gagggcggTGCACCGGGAGGAAGGGATGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAACCGCGGTATGAATCCAACACCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCa                                                                                                                                                                                                                        <  2:39693/248‑1 (MQ=11)

                              ggcggTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGAAGGGTCATcc                                                                                                                                                                                                                     >  2:553088/1‑251 (MQ=14)

                               gcggTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGGGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCg                                                                                                                                                                                                                    <  1:86752/251‑1 (MQ=14)

                                   tGCACCGCGAGGAAGGGCTGAGC‑CCGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGGATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg                                                                                                                                                                                                                >  2:105632/1‑251 (MQ=9)

                                     cACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAc                                                                                                                                                                                                              >  1:70224/1‑251 (MQ=14)

                                     cACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAc                                                                                                                                                                                                              <  1:28567/251‑1 (MQ=14)

                                         gcgAGGAAGGGCTGAGT‑CAGATCCTCGCATACGACAGCCGTGGACAGTTAAGTTCCGTGAAAGACACGAAGGGCCATGAAACGCGGTATGAATACAACAGCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCGGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTTACCAGt                                                                                                                                                                                                          <  2:356483/251‑1 (MQ=11)

                                         gcgAGGAAGGGCTG‑GTGCAGCACCGCGCATACGACAGCCGTGGCCAGTTAATTGCCGTGAAAGACACGTAGGGCCAGGAAACGCGGTATGAATACACCCTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACTGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGt                                                                                                                                                                                                          <  2:122103/251‑1 (MQ=11)

                                                aGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg                                                                                                                                                                                                   >  1:400239/1‑251 (MQ=14)

                                                      gAGT‑CAGTCCCGCGCATACGACAGCCGTGGACAGTTAATGGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGGGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGAGGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAg                                                                                                                                                                                               <  1:438037/249‑1 (MQ=9)

                                                         gccAGTACCGCGCATACGACAGCCGTGGCAAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGCCCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCcac                                                                                                                                                                                           <  2:407603/249‑1 (MQ=9)

                                                            aGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcacc                                                                                                                                                                                        >  1:292958/1‑251 (MQ=14)

                                                                ccGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCt                                                                                                                                                                                    >  1:223859/1‑251 (MQ=14)

                                                                    ccctACGCCAGCCGTGGACAGTTAATTGCCTGGAAAGACCCGCAGGGCCATGAAACGCGGTAGGAATCCAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCg                                                                                                                                                                                <  2:420367/248‑1 (MQ=11)

                                                                      aTACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGtt                                                                                                                                                                              <  1:363438/251‑1 (MQ=14)

                                                                            cAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATg                                                                                                                                                                        >  1:469385/1‑251 (MQ=14)

                                                                             aGCCGTGGCCAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCCGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGt                                                                                                                                                                       >  2:174016/1‑251 (MQ=14)

                                                                              gCCGTGGACAGTTAATGGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTa                                                                                                                                                                      <  2:475489/251‑1 (MQ=14)

                                                                               ccGTGGACAGTTAATTGCCGGGAAAGACACGCGGGTCCAGGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAAGACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTAc                                                                                                                                                                     <  2:194273/251‑1 (MQ=255)

                                                                                cGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAACCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACt                                                                                                                                                                    <  1:404221/251‑1 (MQ=14)

                                                                                      cAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGAACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGAGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcacc                                                                                                                                                                                        >  2:305100/1‑225 (MQ=9)

                                                                                      cAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCACCGGACGGCAGCAGAAACGGGACACAGTACGATGCGAGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGACGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCCCACCACCTTCCGTTACGATGTACTTGACCg                                                                                                                                                              >  1:491817/1‑251 (MQ=9)

                                                                                           aaTTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTg                                                                                                                                                          >  1:564930/1‑250 (MQ=14)

                                                                                                ccGTGAAAGGCACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGGGACCCGACCGCCGTCATTGCCCCGGACGGCAGCAGAAGCGGGACACAGTACGAGGCGTGGGTAAAGGCTGTCCGTATCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAg                                                                                                                                                    <  1:41911/251‑1 (MQ=11)

                                                                                                ccGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACCCCACCTTCCGTTACGATGTACTTGACCGGCTGATAc                                                                                                                                                      >  2:294839/1‑249 (MQ=14)

                                                                                                 cGTGAAAGACACGCAGGGCCATGAAACGCCGTATGAATCCAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACAAAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAgg                                                                                                                                                   >  2:546286/1‑251 (MQ=9)

                                                                                                   ggAAAGACACGCACGGCCCTGAAACTCGGTATGAATACAACAGCGCCGGGGACCGGACCGCCGGCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGAGTCGTGGGGAAAGGCCGTCCGTTCCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATTCTGCCGGACGGGTCATCCGCCTTACCAGGGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACGTTACCGGCTGATACAGGaa                                                                                                                                                 <  2:540540/250‑1 (MQ=11)

                                                                                                     aaaGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAAc                                                                                                                                               <  1:525593/251‑1 (MQ=14)

                                                                                                     aaaGACACGCAGGACCATGAACCGCGTTATGAACACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCCCGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAAc                                                                                                                                               <  2:227797/251‑1 (MQ=255)

                                                                                                      aaGACCCGCAGGGCCAGGAACCGCGTTATGAATACAACATCGGCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAAcc                                                                                                                                              <  1:309721/251‑1 (MQ=255)

                                                                                                      aaGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAAcc                                                                                                                                              <  1:93265/251‑1 (MQ=14)

                                                                                                        gACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCgg                                                                                                                                            >  1:312882/1‑251 (MQ=14)

                                                                                                             gaaggtccaggAAACCCGGTATGAATTCAACATCGCCGGTGACCCGACCGCCGCGACTGCCCCGGACGGCAGCAGCAACGGGACTCAGTACGATGCGTGGGGAAACGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAGGCTGCCCTACGGGTCAGCCGCCGGACCATTGAAAACGGCAGCCACACCACCTTCCGTTAAGATGTTCTTGACCGTCTGATACAGGAAACCGGCTTTg                                                                                                                                       <  1:259198/241‑1 (MQ=255)

                                                                                                             gCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGAAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCt                                                                                                                                                           >  1:173687/1‑231 (MQ=14)

                                                                                                               aGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTgcc                                                                                                                                     >  1:284414/1‑249 (MQ=14)

                                                                                                                gggCCAGGAAACGCGGTATGAATACACCACCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACg                                                                                                                                    <  2:174011/251‑1 (MQ=9)

                                                                                                                  gCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGc                                                                                                                                  >  1:536510/1‑251 (MQ=14)

                                                                                                                          aCGCGGTATGAATACAACATCGCCGGTGCCCGGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcacaca                                                                                                                          <  2:372633/251‑1 (MQ=14)

                                                                                                                           cgcgGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑g                                                                                                                        <  1:10106/251‑1 (MQ=14)

                                                                                                                             cgGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCg                                                                                                                      >  1:105203/1‑251 (MQ=14)

                                                                                                                             cgGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGAAACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGCCGGGTCATCCGCCTGACCAGTGGAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCg                                                                                                                      >  1:384553/1‑251 (MQ=9)

                                                                                                                                 atgaaaaaaacaCCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCTTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTTATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTAt                                                                                                                  <  2:583021/243‑1 (MQ=255)

                                                                                                                                 aTGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTAt                                                                                                                  >  1:114348/1‑251 (MQ=14)

                                                                                                                                   gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATc                                                                                                                 >  1:303697/1‑250 (MQ=14)

                                                                                                                                   gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATc                                                                                                                 >  1:15261/1‑250 (MQ=14)

                                                                                                                                       acaacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATcacca                                                                                                             <  1:474070/250‑1 (MQ=14)

                                                                                                                                       acaacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATcacca                                                                                                             <  1:184531/250‑1 (MQ=14)

                                                                                                                                         aacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATcaccccg                                                                                                           >  2:163153/1‑247 (MQ=14)

                                                                                                                                         aacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCAGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCAAACA‑GCGTTATCACCACg                                                                                                           >  1:341828/1‑250 (MQ=14)

                                                                                                                                              cGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACAAA‑GCGTTATAACCCCGACCTgc                                                                                                     >  2:459463/1‑250 (MQ=9)

                                                                                                                                               gCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGAc                                                                                                    <  1:200759/251‑1 (MQ=14)

                                                                                                                                               gCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCAC‑CACGCGGTATCACCACGACCTGAc                                                                                                    >  2:437459/1‑251 (MQ=9)

                                                                                                                                                ccGGTGCCCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATTGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCTTTATCACCACGACCTGAcc                                                                                                   <  1:210104/251‑1 (MQ=9)

                                                                                                                                                    ggacctgacgGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATTCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTGCCGTTACGATGTACTTGACCGGCTGATACAGTAAACCTGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCa                                                                                               <  1:255963/250‑1 (MQ=255)

                                                                                                                                                     gacctgacaGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGAAACCGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGCATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGACAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCCTTGAGGGCCGCCCACA‑GCGTTATCACCACGACCTGACCGGCa                                                                                               >  2:369146/1‑250 (MQ=11)

                                                                                                                                                          gaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACt                                                                                           <  1:172576/249‑1 (MQ=14)

                                                                                                                                                          gaccGCCGTCATGGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTa                                                                                         <  1:211786/251‑1 (MQ=14)

                                                                                                                                                           accGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATCCAGGAAACCGGCTTTGACGGCCGCACAAC‑TCGTTCTCACCACGACCTGACCGGCAAACTTAt                                                                                        >  2:171952/1‑251 (MQ=255)

                                                                                                                                                             cgccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATcc                                                                                      <  1:496063/251‑1 (MQ=14)

                                                                                                                                                             cgccgTCAGTGCCACGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCGGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATcc                                                                                      <  1:116027/251‑1 (MQ=9)

                                                                                                                                                                 gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGTCCGTCCGTACCACGCAGGGCGGGTTGACGCGTAGTATGGAATGCGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCTCCTTCCGGTACGATGTGCTGGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAg                                                                                  <  1:291816/251‑1 (MQ=11)

                                                                                                                                                                 gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGGAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAg                                                                                  <  2:563739/251‑1 (MQ=14)

                                                                                                                                                                  tCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATCCGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGACGCACACA‑GCGTTATCACCCCGAACAGCCCGGCAAACTTATCCGCAGc                                                                                 >  2:408485/1‑251 (MQ=11)

                                                                                                                                                                              aCGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGt                                                                     >  1:349465/1‑251 (MQ=14)

                                                                                                                                                                              aCGGCAACAGAAACGGGACACATTGCGCTGCGGGGGTAAAGGCCTTTCGTACCACGCAGGGCGGGCTGACGCGTAGTATGTAATACGATGCTGCCGGACGGGTCATCCGCCTGACCATGGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGTAATCCGGCTTTGACGGCCGCACACA‑GCGGTATAACTACGACCTGACCGGCAAACTTATCCGCAGCGAGTATTAGGGt                                                                     <  2:562299/251‑1 (MQ=11)

                                                                                                                                                                                  cagcagAAACGGGACCCAGTACGATGCGTGGGGAAAGGCCGTCCGTACCTCGCAGGGCGGGCTGACGCGTCGTATGGAATACGATGCTGCCGGAAGGGTCATCCGCCTGCCCAGTGAAAACGGCAGCCACACCACCTTCCGTTTCGCTTTTCTTTCCCTCCTTTTCCCTGTCTCCGGCTTTGACGGCCGCACACA‑GCCTTATCTCCACGACCTGACCGGCAAACTTATCCGCTGCgttgttgaccgtctgg                                                                 <  2:432130/251‑16 (MQ=2)

                                                                                                                                                                                    gcagAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTc                                                               <  1:425625/251‑1 (MQ=14)

                                                                                                                                                                                           cGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGTAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTGGACGGCCGCACTCA‑GCGTTATCTCCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACt                                                        <  1:307523/251‑1 (MQ=9)

                                                                                                                                                                                            gggACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCgg                                                                                                                                            >  2:34862/1‑167 (MQ=17)

                                                                                                                                                                                            gggACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCgg                                                                                                                                            <  1:34862/167‑1 (MQ=17)

                                                                                                                                                                                                cacaGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGACATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCCCCCACTGGCAc                                                   >  1:26528/1‑251 (MQ=14)

                                                                                                                                                                                                   aGTACGATGCGGGGGGAAAGGCCGGCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTAt                                                <  2:14460/251‑1 (MQ=14)

                                                                                                                                                                                                      aCGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCGGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCGCCACGACCTGACCGGCAAGCTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCCCGGGCACTATGAc                                             <  2:502356/251‑1 (MQ=11)

                                                                                                                                                                                                       cGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACg                                            <  1:529395/251‑1 (MQ=14)

                                                                                                                                                                                                         aTGCGTGGGGAAAGGCCGTCCGTCCCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGCTATCACCACGACCTGACCGGCCAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGa                                           >  2:473410/1‑250 (MQ=9)

                                                                                                                                                                                                         aTGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCATGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGCCGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAACCTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGa                                           >  1:95254/1‑250 (MQ=9)

                                                                                                                                                                                                                   aaaGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGcc                                >  1:167179/1‑251 (MQ=14)

                                                                                                                                                                                                                     aGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGGAGGATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTGGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACGGGCACTATGACGAAGCAGACCGCCTc                              <  2:241594/251‑1 (MQ=9)

                                                                                                                                                                                                                      ggCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCa                             <  1:174016/251‑1 (MQ=14)

                                                                                                                                                                                                                              gtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgcaccgc                     >  1:409707/1‑251 (MQ=14)

                                                                                                                                                                                                                                   aCGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGt                >  2:184163/1‑251 (MQ=14)

                                                                                                                                                                                                                                      cAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGCCATCCGCCTGACCAGTGAAAAAGGCAGCCACACAACCTTCAGTTACGATGTACTTGACCAGCTGATACAGGAAACCGGCTTTGACGGCCGCAAACA‑GCGTTATCACCAAGACCTGACAGGCAAACTTATACGCCGCGAGGATGAAGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGa              >  1:357806/1‑250 (MQ=11)

                                                                                                                                                                                                                                       aGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGa              >  1:208699/1‑249 (MQ=14)

                                                                                                                                                                                                                                           cgggcTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGAcc                                                                                                   <  2:208233/160‑1 (MQ=17)

                                                                                                                                                                                                                                               cTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAAc    >  1:102112/1‑251 (MQ=14)

                                                                                                                                                                                                                                               cTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGC‑CTATTACGAAGCCGGCCGCCTCCCGCACCGCCCCGTGAAGGGGGAAAcc   >  1:104910/1‑251 (MQ=11)

                                                                                                                                                                                                                                                tGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAAcc   >  1:569359/1‑251 (MQ=11)

                                                                                                                                                                                                                                                 gACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAAc    <  1:128740/249‑1 (MQ=11)

                                                                                                                                                                                                                                                 gACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGGGAAAACGGCAGACACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGCCGGCCGCACCCA‑GCGTTATCACCACGACCTGACCGGCAAACTTCTACGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCAATATGACGAGACAGACCGCCACACGCACCGAACCGTGAAGGGTGaaaacg  >  2:423973/1‑248 (MQ=11)

TTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGT‑CAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA‑GCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCG  >  NC_000913/3620844‑3621334

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