Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	500,988	C→T	I288I (ATC→ATT)	gsk →	inosine/guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	500,988	0	C	T	93.5%	173.6 / 0.3	63	I288I (ATC→ATT)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (20/38);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACC  >  NC_000913/500741‑501231                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    cGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGc                                                                                                                                                                                                                                                  <  1:418212/251‑1 (MQ=255)       cAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGtat                                                                                                                                                                                                                                            <  1:290147/251‑1 (MQ=255)        aaCCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGtata                                                                                                                                                                                                                                           <  1:515778/251‑1 (MQ=255)         aCCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGACATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTatat                                                                                                                                                                                                                                          <  1:34833/251‑1 (MQ=255)          ccATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATg                                                                                                                                                                                                                                         <  1:422566/251‑1 (MQ=255)             tGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGCTCCGTTGTTGGCATCTGACAAGGCGCTGGCCTGGGTAGATCTGGGGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCg                                                                                                                                                                                                                                      >  2:21027/1‑251 (MQ=255)                   ccATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CCGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCttt                                                                                                                                                                                                                                <  1:604179/251‑1 (MQ=255)                    cATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTa                                                                                                                                                                                                                               <  1:20854/251‑1 (MQ=255)                     aTTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTAc                                                                                                                                                                                                                              >  1:294795/1‑251 (MQ=255)                       tGAGTACGCGAAGAAATATAACGTTCCGGGGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTAAAAGATCACGTCGCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTAcc                                                                                                                                                                                                                             <  2:204506/250‑1 (MQ=255)                        gAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCg                                                                                                                                                                                                                            <  1:233867/250‑1 (MQ=255)                         aGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCtagatcggg                                                                                                                                                                                                                          >  1:285095/1‑243 (MQ=255)                         aGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGGGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCt                                                                                                                                                                                                                                  <  2:285095/243‑1 (MQ=255)                                     aaTATAACGTACCGGTGGGGCTGACGCTGGGCACAAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGGTCCCGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCGGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGaaa                                                                                                                                                                                                              <  2:33565/251‑1 (MQ=255)                                      atatAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAAc                                                                                                                                                                                                             <  1:641469/251‑1 (MQ=255)                                       tataACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAACGGAACg                                                                                                                                                                                                            >  1:371674/1‑251 (MQ=255)                                                        gCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCg                                                                                                                                                                                            <  1:555306/250‑1 (MQ=255)                                                          tGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGGCGAAGCGAAACGTAAAACCCAGCATCCgctg                                                                                                                                                                                         >  1:436929/1‑251 (MQ=255)                                                               cTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAAACCAGCATCCGCTGCTGcc                                                                                                                                                                                    >  1:156308/1‑251 (MQ=255)                                                                 gggCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATACGCTGCTGcc                                                                                                                                                                                    <  1:601561/249‑1 (MQ=255)                                                                  ggCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCggg                                                                                                                                                                                 >  1:441725/1‑251 (MQ=255)                                                                   gCACCACGGTTGTAATTGCCGAGAATCCGCAG‑GTTGGGCAGCAATTCCTCAAAGAACACGTCTCTAGTCTGGCAGTGAACGAAGATTAAGCCGAAGCGTGGACCGGAGAAAGCGATACTTTGTTGGCATCTGTCAAGGCGCGGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCATCATCCGCTGCTGCCGGgc                                                                                                                                                                                <  2:431811/251‑1 (MQ=255)                                                                         aggttggtaTTGCCGTGAATCCGCAGGGGTGG‑CAGCAATTCCTAAAAGATCTCGTCTCTATCCTGGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCGGACAAGGCGCTGGCCTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTGACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCtata                                                                                                                                                                          <  2:504008/243‑1 (MQ=255)                                                                                cATTGCCTAGAATCCGCAGGGGTGG‑CAGCAATTCCCCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGtat                                                                                                                                                                                                                                            <  2:532000/178‑1 (MQ=255)                                                                                cATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGtat                                                                                                                                                                                                                                            >  1:532000/1‑178 (MQ=255)                                                                                      cGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAAcc                                                                                                                                                             >  1:312759/1‑251 (MQ=255)                                                                                          aaTCCGCAGTGTTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTa                                                                                                                                                         <  2:134168/251‑1 (MQ=255)                                                                                                               atttcTCAAAGAT‑ACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAGGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATgcg                                                                                                                                    >  1:341543/5‑251 (MQ=255)                                                                                                                   cTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAAACCAGCATCCGCTGCTGCCGGGCGCTATTACGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGccc                                                                                                                                 >  1:420106/1‑250 (MQ=255)                                                                                                                    cccaaGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGACGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGCCTGTGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCcac                                                                                                                                <  1:659668/248‑1 (MQ=255)                                                                                                                          gccACGGCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGGTGACCGGAGAAAGCGATCCGTGGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGAt                                                                                                                          <  2:294795/249‑1 (MQ=255)                                                                                                                          aTCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATTTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTAGGAGTTTAGCCGCTCCATGCGCCACAAGGAt                                                                                                                          <  1:517359/251‑1 (MQ=255)                                                                                                                          aTCACGTCTCTATCCGTGCGATGAACGATGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGAt                                                                                                                          <  1:243803/251‑1 (MQ=255)                                                                                                                               gTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCCCAAGGATTGccc                                                                                                                     >  1:282644/1‑250 (MQ=255)                                                                                                                                   ctATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGa                                                                                                                   <  1:160588/249‑1 (MQ=255)                                                                                                                                             gATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTTCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATTCGACACAAGGATTGCCAGTATCCGcctcgtg                                                                                                       >  1:397786/1‑245 (MQ=255)                                                                                                                                             gATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTTCCAAAATCCGCTGCgtg                                                                                                       >  1:266637/1‑251 (MQ=255)                                                                                                                                                gAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGta                                                                                                     <  1:597922/250‑1 (MQ=255)                                                                                                                                                     aaGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCTTTGTTGGCATCTGACAAGGCGCTTGACTGGGTAGATCTGGTGCTGTGCACTGCCGGGCCAATTTGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCg                                                                                               <  1:410580/251‑1 (MQ=255)                                                                                                                                                           aaGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATTTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACAtt                                                                                         <  1:150314/251‑1 (MQ=255)                                                                                                                                                           aaGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACAtt                                                                                         >  1:482584/1‑251 (MQ=255)                                                                                                                                                             gccgatgcttttaCCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGTCTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGGATAGGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAACCCCAGCATCCGCTGATGCCGTGCGCTATAGCGGAATTCAACCAGTATGAGTTTTGCCGCGCCATGCTCCACAAGGATTGCCAGAGTCCGCTGCGTGTATATTCGCTCATTgc                                                                                       <  2:511953/239‑1 (MQ=255)                                                                                                                                                                  aGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGt                                                                                  <  1:56675/251‑1 (MQ=255)                                                                                                                                                                        gACCGGAGAAAGCGATCCGTTGTTGGCATCTTACAAGGCGCTGGACTGGGTTGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATg                                                                             <  1:15979/250‑1 (MQ=255)                                                                                                                                                                         aCCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATggg                                                                           <  1:553445/251‑1 (MQ=255)                                                                                                                                                                                  aGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTTAACCAGTATGAGTTTAGCCGCGCCCTGCGCCACAAGGATTGCCAGCATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCggg                                                                  >  1:657204/1‑250 (MQ=255)                                                                                                                                                                                   gcgttCCGTTGTTGGTATCATACAAGGCGCTGTATTGGGTATAGCTGTTGCTGTGCACGGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAATACGAAGCGAAACGTAAAACCCAGAATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAGCCAGGATTAGTTTAGTCGCGCCATTCGCCACAATGATTGCCAGAATCCGCGGCGTGTAGATTCGCACATTGCGCCGTACATGGGTGGGTCGGaa                                                                 <  2:629518/247‑1 (MQ=255)                                                                                                                                                                                    cGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGaa                                                                 >  1:637715/1‑250 (MQ=255)                                                                                                                                                                                     gATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAAACCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTTAACCAGTATGAGTTTAGCCGCGCCATGCGCCACCAGGATTGCCAGAATCCGCTGCGTGTCTATTTGCACATTGCGCCGTACAAGGGAGGGCCGGaaa                                                                >  2:448020/1‑250 (MQ=255)                                                                                                                                                                                        ccGTTTTTGGCATCTGACAAGGCGCGGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAAt                                                            <  2:266637/251‑1 (MQ=255)                                                                                                                                                                                            tgttgGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACAt                                                                                          <  1:356965/217‑1 (MQ=255)                                                                                                                                                                                                ggCATCTGACAAGGGGCGGGACTGGGTAGAGCTGGTGCTGTGCACCGCCGGGCCAATGGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAAGCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAac                                                     <  1:669666/250‑1 (MQ=255)                                                                                                                                                                                                 gCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAacac                                                   >  1:105460/1‑251 (MQ=255)                                                                                                                                                                                                       gACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTaa                                                <  1:448020/248‑1 (MQ=255)                                                                                                                                                                                                        aCAAGGCGCTGGCCTGGGTAGGTCTGGTGCTGTGCCCCGCCGGGCCAAACGGCTTGTATATGGCGGGCTTTACCGAGGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGa                                            <  2:341543/251‑1 (MQ=255)                                                                                                                                                                                                         cAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAg                                           <  1:598762/251‑1 (MQ=255)                                                                                                                                                                                                                     tggggTAGATCTGGTGCGGTGCACCGCCGGGCCAATTTGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGgc                                <  2:152569/248‑1 (MQ=255)                                                                                                                                                                                                                      tggttcggttcggtgtggggaaacgccgGGCCAATTGGGTTGTGTAGGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGGATCCGCTGATGCTGGGCGCTATAGCTGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACATGGATTGCCAGAATCCGCTGCGTGTATATTCGCTCATTGCGCCGTACATGGGCGGGCAGGAAAAAATCATGAACACTAATGGAGCGGGGGATGgcg                               <  2:23449/227‑1 (MQ=255)                                                                                                                                                                                                                      tGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGgc                                >  1:152569/1‑249 (MQ=255)                                                                                                                                                                                                                                 ggTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGCAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGAGGGGGAGGGCGCATTGGCAGCGt                   >  1:245466/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggacACCGCCGGGCCAATTGGATTGTATATGGCGGGCGTTACCGAAGAAGAAGCGAAACGTAAAACCCAGCATCCGCGGCTGCCGGGCGCTAGAGCGGAATGCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCTTACATGGGCGGGCCGGAAAAAATCATGAACACAAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgct               <  2:657204/246‑1 (MQ=255)                                                                                                                                                                                                                                          cACCGCCTGGCCAATTGGCTTGTATATGGCTGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGTGCGCTATAGCGGAATTCAACCAGTATGAGTTTATCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCTGAAAAAAGCATGAACGCTGATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATg          <  1:250112/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ggCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGATGCATGACATTAcc  >  1:502665/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    CGGAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGG‑CAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACC  >  NC_000913/500741‑501231

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