Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	501,181	A→C	21.6%	T353P (ACT→CCT)	gsk →	inosine/guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,181	0	A	C	21.6%	150.7 / 8.1	78	T353P (ACT→CCT)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base A (33/26);  new base C (4/14);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGT  >  NC_000913/500949‑501412                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         gggatGGGTAGTTCGGGTGCTGGTCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATg                                                                                                                                                                                                                                   <  2:314403/235‑1 (MQ=255)      gggTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAACGGAATTCAACCAGTAGGAGTTTAGCCGCCCC‑ATTGGCCACAAGG‑GTTGCCCGAATCCACAGCGGGTAAATTTGCCCATTTCGC‑CGTACATGGGGGGGCCCGCAAAAAACATTAAAA‑CAAAAGGAgag                                                                                                                                                                                                                              >  1:542623/1‑237 (MQ=255)       ggTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCCCAAGG‑ATTTCCAGACTCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGCAAAAATCATGAACA‑ATAATGGCGCgg                                                                                                                                                                                                                             >  2:96572/1‑239 (MQ=255)       ggTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCCGAATCCGCTGCGTGTATATTCGCACATTTCGC‑CGTACATGGGCGGGCCGGGAAAAATCATGAACA‑CTAATGGAGCgg                                                                                                                                                                                                                             >  2:875725/1‑239 (MQ=255)            aTCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATTTGGCGGGCTTTTCCGGAGAAGAAGGGAAACGTAAACCCCAGCATCCGCTGCTGCCGGGCGATCTAGCGGCATTCACCCAGTCAGAGTTTAGGCGCGCC‑ATGCCGCCCAACG‑ATTGTCACAATTCTCCGCGGGTATCTTTCCCCCTTGCGG‑CCGACCTCTGCGGGCCGGGAACAA‑CACTAACA‑CCACTGAAGggggggggg                                                                                                                                                                                                                       >  1:18150/1‑236 (MQ=255)            aTCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGGATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGGGTATATTTGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATACTGAACC‑CTAATGGAGCggggggt                                                                                                                                                                                                                        >  2:682431/1‑237 (MQ=255)                 gTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTAACCAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCCGCCGGGGGCTATAGCGGAAATTAACCAGTATGAGTTTAGCCGCCCC‑ATGGGCCCCACGG‑GTTGCCAGAATCCGCTGCGTGGATATTCACACAATGCGC‑CGGACCTGGGCGGGCCGGACAAAAACAATAACA‑Caaaaaaacgggggggaggggg                                                                                                                                                                                                                   >  1:237414/1‑221 (MQ=255)                    cgtggCACCGCCGTGCCAATGGGCTTGTATATGGAGGGCTTTACCGAAGACGAAGCGAAACGTAACTCCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCAt                                                                                                                                                                                                                 <  2:797584/234‑1 (MQ=255)                      tggccccgccgGGCCAATCGGTTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTgg                                                                                                                                                                                                              <  1:753977/234‑1 (MQ=255)                      gtgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAGACGTCAAACCCAGCATTCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCCCAATGCGTTTAGCCGGGCC‑ATGCGCAACCAAG‑ACTGCCCGAATACGCCGCACGGATATTCCAACATTGCCC‑CGTGCATTGGGGGGTCGGCAAAAATACTCGACA‑CTCATGGGGCGGGGGACGGCGGAATgg                                                                                                                                                                                                              >  1:157743/1‑239 (MQ=255)                       tgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGACGCGAAAACTAAAAACCACCCACCCCCCCCGCCGGGGGCCATCGAGGAATTCAACCAGCATGAATTTACCGCCGCC‑AGGCGCCCCAAGG‑ATTGCCCCAATCCGCTCCGGGGAAATTGCCACATCGCCC‑CCGACCGGGGCGGGGCGGCAAAAAACACTAAAA‑CTAAAGGAACGGGGGGGGGGGCATCggg                                                                                                                                                                                                             >  1:207286/1‑238 (MQ=255)                       tgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGGGCTATAGCGGGATTTAACCAGTATGCGTTTAGCCGGGCC‑ATGCGCCCCAAAG‑ATTGCCAGAAGCCGATCCCGGTAATTTTGCCAATTGCGCACGT‑CCTTGGCGGGCCGGAAAAAA‑CAAGAAAACCCAATGGAACggggggagggcaattggc                                                                                                                                                                                                             >  1:384875/1‑226 (MQ=255)                       tgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACC‑CTAATGGAGCGGGGGGTGGGGCATTGGc                                                                                                                                                                                                             >  2:753977/1‑239 (MQ=255)                       tgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAAAC‑CTAATGGTACCGGGGATGGGCCATTGGc                                                                                                                                                                                                             >  2:126607/1‑239 (MQ=255)                                                         gCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCCTGAACC‑CTAATGGAGCGGGGGCTGGCGCATTGGCACGGGTGCTGCATGAAATTaacccccaaaaccac                                                                                                                                                                           >  1:464678/1‑225 (MQ=255)                                                               ccGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCACCAGCTACCATCgg                                                                                                                                                                     >  2:90809/1‑238 (MQ=255)                                                                   aGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGACCA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCa                                                                                                                                                                 >  2:30874/1‑239 (MQ=255)                                                                         aGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTACAAACGTAc                                                                                                                                                           >  2:869960/1‑239 (MQ=255)                                                                               aCGAAAACCCCAGACTCCCCGTCTGCCGGGGCGTATAGGGAAATTCAACCATAATGAGTTTAGCCTGGCC‑AGGCGCCACAAGG‑ATGGCCAGAATCCGCGGCGTTTATATTCGCACATTGCGC‑GGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACt                                                                                                                                                     <  1:74098/239‑1 (MQ=255)                                                                                     accccgtgctcccgtgcgcccGGGGGGTATAGAGGAATTCATCCATTGTGAGTTTGGCGGCGCA‑AGCCGCCACAAGG‑ATTGCCAGCATCCGCTGCGTGTAAATTCGCACGTTGCGC‑GGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CAAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCACCGTACCAAACTCCAGCa                                                                                                                                               <  1:608040/220‑1 (MQ=255)                                                                                     aaCCCAGCATCCGCTGCGGCCGGGCGCTAAAGCGGAATTCAACCAGTATGAGTTTAGCCGCCCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCa                                                                                                                                               <  2:542623/239‑1 (MQ=255)                                                                                          aGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgctg                                                                                                                                                                                                   >  2:681956/1‑182 (MQ=255)                                                                                          aGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgctg                                                                                                                                                                                                   <  1:681956/182‑1 (MQ=255)                                                                                          aGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCCGCAAACAt                                                                                                                                          >  2:169304/1‑239 (MQ=255)                                                                                           gCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAAACCCAGCaaaaaaa                                                                                                                                         >  1:560583/1‑235 (MQ=255)                                                                                                                 taaGCGGAATTCCACCAGTTTGTGTTTACCCGCGCC‑ATG‑GCCCCAAGGTGTTGCCAGAATCCGGTGCGTTTAAATTCGCACATTGCGC‑CGTCCATGGGCGGGCCGGAAAAAATCATGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTc                                                                                                                   <  1:217837/237‑1 (MQ=255)                                                                                                                  taGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCACCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCa                                                                                                                  >  2:701801/1‑239 (MQ=255)                                                                                                                         aTTCAACCAGGATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCCCATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAg                                                                                                           <  2:780533/239‑1 (MQ=255)                                                                                                                           tCAACCAGTTTGGGTTT‑GTCCCGCCAAGGCGCGAAAATG‑ATTGCCAGAATCCGCGGCGGTGATATTCCCACATTGGGT‑CGTTCATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAgcg                                                                                                         <  1:650304/239‑1 (MQ=255)                                                                                                                            tcaccggtttggttttGGCCGCGCC‑ATGCGCCAAAGGG‑ATTGCCAGGACCGGCGGGGTGTATATTCGCGCATTGGGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAgcgc                                                                                                        <  1:682431/226‑1 (MQ=255)                                                                                                                              accattttgggtttgaccggccc‑ATGGGCCACAGGG‑ATTGCCAGAATCCGGGAGGTTTAAAGTCGCACATTGCGC‑CGGTCATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAg                                                                                                      <  1:390006/218‑1 (MQ=255)                                                                                                                              aCCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCAc                                                                                                                                  <  1:683990/211‑1 (MQ=255)                                                                                                                              aCCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCAc                                                                                                                                  >  2:683990/1‑211 (MQ=255)                                                                                                                              aCCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAg                                                                                                      <  2:657082/239‑1 (MQ=255)                                                                                                                                 ggttggAGTTTTGCCGCGCC‑ATGGGCCACAAGG‑ATTGCCAGAACCCGCGGGGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGgtg                                                                                                   <  1:100133/234‑1 (MQ=255)                                                                                                                                   gtttggTTTAGCCCGGCC‑ATGGGCCACAAGG‑ATTGCCAAAACCCGCTGGGTGTATATTCGGACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGgt                                                                                                    <  1:668895/231‑1 (MQ=255)                                                                                                                                     tGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTa                                                                                               <  2:560583/239‑1 (MQ=255)                                                                                                                                       cgTTTGGCCCCCCC‑ATGCCCCACAGAG‑TGTGACAAAAACCGCTGCGTGGAGATGCCCAAATTGCCC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTTCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTaaa                                                                                             <  1:475834/238‑1 (MQ=255)                                                                                                                                        ttttAGCCGGGCC‑CTGCGCCACAAGG‑ATTGCCAGAATCCGCGGCGTGTATGTTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTaaa                                                                                             <  2:368097/237‑1 (MQ=255)                                                                                                                                        gTTTACGCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTGCCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTaaa                                                                                             <  2:207286/238‑1 (MQ=255)                                                                                                                                             gCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCCAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGACGGTGTGTAAATATGCt                                                                                       >  2:644857/1‑239 (MQ=255)                                                                                                                                             gCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCCGGTGTGTAAATATGCt                                                                                       >  2:893647/1‑239 (MQ=255)                                                                                                                                             gCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAatat                                                                                          >  2:668895/1‑236 (MQ=255)                                                                                                                                             gCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACAAAAATTCACCCGGTTAACTTATTCATCGTTAGCGCAggggggaaaaaataca                                                                                       >  1:472633/1‑225 (MQ=255)                                                                                                                                                 cgcC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAAcc                                                                                   >  2:355372/1‑239 (MQ=255)                                                                                                                                                   cc‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGACAACATAAATTCACCCGGTTAACTTATTTATCTTTAGCGCCgggggtgaaaaagactaacccg                                                                                 >  1:873432/1‑219 (MQ=255)                                                                                                                                                         cgcCACAAGG‑ATTGCCAAAATCCGTTGGGTTTATTTTCGCACATTGGGC‑CGGACATGGGGGGGGCGGAAAAAATCATGAACA‑CTAATGGGGCGGGGGATGGCCCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAg                                                                            <  1:869960/239‑1 (MQ=255)                                                                                                                                                          gcCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAg                                                                            <  2:237414/238‑1 (MQ=255)                                                                                                                                                          gcCACAAGG‑ATTGCCAAAATCCGGTGCGTGTATATTCGGAAATTGCGC‑CGTACAGGGGCGGGCCGGAAAAAATCAGGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCAACTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAg                                                                            <  2:745719/238‑1 (MQ=255)                                                                                                                                                           atgtgtataggagacAGAATCCGCTGTGTGTATATTCGCACATTGCGC‑CGTACAGGGGCGGGCCGGAAAAAATCATGAACT‑CTAATGGAGCGGGGGATGGCGCATTGGAAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAg                                                                            <  2:699457/224‑1 (MQ=255)                                                                                                                                                                     ttGCCAGAATCCGGTGCGTGTAAATTCGCACATTGCGC‑CGTACAGGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTa                                                                  <  2:538365/238‑1 (MQ=255)                                                                                                                                                                          aGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCTCAGGGGGGTAACTATGCTAACCGGGTGAGCTATCAGGCACTTaa                                                             >  2:565642/1‑238 (MQ=255)                                                                                                                                                                             acccgggtggtggATTTTCGCACCTTGCGG‑CGGACATGGGCGGGTCGGAAAAAATCATGAACA‑ATAATGGGGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAAc                                                            <  1:96572/227‑1 (MQ=255)                                                                                                                                                                                     cGTTTATATGCGCAACTGTCGC‑CGGACATGGGCGGGCCGGAAAAAATCATGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCAc                                                  <  2:381082/238‑1 (MQ=255)                                                                                                                                                                                       tgtATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCCTCGTAGCAACGTACCAAACACCAGCAAACATAAATTAACCTGGTTAACTTATCCATCGTTATCGCCGGTGTGTAAATATCATAACCGGGTGAGGTATCAGGTa                                                                  >  1:699457/1‑220 (MQ=255)                                                                                                                                                                                          atatTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGAGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTAc                                                                                                                                                                           <  2:597497/112‑1 (MQ=255)                                                                                                                                                                                                 cctaTTGGGC‑CGTTAAGGGGCGGGCCGGAAAAAAACAGGAACC‑CTAAGGGAGCGGGGGATGGGTGATTGGCCGCGTTGCTGCATGCCATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAAcgcg                                     <  1:126607/236‑1 (MQ=255)                                                                                                                                                                                                 cacaTTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTAc                                                                 >  1:11004/1‑211 (MQ=255)                                                                                                                                                                                                 cacaTTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTAc                                                                 <  2:11004/211‑1 (MQ=255)                                                                                                                                                                                                   caTTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGAAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGGGTAAATATGCAAACCGTGTGGTCATTCAGGGATTGAACCACCATTCACCCTCTTTAAcgccg                                    >  1:772753/1‑236 (MQ=255)                                                                                                                                                                                                    aTTGCGC‑CGTACATGGGCGGGCCGGAAAAAAGCATGAACA‑CTAATGGAGCGGGGGCTGGGGCATTGGCAGCGTTGGTGCATGACATTACCGCCAACAGATACCATCGTAGCa                                                                                                                                                                 >  1:597497/1‑112 (MQ=255)                                                                                                                                                                                                        cgc‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCtt                               >  2:356614/1‑238 (MQ=255)                                                                                                                                                                                                          c‑cGTACATGGGCGGGCCGGAAAAAATCATGCCCA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTTCTGCATGACCTTACCGCCAACAGCTACCATCGTCGCAACGTACCACACTCCCGCACACATAAATTCACCTGGTTAACTTATTCATCTTTAGCGCTGGTGCTTCACTTGTCTAACCCTGTGACTTATCCGGTTCTACTCCTGCCTTCCCCCCGTTTCACAGGCGTCTTGcc                            >  1:791433/1‑239 (MQ=255)                                                                                                                                                                                                             gttcctGGGGGGGCCGGAAAAAATCAGGAACC‑CTAATGGGGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCgg                          <  1:644857/234‑1 (MQ=255)                                                                                                                                                                                                             gTACATGGGCGGGCCGGAAAAAAGCATGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCGGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCgg                          <  2:157743/239‑1 (MQ=255)                                                                                                                                                                                                             gTACAGGGGGGGGCCGAAAAAAATCATGAACC‑CAAAGGGAGCGGGGGAGGGCGCATTGGCGGGGTGCCTGCATGCCATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCgg                          <  1:445088/239‑1 (MQ=255)                                                                                                                                                                                                                   gggcgggcCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTg                    >  2:499284/1‑239 (MQ=255)                                                                                                                                                                                                                     gcgggcCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCgggactggaa                  >  2:385762/1‑231 (MQ=255)                                                                                                                                                                                                                       gggcCGGAAAAAATCATGAACA‑TTAAGGGAGCGGGGGATGGCTCATTGGCAGCGTTTCTGCATTACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCCAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGa                <  1:674046/239‑1 (MQ=255)                                                                                                                                                                                                                        ggcCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCCAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGGGAAGAc               >  2:674046/1‑239 (MQ=255)                                                                                                                                                                                                                        ggcCGGAAAAAACCATGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGa                <  2:474602/238‑1 (MQ=255)                                                                                                                                                                                                                          ccGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGGGCTATCAGGAACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTGCCCGGAGCGGAAAGACa              >  1:474602/1‑238 (MQ=255)                                                                                                                                                                                                                               aaaaaTCAAGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTgg        <  2:218080/239‑1 (MQ=255)                                                                                                                                                                                                                                 aaaTCATGAACC‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa      <  2:472633/239‑1 (MQ=255)                                                                                                                                                                                                                                 aaaTCATGAACA‑CTAATGGGGCGGGGGAAGGGTGATTGGAAGCGTTGCTGAATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa      <  2:18150/239‑1 (MQ=255)                                                                                                                                                                                                                                 aaaTCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa      <  2:491788/239‑1 (MQ=255)                                                                                                                                                                                                                                     ctggAACC‑CTACTGGGGCGGGGGTTGGCGCTTTGGCAGCGTTGCTGGATGACATTACCGCCAACAGCTACCAACGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGt  <  1:866228/236‑1 (MQ=255)                                                                                                                                                                                                                                     cATGAACC‑ATAATGGAGGGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGt  <  2:239332/239‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         GGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGC‑CGTACATGGGCGGGCCGGAAAAAATCATGAACA‑CTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGT  >  NC_000913/500949‑501412

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