Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	3,947,782	A→C	29.0%	L129R (CTG→CGG)	CJ019_03732 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	3,947,782	0	A	C	29.0%	41.0 / 3.3	35	L129R (CTG→CGG)	CJ019_03732	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/10);  new base C (4/7);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTCCCCGGCCGGTGCCGACTGCGCCGGGTCGATCCACGCCCCCGGGGCCATACCCGACAGCG  >  KT2440/3947651‑3947908                                                                                                                                    |                                                                                                                               ccATGGGCCGCCCTGGCAGCTCGATGGGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGGCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCgg                                                                                                                            <  3:366788/136‑1 (MQ=255)    tGGGCCGCCCGGGCAGCTCGATGCGCGCCATGGGGCTGGGCACGCCCCCCCCCGACACCGTCGGGTCGATCACTGTCAGGCGGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGcc                                                                                                                          <  5:78819/135‑1 (MQ=255)        ccgccCTGGCAGCTCGCTGCGCGCCATGGTGCTGGTCACGCCCCCCCCCGCCCCCGGCGGGGCGTGGCCTGTCAGGCGGGCGGGCTTGTTGAGCATCAGGAACCAGTAGCGGTTCACGTCGCTCCGCGGCCGGGc                                                                                                                      >  5:50058/1‑135 (MQ=255)             cTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGCCACCGTCGGGTCGGTCCGTGTCAGGGGGGCGGCGTTGTTGGGGCTGAGCCTGGAGGAGGTGGTGATGTCGTGCCGCGGGCGGGCCCCGa                                                                                                                 >  2:236984/1‑135 (MQ=255)             cTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGGCAGGCTGGCGTGATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGCGCAGCGGCCGGGCCCCGa                                                                                                                 >  7:335509/1‑135 (MQ=255)              tGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGCCACCGTCGGGTCGCTCACTGTCAGGCTGGCGTGATTGTTGAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGag                                                                                                                >  2:49482/1‑134 (MQ=255)               ggCAGCTCGAGGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGGCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGAtgc                                                                                                              <  4:107548/136‑1 (MQ=255)                gCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGAtgc                                                                                                              <  8:372288/135‑1 (MQ=255)                gCAGCTCGATGCGCGCCAGGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGAtgc                                                                                                              <  2:85095/135‑1 (MQ=255)                 cAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGGCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGGCGAGCAGCGGCCGGGCCCCGAtgctg                                                                                                            <  3:37761/136‑1 (MQ=255)                         tGCGCGCCTTGGTGCTGGTCACGCCCCCCACCGACACGGTCGGGTCGATCACTGTCTGGGTGGCGTCTTTGTTCAGAATCAGCAAGAATTAGCTTTTCAGGTCGAGAAGCGGCCGGGCCCCGATGCTGCTCTTga                                                                                                     <  7:273195/135‑1 (MQ=255)                           cgcgcCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGGTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGCGCAGCGGCCGGGCCCCGCTGCGGCTCTTgag                                                                                                    >  6:98815/1‑134 (MQ=255)                                 aTGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGCGCAGCGGCCGGGCCCCGCTGCTGCGCTTGAGCAACAgcg                                                                                            >  3:50375/1‑136 (MQ=255)                                       cTGGACCCGCCCCCCACCGACACCGGCGGGGCGATCACTGTCAGGCGGGCGTCATTGTTCAGCATCAGCAAGAAGTCGCGGTTCCGGGCGACCCGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGAccc                                                                                      <  4:72527/136‑1 (MQ=255)                                        tGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGGCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGGCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGAccc                                                                                      <  7:272754/135‑1 (MQ=255)                                             aCGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTTTTATGCATCAGCAAGAAGTACCTGTTCACGTCCAGCCGCGGCCGCGCCCCGCACCTGCACTTGAGAACAAGCCCGACCCGccccg                                                                                >  4:154370/1‑135 (MQ=255)                                             aCGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGcc                                                                                 >  7:41187/1‑135 (MQ=255)                                             aCGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCt                                                                                >  4:38729/1‑136 (MQ=255)                                             aCGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAACAAGTAGCTGTTGAGGTCGCGCAGCCGCCGGGCCCCGCTGCTGCTCGTGAGAAACATCCCGCCCCCCGcc                                                                                 >  6:119117/1‑135 (MQ=255)                                                  ccccACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCTTTGTTGAGCATCAGCAAGAAGTAGCTGTTAAGGTCGAGCAGCGGCCGGGCCCCGAGGCTGCGCTTGTGAAACAGCGCGCCCCGCGCCTTGCa                                                                            >  7:326072/1‑135 (MQ=255)                                                  ccccACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGCGCAGCGGCCGGGCCCCGCTGCTGCTCTTGAGAAACCGCGCGACCCGCGCCGTGCa                                                                            >  8:196001/1‑135 (MQ=255)                                                  ccccACCGACACCGTCGGGGCGATCACTGTCCGGCTGGCGTCCTTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGAGGCGGCTCGTGAGAAACCGCGCGCCCCGCGCCTTGCa                                                                            >  2:33067/1‑135 (MQ=255)                                                               gTCGGGTCGATCACTGTCAGGCTGGCCTCATTGTTCAGCATCAGCAAGACGTCGCTGTTCAGGGCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCgg                                                               <  7:273528/135‑1 (MQ=255)                                                                gcGGGTCGATCCCTGACCGGCTGGCGTCATTGTGCCGCATCCCCAAGAAGTCGCGGTGCCGGTCGACCCGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTc                                                             <  4:71929/135‑1 (MQ=255)                                                                   ggTCGATCACTGTCAGGCTGGCGTCATTGTTACGCATCAGCAAGAAGTAGCTGTTCAGGTCGCGCCGCGGCCGGGCCCCGATGCGGCGCTTGAGAAACAGCGCGACCCGCGCCGTGCCGCGCTGCTCGCGGTccc                                                           >  5:341/1‑135 (MQ=255)                                                                           aCTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTCCCCGGCCGGt                                                   >  8:287178/1‑135 (MQ=255)                                                                             tGTCTGGCTGGCGGCTGTTCTCCGCACCCGCAAGACGTCGATGATCCGCCCGAGCACCGGCCGGCCCCCGATGCTGCTCTTGACAAACAGCGCGACCCGCGCCTGGCCGCGCTCGTCGCGGTCCCCGGCCGGTGc                                                 <  8:233088/135‑1 (MQ=255)                                                                                       cggTCATTGTTCCGCCCCAGCAAGAAGCCGCGGTGCCGGCCGCCCCCCGGCCGGGCCCCGATGCGGCTCTGGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGGCCCCGGCCGGTGCCGACTGCGCCg                                      <  4:136658/134‑1 (MQ=255)                                                                                          tCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGCAGCTGCTCTTGAGAAACAGCGCGCCCCGCGCCTTGCAGCTCTGCTCCCGGTCCCCGGCCGGTGCCGACTGCGCCggg                                    >  1:253435/1‑135 (MQ=255)                                                                                               gTTCAGCAGCAGCAAGACGTAGCGGTTCCGGCCCAGCCGCGGCCGGGCCCCTCTGGGGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTCCCCGGCCGGTGCCGACTGCGCCGGGTCGAt                               <  4:175517/135‑1 (MQ=255)                                                                                                tgcagcatccccaAGAAGTAGCTGTGCCGGGCGCCCCGCGGCCGGCCCCCGCGGCGGCTCTTGAGAACCCCCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGGCCCCGGCCGGTGCCGACTGCGCCGGGTCGATcc                             <  4:35747/134‑1 (MQ=255)                                                                                                           caAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTCCCCGGCCGGTGCCGACTGCGCCGGGTCGATCCACGCCCCCggg                  >  4:335715/1‑136 (MQ=255)                                                                                                             agaagTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTCCCCGGCCGGTGCCGACTGCGCCGGGTCGATCCACGCCCCCgggg                 >  1:101874/1‑135 (MQ=255)                                                                                                                     cTGTGCCGGGCGCGCCGCGGCGGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTCCCCGGCCGGTGCCGACTGCGCCGGGTCGATCCACGCCCCCGGGGCCATAccc         <  1:54794/135‑1 (MQ=255)                                                                                                                           cggccGAGCCGCGGCCGGGCCCCGATGCTGCTCTTCAGAACCCGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTCCCCGGCCGGTGCCGACTGCGCCGGGTCGATCCACGCCCCCGGGGCCATACCCGACAGCg  <  4:112665/132‑1 (MQ=255)                                                                                                                                    |                                                                                                                               CCATGGGCCGCCCTGGCAGCTCGATGCGCGCCATGGTGCTGGTCACGCCCCCCACCGACACCGTCGGGTCGATCACTGTCAGGCTGGCGTCATTGTTCAGCATCAGCAAGAAGTAGCTGTTCAGGTCGAGCAGCGGCCGGGCCCCGATGCTGCTCTTGAGAAACAGCGCGACCCGCGCCTTGCAGCTCTGCTCGCGGTCCCCGGCCGGTGCCGACTGCGCCGGGTCGATCCACGCCCCCGGGGCCATACCCGACAGCG  >  KT2440/3947651‑3947908

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