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	95.8%	136.2 / ‑2.2	48	I288I (ATC→ATT)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (21/25);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

AACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACC  >  NC_000913/500748‑501231                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      aaCCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGtatt                                                                                                                                                                                                                                             >  1:183277/1‑250 (MQ=255)          aGCCATTGAGTACGCGAAGAAATATAAAGAACCGGCGGTGCTGACGCTGGGGACCAAGTTTATCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAGCGAAGATGAAGCCGCAGCGTTGACCGGAGAAAACGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGGGCTGTGCACCGCCGGGCCAATTGACTAGAATATGGCGGGcc                                                                                                                                                                                                                                    >  1:331579/1‑250 (MQ=255)                   gTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGGttttccaca                                                                                                                                                                                                                            >  2:222024/1‑244 (MQ=255)                     aCGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTGGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGa                                                                                                                                                                                                                          <  1:325640/250‑1 (MQ=255)                         gaagaaAAATAACGTGCCTGTGGTGCTGACTGTTGGTGCCAAGCTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATCCCACAAAGATCACGTCTCTATCCTTTCGATGAACGAAGATGAAGCCGAAGCGTTTACCGGAGAAAGCGATCCGTTGTTGGTACCTGACAAGGCGCTGGACTGGGTAGACCTGGCGCCGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAg                                                                                                                                                                                                                     <  1:208185/251‑1 (MQ=255)                            gaaATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGa                                                                                                                                                                                                                  <  1:261838/251‑1 (MQ=255)                                tataACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGCAACg                                                                                                                                                                                                              >  1:65952/1‑251 (MQ=255)                                     cGTACCGGTGGTGATGACGCTGGGCCCCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGGCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTaaaa                                                                                                                                                                                                         <  2:183277/251‑1 (MQ=255)                                                  gtGACGCTGGGCACCAAGTTTGTCATTACCGAGAACCCGCAGTGGTGGCCGCAATGCCTCAAATATCCCGTCTCTCTCTTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGAGCTGGTGCTGTGCACCGCCGGGCCAAGTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAATCCCAGCATCCgct                                                                                                                                                                                            <  2:339531/250‑1 (MQ=255)                                                     aCGCTGGTCACCAAGGTTGTCATTTCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCgctgct                                                                                                                                                                                         <  1:109021/251‑1 (MQ=255)                                                     aCGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCgctgct                                                                                                                                                                                         >  1:270830/1‑251 (MQ=255)                                                               ccAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCt                                                                                                                                                                               <  1:337221/251‑1 (MQ=255)                                                                cAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGGAACGTAAAACCCAGCATCCGCTGCTGCCGGGGGCtt                                                                                                                                                                              >  1:177404/1‑250 (MQ=255)                                                                   gTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTTAAAGATCACGTTTCTATCCTTGCGATGAACGAAGATTAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGCCAAGGCGCTGGACTGGGTAGATCTGTTGCTGTGCACCGCCGGGCCAAGTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTTCCGGGCGCTATAg                                                                                                                                                                           <  2:68957/251‑1 (MQ=255)                                                                   gTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAg                                                                                                                                                                           <  1:106562/251‑1 (MQ=255)                                                                             gccgaaatccctcaggggtggCAGCAATTCCTCAAAGATCACGGCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCaa                                                                                                                                                                <  2:457215/235‑1 (MQ=255)                                                                                  gaATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGt                                                                                                                                                           <  1:222024/251‑1 (MQ=255)                                                                                     cccGCAGTGGGAGAAGCAATTCTTCTAAGAACGCTTCTCTATTCTTGCTAGGATCGAAGATGTAGCCGAAGCGTTGACCGGAGAAATCGAGCCGTTGTTGGCAGCTGAGAAGGCGCTGGACTGGGTATAGATGGTGCTGTGTACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGGTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATg                                                                                                                                                        <  2:94043/250‑1 (MQ=255)                                                                                       cGCAGTGGTGGCAGCAATTCCCCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGAGGAAGCCGAAGCTGTGACCGGAGAAAACGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTTTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAAACCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAg                                                                                                                                                      <  2:298602/251‑1 (MQ=255)                                                                                                 gaatcaATTCCTCAAAGATAACGTCTCTATCCTGGCGAGGAACGAAAATGAAGCCGAAGCGTTGACCGGAGAAAGCTATCCGTTGTTGGCATTGGACATGGCGCGGGACTGGGTAGATCTGGTGCTGTGCGCCGCCGGGCCAAGTGGCTTGTATATGGCGGGCTGTACCGAAGGCGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGT‑AATTCAACCAGTATGAGTTTAGCcggg                                                                                                                                            <  2:20165/247‑3 (MQ=255)                                                                                                      aagtCCTCAAAGACCACGTCTCTATCCTGGCGATGAACGAAGATGAAGCCGAAGCGGTGACCGGAGAAAGCGATCCGTTGGTGGCATCTGACAAGGCGCGGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATg                                                                                                                                       <  2:72247/248‑1 (MQ=255)                                                                                                           cgcAAAGATCACCTCTCTATCCCGTCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAATCGCTCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCCTGTCCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AAGTCAACCAGTATGAGTTTAGCCGCGCCAGGCGCca                                                                                                                                  <  2:331579/249‑1 (MQ=255)                                                                                                                 gATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTAACGGCGACGGAGCGAAACGTAAAACCCAGCATCCGCTGCCGCCGGGCGCGTTAGGGG‑AATTCCACCCGCATGAATTTAACCCGGCCCTTCGCCACAAggg                                                                                                                            >  2:138747/1‑250 (MQ=255)                                                                                                                   taaCGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGCGCTGTGCACCGCCGGGCCAATTGGCTTGTATCTGGCGGGCTTTAACGAAGACGAAGCGAAACGTAAAACCCAGCATCCGATTCTGCCTGGCGCCATAGGGG‑GATTCAACAAGTATGAGATTAGCCGCGCCATGCGCCACAAGGA‑tt                                                                                                                         >  1:20559/3‑251 (MQ=255)                                                                                                                       gTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGcc                                                                                                                      >  1:195123/1‑250 (MQ=255)                                                                                                                       gTCTCGATCCTTGCGATGTACGAAGAGGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCTTTGTTGGCATCTGACAAGGCGCTGGACTGGGTCGATCTGGTTCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCa                                                                                                                     <  2:398760/251‑1 (MQ=255)                                                                                                                               ccTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGc                                                                                                             >  1:298532/1‑251 (MQ=255)                                                                                                                               ccTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGGAAAAACCAGCATCCGCTGCTGCCGGGAGCTATAGCGG‑AATTCAACCAGTATGAGTTTAACCGCGCCATGCGCAGCAAGGA‑TTGCCAGAAACCgg                                                                                                             >  1:156879/1‑250 (MQ=255)                                                                                                                                       tGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAAAGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGAGG‑AATTCAACCAGTATGAGTTTAACCGCGCCATGCGCCCCCAGGA‑TTGACAGAATCCGCTGCgtgt                                                                                                      >  1:466332/1‑250 (MQ=255)                                                                                                                                         aaCGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGGGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGGTATAGCGG‑AATTCAAACAGGATGAGTTTAGCCGCGCCATGCGCCACAAGGG‑ATGCCAGAAACCGCTGCGTGtat                                                                                                    >  1:105577/1‑250 (MQ=255)                                                                                                                                           cGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTAGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACCAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGAGG‑AATTCAACCAGTATGAGGTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATAtt                                                                                                 >  1:238786/1‑251 (MQ=255)                                                                                                                                               gaggAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCGGACAAGGCGCTGGACTGGGGAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGca                                                                                             <  2:298532/248‑1 (MQ=255)                                                                                                                                               gATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGTCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAACCCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGca                                                                                             <  1:96983/251‑1 (MQ=255)                                                                                                                                                 tGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGGAGAGCAACCTAAAAACCAGCAACCGCTGCGGCCTGGGGCTATAGCGG‑AATTCAACCAGGATGGGTTTAGACGCGCCATGCGCCCCAAAGA‑TTTCCAGAATCCCCGTCGTGTATAtttcaccc                                                                                           >  2:408402/1‑245 (MQ=255)                                                                                                                                                    aGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGGAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATACGGCGGGGTTTAACGAAGACGAAACGATACGGAAAACCCCAAAACCGCTGGTGACGGGCGCTATGACGG‑AGTTCACCCAGAATGAGATTAGCCCAGCCATGCGCCCCAAGGAGTTCCCAGA‑GCCGCTGCATGTATATTCGccccatg                                                                                        >  2:232904/1‑245 (MQ=255)                                                                                                                                                             gTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTAcc                                                                               >  1:234256/1‑250 (MQ=255)                                                                                                                                                                    ggAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATgggcgg                                                                        <  1:464461/251‑1 (MQ=255)                                                                                                                                                                      agaAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTTTATTCGCACATTGCGCCGTAAATgggcgggc                                                                      >  1:393967/1‑251 (MQ=255)                                                                                                                                                                       gaAAGCGATCGGTGGGTGGCAGCGGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTAGATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGcc                                                                     <  2:78937/251‑1 (MQ=255)                                                                                                                                                                        aaaGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACAGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGATATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCCTGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATgggcgggacg                                                                    >  1:141879/1‑248 (MQ=255)                                                                                                                                                                                ccGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGACCCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGGAGACGACGCGCAACGTAAAACCCAGCATCCGCTGCCGCCGGGCGCTATTGCGGCATTTCAACC‑GAATGAATTTAGCCACTCCATGCGCCCCAAGGA‑TTGCACGAATACGCTGCATGTATAATCGCACATTGCCCCGTTCATGGGGGGGCCGGaaaaaa                                                             >  2:408603/1‑250 (MQ=255)                                                                                                                                                                                              ggcaaaGGCGCTGGACTGGGTAGATCTGGGGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATg                                              <  2:393967/247‑1 (MQ=255)                                                                                                                                                                                                aCAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCtat                                                                                                                                                                             <  1:127378/124‑1 (MQ=255)                                                                                                                                                                                                aCAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCtat                                                                                                                                                                             >  2:127378/1‑124 (MQ=255)                                                                                                                                                                                                   aGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCg                                         <  1:122081/251‑1 (MQ=255)                                                                                                                                                                                                               gggTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGgcg                               <  1:389901/249‑1 (MQ=255)                                                                                                                                                                                                                          gTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGtt                  <  2:177404/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGGAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAACCGCGCCCTGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCCAACATTGCGCCGTACATGGGCGGGCCGGCAAAAATCCTGAACACTACTGGAGCGGGGGATGGCGAGTGGGCCGCGGTGCTGCATgcaattctc  >  2:408350/1‑243 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      AACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATGAAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGACAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAGCATCCGCTGCTGCCGGGCGCTATAGCGG‑AATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGA‑TTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACC  >  NC_000913/500748‑501231

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