BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,984	2 bp→CG	coding (537‑538/1611 nt)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,984	0	T	C	98.1%	313.8 / ‑2.4	105	S180G (AGC→GGC)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (42/61); total (42/63)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01
*	NC_000913	73,985	0	A	G	98.1%	285.9 / ‑2.3	105	R179R (CGT→CGC)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (42/61); total (42/63)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 105 total aligned reads displayed.

AGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGG  >  NC_000913/73746‑74226

||

     caacccGAGGCAATACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTATTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGGTGCCGGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTggg                                                                                                                                                                                                                                    <  2:127998/248‑1 (MQ=255)

     caacccGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTATCTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGGGGCGAAGTTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAAGGCCAGCCGCGCATCCCAAGCTgtg                                                                                                                                                                                                                                    <  1:790568/248‑3 (MQ=255)

         ccGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGggcggc                                                                                                                                                                                                                                <  1:77157/251‑1 (MQ=255)

           gAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCaa                                                                                                                                                                                                                              <  1:851094/251‑1 (MQ=255)

               catcacaCCATCTGGAGCAGCGCCAGCTTCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAgctg                                                                                                                                                                                                                          <  2:525244/248‑4 (MQ=255)

                  cccacCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACg                                                                                                                                                                                                                       <  1:861800/249‑1 (MQ=255)

                    cacCATCTGGATCAGCGCCAGCACCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACgtt                                                                                                                                                                                                                     <  2:862751/251‑3 (MQ=255)

                     cccATGTGGACCAGCGCCAGCAGCTCCGCGCGGTCAGGGACGTAGTTGTAACTCAGCGCCTGATAGAGTGCCAGCGCTATAGTGGTTGCCTGCGGACCGCCGCCCAGCGAGAGCACGGTGGCGAAGCTGTCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACgttg                                                                                                                                                                                                                    <  2:482847/250‑4 (MQ=255)

                            gggtcagcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGcc                                                                                                                                                                                                              <  2:464379/246‑1 (MQ=255)

                            ggatcagcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCg                                                                                                                                                                                                             <  1:466743/247‑1 (MQ=255)

                            ggatcagcGCCAGCACCGCCGCGCGGGCAGGATCGCAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGTCCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTc                                                                                                                                                                                                                 <  1:344529/243‑1 (MQ=255)

                               tcagcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCgggg                                                                                                                                                                                                          >  1:767680/2‑251 (MQ=255)

                                 agcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCgggaat                                                                                                                                                                                                        >  1:811311/1‑248 (MQ=255)

                                      cAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGtttt                                                                                                                                                                                                   <  1:358905/251‑1 (MQ=255)

                                         catcgCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGCTAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGACGCCGAAACCACGGCCCTTCGACGAAGGGGCAAAACAGCCAGGCGCGCATCCCAAGCTGGGGGGCCCGCTGACGTTGTTCGCCGGGAATGttttac                                                                                                                                                                                                 >  2:574470/5‑248 (MQ=255)

                                            cgCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTGGCCTGCGGACCGCCGCCCAGCGATAGCTCGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTg                                                                                                                                                                                             <  1:420721/250‑1 (MQ=255)

                                                   cgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAg                                                                                                                                                                                      <  1:193537/251‑1 (MQ=255)

                                                   cgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTAACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTgg                                                                                                                                                                                        >  1:439899/1‑249 (MQ=255)

                                                      cttccgatctAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCCCGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGta                                                                                                                                                                                    >  1:773032/10‑250 (MQ=255)

                                                        aGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGtaata                                                                                                                                                                                 <  1:402159/251‑1 (MQ=255)

                                                        aGGATCGTAGTCGGAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGtaata                                                                                                                                                                                 <  2:469714/251‑1 (MQ=255)

                                                         ggATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCCTCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGtaata                                                                                                                                                                                 >  1:919819/1‑250 (MQ=255)

                                                           aTCGTAGTCGTAACTCATCGCCTGATAGATTGCCAGCACGATAGTGGTTGCCTGCGGACCGCAGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAACCAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGc                                                                                                                                                                              <  1:958834/251‑1 (MQ=255)

                                                           aTCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGCCGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGc                                                                                                                                                                              >  1:901332/1‑251 (MQ=255)

                                                             cggaggcGTAACTCAGCGCCTGATAGATTGCCAGCGCGAGAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCgg                                                                                                                                                                            <  2:945537/245‑1 (MQ=255)

                                                              gTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGGTTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGtaatatacggc                                                                                                                                                                           >  1:764700/1‑245 (MQ=255)

                                                                    gTAACTCGTCGCCTGTTTGAGTGCCAGGTAGATAGGGGTTGCCGGCGGACCGCCGCCCAGCGATAGCACGGGGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCa                                                                                                                                                                     <  2:983258/251‑1 (MQ=255)

                                                                        cccAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGAACGGGGGCGAAGCTGGCTAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCgg                                                                                                                                                                 <  2:811289/249‑1 (MQ=255)

                                                                          cAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCggt                                                                                                                                                                <  1:395699/250‑2 (MQ=255)

                                                                           aGCGCCTGATAGAGTGCCAGCGCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTTGAGTAATAAGCGGCTCGCCATCggtag                                                                                                                                                              <  2:905276/251‑4 (MQ=255)

                                                                            gcgcATGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGa                                                                                                                                                             <  1:11328/251‑1 (MQ=255)

                                                                              gcctggaagtatGCCCGCTCGATAGTGGTCTCCTGCGGACCGCCGCCCAGCGACAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCACGCGGCGGGATTGGTCGCCGTAACCACGGCCATTCGACGAAGCGGATAAAATGCCGGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCTCCATCGGTAGAtt                                                                                                                                                           <  2:144775/240‑1 (MQ=255)

                                                                                  gATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGGGCAACCCAAGCTGGGCGGCAAGCTAACGTTGGTCGCCGGGGATGTTTTCCAGTGCCTGGGGCAATAAGCGGGTCCCCATCggtgggttaaa                                                                                                                                                        >  2:160580/1‑241 (MQ=255)

                                                                                  gATAGATGGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaa                                                                                                                                                       <  1:14051/251‑1 (MQ=255)

                                                                                   aTAGATTGCCAGCCCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaa                                                                                                                                                      <  2:247155/251‑1 (MQ=255)

                                                                                    tAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGACGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGGGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCAGGAGTGATAAGCGGCTCGcccacggttggttaaaaa                                                                                                                                                      >  1:69374/1‑234 (MQ=255)

                                                                                    tAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGCAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGGAGATTaaaaa                                                                                                                                                      >  1:125653/1‑250 (MQ=255)

                                                                                    aaGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaa                                                                                                                                                     <  2:265988/250‑1 (MQ=255)

                                                                                       aTTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatac                                                                                                                                                  <  1:121873/251‑4 (MQ=255)

                                                                                          gCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACATAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatacgtg                                                                                                                                               >  1:647700/1‑245 (MQ=255)

                                                                                          gCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGATCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCTGGGAATGTTTTCCAGTGTCAGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatacgt                                                                                                                                                <  1:310966/250‑6 (MQ=255)

                                                                                              gCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatacgtgcgcc                                                                                                                                           <  1:238870/251‑11 (MQ=255)

                                                                                                          gtgCCTGCGGACCCCCGCCCAGCGTTAGCACGGTGGCGAAGCTGGCGAAACAGAGCACAAAGATTAGTGACGCAACCTGCGGGAGTTGTCGCCGTACCCACGGCCATGCGACGAAGCGTACAAAATGCCAGCCGTGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTCTCCAGTGCCTGGAGTAATAAGCGGCTCGCCAGCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCt                                                                                                                               <  2:266175/249‑1 (MQ=255)

                                                                                                            gCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAAAGCCAGCCGCGCATCCCAAGCTGGGCGACAAGCCGGCGTTGTTCGCGGGGGATGTTTTTCCGGGCCTGGAGTAAAAAGGTGCTCGCCATCAGGAGATTaaaaaatcagagcgccagaaaaaaaccttc                                                                                                                             >  2:938461/1‑227 (MQ=255)

                                                                                                                gCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCGAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCa                                                                                                                           <  1:660318/249‑1 (MQ=255)

                                                                                                                gCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAgg                                                                                                                         <  1:586278/251‑1 (MQ=255)

                                                                                                                 cGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCCACGAAACGGAAAAAATGCCAGCCGCCCATCCCAAGCTGGGCGGCAAGCTGACGTTGGTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTCAAAACTACGTGCGCCAGCAAAATACCTTGCAGGc                                                                                                                        >  1:631750/1‑251 (MQ=255)

                                                                                                                    aCCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGt                                                                                                                     <  1:187495/251‑1 (MQ=255)

                                                                                                                      cgccgcCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGGAACCACGGCCATTCGGCGAAGCGGAAAAAATGCCCGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGGTGTTCGCCGGGGATGTTTTCCAGTGCCCGGAGTAATAACCGGCTCGCCCTCGGTCGATTAAAAAATACGTGAGCCAGCAAAATACCTTGCAAGCCGt                                                                                                                     >  2:116783/1‑249 (MQ=255)

                                                                                                                       gccgccaaGCGATAGCACGGGGGCGAAGCTGGCGAAACAGAGCTTAAAGATAAGTGCAGCAACCGGCGGGACTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAAGGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAgg                                                                                                                  <  2:108061/244‑1 (MQ=255)

                                                                                                                          gcccAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTAGACGTAACCACGGCCATTCGGCGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGACAGCTGACGTTGTTCGCCGGGAATGTTTTCCGTTGCCTGGAGTAATAAACGGCTCGCCCTAGGGAGATTAAAAAATACGTGCGCCAGCAAAACACCTTTCAGGcaggaggggg                                                                                                               >  2:421600/2‑242 (MQ=255)

                                                                                                                            ccAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAACGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAACCGACGTTGTACGCCGGGAATGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGCCCTCGGGAGAGTCAAAAATACGTTCCCCCGCCAAAAACCTTGCAGGCCGtgaggcga                                                                                                              >  2:302259/1‑243 (MQ=255)

                                                                                                                                    aGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGAGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGACGGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTAGTTCGCCGGGAAAGTTTTCCAGTGCATGGAGTAAAAAGAGGCTCGCCATCGGTAGATTAAAAAATACGTGAGCCAGCAAAATACCGTGCAGGCCGAAGGGCGAAAAGGTcc                                                                                                      >  1:550669/1‑250 (MQ=255)

                                                                                                                                       ccGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTc                                                                                                  <  1:730117/250‑1 (MQ=255)

                                                                                                                                          gTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCCCCGTAAACAAGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGGTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGAGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTTCAAGCCGGAAGGCGAAAAGGGCCAATCCa                                                                                                >  1:601635/1‑250 (MQ=255)

                                                                                                                                             gCGAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGAcc                                                                                            <  2:811311/251‑1 (MQ=255)

                                                                                                                                               gAAGCTGGCGAAACAGAGCATAAAGATAAGTGCAGCAACCGGCGGAATTTGTCGCCGTAACCACGGCCATTTGACGAAGCGGCAAAAAAGCCAGCCGCGCCTCCCCAGCTGGGGGGCAAGCTGAAGGTGTTCGCCGGGGATGTTTTTCAGTGCCTGGAGTCATAAGCGGCCCGCCACCGGTAGCTTTAAAAATAAGTTCGCCTGTAAAATACCCTGCAGGCCGTAGGGGGAACAGGGCCACTCCAGACCgg                                                                                          >  2:283797/1‑250 (MQ=255)

                                                                                                                                                          aaCAGAGCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCAATTCGACGTAGCGGAAAAAATGCCAGGCGCGCCTCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCAGGGGATGGTTTCACGTGCCTGGAGTAATACGCCGCCCGCCATCGGAAGATTAAAAAATACGTGCGACAGCAAAATACCTTGAAGGACGTAGGGGGAAAAGGTCCACTCCAGACCACGCGATTGGc                                                                                 >  1:67635/1‑249 (MQ=255)

                                                                                                                                                           aCAGAGCATAAAGAGTAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGGGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTGGATTAAAAAATACGTGCGCCATCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGGCCACTCCAGACCGAGCGATTGGCaga                                                                              <  1:440769/251‑1 (MQ=255)

                                                                                                                                                              gagCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCTATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGTCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATCCCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAgtg                                                                           <  1:172318/251‑1 (MQ=255)

                                                                                                                                                              gagCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCagaggg                                                                           >  1:19785/1‑249 (MQ=255)

                                                                                                                                                              gagCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGGAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCagaggg                                                                           >  1:455465/1‑249 (MQ=255)

                                                                                                                                                              gagCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCGGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAgtg                                                                           <  2:693665/251‑1 (MQ=255)

                                                                                                                                                                gCATAAAGATAAGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGCCCGAGAGATTGGCAGAGTGtt                                                                         >  1:883426/1‑251 (MQ=255)

                                                                                                                                                                 cATAAATATAAGTGCAACAACCGGCGGGATTTGTCGCCGGAACCACGGCCATCCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCGTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGtt                                                                         <  2:409018/250‑1 (MQ=255)

                                                                                                                                                                      aGATAAGTGCAACAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGGAAAAATGCCAGCCGCGCAGCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTGCCAGTGCCTGGAGTAATAAGCGGCTCGCCAGCGGTAGATTAAAAAATACGTGCGCCAGCAAAAGCCCGTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagc                                                                   <  2:74493/251‑1 (MQ=255)

                                                                                                                                                                          aagtgcagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGGCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCcagccagc                                                               >  1:734036/5‑251 (MQ=255)

                                                                                                                                                                           agtgcagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCagccagcc                                                              <  1:825788/248‑1 (MQ=255)

                                                                                                                                                                               cagcaACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTgg                                                          <  1:722793/251‑1 (MQ=255)

                                                                                                                                                                                  caACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCg                                                        >  1:841638/1‑250 (MQ=255)

                                                                                                                                                                                    aCCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCATCCGCGCATCCCAAGCTGGTCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGGAGATTAAAAAATACGTGCGCCAGCAAAATACCGTGCAGGCCGTAGGGCGAAAAGGGCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAc‑c                                                    <  2:96920/251‑1 (MQ=255)

                                                                                                                                                                                      cggcggGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGGAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑Ca                                                   >  1:104803/1‑250 (MQ=255)

                                                                                                                                                                                       ggcggTATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATa                                                 <  1:941766/251‑1 (MQ=255)

                                                                                                                                                                                       ggcggGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CAt                                                  >  1:354854/1‑250 (MQ=255)

                                                                                                                                                                                          gggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATaaac                                              >  1:442617/1‑248 (MQ=255)

                                                                                                                                                                                          gggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGGGTTGCCAGCCAGCCCTGGCGAC‑CATaaac                                              >  1:957837/1‑248 (MQ=255)

                                                                                                                                                                                           ggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACg                                             >  1:594716/1‑251 (MQ=255)

                                                                                                                                                                                            gATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGc                                            <  1:845810/251‑1 (MQ=255)

                                                                                                                                                                                              tttggtgcgttacccACGGCCTTGCGACGAAGCGGAAAACAGGCCAGCCGCGCAGCCCAAGCTGGGCGGCAAGATGACGTTGTTCGCCGGGAATGTTTTCCAGGGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCATCCAGCCCTGGCGAC‑CATAGACGCt                                           <  2:835479/237‑1 (MQ=255)

                                                                                                                                                                                                   cGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTA‑TAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAAt                                     <  1:278189/250‑1 (MQ=255)

                                                                                                                                                                                                       gTAACCACGGCCATTCGACGAAGCGAAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTACCTGGAGTAATAAGCGGATCGCCATCGGTAGATTAAAAAAAACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGAGCCACTCCAGCCCGAGAGATTAGCAGCGTGGAGCCAGCCCGCCCAGGCGCC‑Cagagcagcgacgaatgccg                                 >  2:827866/1‑233 (MQ=255)

                                                                                                                                                                                                         aacaaCGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGACGGTCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGGAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGAGCGCAAAGGTCCACTCCAGCCCGAGCGATTGGACGAGTATTGCCAGCCGGCCCGGCCGAC‑CAACGACGCTACGaagaccga                                >  2:711373/5‑244 (MQ=255)

                                                                                                                                                                                                         aaCCACGGGCATTCGACGAAGCGGAAACAATGCCAGCCGCGCATCCCACGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTACAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAGAAATAGGTGCGAAAGAAAAATTCCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCATACCGAGAGATTGGCAGTGTGATGGCAGCCAGCCCTGGCTCC‑CATATACGCTAAGAATACCGa                                >  2:596502/1‑250 (MQ=255)

                                                                                                                                                                                                         aaCCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGAAGGGCGAAAAGGTCCACTCCAGACCGCGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGa                                >  1:617608/1‑250 (MQ=255)

                                                                                                                                                                                                         aaCCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCGGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGGGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGa                                <  2:764371/250‑1 (MQ=255)

                                                                                                                                                                                                                  cATTCGATGAAGCGGAAAAAATGCCAGGCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCAGCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGTCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGcaac                      <  2:678659/251‑4 (MQ=255)

                                                                                                                                                                                                                  cATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGcaac                      <  1:310878/251‑4 (MQ=255)

                                                                                                                                                                                                                    ttCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCCGAGTGTTGCAAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCAAcc                     >  1:44731/1‑250 (MQ=255)

                                                                                                                                                                                                                     tCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCAACCa                    <  2:354854/250‑1 (MQ=255)

                                                                                                                                                                                                                     tCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCAACCa                    <  2:767680/250‑1 (MQ=255)

                                                                                                                                                                                                                        acgaAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAACCAGCAACCAgcac                <  2:757967/251‑1 (MQ=255)

                                                                                                                                                                                                                         cgaAGCGGAAAAAAGGCCAGCCGCGCATCCCACGGGGGGCGGCAAGCTTCCGTGGGTCGTCGGGAATGTTTGCCAGTGCCTGGAGTAATAAGCGGGTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTCGCAGAGTGGTGCCAGCCAGCCCTCGTGAC‑CATAGAGGCTGAGAATGCAGGAATCAGCAACCAGCAcc               <  2:412251/251‑1 (MQ=255)

                                                                                                                                                                                                                             gCGGAAAAAATACCTGCCGCGCATCCCAATCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGCGTAATAAGAGGCTCGCCCTCGGTCGATTAAAAAATACGTGCGCCAGAAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCAAGCGATAGGCAGAGTGTGAGCAGCCAGCCCCGGCGGCACA‑AGTCGCTGAGAATGCCGCAAACAGCAAcccgaaccgg             >  2:417515/1‑241 (MQ=255)

                                                                                                                                                                                                                              cGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAACACAGCAACCAGCACCgggcg          >  1:172981/1‑249 (MQ=255)

                                                                                                                                                                                                                              cGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGACAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCggggg          >  2:416506/1‑249 (MQ=255)

                                                                                                                                                                                                                               ggAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGGTCGCCGGGAATGTTTTCCAGGGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCGTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGggagg         <  2:569713/251‑1 (MQ=255)

                                                                                                                                                                                                                                gAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGggagg         >  1:832957/1‑250 (MQ=255)

                                                                                                                                                                                                                                gAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTGGTTCGCCGGGAATGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGggagga        <  2:150563/251‑1 (MQ=255)

                                                                                                                                                                                                                                  aaaaaGGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGTAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTGGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCAACCAGCACCGGGAGGAtt      <  2:919819/251‑2 (MQ=255)

                                                                                                                                                                                                                                      aTGCCAGCCGCGCCTCCCAAGCTGGGCGGCAAGCTGGCGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTAAAAAATACGTGCGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAAAGCCGAAACCAGCAACCAGCACAgggcgggagaagg  >  2:694755/1‑241 (MQ=255)

||

AGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACATGGGCCAGCAAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGAC‑CATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGG  >  NC_000913/73746‑74226

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