Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	501,181	A→C	23.1%	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	23.1%	136.0 / 8.6	69	T353P (ACT→CCT)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base A (29/23);  new base C (3/14);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGT  >  NC_000913/500949‑501412                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        gggatGGGTAGTTCGGGTGCTGGTCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATg                                                                                                                                                                                                                                   <  1:314403/235‑1 (MQ=255)      gggTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAACGGAATTCAACCAGTAGGAGTTTAGCCGCCCC‑ATTGGCCACAAGG‑GTTGCCCGAATCCACAGCGGGTAAATTTGCCCATTTCGCCGTACATGGGGGGGCCCGCAAAAAACATTAAAACAAAAGGAgag                                                                                                                                                                                                                              >  2:542623/1‑237 (MQ=255)       ggTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCCCAAGG‑ATTTCCAGACTCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGCAAAAATCATGAACAATAATGGCGCgg                                                                                                                                                                                                                             >  1:96572/1‑239 (MQ=255)       ggTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCCGAATCCGCTGCGTGTATATTCGCACATTTCGCCGTACATGGGCGGGCCGGGAAAAATCATGAACACTAATGGAGCgg                                                                                                                                                                                                                             >  1:875725/1‑239 (MQ=255)            aTCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATTTGGCGGGCTTTTCCGGAGAAGAAGGGAAACGTAAACCCCAGCATCCGCTGCTGCCGGGCGATCTAGCGGCATTCACCCAGTCAGAGTTTAGGCGCGCC‑ATGCCGCCCAACG‑ATTGTCACAATTCTCCGCGGGTATCTTTCCCCCTTGCGGCCGACCTCTGCGGGCCGGGAACAA‑CACTAACACCACTGAAGggggggggg                                                                                                                                                                                                                       >  2:18150/1‑236 (MQ=255)            aTCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGGATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGGGTATATTTGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATACTGAACCCTAATGGAGCggggggt                                                                                                                                                                                                                        >  1:682431/1‑237 (MQ=255)                    cgtggCACCGCCGTGCCAATGGGCTTGTATATGGAGGGCTTTACCGAAGACGAAGCGAAACGTAACTCCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCAt                                                                                                                                                                                                                 <  1:797584/234‑1 (MQ=255)                      tggccccgccgGGCCAATCGGTTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTgg                                                                                                                                                                                                              <  2:753977/234‑1 (MQ=255)                      gtgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAGACGTCAAACCCAGCATTCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCCCAATGCGTTTAGCCGGGCC‑ATGCGCAACCAAG‑ACTGCCCGAATACGCCGCACGGATATTCCAACATTGCCCCGTGCATTGGGGGGTCGGCAAAAATACTCGACACTCATGGGGCGGGGGACGGCGGAATgg                                                                                                                                                                                                              >  2:157743/1‑239 (MQ=255)                       tgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGACGCGAAAACTAAAAACCACCCACCCCCCCCGCCGGGGGCCATCGAGGAATTCAACCAGCATGAATTTACCGCCGCC‑AGGCGCCCCAAGG‑ATTGCCCCAATCCGCTCCGGGGAAATTGCCACATCGCCCCCGACCGGGGCGGGGCGGCAAAAAACACTAAAACTAAAGGAACGGGGGGGGGGGCATCggg                                                                                                                                                                                                             >  2:207286/1‑238 (MQ=255)                       tgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACCCTAATGGAGCGGGGGGTGGGGCATTGGc                                                                                                                                                                                                             >  1:753977/1‑239 (MQ=255)                       tgCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAAACCTAATGGTACCGGGGATGGGCCATTGGc                                                                                                                                                                                                             >  1:126607/1‑239 (MQ=255)                                                               ccGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCACCAGCTACCATCgg                                                                                                                                                                     >  1:90809/1‑238 (MQ=255)                                                                   aGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGACCACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCa                                                                                                                                                                 >  1:30874/1‑239 (MQ=255)                                                                         aGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCACCGTACAAACGTAc                                                                                                                                                           >  1:869960/1‑239 (MQ=255)                                                                               aCGAAAACCCCAGACTCCCCGTCTGCCGGGGCGTATAGGGAAATTCAACCATAATGAGTTTAGCCTGGCC‑AGGCGCCACAAGG‑ATGGCCAGAATCCGCGGCGTTTATATTCGCACATTGCGCGGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACt                                                                                                                                                     <  2:74098/239‑1 (MQ=255)                                                                                     aaCCCAGCATCCGCTGCGGCCGGGCGCTAAAGCGGAATTCAACCAGTATGAGTTTAGCCGCCCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCa                                                                                                                                               <  1:542623/239‑1 (MQ=255)                                                                                          aGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgctg                                                                                                                                                                                                   <  2:681956/182‑1 (MQ=255)                                                                                          aGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgctg                                                                                                                                                                                                   >  1:681956/1‑182 (MQ=255)                                                                                          aGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCCGCAAACAt                                                                                                                                          >  1:169304/1‑239 (MQ=255)                                                                                           gCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAAACCCAGCaaaaaaa                                                                                                                                         >  2:560583/1‑235 (MQ=255)                                                                                                                 taaGCGGAATTCCACCAGTTTGTGTTTACCCGCGCC‑ATG‑GCCCCAAGGTGTTGCCAGAATCCGGTGCGTTTAAATTCGCACATTGCGCCGTCCATGGGCGGGCCGGAAAAAATCATGAACCCTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTc                                                                                                                   <  2:217837/237‑1 (MQ=255)                                                                                                                  taGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCACCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCa                                                                                                                  >  1:701801/1‑239 (MQ=255)                                                                                                                         aTTCAACCAGGATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCCCATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAg                                                                                                           <  1:780533/239‑1 (MQ=255)                                                                                                                           tCAACCAGTTTGGGTTT‑GTCCCGCCAAGGCGCGAAAATG‑ATTGCCAGAATCCGCGGCGGTGATATTCCCACATTGGGTCGTTCATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAgcg                                                                                                         <  2:650304/239‑1 (MQ=255)                                                                                                                              aCCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCAc                                                                                                                                  <  2:683990/211‑1 (MQ=255)                                                                                                                              aCCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCAc                                                                                                                                  >  1:683990/1‑211 (MQ=255)                                                                                                                              aCCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAg                                                                                                      <  1:657082/239‑1 (MQ=255)                                                                                                                                 ggttggAGTTTTGCCGCGCC‑ATGGGCCACAAGG‑ATTGCCAGAACCCGCGGGGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGgtg                                                                                                   <  2:100133/234‑1 (MQ=255)                                                                                                                                   gtttggTTTAGCCCGGCC‑ATGGGCCACAAGG‑ATTGCCAAAACCCGCTGGGTGTATATTCGGACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACCCTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGgt                                                                                                    <  2:668895/231‑1 (MQ=255)                                                                                                                                     tGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTa                                                                                               <  1:560583/239‑1 (MQ=255)                                                                                                                                       cgTTTGGCCCCCCC‑ATGCCCCACAGAG‑TGTGACAAAAACCGCTGCGTGGAGATGCCCAAATTGCCCCGTACATGGGCGGGCCGGAAAAAATCATGAACCCTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTTCCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTaaa                                                                                             <  2:475834/238‑1 (MQ=255)                                                                                                                                        ttttAGCCGGGCC‑CTGCGCCACAAGG‑ATTGCCAGAATCCGCGGCGTGTATGTTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTaaa                                                                                             <  1:368097/237‑1 (MQ=255)                                                                                                                                        gTTTACGCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTGCCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTaaa                                                                                             <  1:207286/238‑1 (MQ=255)                                                                                                                                             gCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCCAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGACGGTGTGTAAATATGCt                                                                                       >  1:644857/1‑239 (MQ=255)                                                                                                                                             gCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCCGGTGTGTAAATATGCt                                                                                       >  1:893647/1‑239 (MQ=255)                                                                                                                                             gCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAatat                                                                                          >  1:668895/1‑236 (MQ=255)                                                                                                                                                 cgcC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAAcc                                                                                   >  1:355372/1‑239 (MQ=255)                                                                                                                                                         cgcCACAAGG‑ATTGCCAAAATCCGTTGGGTTTATTTTCGCACATTGGGCCGGACATGGGGGGGGCGGAAAAAATCATGAACACTAATGGGGCGGGGGATGGCCCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAg                                                                            <  2:869960/239‑1 (MQ=255)                                                                                                                                                          gcCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAg                                                                            <  1:237414/238‑1 (MQ=255)                                                                                                                                                          gcCACAAGG‑ATTGCCAAAATCCGGTGCGTGTATATTCGGAAATTGCGCCGTACAGGGGCGGGCCGGAAAAAATCAGGAACCCTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCAACTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAg                                                                            <  1:745719/238‑1 (MQ=255)                                                                                                                                                                     ttGCCAGAATCCGGTGCGTGTAAATTCGCACATTGCGCCGTACAGGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTa                                                                  <  1:538365/238‑1 (MQ=255)                                                                                                                                                                          aGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCTCAGGGGGGTAACTATGCTAACCGGGTGAGCTATCAGGCACTTaa                                                             >  1:565642/1‑238 (MQ=255)                                                                                                                                                                             acccgggtggtggATTTTCGCACCTTGCGGCGGACATGGGCGGGTCGGAAAAAATCATGAACAATAATGGGGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAAc                                                            <  2:96572/227‑1 (MQ=255)                                                                                                                                                                                     cGTTTATATGCGCAACTGTCGCCGGACATGGGCGGGCCGGAAAAAATCATGAACCCTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCAc                                                  <  1:381082/238‑1 (MQ=255)                                                                                                                                                                                       tgtATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCCTCGTAGCAACGTACCAAACACCAGCAAACATAAATTAACCTGGTTAACTTATCCATCGTTATCGCCGGTGTGTAAATATCATAACCGGGTGAGGTATCAGGTa                                                                  >  2:699457/1‑220 (MQ=255)                                                                                                                                                                                          atatTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGAGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTAc                                                                                                                                                                           <  1:597497/112‑1 (MQ=255)                                                                                                                                                                                                 cctaTTGGGCCGTTAAGGGGCGGGCCGGAAAAAAACAGGAACCCTAAGGGAGCGGGGGATGGGTGATTGGCCGCGTTGCTGCATGCCATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAAcgcg                                     <  2:126607/236‑1 (MQ=255)                                                                                                                                                                                                 cacaTTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTAc                                                                 >  2:11004/1‑211 (MQ=255)                                                                                                                                                                                                 cacaTTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTAc                                                                 <  1:11004/211‑1 (MQ=255)                                                                                                                                                                                                   caTTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGAAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGGGTAAATATGCAAACCGTGTGGTCATTCAGGGATTGAACCACCATTCACCCTCTTTAAcgccg                                    >  2:772753/1‑236 (MQ=255)                                                                                                                                                                                                    aTTGCGCCGTACATGGGCGGGCCGGAAAAAAGCATGAACACTAATGGAGCGGGGGCTGGGGCATTGGCAGCGTTGGTGCATGACATTACCGCCAACAGATACCATCGTAGCa                                                                                                                                                                 >  2:597497/1‑112 (MQ=255)                                                                                                                                                                                                        cgcCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCtt                               >  1:356614/1‑238 (MQ=255)                                                                                                                                                                                                          ccGTACATGGGCGGGCCGGAAAAAATCATGCCCACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTTCTGCATGACCTTACCGCCAACAGCTACCATCGTCGCAACGTACCACACTCCCGCACACATAAATTCACCTGGTTAACTTATTCATCTTTAGCGCTGGTGCTTCACTTGTCTAACCCTGTGACTTATCCGGTTCTACTCCTGCCTTCCCCCCGTTTCACAGGCGTCTTGcc                            >  2:791433/1‑239 (MQ=255)                                                                                                                                                                                                            gttcctGGGGGGGCCGGAAAAAATCAGGAACCCTAATGGGGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCgg                          <  2:644857/234‑1 (MQ=255)                                                                                                                                                                                                            gTACATGGGCGGGCCGGAAAAAAGCATGAACCCTAATGGAGCGGGGGATGGCGCATTGGCGGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCgg                          <  1:157743/239‑1 (MQ=255)                                                                                                                                                                                                            gTACAGGGGGGGGCCGAAAAAAATCATGAACCCAAAGGGAGCGGGGGAGGGCGCATTGGCGGGGTGCCTGCATGCCATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCgg                          <  2:445088/239‑1 (MQ=255)                                                                                                                                                                                                                  gggcgggcCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTg                    >  1:499284/1‑239 (MQ=255)                                                                                                                                                                                                                    gcgggcCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCgggactggaa                  >  1:385762/1‑231 (MQ=255)                                                                                                                                                                                                                      gggcCGGAAAAAATCATGAACATTAAGGGAGCGGGGGATGGCTCATTGGCAGCGTTTCTGCATTACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCCAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGa                <  2:674046/239‑1 (MQ=255)                                                                                                                                                                                                                       ggcCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCCAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGGGAAGAc               >  1:674046/1‑239 (MQ=255)                                                                                                                                                                                                                       ggcCGGAAAAAACCATGAACCCTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGa                <  1:474602/238‑1 (MQ=255)                                                                                                                                                                                                                         ccGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGGGCTATCAGGAACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTGCCCGGAGCGGAAAGACa              >  2:474602/1‑238 (MQ=255)                                                                                                                                                                                                                              aaaaaTCAAGAACCCTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTgg        <  1:218080/239‑1 (MQ=255)                                                                                                                                                                                                                                aaaTCATGAACCCTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa      <  1:472633/239‑1 (MQ=255)                                                                                                                                                                                                                                aaaTCATGAACACTAATGGGGCGGGGGAAGGGTGATTGGAAGCGTTGCTGAATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa      <  1:18150/239‑1 (MQ=255)                                                                                                                                                                                                                                aaaTCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa      <  1:491788/239‑1 (MQ=255)                                                                                                                                                                                                                                    ctggAACCCTACTGGGGCGGGGGTTGGCGCTTTGGCAGCGTTGCTGGATGACATTACCGCCAACAGCTACCAACGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGt  <  2:866228/236‑1 (MQ=255)                                                                                                                                                                                                                                    cATGAACCATAATGGAGGGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGt  <  1:239332/239‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        GGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCC‑ATGCGCCACAAGG‑ATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAACGTACCAAACTCCAGCAAACATAAATTCACCTGGTTAACTTATTCATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGT  >  NC_000913/500949‑501412

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