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	96.9%	82.5 / NA	32	I288I (ATC→ATT)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/19);  minor base G (0/1);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGT  >  NC_000913/500743‑501214                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                     gAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCtt                                                                                                                                                                                                                                 <  1:263765/251‑1 (MQ=255)          aTGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGc                                                                                                                                                                                                                        <  1:103075/251‑1 (MQ=255)              aaGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTTACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGc                                                                                                                                                                                                                    <  1:312461/251‑1 (MQ=255)                gCCATTGAGTACGCGAAGAAATATAACGTACCTGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCTCGTCTCTTTCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGTGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCtt                                                                                                                                                                                                                  <  2:107412/251‑1 (MQ=255)                   aTTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTAc                                                                                                                                                                                                               >  1:7921/1‑251 (MQ=255)                         tACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCgcagg                                                                                                                                                                                                         >  1:372189/1‑247 (MQ=255)                              gaagaaATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGCCTGGGTAGATCGGGTGCTGTGCACCGCCGGGCCAATTGGCTTGGATATTGCggggctttccagcgaagaa                                                                                                                                                                                                     >  2:144006/1‑234 (MQ=255)                                     tcaaaCGTACCGGGGGTGCTGCCGCTGGGCACCAAGTTGGGCATTACCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCTCGTCTCTAACCTTGCGATGAACGAAGAGG‑AAGCAGAAGCGTTGACCGGATAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGTTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTTTATCTGGCGGGCTTTACCGAAGACGAAGTGAAAc                                                                                                                                                                                              <  2:356467/247‑1 (MQ=255)                                          agTACCCGTGGTGCTGACGCTGAGCACAAAGTTTGGCATTGACGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGGTGAACGAGGATG‑AAGCCGAAGCGATGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGTCTGGGTAGCACTGGTGCTGTGCACCGCCGGGCCGATTGCCTTGTATATGGCGGACTTTAAcgacgacgagtcggaatatc                                                                                                                                                                                           >  2:130662/2‑237 (MQ=255)                                             atcGGTGGTGCTGACGCGGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAAcc                                                                                                                                                                                      <  2:252237/248‑1 (MQ=255)                                                        tGACGATGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCgctg                                                                                                                                                                         >  1:130572/1‑251 (MQ=255)                                                              tGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑Cg                                                                                                                                                                  >  1:382362/1‑251 (MQ=255)                                                                      aaGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGAc                                                                                                                                                                                                        <  1:141690/207‑1 (MQ=255)                                                                        gTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTAGCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGGAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGTTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAg                                                                                                                                                        <  1:50834/251‑1 (MQ=255)                                                                           tGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACAGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGACGGGCCACTTGGCTTGTATATGGCGGACTTTAACGAAGAAGACACGAAACCTAAAACCCAG‑CCTCCCCCGCTGC‑CGGGCGCtattgggg                                                                                                                                                     >  2:245563/1‑246 (MQ=255)                                                                            gTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCgg                                                                                                                                                     <  1:144006/250‑1 (MQ=255)                                                                              cATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGaa                                                                                                                                                   <  1:6835/250‑1 (MQ=255)                                                                                 tGCCGAGAATCCGAAGTGTTGGCATCAATGCCTCAAGGATTACGTCTCTCTCCTTGCTAGGAACGATG‑TGTAAGCCGGAGCGTTGATCGGAGTAAGCGATCCGTGGGTGGCATCTGA‑CAAGGCGCTGGACTGTGTATATCTGGTGCTGTTCAGCGCCGGGCCAATTAGCTGGTATATGTCGGGCTTCACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCTGCTGTTGC‑CGGCAGTTATAGCGGAATTc                                                                                                                                                <  2:313410/250‑1 (MQ=255)                                                                                              cAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGAGGAATTCAACCAGTATGAGtt                                                                                                                                  >  1:160447/1‑251 (MQ=255)                                                                                                      gcagcaATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCcgcg                                                                                                                          <  1:206751/251‑1 (MQ=255)                                                                                                         gcaATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGcc                                                                                                                        >  1:167823/1‑250 (MQ=255)                                                                                                           aaTTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑ACGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTTTTGGCATCTGA‑CAAGTCGCTGTAATGGGTAGATCTTGAGCTGTGCACCGCCGGGCCAATTGGCTTGGCTATGGCGGGGTTTAACGAAGACGCAACGAAACGTAAAACC‑AGACATCC‑CCGGTGCACGGGCGCTATAGTGGAAATCCACCAGtgagggtttaaccccagcctt                                                                                                                     >  2:40192/1‑231 (MQ=255)                                                                                                                  cAAAGACCACGTCTCTATCCTTGCGATGAACTAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCTGCTGCTGC‑CGGGCGCTATAGCTGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCca                                                                                                                <  1:72290/249‑1 (MQ=255)                                                                                                                            gTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCa                                                                                                    <  1:76929/251‑1 (MQ=255)                                                                                                                            gTCTCTATCCTGGCGAGGAACGAAGAGG‑AAGCCGAAGCGTTGACCGGAGAAAGCGTTCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGAGCTGTTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCGCAAGGATTGCCa                                                                                                    <  2:382362/251‑1 (MQ=255)                                                                                                                                                  aaGATG‑AAGCCGAAGAGTAGACCGGAGAAAGCGATCCGGTGATGCAATCTGA‑AAAGGCGCGGGAATGGGTAGATCTGGGGCTATGCACCGACGGGCCAATTGGCTTGTATATGGCGTGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCg                                                                              <  2:106600/251‑1 (MQ=255)                                                                                                                                                              gAAGCGTTGACCGGAGAAAGCGATCTGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGCCGAAGCGAAACGTAAAACCCGG‑CATCCGCTGCTGC‑TGGGCGCTATAGCGGAATTCa                                                                                                                                               <  2:191602/175‑1 (MQ=255)                                                                                                                                                              gAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGCATTCa                                                                                                                                               >  1:191602/1‑175 (MQ=255)                                                                                                                                                                      gacgggggaaagtgtgccggggttgGCATATGA‑CAAGGGGCAGGACTGTGTATATCTGGTGCTGTGCCCCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTACCGTAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAAGTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATgg                                                           <  2:247080/230‑1 (MQ=255)                                                                                                                                                                       gtcGGAGAAAGCGACGTGTTGTTGGC‑TCCGAACAAGCCGTTGGATTGGGGAGATTTTGTGCTGTGCACCGCCGGTCCAATTGGCTTGTATAGGGCGGGCCTTACCGAAGACGATGCGAATCGTAAAACCCTG‑CCTCCGCTGCTGC‑CGGGCTCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTTTATATTCGCACATTGCGCCGGTCATggg                                                          <  2:124839/249‑1 (MQ=255)                                                                                                                                                                                                     ggagaaGGCGCTGGACGGGGTAGTTCTGGTTCTGTGCACCGCCGGGCCAATGTGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATg                             <  2:143693/247‑1 (MQ=255)                                                                                                                                                                                                                                ggTGCTGTGCACCGCCGGGCCAATTGGCTTGTATATGGCGGGCTTTAACGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGATTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATAATGAAAACTCATGGAGCAGGGGATGGCGCATTGGCAGCGt  >  1:152967/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                     GAAGCAACCATGAAAGCCATTGAGTACGCGAAGAAATATAACGTACCGGTGGTGCTGACGCTGGGCACCAAGTTTGTCATTGCCGAGAATCCGCAGTGGTGGCAGCAATTCCTCAAAGATCACGTCTCTATCCTTGCGATGAACGAAGATG‑AAGCCGAAGCGTTGACCGGAGAAAGCGATCCGTTGTTGGCATCTGA‑CAAGGCGCTGGACTGGGTAGATCTGGTGCTGTGCACCGCCGGGCCAATCGGCTTGTATATGGCGGGCTTTACCGAAGACGAAGCGAAACGTAAAACCCAG‑CATCCGCTGCTGC‑CGGGCGCTATAGCGGAATTCAACCAGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGT  >  NC_000913/500743‑501214

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