Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	4,388,450	A→T	14.9%	V56E (GTG→GAG)	CJ019_04138 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	4,388,450	0	A	T	14.9%	49.9 / 4.7	27	V56E (GTG→GAG)	CJ019_04138	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/9);  new base T (2/2);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

GGTCGCCTTGAGCGGCAGGCACCAGACCTCGCCCACCAGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGACCTCGGCGTAGAACAGCGGCACACGCAGGT  >  KT2440/4388329‑4388585                                                                                                                          |                                                                                                                                        ggTCGCCTTGAGCGGCAGGCACCAGACCTCGCCCACCAGGTCGCTGCCACGCCAGGTCTCGTACATGTCGGCCCGCATCGCACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGcct                                                                                                                            >  3:270722/1‑133 (MQ=255)     gCCTTGAGCGGCAGGCACCAGACCTCGCCCACCAGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCg                                                                                                                        >  4:195001/1‑135 (MQ=255)           cgCGGCAGGCACCAGACCTCGCCCACCAGGGCGAGGCCACGCCAGGTCTCGTACATGTCGGCCAGCAGCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTc                                                                                                                   <  6:259467/133‑1 (MQ=255)                   gcaCCAGACCTCGCCCACCAGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCgcatcg                                                                                                         >  7:172983/1‑136 (MQ=255)                   gcaCCAGACCTCGCCCACCAGGTCGATGCCACGCCAGGTCTCGTACAAGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGTGAGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGAgcatc                                                                                                          >  6:252019/1‑135 (MQ=255)                    caCCAGACCTCGCCCACCAGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCgcatcg                                                                                                         <  3:92388/135‑1 (MQ=255)                    caCCAGACCTCGCCCACCAGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCgcatcg                                                                                                         <  5:119568/135‑1 (MQ=255)                                  caccaGGTCGATGCCACGCCAGGTCTGGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGAc                                                                                           >  2:81975/1‑135 (MQ=255)                                  caccaGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGCACCGACGCCGCCGCTTGTCTTGGCCAGGCGGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGAATCGTCATTGACt                                                                                          >  8:136898/1‑136 (MQ=255)                                  caccaGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGCTTGTCTTGGCCAGGCTGCGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGGCATGGAc                                                                                           >  1:135914/1‑135 (MQ=255)                                  caccaGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGAc                                                                                           >  6:81031/1‑135 (MQ=255)                                               ccACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGCCTTGGCCGACGAt                                                                               >  1:172784/1‑134 (MQ=255)                                                                cATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGcc                                                             >  6:255750/1‑135 (MQ=255)                                                                      ggCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCTCGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGgc                                                       >  2:89836/1‑135 (MQ=255)                                                                      ggCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCTCGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGgc                                                       >  2:93516/1‑135 (MQ=255)                                                                      ggCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGgcg                                                      >  8:260169/1‑136 (MQ=255)                                                                           gCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGAGCAAGCGCAGCAGGCTGCCCGTGGCGCTgg                                                  >  5:12402/1‑134 (MQ=255)                                                                                           gccgATTGTCTTGGCCAGGCGGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTc                                  <  3:59292/135‑1 (MQ=255)                                                                                            ccgATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGa                                <  8:115017/136‑1 (MQ=255)                                                                                             cgATTGTCTTGGCCAGGCTGAGGCTCGGCTCGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGa                                <  1:93516/135‑1 (MQ=255)                                                                                             cgATTGTCTTGGCCAGGCTGAGGCTCGGCTCGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGa                                <  1:89836/135‑1 (MQ=255)                                                                                             cgATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGa                                <  6:99860/135‑1 (MQ=255)                                                                                             cgATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGa                                <  1:81975/135‑1 (MQ=255)                                                                                             cgATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGa                                <  1:172659/135‑1 (MQ=255)                                                                                                   tCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGACCTCg                           <  2:148140/134‑1 (MQ=255)                                                                                                           aGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGACCTCGGCGTAGa                    >  1:36420/1‑133 (MQ=255)                                                                                                               tGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGACCTCGGCGTAGAACAGCgg             >  7:134487/1‑136 (MQ=255)                                                                                                                          aCGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGACCTCGGCGTAGAACAGCGGCACCCGCAGGt  >  8:115748/1‑136 (MQ=255)                                                                                                                          aCGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGACCTCGGCGTAGAACAGCGGCACACGCAgg   >  6:208786/1‑135 (MQ=255)                                                                                                                          |                                                                                                                                        GGTCGCCTTGAGCGGCAGGCACCAGACCTCGCCCACCAGGTCGATGCCACGCCAGGTCTCGTACATGTCGGCCAGCATCGAACCGACGCCGCCGATTGTCTTGGCCAGGCTGAGGCTCGGCACGAGGGTGAGCTTGCCGATTTCCGGCGCATCGGCATCGTCATTGACTTGGCCGACGATCAAGCGCAGCAGGCTGCCCGTGGCGCTGTTGGCCTGGCTGTTGTCGACCTCGGCGTAGAACAGCGGCACACGCAGGT  >  KT2440/4388329‑4388585

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