Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	5,685,825	G→C	36.5%	P122P (CCG→CCC)	CJ019_05336 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	5,685,825	0	G	C	36.5%	30.6 / 5.8	37	P122P (CCG→CCC)	CJ019_05336	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/14);  new base C (9/4);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 8.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTCCAGCCAGAT  >  KT2440/5685693‑5685956                                                                                                                                     |                                                                                                                                    gCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCCCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCACCCCCGa                                                                                                                                    >  5:38445/1‑134 (MQ=255)             cACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCCGGCCCGAGGTGGTCCCCCCCCCGCCCCCCCAAAAGTCGGCact                                                                                                                       >  5:234547/1‑133 (MQ=255)             cACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCCCCCCCCAAAAACTCCGccgg                                                                                                                       >  2:195509/1‑132 (MQ=255)                     aCAAGACAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTGCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTgac                                                                                                                 <  2:32552/133‑1 (MQ=255)                     aCAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCCAAAAAGTCGGCCGGGGTgtcga                                                                                                               >  3:105792/1‑135 (MQ=255)                       aaGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCCCCGCCTGCGGCCACGCCCGCGGTGGTGCCACCCCCGACCCCCCAAAATTCCGGAGTGGTgacgccc                                                                                                             >  8:241856/1‑132 (MQ=255)                       aaGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCCCCCCCCAAAAACTCGGCCGTGGGGtcgacg                                                                                                             >  5:305362/1‑135 (MQ=255)                        aGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACgg                                                                                                            <  8:117571/135‑1 (MQ=255)                             gCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGccc                                                                                                       <  3:108116/135‑1 (MQ=255)                             gCGGGGCCCAAGACCACGGAGCTGGCGCCGTTCCAGCACAGCAAGATCAACAAGATCCCCCCCCCGCCGGCGGCCAGGCCCGCGGGGGTCCCACCCCCACCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGccc                                                                                                       <  2:308070/135‑1 (MQ=255)                                 tGCCCAAGACCACTGCGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGAGGCCCGGCACGAGGTGGTCCCACCCCCTCACCCGAAAAATAATGCAATGGTGACGACGGCGccaca                                                                                                     >  7:52566/1‑130 (MQ=255)                                      aaGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCCAAAAACTAGGCAGTGGTGACGCCCGCGCCCCAGCCGca                                                                                               >  4:125981/1‑134 (MQ=255)                                              tGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCACCCCCCAAAAAGTCGGGAGTGGTGACGACCGCGCCCCAGCCGCCCAAGGGCgg                                                                                      >  3:153958/1‑134 (MQ=255)                                                        ccGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAAccc                                                                            <  7:55268/135‑1 (MQ=255)                                                             ccAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCGCAAAAGTCGGCAGTGGTGGCGGCGGCGCCCCAGCCGCCCCAGGTCGCGCCCCACCCCAAgg                                                                        >  1:325481/1‑134 (MQ=255)                                                             ccAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCACCCCCGAAAAAGGCGGCAGTGGTGACGCAGGCGCCCCAGCCGCCCCAGGTCGAGCCCCAACCCAAgg                                                                        >  5:72910/1‑134 (MQ=255)                                                                    aGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAg                                                                <  2:241219/135‑1 (MQ=255)                                                                         gATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCCAAAAAGTCGGCAGTGGTGACGACCGCGCCCAAGCCGCACCAGGTCGAGCCCAAACCCCACGAAAGCAACGCGCa                                                           >  8:193358/1‑135 (MQ=255)                                                                         gATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCGAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCa                                                           <  7:193089/135‑1 (MQ=255)                                                                         gATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCA‑CCCCTCCCCCGAAAAAGTCGGCAGTGGTGCCGACGGCGCCCAAGCCGCACAAGGTCGCGCCCAAACCCCAGGAAAGCAAGgcgccg                                                          >  8:234946/1‑133 (MQ=255)                                                                                      aCCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCGCAAAAGTCGGCAGTGGGGCCGACGGCGCCCCAGCCGCACAAGGTCGCGCCCCCACCCCAGGGAAGCCAGGCGCCGCCCCagcccgcc                                              >  1:52173/1‑135 (MQ=255)                                                                                        cccACCGCCTGCGGCCAGGCCCGAGGTGGTCCCCCCCCCGCCCCCCCAAAACTCGGCCGTGGGGGCGGCCGCGCCCCCGCCCCCCCACGGCGCGCCCCACCCCAACGGAAGCCACGCGCCGCCCCagcccccccc                                            >  8:134933/1‑129 (MQ=255)                                                                                             ccccgccggcccggccCGCGGGGGCCCCCCCCCCCCCCCCCAAAACGCCGGCCGGGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGa                                       <  3:232408/129‑1 (MQ=255)                                                                                             cGCCTGCGGCCCGGCCGGAGGGGGCCCCACCCCCACCCCCTAAAACGTCGGCTGGGGTGACGACGGCGCCCAAGCCGCACAAGGGCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGa                                       <  2:83016/135‑1 (MQ=255)                                                                                             cGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGa                                       <  7:33349/135‑1 (MQ=255)                                                                                             cGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCCAAAAAGGCGGCAGTGGGGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGa                                       <  5:261032/135‑1 (MQ=255)                                                                                               ccTGCGGCCAGGCCCGCGGTGGACCCACCCCCTCCCCCGAAAAAGGCGGCCGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACg                                     <  3:217656/135‑1 (MQ=255)                                                                                               ccTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACg                                     <  7:175375/135‑1 (MQ=255)                                                                                                     gCCAGGCCCGAGGTGGTCCCCCCCCCTCCCCCGCAAAAGTCGGCAGTGGTGACGCCCGCGCCCCACCCCCCCAACGTCGCGCCCCCACCCCACGAAAGCCAGGCGCCCCCCCCGCCCGCCCCGCCCCCGCCgg                                 >  1:168969/1‑133 (MQ=255)                                                                                                       cgggCCGGTGGGGGCCCCACCCCCCCCCCCCAAAACGGGGCCCGGGGGGACGAGGGCGCCCACGCCGCACAAGGGCGAGCCCAAACCCAAGCAAACCACGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACtt                             <  6:201607/133‑1 (MQ=255)                                                                                                       cAGGCCCTCGGTGGTCCCACCCCCTCCCCCGAAAACGGCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACtt                             <  1:310348/135‑1 (MQ=255)                                                                                                       cAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACtt                             <  3:246597/135‑1 (MQ=255)                                                                                                                ggtgggCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTc                    <  6:305362/135‑1 (MQ=255)                                                                                                                 gtggtCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCg                   <  7:76956/135‑1 (MQ=255)                                                                                                                 gtggaCCCACCCCCCCCCCCAAATAGGCCGCCGGTGGGCCCGCCGGCGCCAAAGCCGCACACGGGCGACCCCAACCCCAAGGAAAGAACGCCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCg                   <  6:38445/135‑1 (MQ=255)                                                                                                                 ggggacccacccccaccCCCCAAAACGGCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCg                   <  1:93748/130‑1 (MQ=255)                                                                                                                            ccccccccccAACACGCCGGCCGGGGGGACGACGGCGCCCAAGCCGCACACGGTCCACCCAAAACCCAAGCAAAGCAAGCCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTccag        <  5:307750/135‑1 (MQ=255)                                                                                                                                  cccGAAAAAGTCGGCAGTGGTGACGCCGGCGCCCAAGCCGCACAAGGTCGCGCCCCAACCCAACGCACTCAACGCGCAGCCCACGCCCGCACCGCCGCCGCCCGCCTTCGACAGCTCGCACCGGTccagcccgca  >  8:301108/1‑131 (MQ=255)                                                                                                                                     |                                                                                                                                    GCAAGGCAGCGGCACCCTGGACAAGAAAGCGGTGCCCAAGACCACTGAGCTGGCGCCGTTCCAGGACAGCAAGATCAACAAGATCACCCCACCGCCTGCGGCCAGGCCCGAGGTGGTCCCACCCCCTCCCCCGAAAAAGTCGGCAGTGGTGACGACGGCGCCCAAGCCGCACAAGGTCGAGCCCAAACCCAAGGAAAGCAAGGCGCAGCCCAAGCCAGCCGCGCCGACGCCGGACTTCGACAGCTCGCAGCTGTCCAGCCAGAT  >  KT2440/5685693‑5685956

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