Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	KT2440	3,440,343	A→C	S18R (AGC→CGC)	CJ019_03266 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	3,440,343	0	A	C	85.7%	10.5 / ‑2.7	22	S18R (AGC→CGC)	CJ019_03266	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (16/2);  total (19/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

GATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGG  >  KT2440/3440211‑3440468                                                                                                                                     |                                                                                                                              gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgc                                                                                                                             >  4:219408/1‑135 (MQ=255)                    ggtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAg                                                                                                          >  1:51305/1‑135 (MQ=255)                         tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCCCCCGGgccgccacg                                                                                                     >  5:145984/1‑132 (MQ=255)                           gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCTGGgt                                                                                                   >  7:152562/1‑135 (MQ=255)                           gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAACCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACACGCCTGGTCGCCGTGGCGCTGCTGCTGGGCCCCCCCCCAGGCGGCCCCCCGGCAGccgcggt                                                                                                   >  6:2222/1‑130 (MQ=255)                                  tCTGGAGGATCCCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAg                                                                                            >  4:134927/1‑135 (MQ=255)                                       aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCACCCGGGGGGTGGCAGGCCac                                                                                       >  8:35338/1‑134 (MQ=255)                                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGTGTGGCAGGCCACCGCCta                                                                                 >  5:100310/1‑133 (MQ=255)                                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCTGGGGGTGGCCGGCCCACGccga                                                                                 >  1:36611/1‑131 (MQ=255)                                               aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGGCGGCCCACGCCgata                                                                               >  8:154002/1‑135 (MQ=255)                                                                agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCCTGAGGCGCCCCGGGAGCCTGGGGGTGGGAGGGAAACCCCGATAATGATCAGCTCGCTgc                                                               >  8:225630/1‑133 (MQ=255)                                                                agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCTGGCGGCGCCCGCGCCGCCTGGGTGTGGCAGGCCAACGCCAATAAGGCGCAGCTCGCGGAt                                                              >  8:178386/1‑135 (MQ=255)                                                                 gcgcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACCCCTATAAGCGGCAGCTCGCGGa                                                               >  4:198883/1‑133 (MQ=255)                                                                   gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCCACGCCGAGAAGGCGCAGCTCGCTGCTCAg                                                           >  8:39070/1‑135 (MQ=255)                                                                                        gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGGGGGGGGGCCCCCCCGAGAAGGCGCAGCGCGCGGATCGGGCCGCTGGCGTCGGCGGGCAg                                      >  1:78053/1‑135 (MQ=255)                                                                                        gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCGCGGGCCGCCGGGGGGTGGGCGGCCCCCGCCCATAAGGCGCCGCACGCGGCTCCGGGCGCTGACTTCGGCCCGCAg                                      >  4:140397/1‑135 (MQ=255)                                                                                        gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGGGCCCCGGCCGCCGGGGGGTGGCGGGCCCCCCCCCATAAGGCGCCGCGCGCGGGTGACGCCGCGGCCATCGGCcagcag                                      >  4:232746/1‑135 (MQ=255)                                                                                          ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGGGTGGGAGGGCCCCGCCGAAAATGCGCCGCTCGCGGATGAGGCCGCTGTCTTAGGCAACCACtt                                    >  6:39838/1‑135 (MQ=255)                                                                                             ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCCCCGCCGATAAGGCGCAGCGCGCGGGTCAGGCCGATGCCGTCGGCCAGCAGTTGGc                                 >  8:39203/1‑135 (MQ=255)                                                                                                     gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGTGTGGGCGGGCCCCCCCGATAATGCGCCGCGCGCGGCTCAGGCCCCGGACGTCGGCCAGCCGTTGGGGGGGCAgg                         >  4:187943/1‑135 (MQ=255)                                                                                                            cgcgGTGGCGGTGCGGCGGGCCGCCTGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                  <  7:35338/135‑1 (MQ=255)                                                                                                                           gcgggCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCg   <  4:139009/132‑1 (MQ=255)                                                                                                                            cgggCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg  <  8:48726/133‑1 (MQ=255)                                                                                                                                     |                                                                                                                              GATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGG  >  KT2440/3440211‑3440468

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: ‑