Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	4,912,863	T→G	39.3%	L480R (CTG→CGG)	CJ019_04608 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	4,912,863	0	T	G	39.3%	18.7 / 3.3	36	L480R (CTG→CGG)	CJ019_04608	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/17);  new base G (7/8);  total (11/25)
Fisher's exact test for biased strand distribution p-value = 6.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGACACC  >  KT2440/4912728‑4912994                                                                                                                                        |                                                                                                                                      tCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGGGCt                                                                                                                                       >  7:329782/1‑136 (MQ=255)      cTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGGGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTgcgc                                                                                                                                   <  3:106231/135‑1 (MQ=255)                 ccTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCGGCGCGCGCGCcggca                                                                                                                        >  4:305843/1‑134 (MQ=255)                    ggccccGCGCGGCCCCCTGCAAGTGCAGGCACAACTGCGCGAGGGCGGCCTGAGCGGGCAGCTAGGGGCGGAGCGCGAGCAGACCGCGCGGTGGATCGACACCCAGGTTGGCGCGCTGCGCGCGCGCCTGC‑TGGc                                                                                                                    <  4:31499/133‑1 (MQ=255)                      ccccGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCTGTTGATCGACACCCAGCTTGGCGCGCGGCGCGCGCGCCGGC‑TGGCCCg                                                                                                                 >  8:214332/1‑136 (MQ=255)                        ccGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGGGCGGGGCCCGCGCCCGC‑GGGgcccg                                                                                                                >  1:251972/1‑133 (MQ=255)                                  cccTGCAAGTGCAGGCACAACTGCGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCCGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAg                                                                                                      <  3:305843/135‑1 (MQ=255)                                  cccTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCTACACCCAGCTGGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAg                                                                                                      <  2:8996/135‑1 (MQ=255)                                  cccTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGGGCGGCGCGCGCCCCAGCTTGGCCCCGGGCC‑GGAAg                                                                                                      >  2:248417/1‑135 (MQ=255)                                     tGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCg                                                                                                   <  4:284960/135‑1 (MQ=255)                                        aaGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGGGCGGCGGGCGCGCCAGC‑GGGGCCGGGGGCGGGGAGGCgggg                                                                                                >  2:366930/1‑133 (MQ=255)                                        aaGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCGGCGCGCCCCCCCGC‑TTGGCCGGGGCCGGGGAGTAGGCgg                                                                                               >  8:255238/1‑135 (MQ=255)                                             cAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGCCCAgg                                                                                          >  8:346744/1‑136 (MQ=255)                                                             cgggCCGCGTGAGCGGGCAGCGTGGGGCGGCGCGCGAGCAGACCGCGCGGGTGAGCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCa                                                                          <  2:332657/134‑1 (MQ=255)                                                             aGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGGGCGGCGCGCGCGCCGGC‑TGGCCCGGGGGCGGGACGTCCGGGCCCCGGGGTG‑CCCCCCcgggc                                                                          >  6:162787/1‑133 (MQ=255)                                                              gggCCGCCTGAGTGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCCACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCAt                                                                         <  5:189912/135‑1 (MQ=255)                                                              gggCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCCCGCCAGC‑AGGCCCCGGGGCGGGTATTCGGGGGCCTTGAGTT‑CCACCCCGGcct                                                                         >  4:329564/1‑133 (MQ=255)                                                              gggCCGCCGGCGCGGGCAGCTAGGGGCGGCGCGCGAGCCGCCCGCGGGGTTGATCCACCCCCCGCTGGGCGCGGGGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCAt                                                                         <  5:330549/135‑1 (MQ=255)                                                                   gccTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCAGCGCGCGCGCCGGC‑AGGCCCCGGGGCTGGGAAGCGGGGGCCCGGTGGGACCACCCGGG‑ATCccgg                                                                    >  6:134643/1‑134 (MQ=255)                                                                           gggCAGCTATGGGCGGAGCGCGAGCAGACGGCGCGGTTGATCGACACCCAGATTGGCGGGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAgggg                                                            <  2:78599/135‑1 (MQ=255)                                                                                 cTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGca                                                     <  8:131412/136‑1 (MQ=255)                                                                                     gggCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCacac                                                  <  5:211256/135‑1 (MQ=255)                                                                                     gggCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCacac                                                  <  6:235371/135‑1 (MQ=255)                                                                                      ggCGGAGCGCGAGCAGCCCGCGCGGTTGCTCGCCACCCAGCTTGGCGCGCGGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCacac                                                  <  7:146869/134‑1 (MQ=255)                                                                                      ggCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCTGCTGGGCGCGCGGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCacac                                                  <  3:327618/134‑1 (MQ=255)                                                                                      ggCGGAGCGCGAGCAGACCGCGCGGTTGAGCGACACCCAGAGTGGCGCGCGGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCacac                                                  <  6:265091/134‑1 (MQ=255)                                                                                       ggggcgcgcGAGCATCCCGCGGGGTTGAGCGACACCCAGGTGGGCGCGCGGCGGGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACgg                                                <  2:200932/130‑1 (MQ=255)                                                                                       gCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCacac                                                  <  3:176672/133‑1 (MQ=255)                                                                                              gcgAGCAGACCGCGCGGTGGATCTCCACCCTGCTTGGCGCGCGGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaaca                                        <  7:185500/136‑1 (MQ=255)                                                                                                gAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaacaac                                      <  7:266076/136‑1 (MQ=255)                                                                                                gAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGaacaac                                      <  8:208231/136‑1 (MQ=255)                                                                                                    aGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCAGG‑AGGCCCCGGGGCGGGCAGTCGGCGCCCGGGAGGG‑CCCCCCGGGCCTCCCCCCCCAGGGGGCCCGCCCCCCGCTGGGACAACGcg                                   >  3:426735/1‑134 (MQ=255)                                                                                                           gcgGTTGATCGACCCCCAGCTGGGCGCGCGGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGAc                           <  7:272752/136‑1 (MQ=255)                                                                                                           gcgGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGa                            <  4:180661/135‑1 (MQ=255)                                                                                                                     cacccccAGCTGGGCGCGCGGCCCGCGCGCCGGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCAt                 <  7:172719/132‑1 (MQ=255)                                                                                                                                 ggCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGa      <  7:385908/135‑1 (MQ=255)                                                                                                                                     cgcTGGGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGAcacc  <  4:34968/135‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      TCGAGCTGGCCTTCGACCTGTCCCCGCTCGGCCCCCTGCAAGTGCAGGCACAACTGAGCGAGGGCCGCCTGAGCGGGCAGCTATGGGCGGAGCGCGAGCAGACCGCGCGGTTGATCGACACCCAGCTTGGCGCGCTGCGCGCGCGCCTGC‑TGGCCCGGGGCCTGGAAGTCGGCGACCTGGAGTG‑CCACCCGGGCATCCCGCCGCAGGGGCCACGCACACGGATGGAACAACGCTGGGTGGACGAAACCGCATGAACGATACGACACC  >  KT2440/4912728‑4912994

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