BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,003	G→A	D225N (GAT→AAT)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,003	0	G	A	96.8%	367.6 / NA	124	D225N (GAT→AAT)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand): ref base G (4/0); new base A (49/70); total (54/70)
Fisher's exact test for biased strand distribution p-value = 3.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

ATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955757‑3956249

aTCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGAGGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCTGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGc                                                                                                                                                                                                                                                     <  1:121759/251‑1 (MQ=255)

     aCTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGGGGAAGACTCCGCCTGACTG‑AAACCAGCGCTGAATTCGGGTc                                                                                                                                                                                                                                                >  1:317837/1‑251 (MQ=255)

     aCTGTCACAGCAGCAAGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGACCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGGAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTc                                                                                                                                                                                                                                                <  2:418851/251‑1 (MQ=255)

            cagcagcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGAGCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGG‑AGAAGCGGAAGACTCCCCCCGCCTG‑CAAGCAGCGCTGAATGCGGGGCATCCggg                                                                                                                                                                                                                                         >  2:155916/1‑250 (MQ=255)

             agcagcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGtt                                                                                                                                                                                                                                        <  1:358022/251‑1 (MQ=255)

               cagcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGa                                                                                                                                                                                                                                      <  1:22586/251‑1 (MQ=255)

               cagcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTACTGATCAAACCACTGATGCCGCAAATCCAAGTGATCGCCCT‑AGAAGCGGAAAACTCCGCCTTCCGG‑AAAGCAACGCTGAATGCGGGTCCTCCGGTTgt                                                                                                                                                                                                                                      >  1:96964/1‑250 (MQ=255)

                   agGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACCACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCCCCTGCCCG‑AAAGCAGCGCTGAATGCGGGGCATCCGGTTGATCTg                                                                                                                                                                                                                                  >  2:496068/1‑251 (MQ=255)

                   agGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTAGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATAAAACAACTGATGCCGCAACTCAAAGGGATTGCCGT‑AGAAGCGGAAAACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGCTGATCtt                                                                                                                                                                                                                                  >  1:534815/1‑250 (MQ=255)

                     gggTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGcc                                                                                                                                                                                                                                <  1:92851/251‑1 (MQ=255)

                         tCACCTGAGAGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCCCGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCCAACCACTGCAGCCGCAAAAACAAGTGATCGCCGC‑AGAAGCGGAAGACTCAGCCCGCCCG‑AAAACACGGCTGACTGCGGGGCATCCGGCTGGTCTGCcgcgc                                                                                                                                                                                                                            >  2:533406/1‑251 (MQ=255)

                          cACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCcgcgc                                                                                                                                                                                                                            <  1:155916/250‑1 (MQ=255)

                                   gccgccGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTTCTTATCAAACAACTGATGCCGCAAATCCAAGTGAGCGCCGTCAAAACGGGCAAAC‑CCGCCCTCCTG‑AAAGCAGCGCTGAAAGCGGGTCATCCGGTTGATCTGcccgccggtgggcgat                                                                                                                                                                                                                  >  2:529482/1‑237 (MQ=255)

                                      gccgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTg                                                                                                                                                                                                              <  1:122436/251‑1 (MQ=255)

                                       ccgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGc                                                                                                                                                                                                             <  1:496068/251‑1 (MQ=255)

                                         gTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGGATGCGGGGCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTg                                                                                                                                                                                                           >  1:175226/1‑251 (MQ=255)

                                            cGACCATCCGATGGGGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACGGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAg                                                                                                                                                                                                        <  1:589287/251‑1 (MQ=255)

                                             gACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCCAAGTGATCGCCGT‑AGAAGCGGCAGACACCGCCTGCCTG‑AAAGCAGCGCTGCATGCGGGTCATCCGGGTGATCTGCCGCGCGGA‑GGGCTCTTTGCTGAAgg                                                                                                                                                                                                       >  2:202199/1‑251 (MQ=255)

                                                   ccGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGTACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAAGCAAAGTGATCGCCGT‑AGAAGCGTAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGc                                                                                                                                                                                                 <  2:314300/251‑1 (MQ=255)

                                                   ccGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGc                                                                                                                                                                                                 <  2:217067/251‑1 (MQ=255)

                                                      aTGGTGATTGGCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTGGGGCCGGTCGGCGGCGGCGGTCGGGCTGCGTGCGTGGCGGTGCTGAACAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGTCTCCGTCTGCCTT‑AAAGGAGCGCTGAAGGCGGGTCAGCCGGTTGAGCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCTTAGCGGt                                                                                                                                                                                              <  2:224322/251‑1 (MQ=255)

                                                           gATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAAc                                                                                                                                                                                         <  1:477399/251‑1 (MQ=255)

                                                           gATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGCATCCGCCTGCCTG‑GAAGCAGCGCTGAATTCGGGTCATCCGGTTGATCTGCCGCGCGGA‑GGGGTGTTTGCTTAAGGCGTAGCGGTaaaa                                                                                                                                                                                          >  2:440279/1‑250 (MQ=255)

                                                               gCCGTGCAAGGCACGCCGGCGCTGGAACTGCTCCATCCGGACGCCCATCTCGACCGCGTATGTGTGCCAGTCGGCGGCGGCGGTCGGGCTGCGGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCAGCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCa                                                                                                                                                                                      <  2:585770/250‑1 (MQ=255)

                                                                  gggCAAGGCACGCTGGCGCTGGACCTGCTCCAGCAGGACGCCCATCGCGCCCGCGTATTTGGGCCAGTCGGCGGCGGCGGTCTGGCTGCGGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCgg                                                                                                                                                                                  <  2:135444/251‑1 (MQ=255)

                                                                         gCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGaa                                                                                                                                                                           <  1:465882/251‑1 (MQ=255)

                                                                              ctggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTTAAACAACTGATTACGCAAATTAAAAGGATTGCCGT‑TGAAGCGGAAGGCTCCGCCTGCCAG‑ACAGCAGCGCTGGATGCGGGTCATCCGGGTGTTCTGTCGCGCGTA‑GGGCTATTTGCTGAAGGCGGAGGGGTTAAACCCCTATGTGACGAAACCtt                                                                                                                                                                      >  2:446505/1‑251 (MQ=255)

                                                                               tggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTTATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTTCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGGA‑GGGCTATTTGTTGTAGGCGTAGCGGTAAAACGCATCGGGGACGGAACCTTc                                                                                                                                                                     >  1:20337/1‑251 (MQ=255)

                                                                               tggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCCAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTc                                                                                                                                                                     >  1:335458/1‑251 (MQ=255)

                                                                                ggcgctggAACTGCTCCAGCAGGACTCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTTCTGATCAAACAACTGAGGCCGCAAATCAAAGTGATCGCCGC‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCGGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATGTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTcc                                                                                                                                                                    <  1:362175/251‑1 (MQ=255)

                                                                                ggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGGA‑GGGCTATTTGCTGGAGGCGTAGCGGTAAAACGCATCGGGGACGAAACCTTcc                                                                                                                                                                    >  1:145945/1‑251 (MQ=255)

                                                                                  cgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGGAGAATCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGt                                                                                                                                                                  >  1:227215/1‑251 (MQ=255)

                                                                                    cgggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGGGACGAAACCTTCCGttt                                                                                                                                                                <  1:479445/249‑1 (MQ=255)

                                                                                        aaCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATAAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAAACTCCGCCTGCCTG‑AAAGCAGCGCTGAAAGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTTCTGAAGGCGGAGCGGTGAAAAGCATCGGTGAAGAAACCTTCCGATTAttc                                                                                                                                                            >  1:410206/1‑249 (MQ=255)

                                                                                          cTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGcc                                                                                                                                                           >  1:262860/1‑250 (MQ=255)

                                                                                              gccAGCGGGGCGCCCCTCTGGACCGCGTATTTGGGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGTTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTTATCTGCCGCGCGTA‑GGGTTATTTGCTGAAGGCGTAGCGGGAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAg                                                                                                                                                      <  2:69114/250‑1 (MQ=255)

                                                                                               ccAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAg                                                                                                                                                      <  1:407475/250‑1 (MQ=255)

                                                                                                cagcagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGGTGATCTGCCGCGCGTA‑GGGCTATTTGCTGTAGGCGTAGCGGGAAAACGCATCGGGGACGAAACCTTCCGGTTagggcagggggt                                                                                                                                                    >  1:137202/1‑240 (MQ=255)

                                                                                                 agcagGACGCCCATCTCGACCGCGTATTTGTGCCAGGCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAt                                                                                                                                                   <  1:353907/251‑1 (MQ=255)

                                                                                                   cagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGCAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑g                                                                                                                                                                                                                       >  1:275289/1‑181 (MQ=255)

                                                                                                   cagGACGCCCATCTCGACCGCGTATTTGTGCCAGGCTGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACCGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTGGATCTGCCGCTCGTA‑g                                                                                                                                                                                                                       <  2:275289/181‑1 (MQ=255)

                                                                                                    agGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGGTGGTCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGGTGCGGTAAAACGCCTCGGTGAAGAAACCTTCAGTTTATGCCAGGAGTAtctc                                                                                                                                                >  2:67298/1‑251 (MQ=255)

                                                                                                        cGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGGCTGATCAAAAAACCGGTGCCGCAAATCAAAGTGATCCGCGG‑AGAAGCGGAAGAAACCGCCTGCCTGGAAAG‑AGCGACGAATGCGGGGCAACCGGGTGATCTGCCGCGCGTA‑GGGCTAATTGCTGGAGGCGTTGCGGTAAAACGCAGGGGTGGCGGAAACTTCCGGTTATGCCAGGAGGATCTcgagg                                                                                                                                            >  2:491488/1‑249 (MQ=255)

                                                                                                              tctcGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGTCTGCTGGCGTGGCGGTGCTGATAAAACCAATGGTGCCGCAAATAGAAGTGATCGGCGT‑AGAAGCGGGAGACTCCGCCTGCCTG‑AAATCCGCGCTGGATGGGGGGACTCCGGTTGAGCTGCCACGCGTA‑GGGCTGTTTGCTTAAGGCTTAGCGGTAAAACGCATCGGTGCCGAAACCTTCCGTTTATGCCCGGAggtacgcggcggaatc                                                                                                                                       >  2:479952/1‑235 (MQ=255)

                                                                                                               ctcGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTTTGCCAGGAGTATCTCGACGCcatcat                                                                                                                                     >  1:284684/1‑251 (MQ=255)

                                                                                                                 cGACCGCGTTTTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatca                                                                                                                                      <  2:226842/248‑1 (MQ=255)

                                                                                                                 cGACCGCGTATTTGTCCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGACGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatcat                                                                                                                                     <  2:28285/249‑1 (MQ=255)

                                                                                                                  gtctgcgTATGTGTGCCAGTCGGAGGCGTGGGTCCGGCTTCTGGCGTGGCGGTGCGGATCAAACAACTGATGCCGCAAATCAAAGTGATCGGCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAAGGCGGGTCATCCGGTTGATCTGCCGCGCTTA‑GGGCTTTTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAGTCGTCCGTTTAGCCCAGGAGGATCTCGACGTCATCATCAc                                                                                                                                  <  2:354089/247‑1 (MQ=255)

                                                                                                                  gACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAAATGATGCCGCCAATCCAAGTGGTCGCCGG‑AGAAGGGGAAGACCCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGGCATCCGGGTGGTCTGCCCCGCGTC‑GGGGTATTTGCTGGAGGCGTAGCGGGAAAACGCATTGGTGGCGAAACCTTCCGTTTATGCCCGGGGTATCTCGACGGcatcataac                                                                                                                                  >  2:555644/1‑248 (MQ=255)

                                                                                                                   aCCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCGCGACGACATCATCAcc                                                                                                                                 <  1:514811/251‑1 (MQ=255)

                                                                                                                     cgcgTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAAGCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGGTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGt                                                                                                                               <  1:24494/251‑1 (MQ=255)

                                                                                                                          atttggtccAGGCGGCGCCGGCGGTCCGGCTGCTGGCGGGGCGGTGCTGATCAAACAACTGATGCCGCAAAGCAACGGGATGGCCGT‑ATAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGGA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACTCATCGGTGACTAATCCTGCCTTTTATGCCAGGAGTATCTCTACGCCATCATCACCGTCggta                                                                                                                          <  2:201284/244‑4 (MQ=255)

                                                                                                                            ttggtccAGTCGGCGGCGGTGGTCTGGCGGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAGCCGTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAg                                                                                                                         <  2:387351/245‑1 (MQ=255)

                                                                                                                             tgggtcAGTTGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAACCAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGTCGCGCGTA‑GGGCTATTTGCTTAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgc                                                                                                                        <  1:479351/245‑1 (MQ=255)

                                                                                                                              gtgCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcg                                                                                                                       <  1:310316/250‑1 (MQ=255)

                                                                                                                               ggCCAGTCGGTGGCGGCGGTCGGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGAGCGCTGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTAGTTGCTGAAGGCGTATCGGTAAAACGCATCGGTGACGAAACCTGCCGGTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgat                                                                                                                     <  1:289986/250‑1 (MQ=255)

                                                                                                                                         gggcggcggTCTGGCTGCTTGCGTGGCGGTGCTGATCAAACAACGGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCtgtg                                                                                                           <  2:242847/250‑1 (MQ=255)

                                                                                                                                             ggcggTCTGGCTGCTGGCGTGGCGGTTCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGAGCTGAATGCTGGTCATCCGGTTGATCTGCCGGGCGTA‑GGGCTATTTGCTGGAGGCGGAGCGGGAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCCGGAGTATATCGTCGACATCATCACCGTCGGTAGCGATGCGATCTGTgcggc                                                                                                       >  1:585708/1‑251 (MQ=255)

                                                                                                                                             ggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGCTCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑TAAGCCGCGCTGAATGCGGGTCATCCGGGTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGGGACGCACCCTTCCGGTTATGCCAGGGGGATTACGACGACCTCATCACCCTCGATCGCGGTGCGATCTTTgcggg                                                                                                       >  2:194848/1‑250 (MQ=255)

                                                                                                                                                ggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGGTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGAt                                                                                                    >  1:195527/1‑251 (MQ=255)

                                                                                                                                                 gTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATg                                                                                                   <  2:560593/251‑1 (MQ=255)

                                                                                                                                                     ggcggctgGCGTGGTGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGGTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTGCCGTGTATGCCAGTAGTAGCTCGACGACATCATCACCGGCGATAGCGATGCGATCTGTGCGGGGATGaa                                                                                                 <  2:391307/245‑1 (MQ=255)

                                                                                                                                                     ggCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGaa                                                                                                 <  1:186513/249‑1 (MQ=255)

                                                                                                                                                     ggCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGAGGAAgg                                                                                               <  1:468241/251‑1 (MQ=255)

                                                                                                                                                            ggcgtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTTTGCGGCGATGAAGGATTTAt                                                                                         <  1:198689/250‑1 (MQ=255)

                                                                                                                                                             gcgtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTc                                                                                       <  2:297175/251‑1 (MQ=255)

                                                                                                                                                             gcgtggcgGTGCTGATAAAACAACGGATGCCGCAAATCAAAGTGAGCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGAGCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCCCCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTc                                                                                       <  2:413617/251‑1 (MQ=255)

                                                                                                                                                                    ggTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCATG‑AAAGCAACGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑CGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGAGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAtg                                                                                >  1:31237/1‑251 (MQ=255)

                                                                                                                                                                     gTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGAGCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAtgt                                                                               <  2:335458/251‑1 (MQ=255)

                                                                                                                                                                           aTCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTgcgcgc                                                                         >  2:64091/1‑251 (MQ=255)

                                                                                                                                                                              aaaCAACTGCTGCCGCAAATCAAAGGGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTACGTTTATGCCAGGAGTAGCACGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAGGGTCGCGCGGt                                                                      <  2:299536/251‑1 (MQ=255)

                                                                                                                                                                              aaaCAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGt                                                                      <  1:67298/251‑1 (MQ=255)

                                                                                                                                                                                acaacTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTgg                                                                    >  2:281297/1‑251 (MQ=255)

                                                                                                                                                                                   acTGACGCCGCCAATCAAACCGATCGCAGT‑AGAAGCGGAAGACGCCGCCTGACTG‑AAAGCCGCGCTGAATGCCGGTCATACGTGTGATCTGCAGCGCGTA‑GGGATTTTTGCTGAAGGCGGAGCGGCAAAACGCATCGGTGACGACACCTTCAGTTTATGCCAGGAGTATCTCGACGACATCATCCCCGTCGATAGAGATGCGATCTGTGCGGCGATGAGGGATTTCTTCGAAGATGTGCGCGCGGTGGCgg                                                                 >  2:585021/1‑251 (MQ=255)

                                                                                                                                                                                       aTGCCGCAAATCAAAGTGAGCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAAGGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGGCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAAcc                                                             <  2:449107/251‑1 (MQ=255)

                                                                                                                                                                                             cAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTgg                                                       <  1:397394/251‑1 (MQ=255)

                                                                                                                                                                                                    aagtgcacGCCGT‑AGAAACGGAAGACTCCGCCTGCCTG‑TAAGCAGCCCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGGTATTTGCTGAAGGCGTAGCGGTAAAACCCATCGGCGACGAAACCTTCCGTTTATGACAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATCCGATCTGTGCGGCGCTGAAGGATTTAATCGAAGATGTGCGCGCGGGGGCGGCACCCTCTGGCgcgctg                                                >  2:270675/8‑251 (MQ=255)

                                                                                                                                                                                                     aGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctgg                                               >  1:416892/1‑251 (MQ=255)

                                                                                                                                                                                                      gTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggc                                              >  1:567484/1‑251 (MQ=255)

                                                                                                                                                                                                      gTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggc                                              >  2:488454/1‑251 (MQ=255)

                                                                                                                                                                                                      gTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggc                                              <  1:315678/251‑1 (MQ=255)

                                                                                                                                                                                                       tGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATATCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGAgcgctggcg                                             >  2:203197/1‑251 (MQ=255)

                                                                                                                                                                                                         aTCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgct                                           <  1:391617/251‑1 (MQ=255)

                                                                                                                                                                                                            gCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCgctggcgctggc                                        >  1:438686/1‑251 (MQ=255)

                                                                                                                                                                                                                t‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGTACCCTCTGGCGCGCTGGCGCTGGCgggg                                    >  2:365881/1‑250 (MQ=255)

                                                                                                                                                                                                                t‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCggg                                     >  2:409075/1‑250 (MQ=255)

                                                                                                                                                                                                                    aaGCGGAAGACTCCGCCTGCATG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGAGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCCGGAGTATCTCGACGACATCATCACCGTCGATCGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGac                                   >  2:403055/1‑248 (MQ=255)

                                                                                                                                                                                                                    aaGCGGAAGACCCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATg                                 <  2:54538/251‑1 (MQ=255)

                                                                                                                                                                                                                     aGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGCAGGCGTAGCGGTAAAACGCAGCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATATCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTTGCGCTGGAGGGAATg                                 >  2:139479/1‑250 (MQ=255)

                                                                                                                                                                                                                        ggAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaa                              >  1:94750/1‑250 (MQ=255)

                                                                                                                                                                                                                           aGACTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTTGCGCGCTGGCGCTGGCGGGAATGAAAAAat                          <  1:524583/251‑1 (MQ=255)

                                                                                                                                                                                                                             aCTCCGCCTGCCTG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAatat                        >  1:57646/1‑251 (MQ=255)

                                                                                                                                                                                                                                  gcctgcctG‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGc                   >  1:454262/1‑251 (MQ=255)

                                                                                                                                                                                                                                      ggctg‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGAGCGGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGGCGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGATGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGTAATGAAAAAATATATCGCCCTg               <  2:446660/249‑1 (MQ=255)

                                                                                                                                                                                                                                         cg‑AAAGGAGCGCTGACTGCGGATCATCCGGATGATCTTCGGCGCGTAAGGGC‑ATATGCTGAAGGAGTAGCGGTAAAACGCATCGCTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGCCATCATCACCGTCGCTAGCGATGCGATCTGGGCGGCGCTGCAGGATTTATTCGAAGATGTGCGCGCGTAGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAGTGAGAAActatgacgacctgca             >  2:65877/2‑239 (MQ=255)

                                                                                                                                                                                                                                          g‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGcaca           <  1:533406/251‑1 (MQ=255)

                                                                                                                                                                                                                                          g‑AAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGcac            <  2:416892/250‑1 (MQ=255)

                                                                                                                                                                                                                                                 agcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCccaaccatc     >  2:529859/1‑247 (MQ=255)

                                                                                                                                                                                                                                                 agcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTc     >  1:88883/1‑251 (MQ=255)

                                                                                                                                                                                                                                                    gcTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTcgcg  <  2:175226/251‑1 (MQ=255)

ATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGT‑AGAAGCGGAAGACTCCGCCTGCCTG‑AAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTA‑GGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955757‑3956249

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