Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,638	+C	coding (410/1323 nt)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,638	1	.	C	97.8%	272.9 / ‑2.8	92	coding (410/1323 nt)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (31/59);  total (32/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCC‑GGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCC  >  NC_000913/4618391‑4618886                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                              gAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑cccc                                                                                                                                                                                                                                                               <  1:680889/249‑2 (MQ=255)      ggtCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑c                                                                                                                                                                                                                                                       <  2:221579/250‑1 (MQ=255)           cGCTGTCCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AAcc                                                                                                                                                                                                                                                  <  2:455037/250‑1 (MQ=255)           cGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑aa                                                                                                                                                                                                                                                    <  1:566130/248‑1 (MQ=255)           cGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AAc                                                                                                                                                                                                                                                   <  2:240443/249‑1 (MQ=255)              tGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATg                                                                                                                                                                                                                                                          <  1:377494/241‑1 (MQ=255)              tGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCCCCACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATgggatcggaa                                                                                                                                                                                                                                                 >  2:377494/1‑241 (MQ=255)               gACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCAGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTc                                                                                                                                                                                                                                             <  2:155313/251‑1 (MQ=255)                aCCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTcc                                                                                                                                                                                                                                            <  1:141951/251‑1 (MQ=255)                   aTGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCcgcg                                                                                                                                                                                                                                         <  2:139059/251‑1 (MQ=255)                             gagaTTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTaa                                                                                                                                                                                                                              <  1:210282/251‑1 (MQ=255)                              agaTTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGTTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTaaa                                                                                                                                                                                                                             <  1:385044/251‑1 (MQ=255)                               gaTTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAg                                                                                                                                                                                                                            <  2:201447/251‑1 (MQ=255)                               gaTTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCCCGACCAACTTGAAACCCTCCCTGG‑CTTCGACATTTT‑CCCCGGATTA‑C‑AACCGTTTCcgcgcactttttata                                                                                                                                                                                                                              >  2:460335/1‑239 (MQ=255)                                 ttCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATTG‑C‑AACCGTTTCCGCGA‑AATTATTTAAG‑Ac                                                                                                                                                                                                                         >  1:296297/1‑251 (MQ=255)                                  tCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑Ac                                                                                                                                                                                                                         <  1:238978/250‑1 (MQ=255)                                  tCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑Ac                                                                                                                                                                                                                         <  2:526875/250‑1 (MQ=255)                                   cAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATGGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGt                                                                                                                                                                                                                       <  2:236111/251‑1 (MQ=255)                                    aGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTc                                                                                                                                                                                                                      >  1:237821/1‑251 (MQ=255)                                     ggAACCGTTCTCGACTGGAAAAGCCAGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGCTCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGC‑AATTATTAAAG‑ACggcg                                                                                                                                                                                                                     >  1:244637/1‑248 (MQ=255)                                           gTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTgg                                                                                                                                                                                                               >  1:106983/1‑251 (MQ=255)                                                    tGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑aa                                                                                                                                                                                                                                                    <  1:608669/207‑1 (MQ=255)                                                     ggAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATc                                                                                                                                                                                                       <  2:472324/249‑1 (MQ=255)                                                     ggAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATc                                                                                                                                                                                                       <  1:132809/249‑1 (MQ=255)                                                     ggAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCg                                                                                                                                                                                                      <  2:116119/250‑1 (MQ=255)                                                       aaaaGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑Aattat                                                                                                                                                                                                                                 >  1:536824/1‑222 (MQ=255)                                                       aaaaGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑Aattat                                                                                                                                                                                                                                 <  2:536824/222‑1 (MQ=255)                                                        aaaGCCTGCATCTGAATGGACCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGGACCCTCGGCAAAAAGGAAACCATCCCTGG‑CTGCGACATTTTACCC‑GGGTGA‑A‑AAACGTTTCC‑CGCCAATTATTTAAG‑AAGTCGGCGTGGCGAttttccggt                                                                                                                                                                                                   >  2:408503/1‑243 (MQ=255)                                                            ccTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGaac                                                                                                                                                                                              >  1:152199/1‑249 (MQ=255)                                                                cATCTGAATGGCCCGATTGTTGATAAACCCTCCACCGGTGGCGGCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGtt                                                                                                                                                                                          <  2:56044/251‑1 (MQ=255)                                                                     gAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCAc                                                                                                                                                                                       <  1:570066/249‑1 (MQ=255)                                                                     gAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTg                                                                                                                                                                                     <  1:460335/251‑1 (MQ=255)                                                                         ggCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCt                                                                                                                                                                                  <  2:25474/250‑1 (MQ=255)                                                                          gCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTcc                                                                                                                                                                                <  2:186260/251‑1 (MQ=255)                                                                           cccGATTGTTGATAAACCCTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCg                                                                                                                                                                               <  1:314381/251‑1 (MQ=255)                                                                                        aaaCACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGGTTATCGGGCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGtt                                                                                                                                                                 >  1:684588/1‑251 (MQ=255)                                                                                        aaaCACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCGGTTCACTGGCTCCGG‑CCGATAAACGtt                                                                                                                                                                 >  1:695349/1‑251 (MQ=255)                                                                                         aaCACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AAttatta                                                                                                                                                                                                                               >  1:603137/1‑190 (MQ=255)                                                                                         aaCACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AAttatta                                                                                                                                                                                                                               <  2:603137/190‑1 (MQ=255)                                                                                          acacTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGCCAAACTCGAATCCACCCCTGG‑CTTCGCCATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑AAGTCGGCGCGGGGATTATCGGTCAAACCCGGTCACTGGCTCCCG‑CTGGTaaaaggttt                                                                                                                                                               >  2:431130/1‑245 (MQ=255)                                                                                          acacTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTc                                                                                                                                                               <  1:431130/251‑1 (MQ=255)                                                                                               ccACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTAc                                                                                                                                                            <  2:378745/249‑1 (MQ=255)                                                                                               ccACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAATGGAAACCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGGTTTCGCGA‑AATTATTAAAG‑AAGTCGGCGGGGGGGTTATCGGGCAGACCAGGTCCATGGCTCCGG‑CTGATAAACGGTTCTaagc                                                                                                                                                          >  2:482368/1‑248 (MQ=255)                                                                                                cACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCAGG‑Ct                                                                                                                                                                           <  1:610589/233‑1 (MQ=255)                                                                                                cACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCCAAATTCCGATGATCTCTGGTCGCGGCCCCGGTCATACTGGCGGTACGCTCGACAAACTTGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGCAGG‑C‑CACCGTTTCCGCGA‑AATTCTTAAAG‑ACGTCCGCGTGGCGGTTAACCGTCCGACCACGTCACTTGCTcccgggtgacccgaagaaccgagg                                                                                                                                                          >  2:610589/1‑229 (MQ=255)                                                                                                  ccGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGAc                                                                                                                                                       <  1:320395/251‑1 (MQ=255)                                                                                                    ggTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGGCAGACCCGTTCACTGGCTCCGG‑GTGATAAACGGTTCTACGCGAccc                                                                                                                                                     >  1:73286/1‑251 (MQ=255)                                                                                                          gTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGat                                                                                                                                                >  1:9117/1‑250 (MQ=255)                                                                                                           gcGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGatat                                                                                                                                              <  2:340091/250‑1 (MQ=255)                                                                                                            cGGCGATGTGACTTCGCTGATGATGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACCGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTTTTTAAA‑ACGTCGGCGTGGGGATTATCGGTCAGACCAGTTTCCCGGGTCCGG‑CTGCTAAACGTTTCTACGCGGCCCgcggtttt                                                                                                                                             >  1:471162/1‑244 (MQ=255)                                                                                                             ggCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCCGACCCGTTCAATGGCTCCGG‑CTGATAAACGGTTCTACGCGACCCGTGATAtt                                                                                                                                             >  1:12878/1‑250 (MQ=255)                                                                                                             ggCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGAAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGGACGCTCGACACACTGGAATCCCTCCCTGG‑CTTCTACATTTT‑CCCCGGATGAAC‑AAC‑GTTTCCGCGA‑AATTATTTACA‑ACGTCGGCGGAGCGATTATAGGTCCGGCCAGGTCACCAGTTCCGG‑CTGATAAAAGTTTCTACGCGACCCggggtattt                                                                                                                                            >  1:574680/1‑243 (MQ=255)                                                                                                                gATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCCCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGGTGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCCGACCAGTTCACTGGCTCCGG‑CTTGTAAACGTTTTTTCGCGACCCGTGGTATTaacc                                                                                                                                         >  2:171251/1‑248 (MQ=255)                                                                                                                   gtgACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCa                                                                                                                                       <  2:106983/250‑1 (MQ=255)                                                                                                                   gtgACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCa                                                                                                                                       <  2:237821/250‑1 (MQ=255)                                                                                                                          cGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGa                                                                                                                               <  2:612335/251‑1 (MQ=255)                                                                                                                             ccgtagcttacggagtcttttcAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTc                                                                                                                            <  2:329675/230‑1 (MQ=255)                                                                                                                              gaggTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTcc                                                                                                                           <  2:453959/248‑1 (MQ=255)                                                                                                                              gATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGCCCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTc                                                                                                                            <  2:296297/250‑1 (MQ=255)                                                                                                                              gATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTc                                                                                                                            <  1:352834/250‑1 (MQ=255)                                                                                                                                   gggggCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATccc                                                                                                                      <  2:205826/250‑1 (MQ=255)                                                                                                                                           aTGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCa                                                                                                              <  2:152199/251‑1 (MQ=255)                                                                                                                                              gTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGGTAAAAGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGGCTCCATCCCCCTGATCAccc                                                                                                           >  1:55682/1‑250 (MQ=255)                                                                                                                                                 gCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGGTTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGGGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCct                                                                                                        >  1:1892/1‑251 (MQ=255)                                                                                                                                                 gCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGTTGATCACCGCct                                                                                                        <  1:482368/251‑1 (MQ=255)                                                                                                                                                        gcggcggttttttgcggatgatCTCGGGGCGCGGCCTCGGTACTCCTGGCGGTACGCTCGACAAACTGGAATCCACCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTCCCGCAACGGTGGACTCCATCCCGCTGATCACCGCCTCTATTCt                                                                                                 <  2:55682/235‑1 (MQ=255)                                                                                                                                                                natatTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAAATGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGGTCACTGGCTCCGG‑CTGATAAAAGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCga                                                                                            >  2:324952/2‑248 (MQ=255)                                                                                                                                                                 atagtCCGATGATCTCTGGCCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCTTCCCTGG‑CTCCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCTTGGCGATTAGCGGTCAGTCCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTCCGCGACCGGTGATATTCCCGCAACCGGGGACCCCACCCCGCTGATCACCGCCTCTATTCTGGCgaagaa                                                                                        <  1:287540/247‑1 (MQ=255)                                                                                                                                                                     cccGATGATCTCGGGCCGCGGCCTCGGTCATCCTGGCGTTACGCTCCACAAAGTGGAATCCATCCCTGG‑TTTCGCCATTTC‑CCCCGGATGA‑C‑ACCCGTTTCCGCGA‑ATTTATTAAAG‑CCGCCGGCGTGGGGATTATCGGTCAGACAAGTTCACTGGCTCCGG‑CTGATAAACTTTTCTACGCGACCCGTGATATTACCGCAACCGTGCACACCTTCCAGCTGATCACCGCCTCTATTCTGGCGAAGAAACtt                                                                                    <  2:432912/250‑1 (MQ=255)                                                                                                                                                                        gatgatCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCg                                                                                                                                                                                                             >  1:465234/1‑128 (MQ=255)                                                                                                                                                                        gatgatCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCg                                                                                                                                                                                                             <  2:465234/128‑1 (MQ=255)                                                                                                                                                                        gatgatCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTTGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCg                                                                                 >  1:482098/1‑251 (MQ=255)                                                                                                                                                                             cccctGGTCGCGCCCTCGGTAATACTGGCGTTACGCTCGACAAACTGGAACCCATCCCTGG‑CTTCCACATTTT‑CCCCGGATGA‑C‑ACCCGTTTCCCCAA‑AATTATTAAAG‑ACGTCGGCGTGGCAATTTTCGGTCAGACCATTCCACTGGCCCCGG‑CTGATAAACTTTTCTACGCGCCCGGTGATATTCCCGCAACCGTGGACCCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAgg                                                                            <  1:222150/248‑1 (MQ=255)                                                                                                                                                                                cgggtcggggccgcgttattccgggCGGTCCGCTCAACAAACGGGAATCC‑CCCCTGGGCTTCGACATTTC‑CCCCGGAGGC‑C‑AACCGTTCCCGCGA‑ATTTTTTAAAG‑ACGCGGGCGTGGGGATTATCGGTCAGACCGGTTCACTGGCCCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACGGTGGACTCCACCCCGCTGATCGCCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTc                                                                          <  2:695349/227‑1 (MQ=255)                                                                                                                                                                                cgggtcgggcccTCGGTCATCCTGGCGGTCCGCTCGACAACCTGGAACCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑ATTTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCCCCGG‑CTGATAAACGTTTCTACGCGCCCCGTGATATTCCCGAACCCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGt                                                                           <  1:379449/239‑1 (MQ=255)                                                                                                                                                                                   gTCGGGGCCCCGTTCATACTGGCGGTCCGACCGAAAAACTGGAATCCACCCCGGT‑CTTCGGCATTTT‑CCCCGGATGA‑C‑ACCCGTTTCCGCGA‑ATTTTTTAAAG‑CCGTCGGCGGGGCGATTATCGGTCATCCCTGTTCACGGGCCCGGG‑CTGATAAACGTTTCTACGCGACCCGTTATATAACCGCAACCGGGGACTCCATCCTGATGCTCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGa                                                                      <  1:295591/251‑1 (MQ=255)                                                                                                                                                                                         gCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCt                                                                >  2:561048/1‑251 (MQ=255)                                                                                                                                                                                           cTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTgg                                                              >  2:139382/1‑251 (MQ=255)                                                                                                                                                                                             gggtcatagtgcggttgcccTCGACAAAGGGAAATCCCTCCCTGG‑TTTCGCCTTTTT‑CCCCGGTTAA‑C‑ACCCGTTTCCGCGA‑AATTATAAAAG‑ACGTCGGGGTGGCGATTATGGGTCAGCCCAGTTCCCGGGCTCCGG‑CTGATAAACTTTTCTACGCGACCCGTGATATTCCCGCACCCGGGGACTCCACCCCGCGGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTg                                                            <  2:684588/232‑1 (MQ=255)                                                                                                                                                                                                 aaTACTGGGGGTCCGCCCGACAAACGGGAATCCATCCCTGA‑CTTCGAAATTTT‑CCCCGGATGA‑C‑CACCTTTTCCAGAA‑AATTATTAAAT‑ACGCCGGCGTGGCGATTTTCGGACAGACCAGTTCGCTGGCTGCGG‑CTGATAAACGTTTCTCCGCGACCCGTGATATTCCCCCACCCGTGGCCTCCACCCCGCTTAACACCGCCGCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATgg                                                        <  2:608194/250‑1 (MQ=255)                                                                                                                                                                                                   tactgtgggtccgcgcgacaacgtGGATTCCATCCCTGG‑TTTCGTAATTTT‑CCCCGGATGA‑C‑AACCGTTTCGGCGA‑AATTATTAAAG‑CGGTTGGCGAGGCAATTTTCGGTCAGACCAGTTCACTGG‑TCCGGGCTGATAAACGTTTCTACGCGACCCGTATTATTCCCGCACCCGTGGGCCCCTTCCCGCTGATCGCCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGa                                                       <  1:621694/227‑1 (MQ=255)                                                                                                                                                                                                         cGGTACGCTCGACAAACTGGAATCCACCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTCCCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGaaa                                                <  1:314083/251‑1 (MQ=255)                                                                                                                                                                                                             ccgcacgaaaaaCGGGACTCCACCCCTTG‑CTTCTCCATTTC‑CCC‑CGAAGA‑CCAACCGTTCCCGAGA‑AATTAT‑AAAGTACGTCGGCGTGGCTATTATCGGTCAGCCCAGTTCACTGGCTCCGG‑CTGATAACCTTTTCTACGCGCCCCGTGATATTACCGCAACCGTGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTg                                             <  2:1892/241‑1 (MQ=255)                                                                                                                                                                                                             aCGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTgg                                            >  1:682416/1‑251 (MQ=255)                                                                                                                                                                                                             aCGCCCGACAAACGGGAATCCACCCCTTG‑CTTCGACATTTC‑CCCCGAATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCGGTTCACTGGCTCCGG‑CTGATAAACGTTTCTCCGCGCCCCGTGATATTACCGCACCCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTg                                             <  1:165274/250‑1 (MQ=255)                                                                                                                                                                                                                cTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGGTTCTACGCGACCCGGGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAAAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGCGAAAGGGGGt                                          >  1:31808/1‑250 (MQ=255)                                                                                                                                                                                                                   gACAAACGGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTGGGCGTGGCGATTATCGGTCAGACCGGTTCCCGGGCCCCGG‑CTGATAAACGTTTCTACGCGCCCCGTGATATTCCCGCACCCGTGGACTCCACCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgcg                                      <  2:12878/251‑1 (MQ=255)                                                                                                                                                                                                                     cAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑CCGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCCCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTCCCGCACCCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAgcggc                                    <  1:139382/251‑1 (MQ=255)                                                                                                                                                                                                                             aacccatccCTGG‑CTTCGACATTTT‑CCCCGGATGA‑C‑ACCCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTCCGCGCCCCGTGATATTCCCGCACCCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTAt                            <  1:129527/248‑1 (MQ=255)                                                                                                                                                                                                                             aaTCCATCCCTGG‑CTCCGAAATTTT‑CCCCGGATGA‑A‑AACGGTTCCCGCGA‑AATTATTAAAG‑AGGGCGGGGTGGGGATTAGCGGCAAGACCGTTTCCCTGGCTCCGG‑CTAATAAACTTTTCTACGCGCCCCGGGATATTACCGCACCCGTGGACTCCATCCCGCGGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTAt                            <  1:650484/251‑1 (MQ=255)                                                                                                                                                                                                                                       cggctttCGACATTTT‑CCCCGGATGA‑C‑ACCCGTTTCCGCGA‑AATTATTAAAG‑ACGGCGGCGTGGCGATTATCGGTCAGCCCAGTTCCCTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTa                   <  2:9117/246‑1 (MQ=255)                                                                                                                                                                                                                                            ttCGACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGGCGACCTACGaa               >  2:389546/1‑250 (MQ=255)                                                                                                                                                                                                                                               gACATTTT‑CCCCGGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAActct           >  1:192043/1‑251 (MQ=255)                                                                                                                                                                                                                                                        ccc‑GGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTCGCGGCGCGTTTATGCCGACCCACGAACCCTCTGAAGccc  >  2:543680/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                              GAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCGCTGATGTTGGGGCCGATGGTCGCAGCCTGCGGCGGCTATATTCCGATGATCTCTGGTCGCGGCCTCGGTCATACTGGCGGTACGCTCGACAAACTGGAATCCATCCCTGG‑CTTCGACATTTT‑CCC‑GGATGA‑C‑AACCGTTTCCGCGA‑AATTATTAAAG‑ACGTCGGCGTGGCGATTATCGGTCAGACCAGTTCACTGGCTCCGG‑CTGATAAACGTTTCTACGCGACCCGTGATATTACCGCAACCGTGGACTCCATCCCGCTGATCACCGCCTCTATTCTGGCGAAGAAACTTGCGGAAGGTCTGGACGCGCTGGTGATGGACGTGAAAGTGGGTAGCGGCGCGTTTATGCCGACCTACGAACTCTCTGAAGCCC  >  NC_000913/4618391‑4618886

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