BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,405	A→G	T311T (ACT→ACC)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,405	0	A	G	98.2%	347.6 / ‑4.3	110	T311T (ACT→ACC)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (45/61); total (47/63)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 110 total aligned reads displayed.

CGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCG  >  NC_000913/67158‑67648

cGACAATCGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGtt                                                                                                                                                                                                                                                  >  2:498733/1‑251 (MQ=255)

     aTCGCTGGACCGGTTTGGGAATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCCTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGGGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTcc                                                                                                                                                                                                                                              >  2:281478/1‑250 (MQ=255)

       cGCTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGGTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGGGCCGCCCTGCATACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGAGGTTTTCCAGt                                                                                                                                                                                                                                           >  2:18783/1‑251 (MQ=255)

        gcttggCCGGTTTGGGTATTTAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGGTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGCGCACCTGGACATTACCTTTCTCGAAGTGATAGGTGTAGGCCTCCATAAAGGAGTTGCCGCCCTGCTGAACGGTTGAAATCACCTTCATGCTTCGAGGCGGGGCGGAGGTTTTACAGTc                                                                                                                                                                                                                                          >  2:22423/5‑251 (MQ=255)

         cTGGACCGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATACGAAGCAGGGCGGCGGGTTTCCAgacg                                                                                                                                                                                                                                         >  2:482941/1‑248 (MQ=255)

              ccGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAAGCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTc                                                                                                                                                                                                                                    <  1:252470/251‑1 (MQ=255)

               cgggttGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTGCTCTGCTACGGCAAGCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACAGCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCg                                                                                                                                                                                                                                   <  2:179960/247‑1 (MQ=255)

               cGGTTTGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCg                                                                                                                                                                                                                                   <  1:521982/251‑1 (MQ=255)

               cGGTTTGGGTAGTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCg                                                                                                                                                                                                                                   <  2:614736/251‑1 (MQ=255)

                   ttGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCg                                                                                                                                                                                                                               >  1:645998/1‑251 (MQ=255)

                    tGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGCCGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGc                                                                                                                                                                                                                              >  1:51571/1‑251 (MQ=255)

                    tGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGc                                                                                                                                                                                                                              >  1:258273/1‑251 (MQ=255)

                    tGGGTATTGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGc                                                                                                                                                                                                                              >  1:373143/1‑251 (MQ=255)

                     gggTATTTAAGCTCAGTCGGGCCGGATTGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAACACCAGGTCATTAACTTTCTCGAGGTGATCGATGTACTACTCCCTAAAGGAGGTGCCGACCTGCAGCCCGGTTGACATCACCTTCATAGTGCGAAACAGCTCTGCAGTTTTCACGTCGCCTTCGCCCGc                                                                                                                                                                                                                              >  2:625593/1‑250 (MQ=255)

                     gggTATTGAAGATCAGTCGGGCAGGCTCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAAACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTCACATTCGCCTTCCCCCGc                                                                                                                                                                                                                              >  2:699692/1‑250 (MQ=255)

                      ggTATTGAAGAGCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCaa                                                                                                                                                                                                                            <  1:73815/251‑1 (MQ=255)

                           tGAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTTCCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGcc                                                                                                                                                                                                                        <  2:432032/250‑1 (MQ=255)

                            gAAGATCAGTCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCGTCGCCCGCAAAGCCg                                                                                                                                                                                                                       <  1:508549/250‑1 (MQ=255)

                              agctCAGTCGGGCAGGTTCGTCCTTTCCACCAATACCGGGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATTGACGGGCAGACTTTCAGCATATTGGAGTCCAGCCCCAGGTCATTTCCTTTCGTGGAGTGATAGGTGTAGACCTCCATAAAGGATGTGCCGCCCGGTAGTCCGGTTGACATCACTTTCATGATGCGATGCAGGGCGGCGGTGCTCCATTCGCCTTCTGCCGCTAAGCCGTaa                                                                                                                                                                                                                    <  2:395886/248‑1 (MQ=255)

                                  cagtcGGGCAGGATCGTCCTTACAACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCGGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGGAAcccc                                                                                                                                                                                                                >  2:146945/5‑250 (MQ=255)

                                    gtcGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCctgc                                                                                                                                                                                                              <  1:157323/249‑1 (MQ=255)

                                    gtcGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCctgc                                                                                                                                                                                                              <  1:190868/249‑1 (MQ=255)

                                      cGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCctgctg                                                                                                                                                                                                            <  1:166779/251‑1 (MQ=255)

                                      cGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATAGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCctgctg                                                                                                                                                                                                            >  1:704344/1‑251 (MQ=255)

                                                   ccgtACCCCCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTTCCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTgg                                                                                                                                                                                               <  2:234604/248‑2 (MQ=255)

                                                      tccccccaAGACCGAGATGCTGAACGTCTAGGATCGGTTTCTTATCTACGGCAATCGACGGGCCGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTa                                                                                                                                                                                              <  2:568041/244‑1 (MQ=255)

                                                          cccaATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCa                                                                                                                                                                                        <  1:266684/250‑1 (MQ=255)

                                                          accaATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGcc                                                                                                                                                                                         >  1:595258/1‑250 (MQ=255)

                                                            caATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCATACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGa                                                                                                                                                                                      <  1:61890/251‑1 (MQ=255)

                                                             aaTACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGAc                                                                                                                                                                                     >  1:652581/1‑251 (MQ=255)

                                                              aTACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGAcc                                                                                                                                                                                    <  1:282663/251‑1 (MQ=255)

                                                                  cGAGATGGTGACCGTCGTGGATCGGTTACTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCTAGCACCAGGTCATTGCCTGTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCGCCATCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGGAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGa                                                                                                                                                                                <  2:656407/251‑1 (MQ=255)

                                                                        gCTGAACGTCGAGGATCGGTTTCTCGTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTg                                                                                                                                                                           <  2:543186/250‑1 (MQ=255)

                                                                           gAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGttt                                                                                                                                                                        <  1:483272/250‑1 (MQ=255)

                                                                           gAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGttt                                                                                                                                                                        <  1:547107/250‑1 (MQ=255)

                                                                           gAACGTCGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGtt                                                                                                                                                                         <  1:486977/249‑1 (MQ=255)

                                                                            aaCGTCGCGGATCGGTTTCTGTTCTACGGCAATCGACGGGCAGACTTCCAGAATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCTCCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGTCCAGACCTTGAAGCTGTTTCa                                                                                                                                                                      <  2:227312/251‑1 (MQ=255)

                                                                              cGTCGAGGATCGGTTTCTCTTCTACGGCAATCTACGGGCAGACTTCCATCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCTAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATTCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTATCCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGa                                                                                                                                                                    <  1:670629/251‑1 (MQ=255)

                                                                                 cGAGGATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCg                                                                                                                                                                 <  1:514711/251‑1 (MQ=255)

                                                                                    ggATCGGTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCTTGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCCGACGCTGTACGGCCAGCCCTGGGAGCTGGTTCAGACCGTGc                                                                                                                                                              >  2:62477/1‑251 (MQ=255)

                                                                                         ggTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATc                                                                                                                                                         >  1:224454/1‑251 (MQ=255)

                                                                                         ggTTTCTCTTCTACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATc                                                                                                                                                         <  1:531685/251‑1 (MQ=255)

                                                                                              gcccactactgcaaTCGACGGGCAGACTGCCAACATATGGGAGCCGAGCACCAGTGCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGATATCCCCTTCAGGATGCGAAGCAGGGCGGCGGTTTTCCAGCCGCCTTCGCCCGCAAAGCCGTATCCCTGCTGCATCATACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAAACTTCa                                                                                                                                                     <  2:288936/237‑1 (MQ=255)

                                                                                               tcttctACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCaa                                                                                                                                                    >  1:165335/1‑250 (MQ=255)

                                                                                                cttctACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCaa                                                                                                                                                    >  1:160168/1‑249 (MQ=255)

                                                                                                cttctACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCa                                                                                                                                                     <  1:116081/248‑1 (MQ=255)

                                                                                                 ttctACGGCAATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATGCCCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTGCAAAg                                                                                                                                                  <  2:46218/250‑1 (MQ=255)

                                                                                                  tctacggcaATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCaaa                                                                                                                                                   <  1:482941/244‑1 (MQ=255)

                                                                                                       ggcaaTCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAggtggtgg                                                                                                                                           >  1:660688/5‑251 (MQ=255)

                                                                                                              gACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCCAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCCGACCGTGCAAATCTTCAAAGGTGGTGGGGAAcgg                                                                                                                                    >  2:560468/1‑250 (MQ=255)

                                                                                                               aCGGGCAGACTTCCAGCAGATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAAcgcg                                                                                                                                   <  1:467533/251‑1 (MQ=255)

                                                                                                                    cagtcTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTGCTCGAAGTGATTGTTGTAGTCCGCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGACGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGGGGTGGTGAACGCGTGGaa                                                                                                                              <  2:595015/247‑1 (MQ=255)

                                                                                                                    cAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGCAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCACAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCCGACCTGGAAGCTGTTTCAGACCGATCAAATCTTCAAAGGGGGTGGTGAACGGGTGGa                                                                                                                               >  2:385982/1‑250 (MQ=255)

                                                                                                                                  atGGGAGCCGTGCACCAGGTCATTTCCATTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGGGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACTCGTGGAAGCCGCCTTGTTcta                                                                                                                <  2:472589/251‑3 (MQ=255)

                                                                                                                                  atGGGAGCCGAGCACCAGGTCATTCCCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCa                                                                                                                <  1:407891/251‑1 (MQ=255)

                                                                                                                                   tGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGCGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTcc                                                                                                                 >  1:430798/1‑249 (MQ=255)

                                                                                                                                       aGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGAAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCAAGTCGCCTTCGCCCGAAAAGCCGTAACCCTGCTGCATAAGACGATGTACGGCCAGACCCGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAAAGCGGGGGAGCCAACGTGATCCAGGaa                                                                                                            >  2:703675/1‑250 (MQ=255)

                                                                                                                                          cGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGc                                                                                                        >  1:139648/1‑251 (MQ=255)

                                                                                                                                          cGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGc                                                                                                        >  1:263147/1‑251 (MQ=255)

                                                                                                                                            ggCACCAGGTCATTCCCTTTCTCGAAGTGAGAGGTGTAGTCCGCCATAAAGGAGGGCCCGACCTGCAGGACGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTAGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTGCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGTAAGCCACCGTGTTCCATGATACGCtt                                                                                                      <  2:205240/250‑1 (MQ=255)

                                                                                                                                                ccatgTCATTACCTTGCACGAAGTGATAGGGGTAGTCCTCCAGAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCTGCGGTTGTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTGCAAAGGTGTTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATc                                                                                                  <  2:125102/247‑1 (MQ=255)

                                                                                                                                                 cAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGGAAAAGCGTAGAAGCCACCTTGTTCAAGGAAACGCTTCATcc                                                                                                 >  2:390110/1‑251 (MQ=255)

                                                                                                                                                  aGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCAGTAACCCTGCTGCATCAGACGCTGTAAGGCCCGACCTGGAAGCTGTTTCAGCCCGTGCAAATCTTCAAAGGAGGGGGGGAACGCGTGGAAGCCACCTTGTTCCAGCAAACGCttattacc                                                                                                >  2:165345/1‑245 (MQ=255)

                                                                                                                                                    gccATTGCCTTTCTCTAAGTGATAGGTGTAGTCCTCCATAAAGGAGTTGCCGCCCTGCAGACCGGTTGACATCACCCTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCATACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTCGTGAACGTGTGGAAGACGCCTTGTTCCAGTAAACGCGTCATCCCCa                                                                                              <  2:126826/249‑1 (MQ=255)

                                                                                                                                                       atttccTTTCTCGAAGGGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGGTGACATCTCCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCTTCCCCAGCt                                                                                           <  2:373143/247‑1 (MQ=255)

                                                                                                                                                        ttACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTc                                                                                          >  1:286682/1‑251 (MQ=255)

                                                                                                                                                        gttccgtgcgcGAAGTGTTAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCCCCTGCATGATGCGAAGCAGGGCGGCTGGTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGGCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTc                                                                                          <  2:331106/241‑1 (MQ=255)

                                                                                                                                                              ttCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATa                                                                                      <  1:411817/249‑1 (MQ=255)

                                                                                                                                                              ttCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACACTTCATCCCCAGCTCAATa                                                                                      <  1:432504/249‑1 (MQ=255)

                                                                                                                                                               tctcGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATAcg                                                                                    <  1:190454/250‑1 (MQ=255)

                                                                                                                                                               tctcGAAGTGAGAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATAcg                                                                                    <  1:79901/250‑1 (MQ=255)

                                                                                                                                                                      gTGATAGGTGTAGTCCGCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCGGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCtt                                                                            <  2:353336/251‑1 (MQ=255)

                                                                                                                                                                        gATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTc                                                                           <  1:361467/250‑1 (MQ=255)

                                                                                                                                                                           aGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGAGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCCATACGCCCCGCTTcccac                                                                       >  2:488486/1‑248 (MQ=255)

                                                                                                                                                                             gtgtAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCTAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCAc                                                                     >  1:399304/1‑251 (MQ=255)

                                                                                                                                                                             gtgtAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCAc                                                                     <  1:21207/251‑1 (MQ=255)

                                                                                                                                                                             gtgtAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGGGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGATTCATCCCCAGCTCAATACGCGCCGCTTCCAGCa                                                                      >  1:421592/1‑250 (MQ=255)

                                                                                                                                                                                tAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGGGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTCCATCCCCAGCTAAATAAGCGCCGCCTACAGCACggt                                                                  >  2:520251/1‑249 (MQ=255)

                                                                                                                                                                                 aGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTc                                                                 >  1:142712/1‑251 (MQ=255)

                                                                                                                                                                                  gTCCTCCATAAAGGAGGTGCCGCCCTGACGACCGGTTGACCTCACCTTCATGATGCGAAGCAGGTCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTGCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCt                                                                <  2:224454/251‑1 (MQ=255)

                                                                                                                                                                                    cctccATAAAGGAGTTGCCGCCCTGCAGACCGGTTGACATCTCCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAATGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACTTTCTg                                                               <  2:588929/250‑1 (MQ=255)

                                                                                                                                                                                      tccATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAACGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGACGCTGTTTCAGAcc                                                                                                                                                                  >  2:493724/1‑149 (MQ=255)

                                                                                                                                                                                       ccATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATAACCTTCATGATGCGAAGCATCGCGGTGGAATTCAAGTCGTCTCCGCCCGCAAAGCCGTAACCCTGCTGCAGCAGACGCTGTACGGCCAGTCCTGGAAGCTGTTTCATACCGTGCAAATCTTCAACGGTGGGGGTGCACTCGTGGAAGCCGCCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTc                                                             <  1:662380/249‑1 (MQ=255)

                                                                                                                                                                                         aTAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGttt                                                         >  1:607654/1‑251 (MQ=255)

                                                                                                                                                                                           aaaGGAGGTGCCGCCCGGACGTCCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTGCCAGTCGCCTTCGCCCGCAATGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAATGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGtttt                                                        <  2:203982/250‑1 (MQ=255)

                                                                                                                                                                                             aggaggTGCCGCCCTGCAGAACGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGGACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCACATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTACTCCCCTGCTCAATACGCGCCGCTTCCAACCAGTTCTGTCGtttttcc                                                     >  2:279041/1‑249 (MQ=255)

                                                                                                                                                                                               gaggTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGcc                                                   >  2:164334/1‑251 (MQ=255)

                                                                                                                                                                                                       gccCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGAttt                                           <  1:360152/251‑1 (MQ=255)

                                                                                                                                                                                                            gCAGACCGGTTGAAAACACCTTCCAGATCGGAAGCAGGGCGGCGTTTTTCCAGTCGCCTGCGCCCGCAAAGCCGTAACCCTGCGGCCTCAGACGCTGTAGGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGTAAGCTACCTTGTTCCAGGAATCGCTTCATCCCCAGCTCAATACGCGCCGCTGCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTtgtgtt                                      <  2:263108/251‑2 (MQ=255)

                                                                                                                                                                                                                       gACATCACCTTAATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGTTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCAt                           <  2:238008/251‑1 (MQ=255)

                                                                                                                                                                                                                        ccATCTCCTTCATGAAGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTGCGCCGTGGATTTGTGTTGCAGGCGTCATg                          <  2:399304/250‑1 (MQ=255)

                                                                                                                                                                                                                            cACCTTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCGTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATTCGCGCCGCTTCCAGCACGTTCTGTCTTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGgtg                       <  1:587519/250‑1 (MQ=255)

                                                                                                                                                                                                                            cACCGTCATGATGCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCCGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGACCGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGgtgt                      <  1:699692/251‑1 (MQ=255)

                                                                                                                                                                                                                                    tgatgCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCtt              >  1:460072/1‑251 (MQ=255)

                                                                                                                                                                                                                                      atgagAAGCATGGCTGCTGTTTTCCCTTCGCCTGCTGCCGCATAGCAGTAACCCTGCTGCATCACACGCTGTACTGCCTGACCTGGAAGCTGTTTCATACCGTGTAAATCTTCAAAGGGTGTGGTGACCGCGTGGAAGCGTCCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCTTTGTTCGCCGTTGATTTGTGTGTCAGGcttcttggggtatcacctgta            <  2:479236/247‑21 (MQ=255)

                                                                                                                                                                                                                                      atgCGAAGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACATGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGAGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGAGTAGCAGCTTTc            >  1:79539/1‑251 (MQ=255)

                                                                                                                                                                                                                                           aaGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACt       <  1:146945/251‑1 (MQ=255)

                                                                                                                                                                                                                                            aGCAGGGCGGCGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTc      <  1:236986/251‑1 (MQ=255)

                                                                                                                                                                                                                                            aGCAGGGCGGCGGGTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGTCCTGGCAGCTGTTTAAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGGGGAAGCCACCTGGTTCCAGGAAACGCTTCATACCCAGCTCAATACGCGCCGCTTCCAGCCAGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTTCAGTCGTCATGGTGTAGCCTCTTTCGTACTc      >  2:161559/1‑251 (MQ=255)

                                                                                                                                                                                                                                                 ggcggcGGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCTGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTTGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCg  >  1:245108/1‑250 (MQ=255)

CGACAATCGCTGGGCCGGTTTGGGTATTGAAGATCAGGCGGGCAGGATCGTCCTTACCACCAATACCGAGATGCTGAACGTCGAGGATCGGTTTCTCTTCTGCGGCGATCGACGGGCAGACTTCCAGCATATGGGAGCCGAGCACCAGGTCATTACCTTTCTCGAAGTGATAGGTGTAGTCCTCCATAAAGGAGGTGCCGCCCTGCAGACCGGTTGACATCACCTTCATGATGCGAAGCAGGGCGGCAGTTTTCCAGTCGCCTTCGCCCGCAAAGCCGTAACCCTGCTGCATCAGACGCTGTACGGCCAGACCAGGAAGCTGTTTCAGACCGTGCAAATCTTCAAAGGTGGTGGTGAACGCGTGGAAGCCACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTGGCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCG  >  NC_000913/67158‑67648

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