Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	125,463	A→C	44.2%	T309P (ACC→CCC)	CJ019_00134 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	125,463	0	A	C	44.2%	4.6 / 4.8	24	T309P (ACC→CCC)	CJ019_00134	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (9/4);  new base C (7/4);  total (16/8)
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

CGCCAGCGAGTCGTTGGCTGAGGCAGAACCCTTACAAGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGA  >  KT2440/125333‑125581                                                                                                                                   |                                                                                                                       cGCCAGCGAGTCGTTGGCTGAGGCAGAACCCTTACAAGCCCTTGAAGAAACCTTTGGTGACGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGAccc                                                                                                                     <  6:30383/134‑1 (MQ=255)            cGTTGGCTGCGGCATAACCCTTACACGCCTTTAAAGAACCCTTTGGTTAAGACGCCCCCTATCCCCTGCCCTACGAAGTGGAGCCGCATAACCGCCAGGTTGCCGCGCACATCTAGCTGCCCGTTATCGTTCTGCt                                                                                                        <  8:336437/136‑1 (MQ=255)                        cAGAACCCTTACAAGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGGTCGCTTACCGCa                                                                                             >  3:69750/1‑135 (MQ=255)                              ccTTACAAGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCCGACCCGGATCGGTCGGCTCCATGACCACAGCCTGc                                                                                       >  2:86506/1‑135 (MQ=255)                              ccTTACAAGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGCCCCTGATCGGTTTGCTCGATGACCTCAGCCTGc                                                                                       >  4:264607/1‑135 (MQ=255)                              ccTTACAAGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGCCCCTGATCGGTCTGCTCGATGACCTCCGCCTGcc                                                                                      >  7:222153/1‑136 (MQ=255)                                    aaGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGCCCCGGATCGGTCTGCTCGATGACCACAGCCAGCCCgcgc                                                                                 >  2:5519/1‑135 (MQ=255)                                    aaGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGAccc                                                                                                                     >  4:80200/1‑99 (MQ=255)                                    aaGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGAccc                                                                                                                     <  3:80200/99‑1 (MQ=255)                                     aGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACAACGCGCCCACCCAGCTCCCGCGGCTCGCTGCCCTCCCCCGGCCCCAGCg                                                                                >  1:133105/1‑135 (MQ=255)                                        ccTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGac                                                                              >  6:206028/1‑134 (MQ=255)                                           tGAAGAAACCTATGGAGAAGACGGCCCCTATGCCCAGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACAAAGACCAAACCCAGCTCGGTCTTCGCGAAGACCAAAGCCAGCCCGCGCCACACaa                                                                          >  3:106043/1‑135 (MQ=255)                                                     tttGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGCCCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGCGCGACACAAGGCTCAGgcg                                                                >  4:355253/1‑135 (MQ=255)                                                     tttGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCCAGCAGACCCTGATCGGTCTGGTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCGCAGgcc                                                                >  1:62413/1‑134 (MQ=255)                                                                                  ctACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATAGAGCAGCCCCGGAACGGGCTGCTCGATGCCCTCCGCCTGCCCCAGCGACACAAGGCTCAGGCGGCGGAAAAGCGCGAGCGGGGTTCTCGCgg                                   >  1:117913/1‑135 (MQ=255)                                                                                     cgctggggAGCCGCCTAACCGCCCGGTGGCCCCCCAAAGCGCCCAGCCCCTGAGCGGTCGGCGCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCtt                                <  2:318658/129‑1 (MQ=255)                                                                                     ccccgTGGAGCCGCCTAACCGCCTGGTGCCCGCCCACATCCCGCACCCCCTTAGGGGGGTGCGCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCtt                                <  5:206028/131‑1 (MQ=255)                                                                                     cccAGGGGAGCCGCCTAACCCCCGGGCCGCCGCCCACCGCGCGCCGCCCCTGCGCGGTGGGCGCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCtt                                <  1:86506/133‑1 (MQ=255)                                                                                      gctgTGGAGCCGCCTTACAGCCTGGTTGCCGCGCACATCCAGCAGACCCTGAGCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCtt                                <  4:69750/131‑1 (MQ=255)                                                                                        aGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATAGCGCAGACCCTGCTCGGGCTGCTCGATGACCTCACCCTGCCCGCGCGCCACAAGGCTCAGGCGCCGGAAATGCGCGAGCGGGTTGCTCGCGGCGTGAAc                            >  7:263272/1‑136 (MQ=255)                                                                                             aGCCGCCTTACAGCCAGGTTGCCGCGCACATCCAGCAGCCCCTGCTCGGTCTGCTAGATGACCGCAGCCTGCCCGCGCGCCCCAAGGCTCAGGCGCAGGCAATGCGCGAGCGGGTTGCTCGCGGCGTGAACgggga                       >  7:68932/1‑134 (MQ=255)                                                                                                       cAGCCAGGTTGCCGCGCACATCGAGCCGACCCGGCTCGGGCGGCTCGCTGGCCGCCGCCAGCCCCCGCGACACCAGGCGCAGGCGCAGGGAAAGCGCGCGCCGGGTGCGGGCGGCGTGAACTGGGATGAAcctgt              >  4:415545/1‑131 (MQ=255)                                                                                                               ttGCCGCGCACATCGCGCAGCCCCTGATCGGGCTGCGCGATGACCTCAGCCTGCCCGCGCGCCACAAGGCTCAGGCGCAGGAAAAGCGCGCGCGGGGTGCTCGCGGCGTGAACTGGTTTGAACTGATCCCGGTGCt     >  8:94972/1‑136 (MQ=255)                                                                                                                  ccGCGCACATCGAGCAGACGCTGATCGGTCTGCTCGTTGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGa  <  8:222153/136‑1 (MQ=255)                                                                                                                                   |                                                                                                                       CGCCAGCGAGTCGTTGGCTGAGGCAGAACCCTTACAAGCCCTTGAAGAAACCTTTGGTGAAGACGGCCCCTATGCCCTGCCCTACGCAGTGGAGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGA  >  KT2440/125333‑125581

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