Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	3,969,776	A→C	35.9%	L246R (CTG→CGG)	CJ019_03750 ←	putative protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	3,969,776	0	A	C	35.9%	35.0 / 2.7	36	L246R (CTG→CGG)	CJ019_03750	putative protein
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base C (9/5);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGC  >  KT2440/3969649‑3969909                                                                                                                                |                                                                                                                                      gcgggCCTGGGGCGTCATACGCCGGCGAGGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGAAGCCCCACGCCAGGCACAGGGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGgc                                                                                                                                <  8:175991/132‑1 (MQ=255)             cGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCCCGCCAGGCACAGTGGGCGCGCCCCGCGCGCCGCCACCCCCGGGCGGCGCgccggccggc                                                                                                                    >  3:106932/1‑133 (MQ=255)                  tACGGCGGCGATGGCGCCGGGCGAATTCGGCCACCGCCTCGACAACCACTTCTTCGATAGCGGCACCCCCCCGCCCGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGc                                                                                                               <  8:285247/135‑1 (MQ=255)                          cgATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCTCGCGCAGCGCGCCGCCCCACCAGGTCGGCGCGCAGGCCGTCCACGCCGgcggcg                                                                                                       >  3:316300/1‑135 (MQ=255)                          cgATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCCCCCCGCCCCCCCACGGGCGCGCGCCGGCCGGCCCCCCCGgcggcg                                                                                                       >  4:160508/1‑135 (MQ=255)                            aTGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGGGCCGCGCGCCGCCCCCCCCCGGCGGCGCGCCGGCCGCCCCCCCccgcgggggt                                                                                                     >  8:40524/1‑131 (MQ=255)                            aTGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCCGCCACACCCCGGGCGCGCGCACGCCGTCCCCGCCGgcggcgga                                                                                                     >  5:296913/1‑133 (MQ=255)                            aTGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTa                                                                                                     >  3:301929/1‑135 (MQ=255)                             tGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCCCGCCGCCCCCCCCCGTCCGCGCGCCCGCCGTCCCCCCCCGCGGCCTag                                                                                                    >  8:240121/1‑135 (MQ=255)                                        gAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGGGCCGCGCGCCCCGCCCCCCCGTGGGCGCGCCGGCCGGGCACGCCGGGGGCGGGGCCGCGCGCGGt                                                                                         >  8:45809/1‑135 (MQ=255)                                                  aCCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCCCCCCCAGCCCCCCCCCGGCGGCGCGCAGGCCCGCCACCCCGGCGGCGTAGCAGCCCGCGGGAATCCACGcc                                                                               >  2:112596/1‑134 (MQ=255)                                                        gccacacCCATTTCTTCGAACGCTCCCCCCCCCCCCCCGGCACGGGGGGCGGGCGCGGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGcc                                                                         <  8:60852/132‑1 (MQ=255)                                                             aTCCACTTCTTCGATCGCGCCCCCCCCCCGCCAGGCACCGGGGGCGGGCGCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTc                                                                    <  2:32501/135‑1 (MQ=255)                                                             aTCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGGGCCGCGCGCCCCCCCCCCCGGTCGGCGCGCCCGCCCGCCCCCCCGCGCGCGGTGCAGCGCGCCGTTATACAAGCAACCGCCggggg                                                                    >  2:238605/1‑133 (MQ=255)                                                               ccACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCCGCCACACCACGTCGGCGCGCCCGCCCTCCACCCCGGCGGCGTAGCAGCGCGCGGTAATACACGCAAGCGCCAGGTcct                                                                  >  4:257821/1‑133 (MQ=255)                                                                         gATCGCTGCCCCCCCCCCCCAGCCACGGGGGGGGGGGGCTGCGGGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCa                                                          <  4:136200/133‑1 (MQ=255)                                                                           tAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCCCAGCCCCCACCCCCACCAAGTACGCGCCCCGGCCCGCCACGCCGGCGGCGTTGCCGCGCACGGTAATACACGGAAGCGCCAGGGGATTCAGGGGAATGt                                                      >  2:262631/1‑135 (MQ=255)                                                                                 cAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCCCCCCGGCGGCGCGCCGGCCCTCCACCCCGGCGGGGGAGCCGCGCGCGGTTATCCCGGCCACCGCCGGGGCCTCCAGGGCCATGTCGGCCa                                                >  6:54004/1‑135 (MQ=255)                                                                                            cAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCCCCAGGTCGGCGCGCAGGCCGGCCCCGCCGGCGGCGGAGCCGCGCGCGGTAATCCCCGCCCGCGCCGGGTCATCCCGGGCCATGTCGGCCCAGGCCGGgcgc                                     >  8:99484/1‑135 (MQ=255)                                                                                            cAGGCACAGTGGGCGCGCGCAGCGCGCCGCCACACCAGGGCGGCGCGCCGGCCCTCCACGCCCGCGGCCGAGCCGCGCTCCGGTATCCCGGCCAGCGCCTGGGCCTCCCGGGCAATGTCCGCCAAGGCCTGgcgc                                     >  2:212852/1‑135 (MQ=255)                                                                                                       ggCGCGCGCCGCGCGCCGCCCCACCCGGTCGGCGCGCCGGCCGTCCCCGCCGGCGGCGTAGCAGCGCTCGGTAATCCCGGCCAGCGCCGGGGCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCCCgggggc                          >  2:213667/1‑133 (MQ=255)                                                                                                            gccccgcgcgccgcCACACCAGGGCGGCGCGCAGGCCGGCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCt                     <  5:170301/130‑1 (MQ=255)                                                                                                             ctccccgcgaagcCACACCTGCGCGGCGGGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGGTATCGAGGCACGCGCCTGGTGATCCAGGTCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTc                    <  5:238180/130‑1 (MQ=255)                                                                                                                  cccgcagcCACACCAGGTCGGCGCGCAGGCCGGCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCa               <  2:228629/133‑1 (MQ=255)                                                                                                                       cgcCACACCCGGGCGGCGCGCAGGCCGCCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTgc          <  5:42387/134‑1 (MQ=255)                                                                                                                         ccccacCCGGTCGGCGGGCAGGCCGGCCCCGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTgcgc        <  6:213519/132‑1 (MQ=255)                                                                                                                         ccACACCCGGGCGGCGCGCAGGCCGGCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCACGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTgcgc        <  5:66830/135‑1 (MQ=255)                                                                                                                         ccACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTgcgc        <  7:17126/135‑1 (MQ=255)                                                                                                                          cacaCCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCgcc       >  4:71216/1‑135 (MQ=255)                                                                                                                             acCCGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCgcctgg    <  5:262735/135‑1 (MQ=255)                                                                                                                             acCCGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCgcctgg    <  7:324637/135‑1 (MQ=255)                                                                                                                             acCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCgcctgg    <  7:240121/135‑1 (MQ=255)                                                                                                                             acCAGGTCGGCGCGCAGGCCGGCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCACGCGCCTGGTCATCCAGGGCAATGGCGGCCAAGGCCTGGCGCGCCGCGGGGCAGCGCTCGCGCAGTTGCgcctgg    <  6:61868/135‑1 (MQ=255)                                                                                                                             acCAGGGCGCCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCgcctgg    <  6:11880/135‑1 (MQ=255)                                                                                                                               cAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCATTGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCgccggggc  >  5:5493/1‑135 (MQ=255)                                                                                                                               cAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCgcctggg   >  1:257992/1‑134 (MQ=255)                                                                                                                                |                                                                                                                                      GCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTCGCGCAGTTGCGCCTGGGC  >  KT2440/3969649‑3969909

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