Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	2,094,817	T→G	66.5%	Q183P (CAG→CCG)	CJ019_01995 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	2,094,817	0	T	G	66.5%	‑3.8 / 4.9	21	Q183P (CAG→CCG)	CJ019_01995	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base G (5/6);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

CGGTCCAGTACCTGGGCATTGGCCCGGGCC‑AACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGG‑CCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  KT2440/2094691‑2094936                                                                                                                                |                                                                                                                         cGGTCCAGTACCTGGGCATTGGACCGGGCC‑AACACGTTCTCAAGATACGTCAGCTGCTCTTGATACACCTTCCGCTCCGGCTGCTTGAGCAAATACTGCACGCCGCGCTCGACGGCCACGCGCCTCGGG‑CCgggg                                                                                                                 >  8:74758/1‑134 (MQ=255)                 cATTGGCCCGGGCC‑AACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGCG‑CCCGGCGGGGCGGTTTGGAggg                                                                                                 >  1:51602/1‑135 (MQ=255)                 cATTGGCCCGGGCC‑AACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCCGGCCCGGGGGGGCCTTTTTggggg                                                                                                  >  1:277300/1‑131 (MQ=255)                 cATTGGCCCGGGCC‑AACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGG‑CCGGGCTGGGCGTTTTGCAggg                                                                                                 >  4:190243/1‑135 (MQ=255)                 cATTGGCCCGGGCC‑AACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCATCGGCTCCGGCTGCTTGCGCAAAAACTGCACGCAGCGCTCGAAGGCCAGGCGCCCCCGG‑CCCGGCTGGCCTGTTTgaatgg                                                                                                 >  7:44119/1‑130 (MQ=255)                  agtGGCCCGGCCA‑AACCGTTTCGCACGAACCGACCGCGGCTCTACAACCACCTGCGGCTCCGGCTGCTTGCCCAAATCCGCCCCGCCGCGCGCGACGGCCCGGCGCGCCGGG‑CCGGGCTGGCCTTTTTGCAGGGc                                                                                                <  8:162785/133‑1 (MQ=255)                   ttGGCCCGGGCCAAACTGTTACTCAACACCCGTCCGCTGCTCTACAAACACTTCCGGCGCCGGCGGCTGGCCCAAAACCTCCCCCCCGCCCGCCACGGCCAGGCGCGCCGGC‑CCGGGCTGGCCTTTTTGCAGGGc                                                                                                <  8:286949/135‑1 (MQ=255)                            gCC‑AACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCCGG‑CCGGGCCGGGCGGTTTGGAGGGCCGGTTGGGCg                                                                                      >  6:260217/1‑135 (MQ=255)                              c‑AACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGG‑CCGGGCTGGCCCTTTTGCAGGGCCGGTTGGCCGAg                                                                                    >  4:220746/1‑135 (MQ=255)                                         tcAAGATCCGCCAGCTGCTCTACAAACACCTGCGGCGCCGGCGGCTGGCCAAAAAACGCCACCCCGCGCGCCACGGCCCGGCGCGCCGGG‑CCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTga                                                                            <  8:72001/133‑1 (MQ=255)                                                     aGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGG‑CCGGGCGGGCCGTTTTGCAGGGCCTGTTGGGCGAGGTTGTTGTAGAACTCGATATgc                                                              >  7:72001/1‑135 (MQ=255)                                                                          ggCGCCGGGTGCTTGCTAACATCCTGCCCGCCGCGCGCGACGGCCCGGCGCGCCGGG‑CCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                           <  7:193560/133‑1 (MQ=255)                                                                              ccGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGG‑CCGGGCTGGGCGGTTTGCAGGGCCGGGTGGCCGAGGGTGTTGTAGAACTCTATATTCAGCAGCCCCAGGATATGGCGGAttt                                     >  4:120731/1‑135 (MQ=255)                                                                                   ggctGGCGCAAATCCGGCACGCCGCGCTCGACGCCCCGGCGCGCCTGG‑CCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTg                                <  2:83580/134‑1 (MQ=255)                                                                                         cgaaaaTACTGCACGCGGCGCTCGACGGCCCGGCGCGCCGGG‑CCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa                          <  8:227808/132‑1 (MQ=255)                                                                                                       gccgcGCTCGAAGGCCCGGCGCGCCTGG‑GCGGGCTGGCCTTTTTTCAGGGGCTGTTGGCCGAGGTTGGTGAAGAACTCTATATGCAGCAGCCCCAGGATATGGCGGATTTACGTGACCCAATGCTTGGCCTCgcc            >  6:203461/1‑133 (MQ=255)                                                                                                       gccgcGCTCGAAGGCCCGGCGCGCCTGG‑CCGGGCGGGGCTTTTTGCAGGGCCTGTTGGGCGAGGTTGTTGAAGAACTCTATATGCAGCAGCCCCACGATATGGGGGCTTTCCTTGACCCAAAGCTTGGCcgcg              >  3:195508/1‑133 (MQ=255)                                                                                                        cggcgcgcGACGGCCGGGCGCGCGGGG‑CCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACa           <  3:220746/133‑1 (MQ=255)                                                                                                            gcTCGAAGGCCCGGCGCGCCTGG‑CCGGGGGGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTGCTTGACCCAATGCTTGGCCTCGCCAGCgg       >  3:139960/1‑135 (MQ=255)                                                                                                             cTCGAAGGCCAGGCGCGCCTGG‑CCGGGCTGGCCTTTTTGCCGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGGCAGCggg      >  4:47761/1‑135 (MQ=255)                                                                                                                 aaGCCCCGGCGCGCCTGG‑CCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  4:93984/135‑1 (MQ=255)                                                                                                                                |                                                                                                                         CGGTCCAGTACCTGGGCATTGGCCCGGGCC‑AACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGG‑CCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  KT2440/2094691‑2094936

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