Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	1,313,950	T→G	35.1%	R612R (AGG→CGG)	CJ019_01235 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	1,313,950	0	T	G	35.1%	32.8 / 2.6	34	R612R (AGG→CGG)	CJ019_01235	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/7);  new base G (4/9);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 7.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCA  >  KT2440/1313823‑1314080                                                                                                                                |                                                                                                                                   tCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGGGGGTAATGCAACTCGAACTCCTCATTGCGCACGGCCCGGGGCAgg                                                                                                                             >  1:393116/1‑135 (MQ=255)    ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGCCCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGGAATGCAACTCGAACTCCTCATTGCGCAGGGCCCGGCGCCGGTCg                                                                                                                          >  1:323280/1‑135 (MQ=255)    ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCg                                                                                                                          >  6:130938/1‑135 (MQ=255)         cAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCGGCGCAGGGCGTtgtc                                                                                                                     >  5:228226/1‑135 (MQ=255)                ccTTCCTGTCGAGGCCCCCCGCCCACCAGGCCCGCGCCCGCGCCTATGCGCCGGGCGCCCCGGGGGTCGGGGGGCGTGTAATCCAACCCCACCTCCTCATGGCGCCGGCCCCGGCGCAGGTCGTTCTCCAGCTgg                                                                                                              <  8:5010/135‑1 (MQ=255)                       gtgGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGt                                                                                                       >  2:12797/1‑135 (MQ=255)                             gccgccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGa                                                                                                 >  3:386369/1‑135 (MQ=255)                               cgccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATc                                                                                               >  7:9544/1‑135 (MQ=255)                                gccAGCGCACCAAGGCGTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCGCCAGCTGGCGGCGGTAGTGGATcg                                                                                              >  2:98893/1‑134 (MQ=255)                                              gCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCg                                                                                >  3:332432/1‑135 (MQ=255)                                                    aCCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACGGGCTCTGCTGGTTCATTTCGGCAAc                                                                           >  7:326666/1‑134 (MQ=255)                                                     ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAAc                                                                           >  7:40045/1‑133 (MQ=255)                                                       gcgACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGTCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGaa                                                                        >  7:110514/1‑134 (MQ=255)                                                       gcgACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGaa                                                                        >  8:358306/1‑134 (MQ=255)                                                           cgaTTCGCCGGGCGTCCCGGCGGTTGCGTGGCTGGTAATCCAACTCGAACTCCTCATTGCGCTGGGCCGTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAgc                                                                   <  6:278013/135‑1 (MQ=255)                                                           cgaTTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAgc                                                                   >  4:153112/1‑135 (MQ=255)                                                                 gCAGGTCGTCCCGGGGGTAGGGGGGCTGGTAATGCACCTCGAACTCCTCATTCCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCAt                                                              <  6:192152/134‑1 (MQ=255)                                                                          ccAGGGGGGAGGGGGGCTGTTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTc                                                    <  1:17517/135‑1 (MQ=255)                                                                                gtgcgggTGGCTGGTATGCCACCTCGAACTCCTCATGGCGCGGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                              <  2:370801/129‑1 (MQ=255)                                                                                ggttgggTGGCTGTTACTGCAAGTCGACCGCCTCATGGCGCAGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                              <  2:262141/129‑1 (MQ=255)                                                                                ggtcgggTGGCTGGTACTGAACCTCGACCTCCTCTTTGCGCTGGGCCCGGCGCAGGTCGTGCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                              <  7:203508/129‑1 (MQ=255)                                                                                ggTGGCGTGGGTTCTACTGAACCTCGCACTCATCTTGGCGGCGGGCCGTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                              <  3:144178/135‑1 (MQ=255)                                                                                ggTAGGGTGGCTGGTAATGCAACTCGAACTCCTCATGGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                              <  7:73263/135‑1 (MQ=255)                                                                                     ggTGGCGGGTACTCCACCTCCACCCCCTCATGGCCCCGGCCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                         <  2:174882/134‑1 (MQ=255)                                                                                     cGTGGCTGGTAATGCACCGCCAACTCCTAATTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                         <  6:298084/135‑1 (MQ=255)                                                                                          ctgtaaaTGCACCTCGAACTCCTCATGCCGCGGGCCCCTGCCCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTg                                    <  6:334962/130‑1 (MQ=255)                                                                                           tGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCCACGCAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTgg                                   >  5:13509/1‑135 (MQ=255)                                                                                            ggTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGt                                  >  4:18095/1‑135 (MQ=255)                                                                                                          acctcctcATTGCGCAGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCg                    <  8:110514/133‑1 (MQ=255)                                                                                                                   gtGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCg            <  4:386369/133‑1 (MQ=255)                                                                                                                     gcgcCGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGAt         <  6:228226/135‑1 (MQ=255)                                                                                                                          gggCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCt    <  1:98893/135‑1 (MQ=255)                                                                                                                            gCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCa  >  3:138951/1‑135 (MQ=255)                                                                                                                            gCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTAGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCa  >  3:183151/1‑135 (MQ=255)                                                                                                                                |                                                                                                                                   TCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCA  >  KT2440/1313823‑1314080

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