Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	2,094,808	A→C	61.6%	L186R (CTG→CGG)	CJ019_01995 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	2,094,808	0	A	C	61.6%	1.9 / 5.6	31	L186R (CTG→CGG)	CJ019_01995	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/2);  new base C (12/7);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCAGGCGCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  KT2440/2094682‑2094936                                                                                                                                |                                                                                                                                   ggggCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCCGCGCGcc                                                                                                                            >  3:347508/1‑135 (MQ=255)                     cTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCGGCGCGCC‑TGGGCGGGCTGGCCGTTTTg                                                                                                       >  4:97375/1‑135 (MQ=255)                      ggggCTTGGGCCGGGGCAAACGGTTCCTCACGATCCGTCCGCTGCTCTACAAACCCCTGCGGCGCCGGCTGCTGGC‑CAAAAACCGCCCCGCCGCGCGCGACGGCCCGGCGCGCC‑GGGCCGGGCTGGCCTTTTTGc                                                                                                      <  4:54005/134‑1 (MQ=255)                          cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCGGCGCGCC‑GGGCCGGGCGGGGCGTTTTGCAggg                                                                                                  >  6:407367/1‑135 (MQ=255)                          cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCGGCGCGCC‑GGGCCGGGCGGGGCGGTTTGGAggg                                                                                                  >  3:145705/1‑135 (MQ=255)                          cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCGGCCCCCC‑AGGCCGGGCGGGCCTGGTTGCCggg                                                                                                  >  1:319022/1‑135 (MQ=255)                          cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCAGGCGCGCC‑AGGCCGGGCGGGCCGTTTTTCAggg                                                                                                  >  6:242626/1‑135 (MQ=255)                                 ccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACAACTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCGGAGCCCC‑TGGACCGGCGGGCCGGTTTTCAGGGCCGGGTgg                                                                                          >  8:233017/1‑136 (MQ=255)                                   gggCAAACTTTTTCTCAACATCGGACCGCTTCTAATCAAACACCTTCGCCTCGGGCTGGTTCCTAAAACACGG‑CCCGCCGCGCCCGCCGGCCTGGCGCGCC‑TCCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGc                                                                                         <  4:175285/135‑1 (MQ=255)                                       cAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCCGCGCGCCGCGGCCGGCGGGGCGTTTTTCCGGGGCCGGTTGGcgag                                                                                      >  6:430659/1‑132 (MQ=255)                                                  cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCCGCGCGCC‑TGGCCGGGCTGGGCTTTTTTCAGGGCCTGTTGGCCGCGGTTGTTgaag                                                                           >  4:194717/1‑134 (MQ=255)                                                              gcggctCTTCATACACTTTCGGCGCCGGCGGATGGCGAAAATCCTC‑CCCGCCGCGCGCGAGGGCCTGGCGCGCG‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgc                                                               <  4:142334/131‑1 (MQ=255)                                                                                gcGGCGCCGGGTGCTTGCAAAAATACTG‑CACGCCGCGCTCGACGCCCCGGCGCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCg                                           <  8:41188/135‑1 (MQ=255)                                                                                  ggCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCAGGCGCGCC‑GCGCCCGGCTGCCCGTTTTTCAGGGCCTCTTGGCCGAGGGTGTTGAATACCTCTATATTCAGCAGCACCAGGATAAGGCgg                                          >  6:160932/1‑135 (MQ=255)                                                                                  ggCGCCGGCTGCTGGCGAAAAAACTG‑CCCGCCGCGCTCGACGGCCCGGCGCGCC‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCg                                           <  4:347508/134‑1 (MQ=255)                                                                                  ggCGCCGGCTGCTGCCCCAAATCCCC‑CCCCCCGCGCGCGCCGCCCCGGCGCGCC‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                            <  1:15165/133‑1 (MQ=255)                                                                                    cTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCCGGCGCGCC‑GGGCCGGGCTGGGCGGTTTGCCGGGCCTGTTGGGCGAGGTTGGTGGAGCACTCGATATTCAGCCGCCCCCCGTGATGGCGG‑At                                       >  4:364532/1‑135 (MQ=255)                                                                                    cTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCACGCGCGCC‑TGGCCGGGCTGGGCTTTTTGCCGGGCCGGTTGGCCGAGGGTGGTGAAGACATCTATATGCAGCAGCCCCAGGATATGGCggct                                        >  6:274254/1‑133 (MQ=255)                                                                                         gctgctTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCAGGCGCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGGCGAGGGTGTTGAAGAACTCTATATGCAGCAGCCCCCGGATATGGCGG‑ATTTCCt                                  >  6:83223/1‑135 (MQ=255)                                                                                          ctgctTGCGCAAATACTG‑CACGCCGCGCTCGAAGGGCCGGCGCCGC‑AGGCCGGGCGGGGCGTTTTTGATGGCCGGGTGGCCGCGGTTTTTGAATTACGATATAAGCACCAGCCCCCCGAAATCGCGGCATTTct                                   >  6:375584/1‑133 (MQ=255)                                                                                           tgctTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCAGGCGCGCC‑TGGCCGGGCTGGGCGTTTTGCAGGGCCGGTTGGGCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTg                                >  2:27311/1‑135 (MQ=255)                                                                                                                gccgcGCTCGAAGGCCCGGCGCGCC‑TGGCCGGGCTGGCCGTTTTGCAGGGCCTGGTGGCCGAGGTTGTTGTAGTACTCAATATGCACCACCCCCAGGCTATGGCGG‑ATTTCCTTGACCCAATGCTTTGGCTCgcc            >  5:167203/1‑133 (MQ=255)                                                                                                                gccgcGCTCGAAGGCCAGGCGCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGCTATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGAc            >  3:280821/1‑135 (MQ=255)                                                                                                                gccgcGCTCGAAGGCCAGGCGCGCC‑TGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGAc            >  5:231414/1‑135 (MQ=255)                                                                                                                     gcTCGAAGGCCCGGCGCGCC‑TGGCCGGGCTGGGCGTTTTGCAGGGCCGGTTGGCCGCGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑CTTTCCTTGACCCAATGCTTGGGCTCGCCAGCgg       >  1:68161/1‑135 (MQ=255)                                                                                                                      cTCGAAGGCCAGGCGCGCC‑TGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCAACACGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTAGACAGCggg      >  8:380504/1‑135 (MQ=255)                                                                                                                       cccaaGGCCAGGCGCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCa    <  8:108362/133‑1 (MQ=255)                                                                                                                         gccggCCCGGCGCGCG‑GGCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  7:393516/133‑1 (MQ=255)                                                                                                                          ccggCCCGGCGCGCC‑GGCCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  2:68161/133‑1 (MQ=255)                                                                                                                          ccggCCAGGCGCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  2:319022/133‑1 (MQ=255)                                                                                                                          acggCCCGGCGCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag  <  4:19630/133‑1 (MQ=255)                                                                                                                                |                                                                                                                                   GGGGCCATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTG‑CACGCCGCGCTCGAAGGCCAGGCGCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGG‑ATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG  >  KT2440/2094682‑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: ‑