BRESEQ :: Evidence

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

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

Only 100 of 104 total aligned reads displayed.

GAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955760‑3956249

gAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGAAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AAAAAC‑TGATGCCGC‑AAAT‑CACAGTGATAGCCGTAGAAGCGGGAGACTGCGCCTGCCTGAAAGCAACGCCGAATGCggg                                                                                                                                                                                                                                                  >  1:399547/1‑251 (MQ=255)

 aaCTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑CAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCgggg                                                                                                                                                                                                                                                 >  1:113103/1‑250 (MQ=255)

     gTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATc                                                                                                                                                                                                                                             <  1:548137/251‑1 (MQ=255)

       cacaGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATcc                                                                                                                                                                                                                                            <  1:499079/250‑1 (MQ=255)

          agcagcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CCAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGGCATCCGG‑tt                                                                                                                                                                                                                                       >  1:201030/1‑251 (MQ=255)

          agcagcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑tt                                                                                                                                                                                                                                       <  1:561798/251‑1 (MQ=255)

             agcagGGGTTCACCGGGGTGCCGCCGGTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCGGGCGCGGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCGGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCGGAAAGCAGCGCTGAATGCGGGGCATCCGG‑TTGAt                                                                                                                                                                                                                                    <  2:359457/251‑1 (MQ=255)

                   ggTTCACCTGGGGGCCGCCGTTCGACCAGCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCTCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCTGAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCGGAATGCGGGTCATCCGG‑TTGATCTGCcg                                                                                                                                                                                                                              <  2:542549/251‑1 (MQ=255)

                              gTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AAAAAC‑TGATGCCGC‑CAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCt                                                                                                                                                                                                                   >  1:9810/1‑251 (MQ=255)

                                gccgccGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTAt                                                                                                                                                                                                                 <  1:582867/251‑1 (MQ=255)

                                gccgccGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTAt                                                                                                                                                                                                                 <  1:489647/251‑1 (MQ=255)

                                     cgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCt                                                                                                                                                                                                            <  1:429562/251‑1 (MQ=255)

                                      gTTCGACCAGCCGAGGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTg                                                                                                                                                                                                           <  2:490770/251‑1 (MQ=255)

                                           aCCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGc                                                                                                                                                                                                      <  1:340130/251‑1 (MQ=255)

                                            ccATCCGATGGTGATTGCCGGGCAAGGCCCGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCg                                                                                                                                                                                                     <  1:164438/251‑1 (MQ=255)

                                                ccGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGc                                                                                                                                                                                                 <  2:240922/251‑1 (MQ=255)

                                                ccGATGGTGAGTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGc                                                                                                                                                                                                 <  2:465895/251‑1 (MQ=255)

                                                    gggTGAGTGCCGGGCAAGGCACGGTGGCGCTGGAACTTCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCGGGCTGCTGGCGTGGCGGTGCTGATCG‑AACAAC‑TGATGCCGC‑AAAG‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGGAGCGGTa                                                                                                                                                                                             <  2:289723/250‑1 (MQ=255)

                                                     ggggATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCA‑AACAACCTG‑TGACGCCAAAT‑‑AAAGTGATCGGCGTAGAAACGGAAGACTTCGCCCTCCTGAAAGACGCGGTGAATGCGGGGCCTACGG‑TTGAGTTGCCGCGCCGAGGGCTATTTTCTGAAAGGGGAGCGGTa                                                                                                                                                                                             >  2:348919/4‑250 (MQ=255)

                                                     ggTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGCAGCCCATCTCGACCGCGTATTTGCGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCCGGAATGGGCA‑CACAAC‑TGTGGCCGC‑AAAT‑CCAAGTGCACGCCGTAGAAAAGGAAGTCTCCTCCTGCCTGAAAGCAGAGCTGAAATCGGGGCCTACAG‑TTTGTATGCCCCGCGTAGGGCGATTTGTTGGAGGCGTAGCggca                                                                                                                                                                                             >  1:584961/1‑248 (MQ=255)

                                                      gTGATTGCCGTGAAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCGCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTaaa                                                                                                                                                                                           <  2:158415/251‑1 (MQ=255)

                                                            gCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGGTGCCGC‑AAAT‑CAAAGGGATCGGCGGAGAAGCGGGAGACTCCGCCCGCCTGAAAGCAGCGCTGAATGCGGGGCATACGG‑TGGATCCGCCCGGCGGACGGGTATTTGCTGGAGGGGTaaggggaaaaagcat                                                                                                                                                                                     >  2:104815/1‑237 (MQ=255)

                                                              cGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCATGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑CAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑GTGATCTTCCGCGGGTAGGGCTATTTGCTGAAGGCGGAGCGGGAAAACGCATCg                                                                                                                                                                                   >  1:430320/1‑251 (MQ=255)

                                                              cGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTACTGATCA‑AACAAC‑TGATGCCCC‑AAAT‑CAAAGTGATCGCCGTTGAAGCGGAAGACTCCGCCTGCCTGAAAGAAGCGCAGAACGCGGGTCCTCCGG‑TTGATCTACCGCGCGTAGGGCTATTTGCTGTAGGCGTAGCGGTAAAAAGCATCg                                                                                                                                                                                   >  1:61251/1‑251 (MQ=255)

                                                                ggCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGCTGCCGC‑AAAT‑CAAAGGGATCGCCGTAGAAGCGGAAAACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGG‑GTGATCTGCCGCGCGGAGGGCTATTTGCTGAAGGCGGGGCGGTAAAAAGGATCggg                                                                                                                                                                                 >  1:296489/1‑250 (MQ=255)

                                                                  cAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGGAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTg                                                                                                                                                                                >  1:124303/1‑250 (MQ=255)

                                                                  cAAGGCACGCTGGCGCTGGAACTGCGCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGa                                                                                                                                                                               <  1:94315/251‑1 (MQ=255)

                                                                   aaGGCGCGCTGGCGCTGGAACTGCTCCAGCAGGACGCACATCTCGACCGCGTCTTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGGGCTGATTC‑AACAAC‑TGATGCCGC‑AAATACAAA‑TGATCGCCGGCGAAGCGGGAGAAGCCGCCCGCATGCAAGCAGCACAGAATGCGGGACCCCAGG‑TTGATCTGCCGCGGCTCGGGCTTCTTGCTGAAGGCGTAGCGGTAGAACGcctcagagcc                                                                                                                                                                              >  2:113822/1‑242 (MQ=255)

                                                                    aGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑CAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑GTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAAAGCATCGGTGACg                                                                                                                                                                             >  1:523081/1‑251 (MQ=255)

                                                                     ggCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCTCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGa                                                                                                                                                                            <  1:574556/251‑1 (MQ=255)

                                                                       cACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTTATTACAACA‑C‑TGGTGCCGC‑AAAT‑CAAAGTGATCGCCGTAGGAGCGGAAGACTTCGCCCACCTGAAAGCAGCGCTGAAAGCGGGGCATCCGGGTTG‑TCTGCCGCGCGTAGGGCTATTTGTTGACGGCGTAGCGGTTAAACGCATCGGTGTCGaa                                                                                                                                                                           >  1:247064/1‑250 (MQ=255)

                                                                             ggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTcc                                                                                                                                                                    <  1:333516/251‑1 (MQ=255)

                                                                             ggcgctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTc                                                                                                                                                                     <  1:324109/250‑1 (MQ=255)

                                                                               ggcgggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGGGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGt                                                                                                                                                                  <  2:479001/247‑1 (MQ=255)

                                                                                 ctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑CGGTGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCCGCCTTGAAGCAGCGCTGAATGCGGGTAATCCGG‑TTGATCTGCCGCGCCTAGGGCGCTTTGCTGAAGGCGGAGCGGTAAAAAGCATCGGGGAAGGAAACTTCCGttt                                                                                                                                                                >  2:505698/1‑251 (MQ=255)

                                                                                     aaCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCACTCGGCGGCGGCGGTCTGGCTACTGGCGTGGCGGTGCTGCCCA‑AACAAC‑CGATACCCC‑ACAT‑CAAAGTGATCGCCGGAGAACCGGAAGACTCCGCCTGCCTAAAAGCAGCGCTGACTGCGGGTCGTCCGG‑GTGATCTCCCGCGCGGTGGGCTCTTTGCTGAAAGCGTAGCGGGAAACCGCATCGGTGTCGAAACCTTCCGTTTATGc                                                                                                                                                            >  1:338028/1‑251 (MQ=255)

                                                                                     aaCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATGTGTGCCAGGCGGCGGCGGCGGTCTGGCGGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGc                                                                                                                                                            <  2:516294/251‑1 (MQ=255)

                                                                                      aCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGcc                                                                                                                                                           <  1:238764/251‑1 (MQ=255)

                                                                                       cTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAA‑TGATGCCGC‑AAAT‑CAAAGTGCTCGCCGTAGAAGCGGAAAACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGGCATCCGG‑TTGGTCTGCCGCGCGTAGGGCTATTTGTTGAAGGCGTAGCGGTAAAAAGCATCGGTGACGAAACCTTCCCTTTATGcc                                                                                                                                                           >  2:82829/1‑250 (MQ=255)

                                                                                           tCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTa                                                                                                                                                               <  1:99901/242‑1 (MQ=255)

                                                                                            ccatcagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AAAAAC‑TGCTGCCGG‑AAAT‑CCAAGGGGACCCCGGAGAAGCGGAAGACTCCGCCCGCCTGAAAGCAGCGCTGAATGAGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTTCTGAAGGCGTAGCGGTAAAACGCCTCGGTGACGAAACCTTCCGGTTATGCCAggggt                                                                                                                                                     >  2:463269/5‑248 (MQ=255)

                                                                                               gcagGACGCCAATCTCGACCGAGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAA‑TGAAGCCGC‑AAAT‑AAAAGGGATCGCCGTAGAAGAGGAAGAACACGCCTGACTGAAAGCAGCGCTGAATGAGGGGAATCCGG‑GTGATCTGGCGCGCGGAGGGCAATTTGCTGAAGGCGGAGCGGGAAAAAGCATCGGGGACGAGACCTTCCGTTTATGCCAGGAggat                                                                                                                                                   >  2:97688/1‑247 (MQ=255)

                                                                                                 agGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTGCCGTTTATGCCAGGAGTAtctc                                                                                                                                                <  1:97688/251‑1 (MQ=255)

                                                                                                            ctcGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGGTGCCGC‑AAAT‑CCAAGGGGTCGCCGGAGAAGCGGGAGACTCCGCCCGCCTGAAAGCGGCGCTGGATGCGGGGCATCCGG‑GTGATCCGCCGCGGGCAGGGGTTTTTGTTGGAGGCGTCGCGGGTAAACGCATCGGTGACGAAACCTTCCGGTTATGCCCGGGGTATCTcgacgaactac                                                                                                                                      >  2:96753/1‑245 (MQ=255)

                                                                                                              cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatc                                                                                                                                       <  1:170271/247‑1 (MQ=255)

                                                                                                                gccGCTTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTTCTGTCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCGGCCGCGCGTAGGGCTATTTGCTTAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTCTATGCCAGGAGTATCTCGACGACATCATCAcc                                                                                                                                 <  1:127112/250‑1 (MQ=255)

                                                                                                                     gTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCCGATCC‑AACCAC‑TGGTGCCGC‑AACT‑CAAAGGGGTCGCCGGAGGAGGGGAAGGCTCCGCCTGCCCGAAAGCAACGCTGAAAGCGGGTCATCCGGGTTG‑TCTGCCGCCCGTCGGGCTATTTGCTGAAGGCGTCGCGGGAAAACGCATCGGTGGCGCAACCTTCCGGTTATGCCCGGGGTATCTCGACGGCATCATCCCCGTCgc                                                                                                                            >  2:113769/1‑250 (MQ=255)

                                                                                                                                   cggcggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGAGCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGAGCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCt                                                                                                              <  1:480491/251‑1 (MQ=255)

                                                                                                                                   cggcggcggcggTCTGGCGGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATc                                                                                                               <  1:485622/250‑1 (MQ=255)

                                                                                                                                      cggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCtgtg                                                                                                           <  1:382490/251‑1 (MQ=255)

                                                                                                                                       ggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATTCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGTGCTGAATGCGGGTCATCCGG‑TTGTTCTGCCGGGCGTAGGGTTATTTGCTGAAGTCGTAGCGGTAAAACGCATTGGTGACGTAACCTTCCGTTTATGCCATGGGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgc                                                                                                          <  1:235955/251‑1 (MQ=255)

                                                                                                                                        gcggcgtTCTGGCTGCTGGCGTGGCGGTCGTTATCA‑ATCAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCTTAGAAGCGGAAGTTTCCGCCTGCCTGAAAGCAGCGCGGAATGCGGGTCATCCGG‑TTGATCTGCCACGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGGGACGAAACCTGCCGTTTATGCCAGGAGGATCTCAACGACATCATCACCGCCGATAGCGATGCGATCTGTgcg                                                                                                         <  2:515012/251‑1 (MQ=255)

                                                                                                                                          ggcgtTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgcggc                                                                                                       <  2:229806/251‑1 (MQ=255)

                                                                                                                                          ggcggTCTGGCTGCTGGCGTGGCGGTGGTGATCA‑AACAAC‑TGACGCCGC‑AAAT‑CAAAGTGATCGCCGGAGAAGCTGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCACCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAGCTCGACGACATCATCACCGTCGATAGCGATGCGATCTTTgcggc                                                                                                       <  2:61385/251‑1 (MQ=255)

                                                                                                                                                 tGGCTGCTGGCGTGGCGGTGATGATCA‑AACAAC‑TGGTGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCCTCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCCTCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGGCGATAGCGGTGCGATCTGTGAGGCGATGa                                                                                                  >  2:421454/1‑249 (MQ=255)

                                                                                                                                                      gctgGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGAt                                                                                             <  2:303071/249‑1 (MQ=255)

                                                                                                                                                      gctgGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACCCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGAGCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGAttt                                                                                           <  1:96753/251‑1 (MQ=255)

                                                                                                                                                         ggcgtggcgGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCg                                                                                                                                                                                                                                           <  2:518354/105‑1 (MQ=255)

                                                                                                                                                         ggcgtggcgGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑GTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGGTTTAtt                                                                                        >  1:81616/1‑251 (MQ=255)

                                                                                                                                                         ggcgtggcgGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAAACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCg                                                                                                                                                                                                                                           >  1:518354/1‑105 (MQ=255)

                                                                                                                                                              ggcgGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGa                                                                                     >  1:405305/1‑249 (MQ=255)

                                                                                                                                                                gggTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGGAGACGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCTCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCACGACGACATCATCACCGTCGATAGCGATTCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAt                                                                                 <  2:394886/250‑1 (MQ=255)

                                                                                                                                                                cgGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAt                                                                                 <  1:369437/251‑1 (MQ=255)

                                                                                                                                                                cgGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAt                                                                                 <  2:252383/251‑1 (MQ=255)

                                                                                                                                                                               caac‑TGATGCAGC‑ACAT‑CAAAGTGCTCGCCGTAGCAGCAGCAGACTCCGCCTGCCTCACAGCAGCGCTGAATGAGGGTAATCCGG‑TTGCTCTGCCGCGCGTCGGGCTATTTGCTGAAGGCGTAGCGGTAAACCGCATCGGTGCCGAAACCTTCCGTTTCTGCCAGGAGTATCTCGACGACATCATCACCGTCGCTAGCGATGAGATCTGTACGGCGTTGAAGGCTTTATTCGTAGATGTGCGCGCggaggg                                                                   >  2:220408/1‑247 (MQ=255)

                                                                                                                                                                                aac‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGCGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCg                                                                  >  1:251581/1‑251 (MQ=255)

                                                                                                                                                                                 ac‑TGATGCAGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGACTGAAAGCAGCGCTGAATGCGGGACATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCgg                                                                 >  2:293417/1‑251 (MQ=255)

                                                                                                                                                                                    tGATGCCGC‑AAAA‑CAACGTGATCGACGAAGACGCGGAAGACTCCGCCTGCCTCAAAGCAGCCCTCAATACGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGTGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGAATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGa                                                                >  2:294735/1‑250 (MQ=255)

                                                                                                                                                                                     gATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTCTTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCTTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGGGCGCGCGGTGGCGGacc                                                              >  2:95087/1‑249 (MQ=255)

                                                                                                                                                                                     gATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGa                                                                <  2:369123/249‑1 (MQ=255)

                                                                                                                                                                                                      aGTGATCGCCGTAGAAGCGGAAGGCGCCGCCTGCCTGAAAGCAGCGCGGAATGCGGGTCGTCCGG‑TTGAGCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGGAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGGCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTGATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctgg                                               <  2:45385/251‑1 (MQ=255)

                                                                                                                                                                                                       gTGATCGCCGTAGAAGAGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggc                                              <  1:392229/251‑1 (MQ=255)

                                                                                                                                                                                                        tgagcGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATGTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAATCCTTCCGGTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGAGGAAGGATTTATTCGAAGATGTgcg                                                                            <  2:401269/216‑1 (MQ=255)

                                                                                                                                                                                                        tGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGCTGAAGGATTTATTCGAAGATGTgcg                                                                            >  1:401269/1‑220 (MQ=255)

                                                                                                                                                                                                        ggATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTAGTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGGGTAGGCCAGGAGTATCTCGACGACATCAACACCGTCGATAGCGATGCGAGCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCGGGCgtgctggcg                                             <  2:584961/250‑1 (MQ=255)

                                                                                                                                                                                                          aTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCCGGGGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGGAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgcc                                           >  1:47084/1‑250 (MQ=255)

                                                                                                                                                                                                          aTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑GTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTTCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgct                                           >  2:541758/1‑251 (MQ=255)

                                                                                                                                                                                                            cGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGCACCCTCTGGCGcgctggcgctgg                                         >  1:332277/1‑251 (MQ=255)

                                                                                                                                                                                                                gTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGACGATGTGCGCGCGGTGGCGGAGCCCTCTGGCGCGCTGGCGCTGGCggg                                     >  2:123470/1‑251 (MQ=255)

                                                                                                                                                                                                                gTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCggg                                     <  1:472937/251‑1 (MQ=255)

                                                                                                                                                                                                                    aaGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATg                                 >  1:521359/1‑251 (MQ=255)

                                                                                                                                                                                                                      gCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaa                               <  1:208962/251‑1 (MQ=255)

                                                                                                                                                                                                                          aaGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCGTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaaaaa                           <  2:400726/251‑1 (MQ=255)

                                                                                                                                                                                                                             tcTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGCGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGGATGAAAAAttat                        >  2:135402/2‑251 (MQ=255)

                                                                                                                                                                                                                              cTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAatat                        >  1:488545/1‑250 (MQ=255)

                                                                                                                                                                                                                                   cctgcctgAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGGATGAAAAAATATATCGcc                  >  1:296812/1‑251 (MQ=255)

                                                                                                                                                                                                                                     tgcctgAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCt                >  1:62778/1‑251 (MQ=255)

                                                                                                                                                                                                                                      gcctgAAAGCAGCGCTGAATGCTGGTCCTCCGT‑TTGAGCTGCCGCGCGTAGGGCTATTTTCTGAAGGCGTAGCTGTAAAACTCATCGGTTACTAAACCTTCCGTTTATGCCAGGAGTATCTCGACAACATCATCACCGTCTATAGCGAGTCGATCGGTGCGGCGATGAAGTATTTATTCGAAGATTTGCGCTCGGTGGCGGAACCCTCTGGCGCGCTGTCGTTGGCGGTAATGAAAAAATATATCGCCCTg               <  1:422493/251‑1 (MQ=255)

                                                                                                                                                                                                                                      gcctgAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTg               <  2:201030/251‑1 (MQ=255)

                                                                                                                                                                                                                                      gcctgAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTg               <  1:439861/251‑1 (MQ=255)

                                                                                                                                                                                                                                      gcctgAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTg               >  1:321840/1‑251 (MQ=255)

                                                                                                                                                                                                                                          gAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGTATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGcac            <  2:332277/250‑1 (MQ=255)

                                                                                                                                                                                                                                          gAAAGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTGCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGca             <  2:491220/249‑1 (MQ=255)

                                                                                                                                                                                                                                          gAAAGCAGCGCTGAATGCGGGTCATCCGG‑CTGATCTGCCGCGCGTAGGGCTATTTGATGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCAGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATAACGCCCTGccc            >  2:402553/1‑248 (MQ=255)

                                                                                                                                                                                                                                           aaaGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGGTTATGCCAGGAGTAGCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGGGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCacaa          <  1:269471/251‑1 (MQ=255)

                                                                                                                                                                                                                                            aaGCAGCGCTGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTCGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCccaac         >  1:459223/1‑251 (MQ=255)

                                                                                                                                                                                                                                               cagcGATGAATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTACGGCTCTTTGCTCAAGGCGTAGCGGTAAAACGCATCGGTGACGAACCCTTCCGTTTCTGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGACTGAAACAATATATCGCCCTGAACAAAAtt      >  2:245513/1‑251 (MQ=255)

                                                                                                                                                                                                                                                 gcgcTGAATGCGGGTCATCCGG‑TTGATCAGCCGCGCGTAGGGCTAGGTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTACGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGTATGAAAAAATATATTTCCCTGCACAACATTcg    <  2:493508/251‑1 (MQ=255)

                                                                                                                                                                                                                                                  cgcTGAATGCGGGTCAGCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCGCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTcgc   <  2:9810/251‑1 (MQ=255)

                                                                                                                                                                                                                                                   gcTGAATGCGGGTCATCCTG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTcgcg  <  1:355415/251‑1 (MQ=255)

GAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCA‑AACAAC‑TGATGCCGC‑AAAT‑CAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGG‑TTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955760‑3956249

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