Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	3,157,749	C→A	9.8%	P68H (CCT→CAT)	CJ019_02972 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	3,157,749	0	C	A	9.8%	110.7 / 4.7	41	P68H (CCT→CAT)	CJ019_02972	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base A (2/2);  total (21/20)
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.68e-01

TGTGCGCCTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTGCTGGGA  >  KT2440/3157616‑3157878                                                                                                                                      |                                                                                                                                  tgtgCGCCTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGcc                                                                                                                                   >  1:64925/1‑134 (MQ=255) tgtgCGCCTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCt                                                                                                                                  >  1:6937/1‑135 (MQ=255)     cgcCTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAgg                                                                                                                              <  5:116815/135‑1 (MQ=255)     cgcCTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAgg                                                                                                                              <  5:86067/135‑1 (MQ=255)       ccTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGgc                                                                                                                            >  3:32619/1‑135 (MQ=255)       ccTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGgc                                                                                                                            >  5:21898/1‑135 (MQ=255)       ccTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGGGTGTACGTCCGCAGCCTGTCACGCAGCCTCGAGCCTCCGGgc                                                                                                                            >  6:200351/1‑135 (MQ=255)          cGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGgcgcg                                                                                                                         >  8:255413/1‑135 (MQ=255)            tGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGgcgcg                                                                                                                         <  4:32619/133‑1 (MQ=255)            tGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGgcgcg                                                                                                                         <  3:219021/133‑1 (MQ=255)             gACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGgcgcgcgc                                                                                                                      <  4:54766/135‑1 (MQ=255)                     aCCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTa                                                                                                              >  5:150583/1‑135 (MQ=255)                     aCCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTa                                                                                                              >  7:87828/1‑135 (MQ=255)                               gcgcAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCgccggc                                                                                                    <  6:21898/135‑1 (MQ=255)                               gcgcAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCgccggc                                                                                                    >  4:192719/1‑135 (MQ=255)                                            aTTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCga                                                                                         <  2:6937/133‑1 (MQ=255)                                                cGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGGATTGCGCCGGCCGTCTCGCCGAGACCtt                                                                                   >  7:116620/1‑135 (MQ=255)                                                 gACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCtt                                                                                   <  1:199945/134‑1 (MQ=255)                                                                        tgtCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCcgacga                                                           >  5:67732/1‑135 (MQ=255)                                                                           cGCAGCGCACGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGc                                                         <  5:149758/134‑1 (MQ=255)                                                                                   ccgccTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCtt                                                                                   <  8:44378/100‑1 (MQ=255)                                                                                   ccgccTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCtt                                                                                   >  7:44378/1‑100 (MQ=255)                                                                                     gccTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTg                                              >  8:47453/1‑135 (MQ=255)                                                                                     gccTTGCCGGCGGGGGTGTACGCCCGCAGCCTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCGGCTGTATTGCGCCGGCCGTCTCGCCGAGCCCTTGCGTGGGCGGGCGCTGCCCGCCGCGCGGCGCGCATTg                                              >  3:157065/1‑135 (MQ=255)                                                                                                  gtgtgtACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTc                                 >  4:189005/1‑135 (MQ=255)                                                                                                    gtgtACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGCGCCCTTGCGTGGGCTGGAGCTGCCCGCCGCGCGGATCGCCTTGCCCGCGGCCCTTCCg                               >  3:163843/1‑135 (MQ=255)                                                                                                    gtgtACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCg                               <  3:192719/135‑1 (MQ=255)                                                                                                          gACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGt                          <  6:67732/134‑1 (MQ=255)                                                                                                          gACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTg                         <  3:189005/135‑1 (MQ=255)                                                                                                                   cTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGcccc                <  8:116620/135‑1 (MQ=255)                                                                                                                         cGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGt                          <  4:7687/119‑1 (MQ=255)                                                                                                                         cGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGt                          >  3:7687/1‑119 (MQ=255)                                                                                                                         cGCAGCCTCGAGCATCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCg          >  7:229310/1‑135 (MQ=255)                                                                                                                         cGCAGCCTCGAGCATCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCg          >  7:233216/1‑135 (MQ=255)                                                                                                                          gCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGt         >  4:20542/1‑135 (MQ=255)                                                                                                                           cAGCCTCGAGCATCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTg        <  8:229310/135‑1 (MQ=255)                                                                                                                           cAGCCTCGAGCATCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTg        <  8:233216/135‑1 (MQ=255)                                                                                                                                tCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTGCTggg   <  8:87828/135‑1 (MQ=255)                                                                                                                                tCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTGCTggg   <  6:150583/135‑1 (MQ=255)                                                                                                                                tCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTGCTggg   <  2:64925/135‑1 (MQ=255)                                                                                                                                 cGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTGCTGGGa  >  4:44093/1‑135 (MQ=255)                                                                                                                                 cGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTGCTGGGa  >  1:277456/1‑135 (MQ=255)                                                                                                                                      |                                                                                                                                  TGTGCGCCTCGTGACAGTCAACCTCAGGCCGCGCAACCGCTGGATTTCGACCCGGCACGTCTGGCCCATTGTGTCGCAGCGCCCGCCTTGCCGGCGGGTGTGTACGACCGCAGACTGTCCCGCAGCCTCGAGCCTCAGGGCGCGCGCCTGCTGTATTGCGCCGGCCGTCTCGCCGAGACCTTGCGTGGGCTGGAGCTGCCCGACGAGCGGATCGCATTGACCGCTGCCATTCCGGAAGTGGATGGCCCCAGTTCGTGCTGGGA  >  KT2440/3157616‑3157878

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