BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,907	T→A	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,907	0	T	A	95.5%	752.5 / ‑4.0	227	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base T (0/10); new base A (113/101); total (114/113)
Fisher's exact test for biased strand distribution p-value = 7.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

ACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑TTTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTC  >  NC_000913/4558663‑4559154

 ccGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACt                                                                                                                                                                                                                                                      <  2:1003311/251‑1 (MQ=255)

   gtggCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGt                                                                                                                                                                                                                                                    <  1:967970/248‑1 (MQ=255)

   gTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGt                                                                                                                                                                                                                                                    <  2:364717/251‑1 (MQ=255)

   gTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGt                                                                                                                                                                                                                                                    <  1:72401/251‑1 (MQ=255)

     tGCCGTCGGAGCTGAGGGCGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCg                                                                                                                                                                                                                                                  >  2:651320/1‑251 (MQ=255)

       ccggcggggCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGcc                                                                                                                                                                                                                                                <  2:784831/243‑1 (MQ=255)

              aGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTgg                                                                                                                                                                                                                                        >  2:69430/1‑251 (MQ=255)

                  gAGGGTGACGCTTGCCAGCGTAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAac                                                                                                                                                                                                                                     <  1:4305/250‑1 (MQ=255)

                   aGGGTGACGCGTGCCAGCGGAATTCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCGCAAACAACGGTACGTTGCGGTTAACGTTGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAacac                                                                                                                                                                                                                                   <  1:336526/251‑1 (MQ=255)

                      gTGACGCGGGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGt                                                                                                                                                                                                                               <  2:735204/251‑1 (MQ=255)

                       tGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTg                                                                                                                                                                                                                              >  2:1125797/1‑251 (MQ=255)

                        gACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTg                                                                                                                                                                                                                              >  2:592789/1‑250 (MQ=255)

                                        aTACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTAccgcc                                                                                                                                                                                                            <  1:462048/251‑1 (MQ=255)

                                           cccGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCACTGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCCGGTTTAccgcc‑g                                                                                                                                                                                                          >  2:904668/1‑249 (MQ=255)

                                                  gAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGTTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑TGGCAAAc                                                                                                                                                                                                   <  1:764579/249‑1 (MQ=255)

                                                     cGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGcc                                                                                                                                                                                               >  1:594581/1‑250 (MQ=255)

                                                       gcgcGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAAc                                                                                                                                                                                            <  1:691106/251‑1 (MQ=255)

                                                           gggCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCgg                                                                                                                                                                                        >  1:950889/1‑251 (MQ=255)

                                                           gggCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCgg                                                                                                                                                                                        >  2:874276/1‑251 (MQ=255)

                                                               aaTACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGAtt                                                                                                                                                                                    <  2:160238/251‑1 (MQ=255)

                                                               aaTACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGAtt                                                                                                                                                                                    <  1:626834/251‑1 (MQ=255)

                                                                  cccTTCGGCCGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGGGCCGCCTTTGCGCGCGACCTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCt                                                                                                                                                                                  <  2:524057/249‑1 (MQ=255)

                                                                  cccTTCGGCAGGGGCTCCCGTTCCGTCAAGGCTGCTGGTGATACCGATGGTGCCGCCTTGGCGCGCGATCGCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑AGCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTACCGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCCGGTTTACCGCC‑GAGCAATCCTCCAACGCGTGATTCTg                                                                                                                                                                                 <  2:856040/250‑1 (MQ=255)

                                                                  aCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCt                                                                                                                                                                                  <  1:375455/250‑1 (MQ=255)

                                                                       cGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTgccgcc                                                                                                                                                                            <  1:497123/251‑1 (MQ=255)

                                                                         gCAGGGGCGGCCGGTTCGGCAATGAGGATGGTGATTTCGATGGTGCCGCCTTGGCGCGCGACCCCCAGC‑GCTGGTTCAAACAACGGTACGTTGCGGTTAACGGGGG‑GCGGCAGCAGTGTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTGATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑AGGGGGAACAC‑GGTA‑CCGGCAGGGTTCCCGCC‑GAGCAATCCGCCAGCGCGGGTTTCtaaggccat                                                                                                                                                                          <  2:329443/251‑9 (MQ=255)

                                                                           aGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCctat                                                                                                                                                                        >  2:554603/1‑251 (MQ=255)

                                                                            ggggCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATAtt                                                                                                                                                                       <  1:608544/251‑1 (MQ=255)

                                                                                cagccGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCTCCGTGGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAGCGTGTG‑TCGGCAGCATCTGGCTGTTCGTCCC‑ATCGCAGTTGTCCAGCAGGTCATAAATATGCTGGAACGCC‑ATTTGCCTGTCGCCC‑ATGTGGAACAC‑G‑TGTACGCCAGGTTT‑CCGCCCGAGCAAACTGCCTACGCGGGATTCTGCCGCCATATTGGcc                                                                                                                                                                   <  2:299008/248‑1 (MQ=255)

                                                                                 gACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAATCCGCCAACGCGGGATTCTGCCGCCATATTGGCCa                                                                                                                                                                  <  1:1040981/251‑1 (MQ=255)

                                                                                 gACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGcc                                                                                                                                                                   >  1:581366/1‑250 (MQ=255)

                                                                                 gACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGcc                                                                                                                                                                   >  1:1167890/1‑250 (MQ=255)

                                                                                       ttCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGAt                                                                                                                                                            >  2:448843/1‑251 (MQ=255)

                                                                                        tCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTGGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGAt                                                                                                                                                            <  1:284762/250‑1 (MQ=255)

                                                                                          gTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATaa                                                                                                                                                          >  1:1011993/1‑250 (MQ=255)

                                                                                                gtggttgGTGATAGCGATGTGGCCGCCTTTGCGCGCGAAC‑CCAGCCGCCGGCTCAAACAACGGTACGTTGCGGTTAACGGGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTCCCAGCAGTTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACCC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTcc                                                                                                                                                    <  2:227168/245‑1 (MQ=255)

                                                                                                   gatgGTGATACCGAGGGTGCCGCCTTGGCGCGCCAACTCCAGC‑GCATGCTCAAACAACGGTCCGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGACCGTCAC‑ATCGCAGTTTTCCAGCAGTTCATAAATAGGCGGTAACGCC‑ATTTTCCTGTCGCCC‑ATGTGGAACCC‑GGTG‑ACGCCAGGTTTCCCGCC‑GAGCCAACCGCCAGCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTg                                                                                                                                                <  1:52272/249‑1 (MQ=255)

                                                                                                          atatCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGc                                                                                                                                            >  2:458896/1‑248 (MQ=255)

                                                                                                              cGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTGGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATGGGCAC‑ATCGCAGTTTTCCAGCAGGTAATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACg                                                                                                                                      <  2:576359/250‑1 (MQ=255)

                                                                                                               gATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGAGGAACAC‑GGTG‑ACGCCAGGTTTCCCGCC‑GAGCAAACCGCCAACGAGGCATTCTGCCGCCCTATTGGCCAGGGGATAAACGTCCGGTGAGGCAGAACGGt                                                                                                                                    >  2:872902/1‑251 (MQ=255)

                                                                                                                aTGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTg                                                                                                                                   >  2:716786/1‑251 (MQ=255)

                                                                                                                  ggTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGAt                                                                                                                                 >  2:497360/1‑251 (MQ=255)

                                                                                                                   gTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                                >  1:531327/1‑251 (MQ=255)

                                                                                                                     gccgccTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                                >  2:818705/1‑249 (MQ=255)

                                                                                                                      ccgccTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                                >  1:869502/1‑248 (MQ=255)

                                                                                                                            ttGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATgg                                                                                                                        <  1:1069851/250‑1 (MQ=255)

                                                                                                                            ttGCGCGCGAACGCCAGC‑GCCTGCTCAACCAACGGTACGTTGCGGTTAACGTGGG‑TCGGCGGCAGCTTGCTGATCGGCAC‑ATAGCAGTTTGCCAGAAGGTATTAAAGAGGCTGTAACGCC‑ATTTTACTGCCGCCC‑ATGTGGAACAC‑GGGG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATg                                                                                                                         <  2:1105404/249‑1 (MQ=255)

                                                                                                                                         ccAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGcc                                                                                                           >  1:539216/1‑250 (MQ=255)

                                                                                                                                                 cggCTCAAACAACGGTACGTGGCGGTTAACGTGGG‑GCGGCAGCAGTTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCacac                                                                                                   <  2:1048253/249‑1 (MQ=255)

                                                                                                                                                 cggCTCAAAAAACGGTGGGTTGCGGTTAACGGGGG‑TCGGCAGAGGCTGGCTGACGGGCCG‑AGCGCAGTTTTCCGGAGGGTCATAAACAGGCTGTAACGCC‑ATTTTACGGCCGCCA‑GTGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCacac                                                                                                   <  2:163492/249‑1 (MQ=255)

                                                                                                                                                     ccaaccaacGGTCCGTTGCGGTTAACGGGGG‑TCGGTAGCCGCTTGCGGACCGCCAC‑ATCGTAGTTTGCCAGCAGGTCATAAATGGGCCGTATCGCC‑ATTTTACTGTCGCCC‑ATGTGGCACAC‑GGTG‑CCGCCAGGTTTCCCGCC‑GAGCAAACCGCCAACGCGGGATTCGGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatc                                                                                               <  2:689764/246‑1 (MQ=255)

                                                                                                                                                      cAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatcg                                                                                              <  1:1058427/251‑1 (MQ=255)

                                                                                                                                                        aacaacGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatcgat                                                                                            >  1:1120303/1‑251 (MQ=255)

                                                                                                                                                           aacGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGAt                                                                                                                                                            <  1:96149/184‑1 (MQ=255)

                                                                                                                                                           aacGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGGGAt                                                                                                                                                            >  2:96149/1‑184 (MQ=255)

                                                                                                                                                           aacGGTACGTTGCGGTTAACGGGGG‑CCGGCAGCAGCTTGCTGATCGGCCA‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGACCAC‑GGTG‑ACGCCGGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGgcgc                                                                                                                     <  2:537796/223‑1 (MQ=255)

                                                                                                                                                             cGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTCCCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCg                                                                                       <  1:716786/251‑1 (MQ=255)

                                                                                                                                                              ggaaCGTTGCGGTTAACGGGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCACACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGGGATAACCGTCCGGTGCGGCAGAACGGTGGTTAGAGATGGCGCATTTCACTCCAATCACACGATCGATAATCGc                                                                                      >  2:23064/4‑251 (MQ=255)

                                                                                                                                                              ggTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGGGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCCCACGATCGATAATCGc                                                                                      >  2:731603/1‑251 (MQ=255)

                                                                                                                                                                 gcgtggCGGTTA‑CGTGGGGTCGGAATCAGCTTGCTGATGGGCAC‑ATGGCTGTGTTCCAGCAGGGAATAAATGGGCTGTACCCCC‑ATTTTACTGTCGCCC‑ATGTGGAACCC‑GGGG‑ACGCCAGGTTTACCGCC‑GAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCTATAATCGCCa                                                                                    <  2:950889/245‑1 (MQ=255)

                                                                                                                                                                    ttGCGGTTAACGTGGG‑CCGGCAGCAGCTTGCTGGTCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTCACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGAATCCTGCCGCCATATTGGCCAGGAGATAAACGGACGGTGCGGAAGAACGGTGATCAGAGATGGCGCATTTCCCGCCAATCACACGATCGATACTCGCCACGTc                                                                                >  1:388687/1‑251 (MQ=255)

                                                                                                                                                                             ccGTGGG‑TGGGCGGCAGCTGGCTGATCG‑CACCATGGAATTTTTCCAGCAGGTAATAAAT‑GGCCGTAACCCCCATTTTACTGTCGCCC‑ATGTGGACCCC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCCCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCAc                                                                       <  2:1069570/250‑1 (MQ=255)

                                                                                                                                                                             aCGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCGGTTTTCCAGCAGGTCATAAAGAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCAc                                                                       <  2:496475/251‑1 (MQ=255)

                                                                                                                                                                             aCGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTCCCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCAc                                                                       <  1:1102711/251‑1 (MQ=255)

                                                                                                                                                                             aCGTGGG‑TCGGCAGCAGCTTGCTGAGCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCAc                                                                       <  1:137313/251‑1 (MQ=255)

                                                                                                                                                                              cGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACg                                                                      >  2:89410/1‑251 (MQ=255)

                                                                                                                                                                                 ggg‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGa                                                                    >  2:181827/1‑250 (MQ=255)

                                                                                                                                                                                 ggg‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGa                                                                    >  2:1100234/1‑250 (MQ=255)

                                                                                                                                                                                     cgggCAGCAGCTGGCTGACCGGCAC‑TTCGCAGTTTCCCGGCGGGTAAAAAATAGGCTTTAACGCC‑ATTTTAC‑GGCGCCCCATGTGGACCAC‑GTGG‑CCGCCGGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGGCCGGTGCGTCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAccc                                                                <  2:312257/249‑1 (MQ=255)

                                                                                                                                                                                       ggCAGCAGCTTGCTGATCGGCAC‑AACGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTCCCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAAc                                                                                                                                       <  1:144735/178‑1 (MQ=255)

                                                                                                                                                                                         cagcagCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTa                                                             >  2:594117/1‑250 (MQ=255)

                                                                                                                                                                                          agcagcTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAAt                                                           >  1:827131/1‑251 (MQ=255)

                                                                                                                                                                                               cTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGc                                                      >  1:1048545/1‑251 (MQ=255)

                                                                                                                                                                                                 tGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTACCGCC‑ATTTTACTGTAGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCA‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCAGGTGCGGCAGAACGGTGATCAGAGATGGCGAATTTCACGCCAATCACACGATCGGTAATCGCCACGTCTTTTTTCCCGGAAACCGGAATGGTGCgg                                                    >  2:981791/1‑251 (MQ=255)

                                                                                                                                                                                                 tGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCgg                                                    >  1:733586/1‑251 (MQ=255)

                                                                                                                                                                                                        cGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCGGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCGGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGTAACCCGTAATGGTGCGGGAAGGGa                                             <  1:75808/251‑1 (MQ=255)

                                                                                                                                                                                                            cacagcgaaGTTTTCCAGCAGGTCAAAAAGGTGCGGAAACGCC‑ATTTTCCTGCCCCCC‑ATGTGGAACAC‑GGGG‑ACGCCGGGTTTCCCGCC‑GAGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACAt                                          <  2:1040561/243‑1 (MQ=255)

                                                                                                                                                                                                            ac‑aTCGCAGTTTACCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGGCGCCC‑ATGTGGAACCC‑GGGG‑ACGCCGGGTTTCCCGCC‑GAGAAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATg                                         <  2:934668/251‑1 (MQ=255)

                                                                                                                                                                                                                  gCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATa                                      >  1:351757/1‑249 (MQ=255)

                                                                                                                                                                                                                    aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAgcg                                  >  1:913591/1‑251 (MQ=255)

                                                                                                                                                                                                                      ttttCCAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGcc                                >  2:376109/1‑251 (MQ=255)

                                                                                                                                                                                                                       tttCCAGCAGGTCATAAATAGGCTGTAACGCC‑TTTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGcc                                <  1:181827/250‑1 (MQ=255)

                                                                                                                                                                                                                       tctCCAGCGGTTCATAAATAGCCCGTACCGCC‑ATTTTCCTGTCGCCC‑GTGGGTAACGC‑GGTG‑ACGCCATGTTTCCCGCC‑GAGCAACCCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGTGAAGGGACATGATAAGCGcc                                <  2:215063/248‑1 (MQ=255)

                                                                                                                                                                                                                          ccAGCAGGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGCACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTc                            >  2:941197/1‑251 (MQ=255)

                                                                                                                                                                                                                           cagcagGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ATGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTc                            >  1:965204/1‑250 (MQ=255)

                                                                                                                                                                                                                            agcagGTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTc                            >  2:1167555/1‑249 (MQ=255)

                                                                                                                                                                                                                                 gTCATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATcc                     >  1:940752/1‑251 (MQ=255)

                                                                                                                                                                                                                                  tCATAAATAGGCTGTAAAGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATcc                     >  1:214121/1‑250 (MQ=255)

                                                                                                                                                                                                                                   cATAAATAGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATcc                     >  2:429017/1‑249 (MQ=255)

                                                                                                                                                                                                                                    aTAAATAGGCTGTAACGCC‑TTTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTCCCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAgg                  <  1:1138002/251‑1 (MQ=255)

                                                                                                                                                                                                                                    aTAAATAGGATGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTACACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAgg                  >  1:528706/1‑251 (MQ=255)

                                                                                                                                                                                                                                     aaaaTAGGCTGTAACGCC‑ATTTTACGGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgc                 <  1:904668/250‑1 (MQ=255)

                                                                                                                                                                                                                                          aGGCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTg             <  1:397270/250‑1 (MQ=255)

                                                                                                                                                                                                                                           ggcgtaaccGCC‑TTTTTCCTGCCGCCC‑ATG‑GGGACACCGTTG‑AGGCAGGGTTTACCGCC‑GAGCAAACCGCCACCGGGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGAt           <  2:351757/243‑1 (MQ=255)

                                                                                                                                                                                                                                           ggCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGGTAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACCTGATAAGCGCCGGTCAGCATCCAGGCGCTGAt           >  2:526045/1‑251 (MQ=255)

                                                                                                                                                                                                                                            gCTGTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGagg          >  1:679531/1‑249 (MQ=255)

                                                                                                                                                                                                                                               gTAACGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCct       >  2:409237/1‑251 (MQ=255)

                                                                                                                                                                                                                                                 aaCGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGTC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCACCCAGGCGCTGATGCctt      <  1:218341/250‑1 (MQ=255)

                                                                                                                                                                                                                                                 aaCGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttc     <  1:554603/251‑1 (MQ=255)

                                                                                                                                                                                                                                                 aaCGCC‑ATTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttc     <  1:1154023/251‑1 (MQ=255)

ACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGC‑GCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGG‑TCGGCAGCAGCTTGCTGATCGGCAC‑ATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCC‑TTTTTACTGTCGCCC‑ATGTGGAACAC‑GGTG‑ACGCCAGGTTTACCGCC‑GAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTC  >  NC_000913/4558663‑4559154

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