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	98.2%	175.1 / NA	56	I288I (ATC→ATT)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (25/30);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACA  >  NC_000913/500743‑501226                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 gAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGGGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCtt                                                                                                                                                                                                                                             >  1:174765/1‑251 (MQ=255)       aCCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTatat                                                                                                                                                                                                                                       <  1:42062/251‑1 (MQ=255)                        gTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTTCCGa                                                                                                                                                                                                                        >  2:115484/1‑249 (MQ=255)                          aCGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGAc                                                                                                                                                                                                                    <  1:131951/251‑1 (MQ=255)                          aCGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAAACCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGAc                                                                                                                                                                                                                    <  1:345316/251‑1 (MQ=255)                            gcgAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGATGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACg                                                                                                                                                                                                                   >  1:168920/1‑250 (MQ=255)                                agaaATATAACGTACCGGTGGTGCTGACTCTGGGCACCAAGGTGGTCATTTCCGAGAAGCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCg                                                                                                                                                                                                              <  2:371676/251‑1 (MQ=255)                                  aaaTATAACGTACCGGTGGTGCTGACGTTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGaa                                                                                                                                                                                                            <  2:324832/251‑1 (MQ=255)                                      ataACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAGGGCGCTGGACTGGGTAGATCTGGTGCTGTGCAGCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACg                                                                                                                                                                                                         <  1:156222/250‑1 (MQ=255)                                        aaCGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑a                                                                                                                                                                                                      >  1:225581/1‑250 (MQ=255)                                             aCCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGAAGATGCGAAAAGT‑AAAA‑ccc                                                                                                                                                                                               >  1:132272/1‑251 (MQ=255)                                              ccGGTGGTGCTGACGCTGGGCACCAAGTTTGGCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCGTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAATGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGTCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCa                                                                                                                                                                                              <  1:205197/251‑1 (MQ=255)                                               cGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGGAGACGAAGCGCAACGT‑AAAA‑ccc                                                                                                                                                                                               >  1:55330/1‑249 (MQ=255)                                                ggtggaGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGCCTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTGTATGGCGGGCTTTACCTGAGCCGAAGGGAAACGTGAAAA‑ccccc                                                                                                                                                                                             >  2:361616/1‑249 (MQ=255)                                                         gACGCTGGGCACCAAGTTTGTCAGTGCCGAGAACCCGGAGTGTGGGTAGCAATTACTCAAAGACCACGTCGCGATCCTTTTGATGACAGAAGATGAAGCCGATGCGTTGACCGGAGAAAGCGACCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTTCTGTGCACCGCCGGGCTAATTGGCTTGTATATTGCGGGATTTACCGAAGGCTAAGCGAATCGT‑AAAA‑CCCAGCATCCgct                                                                                                                                                                                     <  2:331055/249‑1 (MQ=255)                                                            gCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGCAGCGAAACGT‑AAAA‑CCCCGCATCCGCtgctgc                                                                                                                                                                                >  1:111679/1‑251 (MQ=255)                                                               gggCCCCAAGTTTGTCTTTGCCGAGAATCCGCATTGGTGGCAGCAATTCCGCAAAGATCCCGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTGACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCg                                                                                                                                                                              <  2:55330/250‑1 (MQ=255)                                                                   aCCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCACTTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGGAGACGAAGCGAAACGT‑AAAA‑CCCAGCATcc                                                                                                                                                                                        >  1:400718/1‑236 (MQ=255)                                                                   aCCAAGTGTGTCATGGCCGAGAATCCACAGTGGGGGCAGCAATTCCACAAAGATCACGTCTCTAACCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATcc                                                                                                                                                                                        <  2:400718/236‑1 (MQ=255)                                                                    ccAAGTTTGTCATTGCCGAGAACCCGCAGGGGTGGCAGCAATTACTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCt                                                                                                                                                                        <  2:135149/251‑1 (MQ=255)                                                                       aGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAAACGAAGCGAAACGT‑AAAA‑CCCCACAACCCCTGCTGCCGGGCGCtatt                                                                                                                                                                     >  2:255444/1‑250 (MQ=255)                                                                                 tGCCGAGAAGCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCGGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTTTAGCGGAAtgc                                                                                                                                                            <  2:90934/250‑3 (MQ=255)                                                                                         aTCCGAAGGGGTGGCTGCAATTCCTCAAATATCACGTCTCTATCCTGGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTAt                                                                                                                                                   <  2:481583/251‑1 (MQ=255)                                                                                            cGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGCAGCAGCGAAACCT‑AAAA‑CCCAGACTCCGCTGCTGCCGGGCGCTATAACGGAATTCAACCAGTATg                                                                                                                                                  >  2:402819/1‑249 (MQ=255)                                                                                            cGCAGTGGTGGCAGCAATTCCGCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAg                                                                                                                                                <  2:106495/251‑1 (MQ=255)                                                                                                          caATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTTCCGAAGACGAAGCGAAACGT‑AAAAACCCA‑CATCCGCTGCTTCCGGGCGCCTTTGCGGAATTCAACCAGTATGCGttttgattggaagg                                                                                                                                  >  2:159255/1‑240 (MQ=255)                                                                                                          caATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCAGCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGttt                                                                                                                                             <  1:159255/240‑1 (MQ=255)                                                                                                             ttCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGA‑AAAA‑ACCAGCATCCGCCGCTGCCGGGCGCTATAGCGGGATTCACCCAGTATGAGTTTAACCGCGCCCTgcg                                                                                                                               >  2:449261/1‑251 (MQ=255)                                                                                                                               tctATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGa                                                                                                              >  1:348282/1‑250 (MQ=255)                                                                                                                               tctATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGa                                                                                                              >  1:258139/1‑250 (MQ=255)                                                                                                                                ctATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGGTTTACCGTAGACGAAGCGAAACGT‑AAAA‑CCCAGCCTCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGGGGTTTGACGCGCCCTGGCCCACAAGGGTTGCCAGaa                                                                                                             >  2:470104/1‑250 (MQ=255)                                                                                                                                        gCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAATCGTGCCGTTTTTGGCATCTGACTAGGCGCTGGACTGGGGAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGc                                                                                                     <  1:56107/250‑1 (MQ=255)                                                                                                                                                  aaGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCg                                                                                          <  1:77020/251‑1 (MQ=255)                                                                                                                                                       gAAGACGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCCTCCGCAGCTGCCGGGCGCTACAGCCGAATTTAACCAGTATGAGTTTAGGCGCGCCAGGCGCCCCAAGGCTTGCCCGAATCCGccgggtgttgtttccccact                                                                                     >  2:365737/1‑232 (MQ=255)                                                                                                                                                        aaGCCGAAGCGTTGGCCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGTCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTTCGTGTATATTCGCACAtt                                                                                    <  1:106186/251‑1 (MQ=255)                                                                                                                                                        aaGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACAtt                                                                                    >  1:215877/1‑251 (MQ=255)                                                                                                                                                           ccGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTgcg                                                                                 >  1:149842/1‑251 (MQ=255)                                                                                                                                                               ggCGTTGCCCGGAGAAGACTATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTTTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGACGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCTCATTGCGCCGt                                                                             <  2:435233/250‑1 (MQ=255)                                                                                                                                                               aGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGCAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGt                                                                             >  1:70159/1‑251 (MQ=255)                                                                                                                                                                  gTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCCCAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTAcc                                                                          >  1:52944/1‑250 (MQ=255)                                                                                                                                                                        cGGAGAAAGCGATCCGGTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATgggcg                                                                    <  1:88407/251‑1 (MQ=255)                                                                                                                                                                         ggAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGAtt                                                                                                                    <  1:289331/202‑1 (MQ=255)                                                                                                                                                                           agaAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATgggcgggc                                                                 <  1:200273/251‑1 (MQ=255)                                                                                                                                                                             aaaGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCg                                                               >  1:146728/1‑251 (MQ=255)                                                                                                                                                                                             ggCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAAAGT‑AAAA‑ACCAGCATCCGCTGCCGCCGGGCGCTATAGCGGAATTCAAACAGTATTAGGTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCCCCGCGTGTATATTCCCACCTTGCGCCGGACAAGGGAGGGCCCGAAAACATCATGAac                                                >  2:2839/1‑250 (MQ=255)                                                                                                                                                                                                agcTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTa                                            <  2:149842/249‑1 (MQ=255)                                                                                                                                                                                                    gACAAGGCGCTGGACTGGGTAGAGCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTa                                            <  2:70159/247‑1 (MQ=255)                                                                                                                                                                                                      cAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGGAACGT‑AAAA‑ACCAGCATCCGCTGCTGCCGGGCGCTATAGCGGCATTCCACCAGTATGAGTTTAGCCGCGCCATGCGCCCCAAAGATTGCCAGAAACCGCCGCGTGTATATTCGCACATTGCGCCGTCAATGGGCGGGCCGGGAAAAATCAAGAAACCTAATgggg                                      >  2:484839/1‑249 (MQ=255)                                                                                                                                                                                                        aGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCg                                    <  1:131855/251‑1 (MQ=255)                                                                                                                                                                                                                 aCTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGgc                           <  1:166023/251‑1 (MQ=255)                                                                                                                                                                                                                    gggTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGgcgc                         >  1:148507/1‑250 (MQ=255)                                                                                                                                                                                                                          aTCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTAACGGAGACGAAGCGAAACGT‑ACAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTAAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGGTTGCCAGAATCCGCCGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGCACACTAATGGCGCGGGGGAGGGCGCATTGGcc                  >  2:37313/1‑250 (MQ=255)                                                                                                                                                                                                                              ggggCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGc                <  1:362531/246‑1 (MQ=255)                                                                                                                                                                                                                              ggTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGc                   <  1:449261/246‑1 (MQ=255)                                                                                                                                                                                                                                          cgccgGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACa  <  1:115484/251‑1 (MQ=255)                                                                                                                                                                                                                                          cgccgGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACa  <  1:2839/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 GAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGT‑AAAA‑CCCAGCATCCGCTGCTGCCGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACA  >  NC_000913/500743‑501226

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