BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	793,306	C→T	E161E (GAG→GAA)	modF ←	ABC family protein ModF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	793,306	0	C	T	98.1%	340.8 / NA	106	E161E (GAG→GAA)	modF	ABC family protein ModF
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (41/61); minor base G (1/1); total (44/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 106 total aligned reads displayed.

GCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGC  >  NC_000913/793062‑793550

 cAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCaa                                                                                                                                                                                                                                               <  1:312964/251‑1 (MQ=255)

  ccttGCACACCCTCAAGCTGTTCACTCTGCTTCAGTTGCGTGACTAGTGCTTGTTGGAGCAGTTCCTCCTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAATCTGGACAATCTCCGGGATTTCATCGAAGCGACTGAGCACCAGTCCCAGATTAATACCGGACTGATGTATCGATGCGAGCCGCTCCGCCAGCTGCTGACGTGAGTCATCATCCAGCCCATCGAACGGTTCATCAAg                                                                                                                                                                                                                                              <  2:195459/249‑1 (MQ=255)

      gCNCACCTTCAAGCTTGTCACTATGCGCCAGTTTCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACTCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGAGTGAGCACCAGTGCCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCTAGCGGCGGACGTGAGGCAACATTCAGGCCATCGAACCGTTCATCAAGAAt                                                                                                                                                                                                                                           <  2:586112/250‑1 (MQ=255)

         cacCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGGTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATcaac                                                                                                                                                                                                                                       <  2:515575/251‑1 (MQ=255)

             ttCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGt                                                                                                                                                                                                                                   >  2:30657/1‑251 (MQ=255)

              tCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAg                                                                                                                                                                                                                                    <  2:28541/249‑1 (MQ=255)

                     gTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTcc                                                                                                                                                                                                                           >  1:44231/1‑251 (MQ=255)

                         aCTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGAc                                                                                                                                                                                                                        >  1:474425/1‑250 (MQ=255)

                               cgcCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAgcg                                                                                                                                                                                                                 >  2:503208/1‑251 (MQ=255)

                                      tGCGCGACGAGTGCATGTTGGGGCAGTTCCTCTTTAGCGCCAGTTTCCGCGAACGTGCAATCCGCCATCACGCCAGCAAACTGGACGAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCTAACGGTGCATCAAGAATCAACAAGTTAGGCTCCGACAGCAGCGCCTGac                                                                                                                                                                                                           <  2:426708/250‑1 (MQ=255)

                                        cgcgACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGacac                                                                                                                                                                                                         >  1:181441/1‑250 (MQ=255)

                                        cgcgACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGacac                                                                                                                                                                                                         >  2:536970/1‑250 (MQ=255)

                                        cgcgACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGacac                                                                                                                                                                                                         >  2:368425/1‑250 (MQ=255)

                                        cgcgACGAGGGCTTGTTGGAGCAGTTCCTCCTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGAcaca                                                                                                                                                                                                        <  1:466173/251‑1 (MQ=255)

                                         gcgACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGATGCTGACGTGAGGCAACATCCAGGCCATCGAACGGATAATCAAGAATCAACAAGTCAGACTCCGACATCAGCGCCTGacac                                                                                                                                                                                                         >  2:199148/1‑249 (MQ=255)

                                             cgaGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAgg                                                                                                                                                                                                   <  1:460147/251‑1 (MQ=255)

                                               aGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCACTCAGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGt                                                                                                                                                                                                 >  1:77095/1‑251 (MQ=255)

                                                gTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGtt                                                                                                                                                                                                <  1:598220/251‑1 (MQ=255)

                                                           aGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGtctc                                                                                                                                                                                     <  2:418236/251‑1 (MQ=255)

                                                            gCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCg                                                                                                                                                                                    <  1:318054/251‑1 (MQ=255)

                                                               gTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAGCTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAAAGTAATACCGGACTGATGTAACGAGGCGAGTTGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCa                                                                                                                                                                                 <  1:230301/251‑1 (MQ=255)

                                                                  ccTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGt                                                                                                                                                                               <  1:236488/250‑1 (MQ=255)

                                                                  ccTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTg                                                                                                                                                                              <  1:99916/251‑1 (MQ=255)

                                                                     ctTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGa                                                                                                                                                                            >  1:23093/1‑250 (MQ=255)

                                                                       ttAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGaaa                                                                                                                                                                          <  2:474425/250‑1 (MQ=255)

                                                                       ttAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGa                                                                                                                                                                            <  1:480622/248‑1 (MQ=255)

                                                                           cgcCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGACAGTGGAAAGGTAt                                                                                                                                                                     >  1:59057/1‑251 (MQ=255)

                                                                               aGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATcaac                                                                                                                                                                                                                                       <  1:305373/181‑1 (MQ=255)

                                                                               aGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATcaac                                                                                                                                                                                                                                       >  2:305373/1‑181 (MQ=255)

                                                                                    ccGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgt                                                                                                                                                            <  2:503230/251‑1 (MQ=255)

                                                                                     cGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtc                                                                                                                                                           <  1:30657/251‑1 (MQ=255)

                                                                                      gCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtcg                                                                                                                                                          >  1:358803/1‑251 (MQ=255)

                                                                                          aCGTGCAATCCGCCAGCACGCCAGCAAACTGGTCAAACTCCGTGAACTCAGCGTAGCGATTGAGCCCCAGTCCCAGAGTAATTCCGGACCTATGTAATGAGGCGAGTCGCTCAGCCAGCTGCTTACGTTAGGCAACATCCATGCCATCGAACGGTTCTTCAATAATCAACAAGCCAGGCTCCCACATCAGCGCTTGACACCGCATTGTTTTTCTTGTCTCGCCATTGGTAAGGTTTTTAATGCgtccgtcg                                                                                                                                                      <  2:623868/251‑1 (MQ=255)

                                                                                                   ccGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCg                                                                                                                                             >  1:462007/1‑251 (MQ=255)

                                                                                                      ccagcacgccagcaAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGCGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGCACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGACAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGt                                                                                                                                           >  2:432347/1‑250 (MQ=255)

                                                                                                            cgccagcaAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGAACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGGAATNCCg                                                                                                                                    >  2:599647/1‑251 (MQ=255)

                                                                                                            cgccagcaAACTGGACAAACTCCGGGATATCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGATCGAGGAGGGCGGTAATACCg                                                                                                                                    >  1:552126/1‑251 (MQ=255)

                                                                                                                    aaCTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAgg                                                                                                                                                                        <  1:220759/207‑1 (MQ=255)

                                                                                                                    aaCTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAgg                                                                                                                                                                        >  2:220759/1‑207 (MQ=255)

                                                                                                                    aaCTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgctg                                                                                                                            <  1:22406/251‑1 (MQ=255)

                                                                                                                        ggACAACCTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTgcg                                                                                                                          <  2:77095/249‑1 (MQ=255)

                                                                                                                        ggACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGcc                                                                                                                        <  1:482562/251‑1 (MQ=255)

                                                                                                                        ggACAAACTCCAGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCATGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGACCTGCTGCGcc                                                                                                                        >  2:73479/1‑251 (MQ=255)

                                                                                                                            aaaCTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGATGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGGTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGACTGACACAGCAGGGTTTTTAGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTAGAGGAGGGCGGTAATACCGAACTGATGCGCCAGtt                                                                                                                    >  2:354323/1‑251 (MQ=255)

                                                                                                                             aaCTCCGGGATCTCATCGCAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGAAGTGAGGAAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtg                                                                                                                   >  2:422143/1‑251 (MQ=255)

                                                                                                                              aCTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgc                                                                                                                  <  2:516470/251‑1 (MQ=255)

                                                                                                                               cgccGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCGCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgca                                                                                                                 <  2:497968/249‑1 (MQ=255)

                                                                                                                               cTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgc                                                                                                                  >  1:41535/1‑250 (MQ=255)

                                                                                                                                tCCGGGATCTCATCGACGCGATTGAGCAGCAATACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCAAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGCGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcat                                                                                                                >  1:249493/1‑251 (MQ=255)

                                                                                                                                 ccGGGATCTCATCGACGCGATTGAGCACCAGTACCAGAGTAATACCGGACTTATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCCTCGAACGGTTCATCAAACATCAACAAGTCAGGCTCCGACATCAGCGCCTGGCACAGCAGGGTTTTTCGCTTCTCGCCAGCGGCAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGGACTGATGCGCCATTtgcatg                                                                                                               >  2:481320/1‑251 (MQ=255)

                                                                                                                                 ccGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcatg                                                                                                               <  2:39006/251‑1 (MQ=255)

                                                                                                                                 ccGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcat                                                                                                                <  1:625281/250‑1 (MQ=255)

                                                                                                                                   gggATCTAATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGCCTGATGAAACGAGGCGAGTCCCTCAGCCAGCTGCTGACGTGAGGCCACATCCAGGCCTTCGAACGGTTCATCACGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACCGCAGGGTTTTTCGCGTCTCGCCAGTGGAGAGGTATTTAAGGCATCGGTCGCGGAGGCCGGTAATACCGAACTGCTGCGCCAGTtgcatgc                                                                                                              >  2:52341/1‑250 (MQ=255)

                                                                                                                                    ggATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcatnca                                                                                                             >  2:598894/1‑250 (MQ=255)

                                                                                                                                          catagAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGATCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAAGTGCATGCCACGCggg                                                                                                      >  2:520931/5‑250 (MQ=255)

                                                                                                                                          cATCGAAACGATTGAGCACCAGCACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCACAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGt                                                                                                      <  1:199148/251‑1 (MQ=255)

                                                                                                                                           aTCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTg                                                                                                     >  1:105257/1‑251 (MQ=255)

                                                                                                                                                   gATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCCGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTAATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGCATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCttt                                                                                             >  1:401340/1‑251 (MQ=255)

                                                                                                                                                       gAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACtt                                                                                         <  1:516623/251‑1 (MQ=255)

                                                                                                                                                                 aCCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTcg                                                                                                                                                                                           >  2:177041/1‑143 (MQ=255)

                                                                                                                                                                 aCCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTcg                                                                                                                                                                                           <  1:177041/143‑1 (MQ=255)

                                                                                                                                                                  ccAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATg                                                                              >  1:506527/1‑251 (MQ=255)

                                                                                                                                                                  ccAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCGCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATg                                                                              <  2:455398/251‑1 (MQ=255)

                                                                                                                                                                   caggatgaTACCGGTCTGATGTAACGAGGCGAGTCGCTCAGCCCGCTGCTGACGTGAGGCACCATCCAGGCCATCGAACCGTTCATCAAGAATCAACAAGTCAGACTTCGTCATCAGCGCCTGTCGCAGCAGGGTTTTTCGCGTCGCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGATCTGCTGCGCCAGTTGCATGCACCTCGGTGCATCCTTTACTTCAGCCTGAATGa                                                                             <  2:47156/244‑1 (MQ=255)

                                                                                                                                                                   cAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGa                                                                             <  2:283518/251‑1 (MQ=255)

                                                                                                                                                                        tAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATAc                                                                                                                                      <  1:540906/189‑1 (MQ=255)

                                                                                                                                                                        tAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATAc                                                                                                                                      >  2:540906/1‑189 (MQ=255)

                                                                                                                                                                         aaTACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAg                                                                       <  2:279067/251‑1 (MQ=255)

                                                                                                                                                                            aCCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCg                                                                    <  1:315056/251‑1 (MQ=255)

                                                                                                                                                                             ccGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGt                                                                   >  1:11756/1‑251 (MQ=255)

                                                                                                                                                                               ggACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAg                                                                 <  1:332887/251‑1 (MQ=255)

                                                                                                                                                                               ggACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAg                                                                 <  2:105257/251‑1 (MQ=255)

                                                                                                                                                                                 aCTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                               <  1:164883/251‑1 (MQ=255)

                                                                                                                                                                                       gTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACt                                                                                          >  2:150659/1‑218 (MQ=255)

                                                                                                                                                                                       gTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACt                                                                                          <  1:150659/218‑1 (MQ=255)

                                                                                                                                                                                         aaCGAGGCGAGTCGCGTAGACAGCTGCTGACGTGAGGCAACATCCAGGACATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCTAGTTGCATGCAACGCGGTGCATCCTTGACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCTGgt                                                       <  2:181441/251‑1 (MQ=255)

                                                                                                                                                                                            gAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTc                                                    <  1:93121/251‑1 (MQ=255)

                                                                                                                                                                                            gAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGGCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTc                                                    <  2:491467/251‑1 (MQ=255)

                                                                                                                                                                                             aGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGAGCGTCCGGTGTc                                                    >  2:280538/1‑250 (MQ=255)

                                                                                                                                                                                                cGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCt                                                <  2:640584/251‑1 (MQ=255)

                                                                                                                                                                                                  aGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTc                                              <  2:23093/251‑1 (MQ=255)

                                                                                                                                                                                                     cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGcc                                           <  2:171817/251‑1 (MQ=255)

                                                                                                                                                                                                       cTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGcc                                           >  2:309551/1‑249 (MQ=255)

                                                                                                                                                                                                         cagccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAg                                         <  2:377314/249‑1 (MQ=255)

                                                                                                                                                                                                           gccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGcc                                     <  1:240024/251‑1 (MQ=255)

                                                                                                                                                                                                               gatgctgACGGGAGGCACCATGCAGGCAAACGAAGGGGTAATCAATAAAAAAAAAGTCAGGATCCGACAGCAGCGCCTGACACAGCATGGTGTTTCGCGTCTCGCCAGTGGAAAGGTATTTCAAGCGGCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGc                                 <  2:332350/249‑1 (MQ=255)

                                                                                                                                                                                                                 ggctgACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCAGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTACGTCCGGTGTCATCTTCGCCAGGGCCGAGCat                               <  2:197693/250‑1 (MQ=255)

                                                                                                                                                                                                                  gctgACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCat                               >  1:246541/1‑250 (MQ=255)

                                                                                                                                                                                                                  gctgACGTGAGGCAACAGCCAGGCCATCGAACGGTTCAGCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAGGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCata                              <  2:103514/251‑1 (MQ=255)

                                                                                                                                                                                                                     gACGTGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATc                            <  1:598211/250‑1 (MQ=255)

                                                                                                                                                                                                                         tGAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCggcgt                       >  1:607621/1‑248 (MQ=255)

                                                                                                                                                                                                                          gAGGCAACATCCAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTTCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtgt                       <  1:563070/250‑1 (MQ=255)

                                                                                                                                                                                                                                    ccAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCCGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTaa                                                                                                                                                                 >  1:593288/1‑102 (MQ=255)

                                                                                                                                                                                                                                    ccAGGCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTaa                                                                                                                                                                 <  2:593288/102‑1 (MQ=255)

                                                                                                                                                                                                                                       ggCCATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGGCGGTCGAGGAGGGCGGTAAGACCGAACGGCTGCGCCAGGTGCATGCAACGCGGTGTATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGATCATATCGGTGTTAGTCTGCTGCCAt         <  2:41535/251‑1 (MQ=255)

                                                                                                                                                                                                                                         ccATCGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtc       <  1:133773/251‑1 (MQ=255)

                                                                                                                                                                                                                                             cGAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGGGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtc    <  1:109522/250‑1 (MQ=255)

                                                                                                                                                                                                                                              gAACGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGc  >  2:464942/1‑251 (MQ=255)

                                                                                                                                                                                                                                               aaCGGTTCATCAAGCCTCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGCACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGTGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAt         >  2:254392/1‑243 (MQ=255)

                                                                                                                                                                                                                                               aaCGGTTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATGTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGGGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAt         <  1:254392/243‑1 (MQ=255)

GCAGTTGCACACCTTCAAGCTGTTCACTATGCGCCAGTTGCGCGACGAGTGCTTGTTGGAGCAGTTCCTCTTTAGCGCCAGTTTCCGCTAACGTGCAATCCGCCAGCACGCCAGCAAACTGGACAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGC  >  NC_000913/793062‑793550

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