Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	1,744,031	G→C	36.8%	A92P (GCC→CCC)	CJ019_01649 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	1,744,031	0	G	C	36.8%	30.4 / 7.8	32	A92P (GCC→CCC)	CJ019_01649	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/3);  new base C (2/10);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 4.52e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCC  >  KT2440/1743909‑1744159                                                                                                                           |                                                                                                                                 ccAATGCGGTTTCCAGCCTGCTGCGTCGCCGAGTCCTGTTCCCCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAAATACACGCACGCCATGATTTACGTTTATGCCCAACAGACCAAAACACCCCACACCGCTg                                                                                                                      >  3:245822/1‑135 (MQ=255) ccAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGCATTACCCGAAGGACTGGATTTAAGTTGATCCGAAACAGACCAAAACATCCGACtacgcg                                                                                                                       >  2:221016/1‑129 (MQ=255) ccAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTTACTTGATCCGAAACATACCAAAAAAGCCCACACCGCt                                                                                                                       >  7:267809/1‑134 (MQ=255) ccAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATGCGAAACAGACCAAAACAGCCGACTCCGCt                                                                                                                       >  2:226628/1‑134 (MQ=255) ccAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGCCTCGGCTg                                                                                                                      >  4:213805/1‑135 (MQ=255) ccAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTACGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGGTGATCCGCAACAGACCAAAACACCCCACTAGGCt                                                                                                                       >  1:215010/1‑134 (MQ=255)          tttCCAGCCGGGTGGGGGGCGGTGTCGTGTGCGGGGACGGGGGGGCCGGGGGTGACTGCGGCGTCAATCCCCCTAAGTACTGTATTTAGGTTGAGCCCAAACACACAAAAACCCCCGACTCGGCTGAAGTGGTcc                                                                                                             <  6:53080/135‑1 (MQ=255)                 ccTGCTGCGTCGCCGTGTCCTGTTCCGCGATGGTGGGGCCCGGGGTGACAGCGGCGTGAATGACCCGAAGGACTGTATTTACGTTTATCCGAAACAGACAAAAACACCCGACTCGGCTGAAGTGGTCCGAATCgg                                                                                                      <  1:221016/135‑1 (MQ=255)                 ccTGCTGCGTCGCCGTGTCCTGTGCCGCGATGGGGGGGCCCGGGGTGACATCGGCGTCAATGACCCGAAGGACTGGATTTAGGTTGATCCGAAACAGACCAAAACACCCGACTCGGCTGAAGTGGTCCGAATCgg                                                                                                      <  1:226628/135‑1 (MQ=255)                      tgcGGCGCGGTGGTCTGTGCCGGCAAGGTGGGGCCGGGTGTGCCTTGGGCGTGAATGCCCCTAAGTACTGGATTTAGGTTTATCCCAAACACACCAAAACACACGACTCGGCTGAAGTGGTCCGAATCGGCTCag                                                                                                 <  8:290323/135‑1 (MQ=255)                           cGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCg                                                                                            <  1:85317/135‑1 (MQ=255)                                          cgcgATGGGGGGGCCGGGGGGGACCTCGGCGTGAATCCCCGGAGGCAGGGTATTTCCGTTGAGCCGAAACAGACAAAACCCCCCGACTCGGCTTAAGTGGTCCGACTCGGCTCAGAGTCGAAACAGACCaaacc                                                                              <  8:110958/134‑3 (MQ=255)                                          cgcgAAGGTGGTGCCCGGGGTGACATCGGCGTGCATGACCCGAAGGACAGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGCCTCGGCGGAACTGGTCCGAATCGGCTCACAGTCGAAACACCCCaaaaac                                                                             >  6:85067/1‑133 (MQ=255)                                          cgcgAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCTGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAAAATAACaaaaac                                                                             >  6:270/1‑133 (MQ=255)                                            cgATGGGGGGGCCGGGGGTCACAGGGGGGTTAATGCCCCTACGGACGGTATTTCCGTTTACCCAAAACAGACAAAAACCCCCGACTCGGCTGAAGTGGTCCTAATCGGCTCAGAGTCGAAACAGACCAAAAccgc                                                                           <  6:335194/135‑1 (MQ=255)                                                          ggggTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACCGCCGCCTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCCAAACCGCCGACTcccggcgtt                                                             >  5:180374/1‑129 (MQ=255)                                                          ggggTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGCAGTGGTCCGAATCGGCTAAGAGTCGAAACAGACCAAAACCGCCGACTCCCGTCttt                                                             >  7:162953/1‑135 (MQ=255)                                                              tGACATCGGCGTGAATGACCCGAAGGACTGGATTTAAGTTGATCCGAAACAGCCCAAAACAGCCGACTCGGCTGAATTGGTCCGCATAGGCGCAGAGTCGAAACCTACCAAAACCCCCCACACCCTTCTTTAttg                                                         >  4:20502/1‑134 (MQ=255)                                                                   gcGGCGTCAACGACCCGACGGACTGGATTTACGTTGAGCCAAAACCCCCAAAACCCGCGGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaagaa                                                   <  3:64819/135‑1 (MQ=255)                                                                    cGGCTTGAAGCACCCTAAGCACGGTATTTACGTTGACCCTAAACACCCAAAACCACCCGACTCGGCTGAAGTGGTCCGAATCGGCTCATAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTaagaaa                                                  <  4:245822/136‑1 (MQ=255)                                                                         tGAATGACCCGAAGGCCTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGCATCGGCTCAGAGTCGAAACAGACCAAACCCGCCGCCGCCCTTCTTTATTCTaaggaagacgc                                             >  3:93795/1‑133 (MQ=255)                                                                         tGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAagaaagaac                                              >  7:289182/1‑134 (MQ=255)                                                                         tGAATGACCCGAAGGACTGGATTTAAGTTGATCCGAAACAGACCAAAACAGCCCACTCGGCTGAATTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCCCCTCCCGTGGTTTTTCTaagaaatccg                                              >  6:328700/1‑131 (MQ=255)                                                                             tcccccGAAGTACTGCATTAACTTTCACCCGAAACACACAAAACCCCCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAAccccc                                         <  3:257870/133‑1 (MQ=255)                                                                                    aaGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTAGGt                                    >  8:139504/1‑134 (MQ=255)                                                                                         atgtaTTTACGTGTATCCCAAACACACCAAACCACCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGCCCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATct                              <  5:126564/131‑1 (MQ=255)                                                                                             aTTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAATCCCCCTATGTAAATcttcgt                          >  6:372619/1‑133 (MQ=255)                                                                                                  cGTGGACCCGAAACAGACCAAAACCCCCGACTCGGCGGAAGTGGTCCGAATCGGCTCAGAGGCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGa                     <  6:270838/135‑1 (MQ=255)                                                                                                         cccaaaCCGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGtt              <  8:96058/132‑1 (MQ=255)                                                                                                         ccGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAGCGGCCTATGTTAATATTCGGTCGGAAGATGtg              >  3:46403/1‑134 (MQ=255)                                                                                                                     aaaaCAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGCCTCCCTTCTTTATTCTAAGAAAGACCGCGCCTCTGTAAATCTTCCGTGGGAAGGTGTTACGTGcaccccc  >  7:309007/1‑135 (MQ=255)                                                                                                                     aaaaCAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAAGCCCCTAGGTAAATCTTCCTTCGGAAGCTGTTACATGaccccac  >  2:298402/1‑133 (MQ=255)                                                                                                                           |                                                                                                                                 CCAAGGCGGTTTCCAGCCTGCTGCGTCGCCGTGTCCTGTTCCGCGAAGGTGGTGCCCGGGGTGACATCGGCGTGAATGACCCGAAGGACTGGATTTACGTTGATCCGAAACAGACCAAAACAGCCGACTCGGCTGAAGTGGTCCGAATCGGCTCAGAGTCGAAACAGACCAAAACCGCCGACTCCCTTCTTTATTCTAAGAAAGAAACCCCCTATGTAAATCTTCCTTCGGAAGATGTTACATGCCCCCCC  >  KT2440/1743909‑1744159

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