BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,702	T→C	S158S (AGT→AGC)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,702	0	T	C	97.4%	197.7 / ‑6.8	79	S158S (AGT→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand): ref base T (1/0); major base C (42/33); minor base A (2/0); total (46/33)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

TGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACC  >  NC_000913/4618458‑4618948

tGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGGAAGCTCGACAAACTG‑GAATCCATCCCCGGCTTCGACATTTTCCCGGATGACAAACGTTTCCGCGAAATTATTAAAGACGGCGGCGTGGCGATTATCGG‑TTCA‑ACCAGCTCACTg                                                                                                                                                                                                                                                       >  1:111713/1‑251 (MQ=255)

 gAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCAc                                                                                                                                                                                                                                                         <  1:526620/248‑1 (MQ=255)

  aaTGGCCAGATTGTTGAAAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCATCGGTCATACTGGCGGGACGCTCGACCAACTG‑GAATCCATCCCCGGCATCGCCATTTTCCCGGGTGACAACAGTTTCCGCGAAATTATTTAAGAAGTCGGCGGGGGGATTATCGG‑C‑CAGAACAGCTCACTGG‑c                                                                                                                                                                                                                                                    >  1:447965/1‑251 (MQ=255)

     ggCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTcc                                                                                                                                                                                                                                                 <  2:178080/251‑1 (MQ=255)

           aTTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAATG‑GAATCCATCCCTGGCTTCGACATTTTCCCCGGTGACAACCGTTTCCCCGAAATTATTAAAGACGTCGGCGTTGCGGTTATCGG‑T‑CAGACCAGCCCACTgggcccggcggg                                                                                                                                                                                                                                            >  2:493954/1‑241 (MQ=255)

           aTTGTTGATAAACAGTCCACTGGTGGCGTCGGCGATGTGACTGCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCGCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGcgga                                                                                                                                                                                                                                           <  2:204571/251‑4 (MQ=255)

                tgATAAACATTCCACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGGAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑G‑CAGACCCGCTCACTGG‑CTcccgcgggtaaac                                                                                                                                                                                                                                      >  1:155832/1‑240 (MQ=255)

                  aTAAACATTCCACGGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCGTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGt                                                                                                                                                                                                                                    <  1:52124/251‑1 (MQ=255)

                           ccACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATCTTCCGATGATCTCTGGTCGCGGCCTCGGGCATACTGGCGGTACGCCCGACAAACTG‑GAATCCCTCCCTGGCTTCGGCATTTTTCCGGGAGACAACCGTTTCCCCGAAATTATTAAAGACGCCGGCGGGGCGATTTTCGG‑G‑CCGTCCACATCACCGG‑CTCCGGCGGATAAACGTTTCTAcgcg                                                                                                                                                                                                                           >  2:156873/1‑251 (MQ=255)

                           ccACTGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTAcgc                                                                                                                                                                                                                            <  1:545736/250‑1 (MQ=255)

                               gggTGGGGTCGTCGATGTGACTTCGCTGATTTTGGGGCCGATGGTCGCAGCCTGCGGCGTCTAGATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCAGCCCTTGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAATACGTCTGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACTTTTCTACTCGA‑cct                                                                                                                                                                                                                      <  2:258595/250‑2 (MQ=255)

                                ggTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTAAGCTCGACAAAATG‑GGAACCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGAACAGCTCACTGG‑GCCCGGCGGAAAAACGTTTCTAAGCGC‑CCCg                                                                                                                                                                                                                     >  1:591890/1‑251 (MQ=255)

                                ggTGGCGTCGGCGATGTGACTGCGCTGATGTTGGGGCCGATGGTCGCATCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCTCTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑cc                                                                                                                                                                                                                       <  1:223892/249‑1 (MQ=255)

                                   ggCGTCGGCGATGGGACTTCGCTGATGTTGGGGCCGATGGTTGCAGCCTGCGGCGGCTATATTCCGAGGATCTCTGGTCGCGGCCGCGGTCATACTGGCGGTACGCTCGACAACCTG‑GAATCCATCCCTGGCTTCGACATTTTCCCTGATGAGAACCGTTTCCGCGAAATTAGTAAAGGCGTCGGCGTGGCGATTATCGG‑T‑CAGCCCAGCTCACTGG‑CTCGGGCGGATAAACGTTTTTACTCGG‑CCCGTg                                                                                                                                                                                                                   <  2:410833/250‑1 (MQ=255)

                                    gCGTCGGCGATGTGACTTCGCTGATGTTGGGGTCGATGGTCGCCGCCTGCGGCGGCTAGATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAAGCCATCCCTGGCTTCGCCATTGTCCCGGATGACAACCGTTTCGGCGAAATTATTAAAGTCGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGa‑t                                                                                                                                                                                                                <  1:102959/251‑1 (MQ=255)

                                     cGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑ccc                                                                                                                                                                                                                      <  1:16951/245‑1 (MQ=255)

                                      gTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACCAACTG‑GAATCCATCCCTGGCTTCGACATTTTTCCGGATGACAACCGTTTCCGCGAAATTTTTAAAGACGTCGGGGTGGCGATTATCGG‑T‑CCGACCAGCTCACTGG‑CTACGGCGGATAAACGGTTCTTCGCGG‑CCCggggtat                                                                                                                                                                                                               >  1:288583/1‑245 (MQ=255)

                                           cGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAAAGTTTCTACGGGA‑CCCGTGA‑TA‑TTAccc                                                                                                                                                                                                        >  1:216157/1‑250 (MQ=255)

                                           cGATGTGACTTCGCTGATGTTGGGGCCGATGGGCGCAGCCTGCGGCGGCTATATTCCGAGGATCGCTTGTCGAGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTGTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGC‑TCCGTGA‑TA‑TTACCg                                                                                                                                                                                                        <  2:272316/251‑1 (MQ=255)

                                             aTGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCCGGCTTCCACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTCTCGG‑T‑CAGACCCGCTCACTGG‑ATCCGGCGGGTAAAAGTTTTTACGCGAACCCG‑GA‑TA‑tttacccc                                                                                                                                                                                                      >  2:6125/1‑245 (MQ=255)

                                                          gATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGGCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATACCTGGCTTCGCCCTTTTCCCGGATGAAAACCGTTTCCGCGAAATTATTTAAGACGTCGGCGTGGGGATTATCGG‑G‑CCGACCCGCTCACTGG‑CCCCCGCGGGTAAACGTTTTTACACGG‑CCCGTggcatttccgcacacgtgggctccc                                                                                                                                                                                           >  2:462830/1‑227 (MQ=255)

                                                             gTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATc                                                                                                                                                                                       <  1:264980/250‑1 (MQ=255)

                                                               tGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGCCAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCCTTTCCGCGAAATTATTAAAGACGTCCGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CACGTGA‑TA‑TTTCCGCACCCGTGGACTCCATcccc                                                                                                                                                                                    >  1:48214/1‑250 (MQ=255)

                                                                ggggCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATccct                                                                                                                                                                                    <  2:544011/250‑2 (MQ=255)

                                                                ggggACGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACGGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCg                                                                                                                                                                                    <  1:355238/250‑1 (MQ=255)

                                                                           tCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACCAACTG‑GAATCCCTCCCTGGCTTCGACCTTTTCCCGGGTGGCCACCGTTTTCGCGAAATTATTAAAGACGGCCGCGGGGCGATTCTCGG‑T‑CAGAACCGCTCACTGG‑CTCCGGCGGAGAAACGTTTTTACGCGC‑CCCG‑GACTT‑TTACCCCCACCGCGGGCTCCATCCCG‑ATGaaccacgcc                                                                                                                                                                       >  2:368005/1‑243 (MQ=255)

                                                                             gCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTAGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTCCCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCctc                                                                                                                                                                     <  1:227531/251‑1 (MQ=255)

                                                                             gCAGCCTGCGGCGGCTATATTCCGATGATATCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGGCAAACTG‑GAATCCATCCCTCGCTTCGACATTTTCCCGGAGGACAACCGTTTCCGCGAAAATATTTAACACGTTGGCGTGGCGGTTATCGG‑T‑CAGACCAGGTCACTGG‑CCACGGCGGATAAACGTTTCTACGCGG‑CCCGGGG‑TAATTCCCG‑AACCGCGGAATCCATCCCG‑ccggtccccccctc                                                                                                                                                                     >  2:326361/1‑238 (MQ=255)

                                                                              cAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAAACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACCCCCTCCCC‑CTGATCAccccccct                                                                                                                                                                    >  1:339791/1‑245 (MQ=255)

                                                                                gCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTTCGCGAAAATATTAAAGACGGCGGCGTGGGGATTATCGG‑T‑CCGACCAGCTCACTGG‑CTCCGGCGGGTAAACGTTTCCACGCGG‑CCCGGGA‑TA‑TTTCCGCCACCGTGGACTCCCTCCCg‑atgatttcaccctcttt                                                                                                                                                                  >  1:3506/1‑239 (MQ=255)

                                                                                     cggcggtattatTCCGATGACTTCTGGTCGGGGCCTCGGTCATACGGGCGGTCGGCCCGA‑AAACCGTGAACCCACCCCGGGCTTCGACTTTTTCCCGGATGACAACCTTTTCCGCGAAATTATTAAAGCCGCCGGCTTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCGGGCGGATAAACGTTTCTACGCGC‑CCCGTGA‑TA‑TTCCCGCACCCGTGGACGCCATCCCG‑CTGATCACCGCCTCTATTCTgg                                                                                                                                                             <  2:71654/242‑1 (MQ=255)

                                                                                      ggcggcTCTATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACACACTG‑GAAACCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTTTACGCGC‑CCCGGGA‑TA‑TTACCGCAACCCTGGACTCCAACCCG‑CTGAACACCGCCTCTTTTCTGGc                                                                                                                                                            >  2:419566/1‑251 (MQ=255)

                                                                                         ggcTATATCCCGATGTTCTCTGGTCGCGGCCCCGGTCATCCGGGCGGTACGCTCGACAAACTG‑GAACCCATCCCTGGCTTCGACATTTTCCCGGAGGACAACCGTTTCCGCGAAATTATAAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGCCCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACCCCACCCCG‑CTGATCACCGCCTCTATTCTGGCgaa                                                                                                                                                         <  1:8066/251‑1 (MQ=255)

                                                                                               atTCCGATGATCTCTGGTCGCGGCCTCGGTCATCCGGGCGGTACGCTCGACAATCTG‑GAACCCATCCCTGGCTGCGAAATTTTCCCGGATGACACCCGTTTCCGCAAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACGGG‑CTCCGGCGGATAAACGTTTCTACGCGC‑CCCGTGA‑TA‑TTACCGAACCCGTGGACTCCACCCCG‑CGGATCACCGCCTCTATTCTGGCGAAGAAAc                                                                                                                                                    <  1:160915/250‑1 (MQ=255)

                                                                                                ttCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGCCAAACTG‑GAAACCATCCCTGGCTTCGACCTTTTCCCGGATGACAACCATTTACGCGAAATTATTAAAGAACTCGGCGTGTAGCATATCGG‑T‑CAGACCAGATCACAGG‑ACCCCGCGGATACACGTTTCTACGCGA‑ACCGGGA‑TC‑TTACAGCAACCGTGGGATCCCTCCCG‑CTGATCAACGCCTCTATTCTGGAGAAGAAACtg                                                                                                                                                  >  1:368932/1‑250 (MQ=255)

                                                                                                       gacttcgggTCGCGGCCTCGGTCATACTGGCTGTACGCCCAACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATCACAACCTTTCCCGCGAAATTATTAAAGCCTTCGGCGTGGCTATTACCGG‑T‑CAGCCCAGCTCACTGG‑CTCCGGCGGAAAACCGTTTCTAGGCGA‑CCCGTGA‑TA‑TTACCGCACCCGGGGACTCCACCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAg                                                                                                                                           <  1:233998/244‑1 (MQ=255)

                                                                                                        atCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACCAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAAACGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑AAGACCAACTCACTGG‑CTCCGGCGGAAAAACGGTTCTACGCGA‑CACGCGA‑TA‑TTAACGCAAACGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTACGGCGCAGAAACTGGCGGCAgg                                                                                                                                          >  1:123000/1‑251 (MQ=255)

                                                                                                               gccgcgGCCTCGGTCATACGGGCGGTCCGCTCGAAAAACGG‑GAATCCACCCCTGGCTTCGACATTTCCCCGGATGACACCCGTTTCCGCGAAATTATAAAAGCCGTCGGCTTGGCGATTACCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGAATAAACTTTTCTACGCGA‑CCCGTGA‑TA‑TTCCCGCACCCGTGGACTCCACCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGAc                                                                                                                                   <  1:286996/249‑1 (MQ=255)

                                                                                                                  gcgGCCTCGGTCATACTGGCGGTACGCTCCCCAAACTG‑GAGTCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑C‑CAGACCAGCTCACCGG‑CCCCGGCGGGTAAACGTTTCTACGCGA‑CCCGCGC‑TC‑TTACCGCACCCGTGGACTCCATCCCG‑ACGATCCCCGCCTCTATTCTGGCGGAGAAACtggacggggggttggagacg                                                                                                                                >  2:589056/1‑232 (MQ=255)

                                                                                                                              atccTGGCGGTACGCTCAAAAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTCCCGCGAAATTATTAAAGACGGCGGGGTGGCGATTACCGG‑T‑CAGACCAGCTCACGGG‑CTCCGGCGGATAAACGTTTCTCCGCGC‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGAc                                                                                                                    <  1:214960/248‑1 (MQ=255)

                                                                                                                                      ggTACGCCCGACAAACTG‑GAAGCCATCCCTGGCTTCGATATTTCCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGGGGGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAAGGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGt                                                                                                            <  1:434005/251‑1 (MQ=255)

                                                                                                                                        aaCGCTGGACAAACTG‑GAATCCATCCGGGGCTTCGAAATTTTCCCGGAGTACAACCGTTTCCGCGAAATGATTAAAGCCGTGGGCGTGGCGATAATCGG‑G‑CAGACCAGCTCACGGG‑CTCCGGCGGATAAACGTTTCGACGCGC‑CCCGGTA‑TA‑TTCCGGCAACCGGGGACTCCTCCCCG‑GTGATCCCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTgg                                                                                                          <  1:589056/250‑1 (MQ=255)

                                                                                                                                          cGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGCCTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGCCGCGCTGGTGATGGACGTGAAAGTgggg                                                                                                        >  2:190327/1‑250 (MQ=255)

                                                                                                                                           gCTCGACACCCTG‑GACTCCATCCCTGGATTCGCCATTTTCCCGGATGACAACCGTTTCCGCGAACCTATTAACGACGACGCCGTGGCGATTATCGG‑T‑CAGACCAGTTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCTAGCACTCCATCCCG‑CTGCTCCCCGCCTCTATTCTGGCTAAGAAACTGGCGGAAGGTCTGGACGCGCAGGTGATGGACGGGAAAGAGGGtc                                                                                                       >  2:359188/1‑250 (MQ=255)

                                                                                                                                              cGACAACCTG‑GAACCCACCCCTGGCTCCGTCTTTTCCCCGGATGACACCCGTTTCCGCAAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAACCGTTTCTCCGCGC‑CCCGTGA‑TA‑TTCCCCCACCCGTGGACTCCTCCCCG‑CTGATCTCCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAg                                                                                                      <  1:62263/249‑1 (MQ=255)

                                                                                                                                              cGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgcg                                                                                                    >  1:135312/1‑251 (MQ=255)

                                                                                                                                                aCAAAATG‑GATTCCACCCCTGGCTTCGCTTTTTTCCCGAATGACAACCGTTCCCGCAAAATAATAAAAGACGTCGGCGTGGGGATTATCGG‑T‑CAGACCAGCCCTCGGG‑CCCCGGCGGATAAACGTTTCTACGCGC‑CCGGTGA‑TA‑TTACCGCAACCGTGGCTTCCATCCCG‑CGGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgcggc                                                                                                  <  1:279164/251‑1 (MQ=255)

                                                                                                                                                    aCTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGAATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGt                                                                                              >  1:19679/1‑251 (MQ=255)

                                                                                                                                                        gaaaTCCATCCCTGCCTTCTAAATTTCCCCGGATGACAACCGTTCCCCCAAAATTATTAAACCCGCCGCCGGGGCGATT‑TCGGGT‑CAGCCCACCTC‑CTGGGCTCCGGCGGATAACCGTTTCTACGGGT‑CCCGGGA‑TT‑TTACCTCACCGGGGGCCTCCAACCCC‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTa                                                                                           <  1:204550/249‑1 (MQ=255)

                                                                                                                                                           accccactctgtttccGACATTTCCCCGGATGACACCCTTTGCCGCAAAATTATTAAAGCCGGCGGCGTGGAGTTTATGGG‑T‑CAGACCAGCCCCCGGG‑CCCCGGCGGATAACGGTTTCTCCGCGC‑CCTGTGA‑TA‑TTCCCCCAACCGTGGACTCCAGCCCG‑TGGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGc                                                                                        <  2:556260/236‑1 (MQ=255)

                                                                                                                                                            gccatccCTGTCTGCGACATTTTCCCGTATGACATCCGTTTCCGCGAAATTATTAAAGACGCCGGCTTGGCGATTAGCGT‑T‑CAGCCCAGCTCACTGT‑CCCGGGCGGATAAACTTTTCTACGCGA‑CCCGGGA‑TT‑TTTCCGCAACCGGGGACCCCA‑CCCGTGGGACCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGcc                                                                                       <  1:61376/250‑1 (MQ=255)

                                                                                                                                                              catccTTGGCTACGACATTTTCCCGGATCCCAACCGTTTCCGCGCAATTATTAAAGCCGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCgc                                                                                     >  2:618924/1‑250 (MQ=255)

                                                                                                                                                                  ccTGGCTTCGACATTTTGCCGGATGACAAACGTTTCCGCGAAATTATTCAAGACGTCGGAGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTAAGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACCCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGCGGAAGAGGGTAGCGGCGCGGTTATGCCGCCCt                                                                                  >  2:307707/1‑250 (MQ=255)

                                                                                                                                                                    tGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATAACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACg                                                                               >  1:513870/1‑251 (MQ=255)

                                                                                                                                                                     ggCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGAGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGAAAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTAc                                                                                <  1:179503/249‑1 (MQ=255)

                                                                                                                                                                          cGACTTTTTCCCGGATAAAAACCTTTGCCGCGAAATTATAAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACGGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TT‑TTCCCGCACCCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAActct                                                                         <  1:340089/251‑1 (MQ=255)

                                                                                                                                                                          cGACATTTTCCAGGCTGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTCTCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCCACCGTGGACTCCTTCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGCCTGGACGCGCTGGTGATGGACGTGGAAGTGGGTAGCGGCGCCTTTATGCCGCCCTACGAAccct                                                                         >  2:557999/1‑248 (MQ=255)

                                                                                                                                                                                    acGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAAACGTGGAATCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGGCTGGACGCGCTGGTGATGGACGGGAAAAGGGGTAGCGGCGCGGTTCGTCCGGCCTAGGGACTCTCTGAGGCCCt                                                               >  2:162722/2‑251 (MQ=255)

                                                                                                                                                                                        aTGACAACCGTTTCAGCGAACTTATTAAAGACGTCGACGTGGCGCTTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAACCGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAAAGTGGGTAGGGGCGCGTTTCAGCCGACCTACGAACTCTCTGAAGCCCTTGcc                                                           >  2:244263/1‑251 (MQ=255)

                                                                                                                                                                                         tGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑C‑AAGACCAGCTCACTGG‑CTCCGGCCGAAACACGTTTCTACGCGA‑CCCGGGA‑TA‑TTACCGCAAACGTGGACTCCATCCCG‑CCGATCCCCGCCTCAATTCTGGCGAAGAAAATGGCGGCAGGACTGGACGCGATGGTGATGGACGTGACAGAGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCCTGCCg                                                          >  1:372847/1‑251 (MQ=255)

                                                                                                                                                                                             aaCCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGACACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGTGAACGTGGGTAGCGGCGCGGTTATGCCGACCTACGAACTCTCTGCAGCCCTTGCCGAAGc                                                      >  2:459372/1‑251 (MQ=255)

                                                                                                                                                                                             aaCCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTATGGCCGCGCTGGTGATGGCCGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGc                                                      >  1:627215/1‑251 (MQ=255)

                                                                                                                                                                                                    tCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGACGGTCTGGGAGCGCTGGTGCTGGACGGGGAAGTGGGTAGGGGCGCGTTTATGCCGGCCTACGAACTCTCTTAAGCCCTTGCCGAAGCGgttgtt                                               >  2:332425/1‑251 (MQ=255)

                                                                                                                                                                                                         aaaaTTATAAAAGCCGGCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGC‑CTCCGGCGGAAAAACGTTCCTACGCGC‑CCCGTGA‑TA‑TTCCCGCACCCGTGGACCCCACCCCG‑CGGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGt                                          <  1:360611/250‑1 (MQ=255)

                                                                                                                                                                                                          aaaTTATTAAAGGCGCCGGCGGGGCGATTACCGG‑C‑CAGCCCGGCTCACTGG‑CTCCGGCGGATAAACGTTTCTCCGCGT‑CCCGTGA‑TA‑TTCCCGCAACCGTGGCCTCCATCCGG‑CTGATCACCGCCTCTAGTCGGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTg                                         <  1:414425/251‑1 (MQ=255)

                                                                                                                                                                                                              tattaAAACCGTCGGCGTGGCGATTGTCGG‑T‑AAGCCCGGCCCGCTGG‑CCCCGGCGGATAAACTTTTCTCCTCGA‑CCCGTTA‑TA‑TCCCCGCACCCGTGGACCCCCCCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCt                                      <  1:307707/250‑1 (MQ=255)

                                                                                                                                                                                                                ttaAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGGTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAAc                                   >  1:430499/1‑251 (MQ=255)

                                                                                                                                                                                                                   aaGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCCGGGCTAACGgc                                >  1:434436/1‑251 (MQ=255)

                                                                                                                                                                                                                   aaGACGTCGGCGTGAAGCATATCGG‑T‑CAGACCAGCTAACTGG‑CTCCGGGGGATAAACGTTTCTACGCGA‑CCCGGGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGGGATGGACGGGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCGGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACggg                                >  2:391794/1‑250 (MQ=255)

                                                                                                                                                                                                                         tCGGCGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGCCCTACGGACTCTCTGAAGCCCTTGCCGGAGCGATTGTTGGCGTGGCCAAAGGCGCTGGc                          >  2:174130/1‑251 (MQ=255)

                                                                                                                                                                                                                          cGGCGTGGCGATTATCGG‑T‑CCGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGGACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCg                         >  2:109441/1‑251 (MQ=255)

                                                                                                                                                                                                                             cGTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTACCGGCGCTGGCGTgc                      >  2:593540/1‑251 (MQ=255)

                                                                                                                                                                                                                             cGTGGCGATTATCGG‑T‑CAGACCACCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGCAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGCAATCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGAGCTGGCGTgg                      >  1:153408/1‑250 (MQ=255)

                                                                                                                                                                                                                              gTGGCGATTATCGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGCAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTgcg                     >  1:492749/1‑251 (MQ=255)

                                                                                                                                                                                                                              gTGGCGATTATCGG‑T‑AAGACCAGCTAACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCAG‑CTGATCACCGCATCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTgcg                     >  2:368600/1‑251 (MQ=255)

                                                                                                                                                                                                                                         cGG‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAAAGGCGCTGGCGTGCGCACCACCgcgc          >  2:509016/1‑251 (MQ=255)

                                                                                                                                                                                                                                           g‑T‑CAGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGATGGCGTGCGCACCACCGCgctg        >  2:347686/1‑251 (MQ=255)

                                                                                                                                                                                                                                                aGACCAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGA‑CCCGGGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGACGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGGTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCCGCTc     >  2:356066/1‑251 (MQ=255)

                                                                                                                                                                                                                                                   ccAGCTCACTGG‑CTCCGGCGGATAAACGTTTCTACGCGC‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTGGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGGGGGTAGCGGCGCGTTTATGCCGCCCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCAcc  >  2:548592/1‑251 (MQ=255)

TGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTG‑GAATCCATCCCTGGCTTCGACATTTTCCCGGATGACAACCGTTTCCGCGAAATTATTAAAGACGTCGGCGTGGCGATTATCGG‑T‑CAGACCAGTTCACTGG‑CTCCGGCTGATAAACGTTTCTACGCGA‑CCCGTGA‑TA‑TTACCGCAACCGTGGACTCCATCCCG‑CTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCCTTGCCGAAGCGATTGTTGGCGTGGCTAACGGCGCTGGCGTGCGCACCACCGCGCTGCTCACC  >  NC_000913/4618458‑4618948

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