Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,722,270	T→A	I773N (ATC→AAC)	pka →	protein lysine acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,722,270	0	T	A	98.2%	184.1 / NA	54	I773N (ATC→AAC)	pka	protein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (26/27);  minor base C (1/0);  total (27/27)
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

CGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGT  >  NC_000913/2722125‑2722413                                                                                                                                                  |                                                                                                                                                cGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGa                                                                                                                                              <  7:113661/149‑1 (MQ=255) cGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGa                                                                                                                                              <  3:45405/149‑1 (MQ=255)  gctgctTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGAc                                                                                                                                             <  8:122183/149‑1 (MQ=255)     gctTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTAc                                                                                                                                          >  8:54113/1‑149 (MQ=255)      ctTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTAc                                                                                                                                          <  5:117109/148‑1 (MQ=255)                      ttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGt                                                                                                                          <  1:11872/148‑1 (MQ=255)                      ttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGAGCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGt                                                                                                                          <  4:34858/148‑1 (MQ=255)                          cAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTg                                                                                                                       >  2:17085/1‑147 (MQ=255)                            gaagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGc                                                                                                                   <  7:73693/149‑1 (MQ=255)                             aagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTAATACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCg                                                                                                                  <  1:63466/149‑1 (MQ=255)                              agaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCgg                                                                                                                 <  3:63282/149‑1 (MQ=255)                                    gCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTacgac                                                                                                           >  6:167363/1‑149 (MQ=255)                                       acaACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTa                                                                                                        <  2:91763/149‑1 (MQ=255)                                       acaACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTa                                                                                                        >  1:90796/1‑149 (MQ=255)                                            ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAg                                                                                                                    >  4:4873/1‑132 (MQ=255)                                            ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAg                                                                                                                    <  3:4873/132‑1 (MQ=255)                                            ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGAt                                                                                                   >  7:92919/1‑149 (MQ=255)                                            ttCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGAt                                                                                                   >  6:131014/1‑149 (MQ=255)                                                 cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTAt                                                                                                       >  1:76438/1‑140 (MQ=255)                                                 cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTAt                                                                                                       <  2:76438/140‑1 (MQ=255)                                                 cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAc                                                                                              >  6:13072/1‑149 (MQ=255)                                                 cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAc                                                                                              <  5:131014/149‑1 (MQ=255)                                                 cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAc                                                                                              >  7:116658/1‑149 (MQ=255)                                                 cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAc                                                                                              <  7:26514/149‑1 (MQ=255)                                                   aTTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACg                                                                                             >  3:61050/1‑148 (MQ=255)                                                          tCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAgaga                                                                                      <  4:120975/148‑1 (MQ=255)                                                              gAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTc                                                                                 >  3:110768/1‑149 (MQ=255)                                                                   aCCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGt                                                                            >  6:112321/1‑149 (MQ=255)                                                                   aCCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTCGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGt                                                                            >  7:5414/1‑149 (MQ=255)                                                                   aCCAAAGAAGATATTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGCGAAATCGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGt                                                                            >  7:122858/1‑149 (MQ=255)                                                                     cAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa                                                                          >  7:53721/1‑149 (MQ=255)                                                                      aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAc                                                                         <  4:110768/149‑1 (MQ=255)                                                                       aagaagATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACCCAGACCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                        >  6:42809/1‑149 (MQ=255)                                                                           agaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTg                                                                    <  5:13072/149‑1 (MQ=255)                                                                              tCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCg                                                                             <  5:102056/137‑1 (MQ=255)                                                                              tCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCg                                                                             >  6:102056/1‑137 (MQ=255)                                                                                   attaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTcc                                                            >  8:123386/1‑149 (MQ=255)                                                                                        cGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGACAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTAGTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATcc                                                       >  8:124948/1‑149 (MQ=255)                                                                                             cTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGa                                                    <  8:116658/147‑1 (MQ=255)                                                                                             cTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGAt                                                   >  2:147588/1‑148 (MQ=255)                                                                                                      gaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATg                                         >  4:17316/1‑149 (MQ=255)                                                                                                               aTTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                <  8:92919/149‑1 (MQ=255)                                                                                                               aTTTACCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                <  1:147588/149‑1 (MQ=255)                                                                                                                   aCCCATGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTg                             <  5:112321/148‑1 (MQ=255)                                                                                                                       aTGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAt                                   <  3:114273/138‑1 (MQ=255)                                                                                                                       aTGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTg                        >  6:100885/1‑149 (MQ=255)                                                                                                                       aTGAAGATTTAGCCAACATGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATAGATGCCGAAt                                   >  4:114273/1‑138 (MQ=255)                                                                                                                                         tGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa      <  5:100885/149‑1 (MQ=255)                                                                                                                                         tGACACAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa      <  4:61050/149‑1 (MQ=255)                                                                                                                                         tGACACAGAACGACTACGATCGGGAAATGGCGTTCGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCaaa      <  8:5414/149‑1 (MQ=255)                                                                                                                                           acacAGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAgg    >  6:26505/1‑149 (MQ=255)                                                                                                                                            cacaGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAggg   <  7:54113/149‑1 (MQ=255)                                                                                                                                             acaGAACGACTACGATC‑GGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGt  >  2:27125/1‑148 (MQ=255)                                                                                                                                               aGAACGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGt  >  2:4973/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                CGCTGCTTGTTCCGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGT  >  NC_000913/2722125‑2722413

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