Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,847,969	C→T	G14R (GGA→AGA)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,847,969	0	C	T	96.5%	88.1 / ‑5.9	29	G14R (GGA→AGA)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/11);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGA  >  NC_000913/1847732‑1848206                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              aGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAGGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGGTTCCAACACTTTTGGGGAAATTTTACAGCCTCAACCAGCTCCTTGTCTGTATTgg                                                                                                                                                                                                                                  >  1:462331/1‑251 (MQ=255)             tGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATc                                                                                                                                                                                                                      >  1:93577/1‑251 (MQ=255)                                   gcaaaaaACGAAGGGATTATCAACGAGCGGCAAAAAGAAGGCGGTGGTACTGAAAACGCTGGTGCCTTTGCCCAGATCGTCCACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACAtct                                                                                                                                                                                                <  1:126235/248‑1 (MQ=255)                                             aaaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATg                                                                                                                                                                                      >  1:124655/1‑251 (MQ=255)                                                   ttATAAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATAGTCGCGGGTTTCATTCCCCACAAGCAAATTCGGATCCACAAACTTCGCCTGCTCACTATGCAGGCTGGTTTCCAACACTTTTGGGGACATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTTGGTAAAAAGAATCGGGTTCTGGCTCATTGACATCTCCTTACAATGCAATcc                                                                                                                                                                                >  2:28192/1‑251 (MQ=255)                                                               ggCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGGATTGGGTCAACCGAAACGAGTTCACGCTCATGGACATCTCCTTTCAATGCAATCAGGGTATGGTAAc                                                                                                                                                                    >  2:356432/1‑251 (MQ=255)                                                                      agaagTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAACTTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCAtt                                                                                                                                                             >  1:518519/1‑251 (MQ=255)                                                                       gaagTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTc                                                                                                                                                            >  1:157004/1‑251 (MQ=255)                                                                        aaagcGGTGGGACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCACGGGTTTCATTACCCACAAGAAAATTCGAATCAACAAACTTCGCCTGCACACTATGCAGGATGGTTTCCAACACTTTTGGGGACATTTTAAAGCAGCAACCAGCTCTGTGGCTGTAATGGGTCAAAAGAATCGAGTGACCGCCCATAGACATATACTGTCAAAGCAAACCGGGTATGGTAAACCTCAttgc                                                                                                                                                           >  2:15621/5‑249 (MQ=255)                                                                         agTCGGTGGTACTGATAACGTGGGGGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCg                                                                                                                                                          <  1:365191/251‑1 (MQ=255)                                                                              gtggtACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGACATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGGATGGTAACCCTCATTCCGTGa                                                                                                                                                       >  2:194413/1‑249 (MQ=255)                                                                                        tAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGCAGTGATAAGTg                                                                                                                                            >  1:513407/1‑250 (MQ=255)                                                                                              tGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAACCTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGACTCGAGTTCTCGCTCATGGCAATCTCCTGGCAATGCAATCAGGGTATGGTAACCCACATTCAGAGAAGTGATAAGAGAGAgtgt                                                                                                                                     >  2:78363/1‑251 (MQ=255)                                                                                                                cGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCg                                                                                                                                                                               <  1:331166/191‑1 (MQ=255)                                                                                                                cGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCg                                                                                                                                                                               >  2:331166/1‑191 (MQ=255)                                                                                                                                    cgGGATTCATTACCCACAAGCACATTCGTAGCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCATAACTTTTGGGGAAATTTTACAGCCGCAACCCGCTCTGTGGCTGAATTCGGTCAACCGAATCGAGTTCTCGCTCATGGAACTCTCCTGTCAATGCAATCCGGGCATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGAcgtttt                                                                                               >  2:89085/1‑246 (MQ=255)                                                                                                                                      ggTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTcg                                                                                             <  1:194413/251‑1 (MQ=255)                                                                                                                                               ccccACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTa                                                                                    <  1:384566/250‑1 (MQ=255)                                                                                                                                                              cggctCAACAAACTGCTCCCGTTCACTATGCAGGGTGGTTTCCAACTCTTTAGTGAAAATTGTACATCCGCAGCCAGCGCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCGCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCAGTCCGTGAAGTGATAAGTGAGAGTGTCTGAATGCCTGCGCTCTTGCTCACAAGCCAGTCATTTTCGCAGCACTTAACAGAAATAAGTTAc                                                                     <  2:350009/247‑1 (MQ=255)                                                                                                                                                                   aaacaaACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAAGCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACg                                                                <  2:504072/250‑1 (MQ=255)                                                                                                                                                                    aacaaACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACgg                                                               <  1:356432/251‑1 (MQ=255)                                                                                                                                                                        aaCTTCGCCTGCACACTATGCAGGCTGGTGTCAAAAACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGTGCAAACGAATCGAGTTACCCCGCCTGAACATCGTCAGTCAATGCAATCCGGGCAAGGTAACACTCATTCAGTGAAGTGATAAGCGAGAGTGTCTGAGTTCCTTCGCCTATGCTCCCACTCCAGCCAGATTCGATACAATTATCAGAAATAAGTGACAAACgactac                                                           >  2:410529/1‑246 (MQ=255)                                                                                                                                                                                cTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTCTCAGAAATAAGTTACAAACGGCGTCGGGTCCgg                                                   >  1:257995/1‑251 (MQ=255)                                                                                                                                                                                            cAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTCTCAGAAATCAGTTACAAACTGCGTCGGGTCCGGGACGTTAATCgc                                       >  2:448898/1‑250 (MQ=255)                                                                                                                                                                                                          acacTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTTTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCa                         <  1:332935/251‑1 (MQ=255)                                                                                                                                                                                                                  ggggAAATTTTTCAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTgcggc                 <  2:157004/251‑1 (MQ=255)                                                                                                                                                                                                                            tACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGATATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTCTCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGGCGTTAATCGAAGTAGATGCTTTCAGCTGCGGCGtcccaaggg       >  2:87930/1‑243 (MQ=255)                                                                                                                                                                                                                                gCCGCAACCAGCTCTGTGGATGGCTTGCGTTAAACGAATCGCGTTATCGCTCATGGACATCTCCTGTAAATGCAATCGGGGTAAGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCATGCGCCTTTCCTCACAATACAGACAGTTTCGCGACCATTATCAGAATTCAGTTACAAACGGAGTAGGGCACGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGGGGTaa     >  2:67408/1‑248 (MQ=255)                                                                                                                                                                                                                                 ccGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTGCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAgaga  <  1:87930/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              AGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGA  >  NC_000913/1847732‑1848206

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