Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	3,905,078	A→C	29.2%	E248A (GAG→GCG)	CJ019_03692 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	3,905,078	0	A	C	29.2%	51.5 / 4.0	48	E248A (GAG→GCG)	CJ019_03692	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/7);  new base C (3/11);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 7.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATGGCGCCGGCACCGGGCGTTTGCCGGCGTTGGTGCGTCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCA  >  KT2440/3904946‑3905209                                                                                                                                     |                                                                                                                                    aGGATGGCGCCGGCACCGGGCGTTTGCCGGCGTTGGTGCGTGAGGGTTAAGAGCCCGCGTTGCCGTGGCGGGTGGGTGTGGTGGGCGAGGCGTTCCCCCAGGTGAGCGCGCCGGCGCCGGCGTTTAGGGCGgcgc                                                                                                                                   >  3:204149/1‑135 (MQ=255)     tGGCGCCGGCACCGGGCGTTTGGCGGGGTTGGTGCGTCAGGGGTACGAGCAGGCGTTGCCGTGCCGGGTGGGGGTGGTGCGCGCGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGGGCccgc                                                                                                                               >  5:100926/1‑135 (MQ=255)     tGGCGCCGGCACCGGGCGTTTGCCGGCGTTGGTGCGTCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGCGCAGGCGTTTCGGGCGGAGCccgc                                                                                                                               >  7:243558/1‑135 (MQ=255)     tGGCGCCGGCACCGGGCGTTTGCCGGCGTTGGTGCGTCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTAGGGCGGAGCccgc                                                                                                                               >  7:282106/1‑135 (MQ=255)     tGGCGCCGGCACCGGGCGTTTGCCGGCGTTGGGGCGTCAGGGTTACGCGCAGGCGTTGCCGTGCCGGGTGGGGGTGGTGGGCGAGGCGTTGCGCCAGGTGAGCGCGCCGGGGCAGGCGGTTGGGGCGGGGccccc                                                                                                                               >  1:279110/1‑133 (MQ=255)     tGGCGCCGGCACCGGGCGTTTGCCGGCGTTGGGGCGGCAGGGTTACGAGCAGGCGTTGCCGTGGCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGCGCGCGCCGGCGCAGGCGTTTCGGGCGGAGccccc                                                                                                                               >  4:282670/1‑133 (MQ=255)                    gCGTTTGCCGGCGTTGGTGCGTCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGGGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGCGCCCCCCGCGCCCCggggggg                                                                                                                >  2:67492/1‑135 (MQ=255)                     cGTTTGGCGGCGTTGGTGGGGCAGGGTTTCGAGCAGGGGGTGCCGGGCCGGGTGGGTGTGGTGGGCGAGGCGTGCCGCCAGGGGAGCGCGCCGGCGCAGGCGTTTCGGGCGGAGCCCGCCGCGCCCCgggtgg                                                                                                                 >  1:146056/1‑133 (MQ=255)                     cGTTTGCCGGCGGTGGGGGGGGAGGGGTTCGCGCCGGCGGTGGCGGGGCGGGGGGGTGTGGGGCGCGCGGCGGTCCGCCCGGGGCGCGCGCCGGAGCAGGCGTTTGCGGCGGCGGCCGCCGCCCCCCgggtgggg                                                                                                               >  2:351195/1‑134 (MQ=255)                           ccGGCGTTGGTGCGTCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGACGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGCGCCCGCCCCGCCCACGGTGGGTAGCggg                                                                                                         >  3:179020/1‑135 (MQ=255)                                                     gCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGCGCCCAGGGTGGGTAGCGGGTGCCAAGCCGAGGCAGTGTTTGCCtt                                                                               >  3:377056/1‑135 (MQ=255)                                                         ggggggccgggcggGGGGGGGGGGGGGCGCCCGGCGGCCCCCCTGGGGCGCGCGCCGCCGCCGGGGTGGCGGGCGCCCCCCGCCGAGCCCAGGGGGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTgcgc                                                                           <  7:91577/130‑1 (MQ=255)                                                          cGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCt                                                                          >  3:86922/1‑135 (MQ=255)                                                          cGTTGCCATGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCCGTGTTTGCCTTGCGCt                                                                          >  4:252059/1‑135 (MQ=255)                                                            ttGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGGGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTAc                                                                        >  1:245234/1‑135 (MQ=255)                                                                 gtgccgGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCCGGCGGAGCCCGCCGAGCCCAGGGGGGGTTGCGGGTGGCCTGCCGAGGGCGTGTTTGCCTTGGGGTACCgcgc                                                                   >  1:197677/1‑132 (MQ=255)                                                                    cgggggggggggggggCGCGCGGCGTGCCGCCTGGTGCGCGCGCCGGAGCCGGCGTTGGGGGCGCCGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTg                                                                <  4:86922/133‑1 (MQ=255)                                                                      gggtgggtGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTggg                                                              >  7:372881/1‑135 (MQ=255)                                                                       ggtgggtgTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGCGCCCAGGGTGGGTAGCGGGTGCCCTGCCGAGGCAGTGTTTGACTTGCGCTACCGGCCCTGGGc                                                             >  7:190157/1‑135 (MQ=255)                                                                        gggggggTGGGGCGCGCGGGGTGCCGCCTGGTCCGCGCGCGGGAGCTGGCGTGGCGGGCGCCGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGcc                                                            <  8:110391/133‑1 (MQ=255)                                                                            gtgtGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGACGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCggg                                                        >  8:69201/1‑135 (MQ=255)                                                                             gggggggcgcgcGGCGTCCCGCACGGTGGGCGCGCCGGCGCCGGCGGTGCGGCCGCACCCGGCCGAGCCCAGGGTGGGGAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCgggg                                                       <  4:133477/129‑1 (MQ=255)                                                                              gtggtgCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGCGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGt                                                      >  5:386998/1‑135 (MQ=255)                                                                               gggggcgcgcGGCGCGCCGCCTGGTGCGAGCGCCGCAGCCGGCGTTGCGGGCGGTGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCACGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTa                                                     <  6:186445/131‑1 (MQ=255)                                                                                  tgCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCCAGCCCCGGGTGGGTTGCGGGTGCCCTGCCGCTGCAGGGTTTGCCGTTCGCTAACCGccccgggccccggggttag                                                  >  6:216563/1‑130 (MQ=255)                                                                                         ggCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGCGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCcaca                                           >  3:127127/1‑135 (MQ=255)                                                                                         ggCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCTGGGTGGGTAGGGGGTGGCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATGTTGCgact                                           >  3:112820/1‑134 (MQ=255)                                                                                         ggCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCcaca                                           <  8:372881/135‑1 (MQ=255)                                                                                         ggCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTGTCGGGCGGAGCCCGCCGAGCCCCCGGTGGGGTGCGGGTGGCCTGCGGACGGAGTGGTTGCCTTGGGCTACCGGCCCGGGGCCCGGGGGTTCTTGCcact                                           >  5:229156/1‑134 (MQ=255)                                                                                         ggCGTGCGGCCTGGTGCGCGCGCCGTCGCTGGGGGTGCGGGGGGAGCCGGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCcaca                                           <  7:185360/135‑1 (MQ=255)                                                                                             tgccGCCCGGTGCCCGCGCCGGCGCCGGCGGTGCGGGCGGCGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCaa                                       <  7:366847/133‑1 (MQ=255)                                                                                              gccGCCTGGTGCGCGGGCGGGAGCAGGCGTGTCGGGGGGCGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAAt                                      <  8:126847/134‑1 (MQ=255)                                                                                                 gCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGaa                                   >  7:320483/1‑135 (MQ=255)                                                                                                            gcgcGGGCGCCGGGGTTGCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGggcc                        <  4:39398/135‑1 (MQ=255)                                                                                                              gcCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGCGCCCAGGGGGGGTTGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGGATCTTGGCCCAGCAAAGAATTACCAGggcccg                      >  1:354225/1‑133 (MQ=255)                                                                                                              gcCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGGAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGGCACAGCAATGAATTACCCGggccgg                      >  3:154941/1‑135 (MQ=255)                                                                                                                   aGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCctct                 >  4:182755/1‑135 (MQ=255)                                                                                                                       ggCGTTGCGGGCGGCGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTcgc             <  8:313305/135‑1 (MQ=255)                                                                                                                        gggtggcGGGCGGCCCCGGCCGAGCCCAGGGTGGGTAGCGGGGGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTcgcg            <  4:132293/129‑1 (MQ=255)                                                                                                                        gggTTGCGGGCGGCGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTcgcg            <  3:256798/133‑1 (MQ=255)                                                                                                                        gcggtgcGGGCGGCCCCCGCCGAGCCCAGGGTGGGTAGCGGGGGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTcgcg            <  7:129974/129‑1 (MQ=255)                                                                                                                         ggttggggggggCGCCGGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTcgc             <  8:824/127‑1 (MQ=255)                                                                                                                           tttCGGGCCGAGCCCGCCGCGCCCAGGGTGGGGAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGGCACCGCCATGAATTACCAGGGCCCGCCTCTTCGCGGa          >  1:270947/1‑134 (MQ=255)                                                                                                                                ggCGGAGCCCGCCGCGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAAc     >  8:7725/1‑134 (MQ=255)                                                                                                                                ggCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGCATAAAcc    >  8:132356/1‑135 (MQ=255)                                                                                                                                  cGGCGCCGGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTGTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCca  <  2:354225/135‑1 (MQ=255)                                                                                                                                  cGGATCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCca  <  3:275580/135‑1 (MQ=255)                                                                                                                                  cGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCca  <  7:104409/135‑1 (MQ=255)                                                                                                                                  cGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCca  <  4:204149/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    AGGATGGCGCCGGCACCGGGCGTTTGCCGGCGTTGGTGCGTCAGGGTTACGAGCAGGCGTTGCCGTGCCGGGTGGGTGTGGTGCGCGAGGCGTTCCGCCAGGTGAGCGCGCCGGAGCAGGCGTTTCGGGCGGAGCCCGCCGAGCCCAGGGTGGGTAGCGGGTGCCATGCCGAGGCAGTGTTTGCCTTGCGCTACCGGCCCTGGGCCCGGGGTATCTTGCCACAGCAATGAATTACCAGGGCCGGCCTCTTCGCGGATAAACCCA  >  KT2440/3904946‑3905209

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