Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,031	A→G	N250N (AAT→AAC)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,031	0	A	G	97.9%	318.7 / ‑2.1	94	N250N (AAT→AAC)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (47/45);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

TAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCGCTATTATCGC‑GCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCA  >  NC_000913/64789‑65272                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   aacggcgtcAGCGCAGGGTGATCTTCGGAGAACCTGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGAGTTTCCTTCGCCTAATAACTCCTGGGCGACGGTT‑GCCAGCGAGAATGAATAGAAATTCCATAACTCGGATTTCAGCGCCTCGATACCGTCGATAATGAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGGCCTCGCGCC‑ACTCCAGCTCACTGTTATcgc‑g                                                                                                                                                                                                                                           <  2:470557/243‑1 (MQ=255)        ccAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑Gccccagc                                                                                                                                                                                                                                    >  1:611592/1‑247 (MQ=255)            cgcAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTAGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCa                                                                                                                                                                                                                                <  1:887563/251‑1 (MQ=255)             gcAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑Gccccagccgc                                                                                                                                                                                                                                 >  1:462254/1‑242 (MQ=255)             gcAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAg                                                                                                                                                                                                                               <  1:935608/251‑1 (MQ=255)               aGGTTTATCTTCGGCGAAACGGCGGTCAAATTCGCCCATTCGATACCACGGGTTATCGATAGATTTTCCTTCGCCTAATAAATCCTGCGCGACGGTT‑TCCCGCGGGGATGAAAAGAAATTCCAGAACGCAGATTTTAGCGCCTCGATACCGTCGATAACTAGCCGACCTTTAAACTGGGCAAAAAA‑GCCG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑CCTCCAGCTAACTGTTATCGC‑Gccccagccgccgcg                                                                                                                                                                                                                             >  2:587659/1‑240 (MQ=255)                 gTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGa                                                                                                                                                                                                                           <  1:351493/251‑1 (MQ=255)                     aTCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGaata                                                                                                                                                                                                                        <  1:339578/250‑3 (MQ=255)                     aTCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACCCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCgggata                                                                                                                                                                                                                        >  1:120250/1‑246 (MQ=255)                               aaaCGTCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTc                                                                                                                                                                                                             <  2:171737/251‑1 (MQ=255)                               aaaCGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTc                                                                                                                                                                                                             <  1:103569/251‑1 (MQ=255)                               aaaCGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTc                                                                                                                                                                                                             >  1:305255/1‑251 (MQ=255)                                   ggcggTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGATCGGc                                                                                                                                                                                                          >  2:798476/1‑250 (MQ=255)                                    gcggTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCAt                                                                                                                                                                                                        <  1:930517/251‑1 (MQ=255)                                      ggTCAATTTCGTCCATTCGATCCCACGGGTTATCTATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAAGTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACgg                                                                                                                                                                                                                     <  2:548157/236‑1 (MQ=255)                                      ggTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACgg                                                                                                                                                                                                                     >  1:548157/1‑236 (MQ=255)                                      agTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCTCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGACGATAATTAGCCGACCTTTAGCCTGGGCAAAACA‑GACG‑CCGTTTTTAAAGCCGGGCTCGCGCC‑ACACCAACTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGt                                                                                                                                                                                                      >  2:787095/2‑251 (MQ=255)                                         cAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGttt                                                                                                                                                                                                    <  1:660940/250‑1 (MQ=255)                                               cGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTgcag                                                                                                                                                                                              <  1:389720/250‑1 (MQ=255)                                                  ccATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCAtt                                                                                                                                                                                          >  2:297300/1‑251 (MQ=255)                                                  ccATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCAtt                                                                                                                                                                                          >  1:72330/1‑251 (MQ=255)                                                   catccgaccccACGGGTTACCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGCCTATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTc                                                                                                                                                                                         <  2:682494/243‑1 (MQ=255)                                                    aTTCGATACCACGGGATATAGGTAGCTTTACCTTCGCCTAATAACTCCTGAGCGACGGTC‑TCCAGCGAGAATGAAGAGCAATTCCAGAACGCGGATTTCAGCGACTCGATACCATCGCTAATTAGCCGACCTTTCGCCTGGGCAAACAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGCCACCCATCA‑GCAGCGGCATACGGTACCGCTCGGAATGtttttccggcatacg                                                                                                                                                                                        >  2:410925/1‑241 (MQ=255)                                                       cgaccccACGGGTTATCGTTAGATGTTCCTGCGCCTAAGAACTCCTGTGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGACCGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAg                                                                                                                                                                                     <  2:954131/247‑1 (MQ=255)                                                       cGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCa                                                                                                                                                                                      <  1:840889/250‑1 (MQ=255)                                                        gATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGAGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGc                                                                                                                                                                                       >  1:509447/1‑248 (MQ=255)                                                         aTCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGAC‑GTTCTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATCCCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GTCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGAt                                                                                                                                                                                   <  1:406753/251‑1 (MQ=255)                                                           cccACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGc                                                                                                                                                                                       >  1:327960/1‑245 (MQ=255)                                                             cACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑ACCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTATTGCAGCATTCgacgaaatg                                                                                                                                                                                >  2:153279/1‑242 (MQ=255)                                                                    tATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTGTGCAGCATTCGCAGATCGAACTGcac                                                                                                                                                                         <  1:788722/250‑1 (MQ=255)                                                                    tATCGATAGATTTTCCTGCGCCGAATAACCCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGcaa                                                                                                                                                                         <  2:228135/250‑2 (MQ=255)                                                                       cGATAGATTTTCATTCGCCTAATAACTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACg                                                                                                                                                                     >  1:741151/1‑251 (MQ=255)                                                                                   cTTCGCCTAATAGCTCCTGAGCGACAGTT‑TCCAGCGAGAATGAAGAGAAATTCTAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCGCTATTATCGC‑GCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGt                                                                                                                                                                    >  1:682542/1‑240 (MQ=255)                                                                                   cTTCGCCTAATAGCCCGTGAGCGACAGTT‑TCCAGCGAGAATGAAGAGAAATTCTAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGCCCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCGCTATTATCGC‑GCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGt                                                                                                                                                                    <  2:682542/240‑1 (MQ=255)                                                                                      cGCCTAAAAACTCCTGAGCTACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACTCGGATTTCAGCGCCCCGATACCGTCGATAATTAGCCGACCTCTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATc                                                                                                                                                      <  2:183517/251‑1 (MQ=255)                                                                                             taacTCCTGAGCGACGGTT‑TCCATCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAACGCCTCGATACCGTCGATACTTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTCGCAGCATTCGCAGATCGCACTGCACCACTTTCCACCCgcgtatcacatc                                                                                                                                                 >  2:611922/4‑238 (MQ=255)                                                                                               acTCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAgg                                                                                                                                              >  1:316078/2‑250 (MQ=255)                                                                                                 tCCTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GACCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCGTTCGCAGATCGAACTGCACCCCGTTCCAACCGAGGATAACATCAggccc                                                                                                                                           >  2:282778/1‑248 (MQ=255)                                                                                                  ccTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCg                                                                                                                                          <  1:198777/251‑1 (MQ=255)                                                                                                   cTGAGCGACGGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAC‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCCGGATCGt                                                                                                                                         >  2:229483/1‑251 (MQ=255)                                                                                                     gAGCGACGGTT‑TACAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAg                                                                                                                                       <  1:676429/251‑1 (MQ=255)                                                                                                         gacggTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGCTTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAC‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGACCTGCACCACGTTCCAACCGATGATCACATCAGGCTCGTAGTTgg                                                                                                                                   >  2:651521/5‑251 (MQ=255)                                                                                                                      cGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATCCCGTCGATAATTAGCCGCCCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCGTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGt                                                                                                                       <  1:695633/251‑1 (MQ=255)                                                                                                                        agaATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTg                                                                                                                     <  1:825712/251‑1 (MQ=255)                                                                                                                         gaATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCCCTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCCCGTTCCAACCGATGATCCCATCAGGCTCGTGGTTGGCAAACCAGGCGTTgg                                                                                                                    >  2:370960/1‑250 (MQ=255)                                                                                                                               agagAAATTCCAGAACGAGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGCGttttt                                                                                                              >  2:919965/1‑251 (MQ=255)                                                                                                                                gagaAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGCATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCCCCAAGTTCCAACCGATGATCACATCAGGTTCGTAGTTGGCAAACCAGGCGTTGAGTTTTTc                                                                                                             >  2:935836/1‑251 (MQ=255)                                                                                                                                gagaAAGTCCAGACCGCGGATTTCAGCGCCTCGATTCCGCTGATAATTAGCCGACCCTTATCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGAGCGATCTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTc                                                                                                             <  2:885295/251‑1 (MQ=255)                                                                                                                                    aataccAGAACCCTGATTTCAGCGCCTCGATTCCGGAGATAATTATCCGCCCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTGTTTAAAGCCGTGCTCGCGCC‑TCGCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAgt                                                                                                         <  2:361865/247‑2 (MQ=255)                                                                                                                                     aTTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAgt                                                                                                         >  1:246798/1‑249 (MQ=255)                                                                                                                                       tCCAGAACGCGGATCTCAGCGCCCCGATACCGTCGATAATTAGCCCGCCTTGAGCCTGGGCAAAAAA‑GACT‑CCGATTTTAAAGCCGTGCTCGCGCC‑TCTCCAGCGCACTGTTATCGC‑GCCTCAGCCGCAGCGGAATACGGTAACGCTCGTCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAgtag                                                                                                       <  2:789825/250‑4 (MQ=255)                                                                                                                                              cgcgGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTgcgggc                                                                                               <  1:767449/251‑1 (MQ=255)                                                                                                                                                cgGATTTCAGCGCCTCGATACCGTCGATAATTAGCCTACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACAACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTgcgggcgg                                                                                             <  1:618221/251‑1 (MQ=255)                                                                                                                                                cgGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTgcgggcgg                                                                                             >  1:936929/1‑251 (MQ=255)                                                                                                                                                  gttttCAGCGCCACGATACCGTCGATAATTAGCCGCCCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCt                                                                                           <  2:489983/249‑1 (MQ=255)                                                                                                                                                      tCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCg                                                                                       >  2:229987/1‑251 (MQ=255)                                                                                                                                                        cgcgcCTCGTTACCGTCGATAATTAGCCGACCTTTCGCCTGGGCAAAAAA‑GACC‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGAc                                                                                     <  2:313509/250‑1 (MQ=255)                                                                                                                                                        aGCGCCTCGATACCGACGATATTTAGCCGTCTTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGTTGTGCTCGCGCC‑ACTCCAGCTCGCTGTTATCGC‑GCTCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGACCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGTCGAc                                                                                     <  1:611922/251‑1 (MQ=255)                                                                                                                                                           gcCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGt                                                                                   >  1:293319/1‑250 (MQ=255)                                                                                                                                                              tCGATACAGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCACCCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTc                                                                               >  2:908327/1‑251 (MQ=255)                                                                                                                                                               cGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑TACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTc                                                                               <  1:387540/250‑1 (MQ=255)                                                                                                                                                               cGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTcc                                                                              <  1:798476/251‑1 (MQ=255)                                                                                                                                                                gATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCa                                                                             <  1:556249/251‑1 (MQ=255)                                                                                                                                                                     cGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTc                                                                         >  1:223311/1‑250 (MQ=255)                                                                                                                                                                       ccGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTcaaa                                                                      <  1:587659/250‑4 (MQ=255)                                                                                                                                                                        cGATAATTAGCCGACCTTTAGCCTGGGAAAAAAA‑GACG‑CCGTTTTTAAGGCCGTGCTCGCGCC‑ACTCCAGCTCCCTGTTATCGC‑GCCACAGCCGCAGCGGAATACGGTAACGCTGGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACAACGTTACAACCGATGCTCACATCAGGGTCGTAGTTGGCAAACCAGGCGATGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGttaaaa                                                                      >  2:581017/1‑246 (MQ=255)                                                                                                                                                                         gATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAAt                                                                    >  1:799297/1‑251 (MQ=255)                                                                                                                                                                             aTTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATcg                                                                  >  1:545960/1‑248 (MQ=255)                                                                                                                                                                                aGCCGACCTTTAGCCTGGGCAAAAAA‑GCCG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCgc                                                             >  1:330562/1‑250 (MQ=255)                                                                                                                                                                                  ccGACCTTTGGCCGGGGCAAAAAATGA‑G‑CCGTTTTTAAAGCAGTGCCCGCGCC‑ACGCCAGCTCACTGTTAACGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTGCCAACCGATGATCACAGCAGGATCTTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGGAGCTGCGGGCGTCTGGCGACGTATTCCAGTTTAAAATCGAGCgcgg                                                           <  2:97478/251‑4 (MQ=255)                                                                                                                                                                                   cGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCgcgga                                                          <  1:228499/251‑1 (MQ=255)                                                                                                                                                                                         ttaaccTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCCCAC‑CCAGCTCACGGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCg                                                      <  2:462254/245‑1 (MQ=255)                                                                                                                                                                                            cccTGGGAAAAAAA‑G‑CGTCCGTTTTTAAAGCCTTGCTCGCGCC‑CCTCCATCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGGTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTGCAAAATCGAGCGCGGATGCGTCGcc                                                 <  2:478131/250‑1 (MQ=255)                                                                                                                                                                                                gggCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAAACAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGACGCCATtc                                             >  1:731196/1‑251 (MQ=255)                                                                                                                                                                                                  gCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATtctc                                           >  2:515522/1‑251 (MQ=255)                                                                                                                                                                                                     aaaaa‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATtctctgg                                        >  1:599190/1‑248 (MQ=255)                                                                                                                                                                                                        aa‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGccc                                     >  2:227137/1‑251 (MQ=255)                                                                                                                                                                                                                ccGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCCTGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGCGGCGGCCCCATTATCTGGCCCCAGcct                               >  2:672799/1‑249 (MQ=255)                                                                                                                                                                                                                   tttttAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCatat                             >  1:516656/1‑250 (MQ=255)                                                                                                                                                                                                                   tttttAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGCAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCatat                             >  1:447486/1‑250 (MQ=255)                                                                                                                                                                                                                     tttAAAGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATaa                           >  1:288511/1‑250 (MQ=255)                                                                                                                                                                                                                        aaaGCCGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGAATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCCCTGGCCCCAGCCTACAAACg                        >  2:948984/1‑250 (MQ=255)                                                                                                                                                                                                                            ccGTGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGACGCCCTTCTCTGGCCCCAGCATAAAAACGATggg                   >  2:542203/1‑249 (MQ=255)                                                                                                                                                                                                                               tGCTCGCGCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTg                >  1:213774/1‑251 (MQ=255)                                                                                                                                                                                                                                gcactcgcC‑CCTCCATCTCACTGTTATCTC‑GCCCCAGCCGCAGCTGAATACTGTAACGCTCGGCATGTTTTTGCAGCATTCGCATATCGAACTGCACCCCGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAGCGATGCGCTg                <  1:212058/245‑1 (MQ=255)                                                                                                                                                                                                                                  ccgcgCC‑CCTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGccc             <  1:104840/250‑1 (MQ=255)                                                                                                                                                                                                                                   cgcgCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCg            >  1:868290/1‑251 (MQ=255)                                                                                                                                                                                                                                   cgcgCC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCg            <  1:746808/251‑1 (MQ=255)                                                                                                                                                                                                                                      gcC‑ACTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAg         <  1:892895/251‑1 (MQ=255)                                                                                                                                                                                                                                        c‑ACTCCAGCTCACTGTTATAGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCAGATGGCGACGTAATCCAGTTCAAAAACGAGCGCGGAGGCGTCGCCATTCGCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGcc       >  2:570391/1‑251 (MQ=255)                                                                                                                                                                                                                                           cTCCAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATACACGATGCGCAGCCCGCATCCtt     >  2:820948/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCTGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCtt     <  1:429466/249‑1 (MQ=255)                                                                                                                                                                                                                                             ccAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGGAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCtt     <  1:86245/249‑1 (MQ=255)                                                                                                                                                                                                                                              cAGCTCACTGTTATCGC‑GCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCa  <  1:908327/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   TAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTT‑TCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCCTGGGCAAAAAA‑GACG‑CCGTTTTTAAAGCCGTGCTCGCGCC‑ACTCCAGCTCGCTATTATCGC‑GCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCA  >  NC_000913/64789‑65272

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