Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,454	A→G	N409S (AAC→AGC)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,454	0	A	G	97.7%	219.0 / NA	86	N409S (AAC→AGC)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (38/45);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAACTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGG  >  NC_000913/4619210‑4619696                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         gagaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑c                                                                                                                                                                                                                                                    <  2:450275/250‑1 (MQ=255) gagaACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTTTCCACGCGCAAGACA‑AAAACAGCTGG‑cg                                                                                                                                                                                                                                                   >  2:522965/1‑250 (MQ=255)    aaCTACGCGAAGTATCTGCCCACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCCGG‑CAgg                                                                                                                                                                                                                                                 >  1:195096/1‑250 (MQ=255)      cTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCCGCGAAAGACG‑AAAACAGCGGG‑CAGGa                                                                                                                                                                                                                                                >  1:49819/1‑249 (MQ=255)      atACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTAAGTGCAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTCGACGGTCAGCGTCCGCTGGCGGGTCTCCCAGGGAAAGGCG‑AAAACAGCTGG‑CAggg                                                                                                                                                                                                                                                >  2:386883/2‑248 (MQ=255)        aCGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAg                                                                                                                                                                                                                                              <  1:222720/249‑1 (MQ=255)        aCGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGAGTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAg‑c                                                                                                                                                                                                                                            <  2:363772/250‑1 (MQ=255)         cgcgAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTATTCACCCGGAAGACG‑AAAACCGCGGG‑GAGGAAg‑cgg                                                                                                                                                                                                                                          >  2:544205/1‑251 (MQ=255)             aaGTATCTGCCGACAGCGATGCCGACGAAAGCAGTCTGTGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACGCCATGGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAAG‑AAAACAGCTGG‑CAGGAAG‑CGGAGaa                                                                                                                                                                                                                                      <  2:173445/251‑1 (MQ=255)                gatcggCCGACAGCGGTGCTGGGGAAAGCAGGCTATGGTGATACCGAAGGTGTGGTCAGTGAAGGGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGGCGTCAGGCATCTGACGCCGTTGATTACAGCGGCGGCTGTACTGATATGGGGCGTTTGGGCGACCAGGGAGGCGGTCAGCGTCCGCTGGCGTTTATGCATGCGAAAGACG‑AAAACAGCTGG‑C‑GGTAGTCGGCGAAAGc                                                                                                                                                                                                                                   <  2:90751/246‑1 (MQ=255)                     gCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTg                                                                                                                                                                                                                               <  1:293573/250‑1 (MQ=255)                     gCCGAAAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGt                                                                                                                                                                                                                                <  2:26591/249‑1 (MQ=255)                       cGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGCAAGAAG‑GAAACAGCTGG‑CAGGAAA‑CGGCGAAAGCGGGGaa                                                                                                                                                                                                                             >  1:183051/1‑250 (MQ=255)                       cGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTTTCCCCGCGGAAGACGGAAAACG‑CTGG‑CCGGAAg‑cggcgcaaacggggaa                                                                                                                                                                                                                             >  2:220253/1‑239 (MQ=255)                           aGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGAACAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCACGCGAAAGACG‑GAAACAGCTGG‑CAGGGAG‑CGGCGAAAGAGGGGAAAACgg                                                                                                                                                                                                                        >  1:382379/1‑251 (MQ=255)                                 gCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAAtt                                                                                                                                                                                                                   <  2:212345/250‑1 (MQ=255)                                  cTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGGGAAAGCGGCAATTa                                                                                                                                                                                                                  >  1:503334/1‑250 (MQ=255)                                     aCGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑Ct                                                                                                                                                                                                             <  1:48800/251‑1 (MQ=255)                                       gAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGTTGAAAGCGGCAATTAAA‑CTTg                                                                                                                                                                                                           <  1:487044/251‑1 (MQ=255)                                        aaaGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGGAAGCGGTGAAAACGGAAATTAAA‑CTTGc                                                                                                                                                                                                          >  1:613023/1‑251 (MQ=255)                                           gCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCCATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGCTTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTCTCCCCGCGGCAGACG‑AAAAACGCTGG‑CAGGAAG‑CGGGGGAAGCGGTGAAAACGGCacatta                                                                                                                                                                                                                 >  2:209318/1‑237 (MQ=255)                                           gCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTaa                                                                                                                                                                                                                 <  1:209318/242‑1 (MQ=255)                                            cAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCCTCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGAGG‑GAAAAAGCTGG‑CAGGAAG‑CGGCGGAAGCGGTGAAAACGGCAATTAAA‑ATTGCCGa                                                                                                                                                                                                       >  1:112672/1‑250 (MQ=255)                                            cAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCACGCGAAAGACG‑AAAACAGCCGG‑CCGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGAt                                                                                                                                                                                                      >  1:3793/1‑251 (MQ=255)                                             aGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATa                                                                                                                                                                                                     <  2:88941/251‑1 (MQ=255)                                              gTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATaa                                                                                                                                                                                                    <  1:265315/251‑1 (MQ=255)                                                  aTGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGTCGTCAGGCATCTGAAACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGGCCGCTGGCGGGTATCCACCCGGAAGACG‑AAAACGGCTGG‑GCGGAAA‑ACGCGAAAGCGGAGGAAGGGGCCATTAAA‑Ctgtccgcacaata                                                                                                                                                                                                 >  2:147101/1‑238 (MQ=255)                                                    gCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCACGCGAAAAACG‑AAAACAGCGGG‑CAGGAAG‑CGGCGGAGACGGGGGAAACGGCCATTAAA‑CTTGCCGATAAAGAAc‑c                                                                                                                                                                                             >  2:538598/1‑251 (MQ=255)                                                        aTACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGCTGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGAAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCGGGAAAGACG‑ACAACAGCTGG‑CCGGCAG‑CGGCGTAATCGGGGCAAGCTGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGa                                                                                                                                                                                          >  2:140167/1‑250 (MQ=255)                                                        aTACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCACGCGCAAGACG‑AAAACAGCTGG‑CAGGGAG‑CGGCGAAAGCGGGGAAAGGGGCCATTAAC‑CTTTCCGATAAAGCAC‑Cggg                                                                                                                                                                                          >  2:117941/1‑249 (MQ=255)                                                         tACCGTAGGATTGGTCAGTGAAATGGATACCCGCGCGCTGAGGCTGGCAGTGGTTGCACTGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCAACCCGGTAGACGGTCAGCGCCCGCTGGCGGTTTTCCCCACGCAAGACG‑AAAACAGCTGG‑CCGGGAG‑CCGCGAAAGAGGGGAAAAGGGCAATTAAA‑CTTGCCGGTTAAGCAA‑CGGaa                                                                                                                                                                                         >  2:368570/1‑250 (MQ=255)                                                           ccGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGc                                                                                                                                                                                      >  1:106884/1‑251 (MQ=255)                                                              aaGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CCGGAAG‑CGGCGAAAGCGGTGAAAACGGCAATTAAA‑CTTGCCGGTAAAGCAC‑CGGAAAGcac                                                                                                                                                                                    >  1:246822/1‑250 (MQ=255)                                                              aaGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGGGAAAGCGGCAATTAAA‑ATTGCCGGTAAAGCAC‑CGGGAAGcac                                                                                                                                                                                    >  1:60900/1‑250 (MQ=255)                                                                   tgtGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCTTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑aa                                                                                                                                                                              <  1:391794/248‑1 (MQ=255)                                                                     tGTCAGTGAAATGGATCCCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCGTCGTCAGGCATCGGACACCATTGATTACAGCGTCGGTGTTACTGATATGGCGCGTTTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGTAAGCCACGCGAAAGACA‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACCCC‑AACt                                                                                                                                                                            <  2:106884/250‑1 (MQ=255)                                                                       tCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCACGCGAAAGGCG‑AAAACAGCTGG‑CAGGGAG‑CGGCGAAAGCGGGGAAAGCGGCCATTAAA‑CTTGCCGATAAAGCAA‑CGGAAAGAACACC‑AACTGTc                                                                                                                                                                         >  1:343345/1‑251 (MQ=255)                                                                             aaaTGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGc                                                                                                                                                                   <  1:410557/251‑1 (MQ=255)                                                                                 ggATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGGGGTCGTCGTCAGGCATCTGACCCCATGGATTACAGCGTCGTCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCggt                                                                                                                                                               <  2:108740/251‑3 (MQ=255)                                                                                   aTACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGGTATCCACGCGAAAGACG‑AAAACAGCCGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTCAA‑CTTTCCGCAAAAACAC‑CGGCAAGCAAACC‑AACTGTCTATCCCCGCATc                                                                                                                                                             >  1:630224/1‑251 (MQ=255)                                                                                      cccGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGc                                                                                                                                                          <  1:42293/251‑1 (MQ=255)                                                                                            cgcTGGGGAAGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTTTCCCCGCGCAAGAAG‑GCAACAGCGGG‑CAGGAAG‑CGGGGGAAGGGGTGGAAGCGGGAATTAAAATTTGCGGACAAA‑CAC‑CCGAAAGCACCCCCAACTGTC‑ATCGCCGGAGCAGCGaaaa                                                                                                                                                     >  2:183349/1‑248 (MQ=255)                                                                                             gcTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAAc                                                                                                                                                   <  2:593931/251‑1 (MQ=255)                                                                                               tGGGGATGGCCGTGGTTGCACTGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGGAG‑CGGCGAAAGCGGTGAAAACGGCAATTAAA‑CTTGCCGGTAACACAC‑CGGAAAGCAAACC‑AACTGTCTATCGCCGTATCAGCGAATAACgg                                                                                                                                                 >  1:193068/1‑251 (MQ=255)                                                                                                ggggATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGt                                                                                                                                                <  1:558002/251‑1 (MQ=255)                                                                                                 ggaaTGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGATATCCACGCGGAAGACG‑AAAACAGCTGGACAG‑AAG‑CGGCGAAAGCGGTGAAAGCGGCAAATAAA‑CCTACCGATTAAG‑ACACAGAAAGCCCACC‑AACTGTCTATCGCCGTAtacgcggaaaacgga                                                                                                                                                >  2:120252/4‑236 (MQ=255)                                                                                                         aGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatc                                                                                                                                       <  2:49819/251‑1 (MQ=255)                                                                                                             gTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGAAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGGGAAAGCGGCAATTAAA‑CTTGCCGGTAAAGCAC‑CGGAAAGCACACC‑AACCGTCTATCGCCGGCACAGCGAATAAAGGTATActtatctgat                                                                                                                                   >  1:359375/1‑251 (MQ=255)                                                                                                                 cAATGGGGGGCGGGCGGCGTCAGGCATCTGACACCAGTGATTACAGCGGCGGCTTTAGTGATATGGAGCGTTTGGGCGACCAGGTAGACGGTCAGCGTGCGCTGGCGGTTATGCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACCAC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAtt                                                                                                                               <  1:107936/251‑1 (MQ=255)                                                                                                                 cAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAtt                                                                                                                               <  1:107911/251‑1 (MQ=255)                                                                                                                     gggCGGCGGTTGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTaaa                                                                                                                           <  2:351969/251‑1 (MQ=255)                                                                                                                     gggCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGAGATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTaa                                                                                                                            <  1:380991/250‑1 (MQ=255)                                                                                                                      ggcggcGGTCGTCGTCAGGCAACTGACACCATTGATTACAGCGTGGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATTCACGCGAAAGACG‑GAAAAAACTGG‑CGGGGAG‑CGGGGAAAGCGGTGAACACGGCAATTTAA‑CCTGGCGCTAAAACAC‑CGGAAAAGACCCC‑AAATTTCTATCGCCGTATCAACGAGTAACGGGATActgatctgcacattttaa                                                                                                                           >  2:171003/1‑240 (MQ=255)                                                                                                                        cggcggGCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCGTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTAGCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAAttt                                                                                                                        <  2:600478/251‑1 (MQ=255)                                                                                                                           gggtcgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGaa                                                                                                                     <  1:581872/244‑1 (MQ=255)                                                                                                                             gtcgtcGTCAGGCATCTGACACCAGTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGc                                                                                                                   <  1:332586/246‑1 (MQ=255)                                                                                                                               cgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGGAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑CACTGTCt                                                                                                                                                                        >  2:224561/4‑196 (MQ=255)                                                                                                                               cgtcGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCt                                                                                                                                                                        <  1:224561/193‑1 (MQ=255)                                                                                                                                         cATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGga                                                                                                      <  1:117941/251‑1 (MQ=255)                                                                                                                                         cATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGga                                                                                                      <  1:220253/251‑1 (MQ=255)                                                                                                                                                cacCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACata                                                                                               <  1:616974/251‑1 (MQ=255)                                                                                                                                                 acCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACatat                                                                                              <  2:193068/251‑1 (MQ=255)                                                                                                                                                  ccATTGATTACAGCGTCGGCTTTACTGATATGGCGCGGTTGGGCGACCAGGTAGACGGTCAGCGTCCGCGGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATg                                                                                             <  2:215275/251‑1 (MQ=255)                                                                                                                                                     ttgATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGGGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGaa                                                                                           >  1:476659/3‑250 (MQ=255)                                                                                                                                                         ttACAGCGTCGGATTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTAGCGGTTATTCAAGGGAGAGAAG‑AACACAACTGG‑CCGGAAG‑CGGCGAAAGCGGTGGAAGCGGCAATTAAA‑ATTGCCGATAAAGCAC‑CGGGAAACACACC‑AACTGTCTATCGCCGGATCAGCGAATAACGGAATACTAATCTGATCATTTAAATTTTAAGCACTGAG‑GACGGAGGACATATGaaaactg                                                                                      >  2:623874/1‑247 (MQ=255)                                                                                                                                                            cAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGGGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGGAAGAACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAACACTGAG‑TACGGAGGACATATGGAACGTggat                                                                                   >  1:427127/1‑248 (MQ=255)                                                                                                                                                               cGTCGGCTTTACGGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCTCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATtta                                                                                <  1:400796/251‑1 (MQ=255)                                                                                                                                                                       ttACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATgg                                                                          <  1:503947/249‑1 (MQ=255)                                                                                                                                                                        tACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCCCGCGAAAGACG‑GAAAAAGCTGG‑CAGGAAG‑CGGCGCAAACGGGGAAAGCGGCCATTAAA‑CTTGCCGATAAAACAC‑CGGAAAGCAAAAC‑AAATGTCTATCGCCCTATCCGCGAATAACGGTATACCGATTTGATCATTTAAATTTGAAGAACT‑AGATACGGAGAACATAATAAACGTGCATTTTTTATtgtgg                                                                       >  2:568047/1‑250 (MQ=255)                                                                                                                                                                               atGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAAtt                                                                                                                                                                                                                   <  1:423100/108‑1 (MQ=255)                                                                                                                                                                               atGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACCTTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTc                                                                >  1:224967/1‑251 (MQ=255)                                                                                                                                                                                 ggcgctgttGGGCGAACAGGTAGACGGGCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAACACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAGTTAAAGTTGAAGCACTGAG‑TACGTAGAACATATGAAACGTGCATTTAGTATGGTGCTGCACTCa                                                               <  2:583067/242‑1 (MQ=255)                                                                                                                                                                                   cgcgTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTc                                                            <  2:503334/251‑1 (MQ=255)                                                                                                                                                                                        ttGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGGCTCATTCGGCa                                                        <  1:197471/249‑1 (MQ=255)                                                                                                                                                                                                 cAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTAc                                              <  1:522965/251‑1 (MQ=255)                                                                                                                                                                                                              aGCGTCCGCTGGCGGTTATCCACCCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCCC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGCG‑TACGGAGAACATATGAAAAGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGAAGa                                  >  1:86007/1‑250 (MQ=255)                                                                                                                                                                                                                   ccGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCtt                            >  1:611680/1‑251 (MQ=255)                                                                                                                                                                                                                   ccGCTGGCGGTTATCCACGCGAAAGACG‑AAAAAAGCTCG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGGATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑GACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATTAAGAACGCtt                            >  2:138513/1‑251 (MQ=255)                                                                                                                                                                                                                       tGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCa                                                                                                                  >  2:503888/1‑162 (MQ=255)                                                                                                                                                                                                                       tGGCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCa                                                                                                                  <  1:503888/162‑1 (MQ=255)                                                                                                                                                                                                                        ggCGGTTATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGt                        <  1:538598/250‑1 (MQ=255)                                                                                                                                                                                                                              tATCCACGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCg                 >  1:204855/1‑251 (MQ=255)                                                                                                                                                                                                                                   aCGCGAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGACTAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGCAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCt            >  1:387385/1‑251 (MQ=255)                                                                                                                                                                                                                                      cgAAAGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGac         >  2:331942/1‑251 (MQ=255)                                                                                                                                                                                                                                          aGACG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACAccc     <  1:604624/251‑1 (MQ=255)                                                                                                                                                                                                                                             cG‑AAAACAGCTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTgg  <  1:39987/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         GAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACG‑AAAACAACTGG‑CAGGAAG‑CGGCGAAAGCGGTGAAAGCGGCAATTAAA‑CTTGCCGATAAAGCAC‑CGGAAAGCACACC‑AACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAG‑TACGGAGAACATATGAAACGTGCATTTATTATGGTGCTGGACTCATTCGGCATCGGCGCTACAGAAGATGCAGAACGCTTTGGTGACGTCGGGGCTGACACCCTGG  >  NC_000913/4619210‑4619696

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: ‑