Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	4,825,420	T→A	7.9%	C35S (TGC→AGC)	CJ019_04521 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	4,825,420	0	T	A	7.9%	155.0 / 4.7	51	C35S (TGC→AGC)	CJ019_04521	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/28);  new base A (2/2);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CTCTCCCCTCTGTGACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGCTACCTGTTCCAGATCGTGCGCAACCTGGCG  >  KT2440/4825290‑4825542                                                                                                                                   |                                                                                                                           ctctCCCCTCTGTGACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCg                                                                                                                        >  3:150358/1‑135 (MQ=255)  tctcCCCTCTGTGACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGt                                                                                                                       <  8:294720/135‑1 (MQ=255)  tctcCCCTCTGTGACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATAGCCGCCCGTATCACAGGCTGccct                                                                                                                       >  4:208882/1‑133 (MQ=255)       ccTCTGTGACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTc                                                                                                                     <  2:157496/132‑1 (MQ=255)       ccTCTGTGACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTc                                                                                                                     <  3:234889/132‑1 (MQ=255)              gACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTg                                                                                                             <  8:291643/133‑1 (MQ=255)              gACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGa                                                                                                            <  8:84447/134‑1 (MQ=255)              gACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGa                                                                                                            >  7:291643/1‑134 (MQ=255)              gACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGa                                                                                                            <  3:258420/134‑1 (MQ=255)              gACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGa                                                                                                            <  8:34901/134‑1 (MQ=255)              gACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAg                                                                                                           <  7:5852/135‑1 (MQ=255)               ccacggTGCTTTCCATGGCGGAACAACTATCCACACGTACGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAgg                                                                                                          <  3:204318/131‑1 (MQ=255)                cAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAgg                                                                                                          <  3:238823/134‑1 (MQ=255)                cAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAgg                                                                                                          <  4:176406/134‑1 (MQ=255)                cAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAgg                                                                                                          <  4:292220/134‑1 (MQ=255)                cAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGAGATTCACCATTACTGCAGGCCTTCGGCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAgg                                                                                                          <  5:37038/134‑1 (MQ=255)                   ggTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGt                                                                                                      >  6:88487/1‑135 (MQ=255)                       aTTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGa                                                                                                            >  7:103181/1‑125 (MQ=255)                       aTTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGa                                                                                                            <  8:103181/125‑1 (MQ=255)                         ttCCATGGCTGAACAACTATCCACAAGTAAGTGCGATTCTCCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGTCTGCCGTTCACGCGCTGAGGACGTTGTAc                                                                                                 <  1:120159/134‑1 (MQ=255)                                 cGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCTACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCt                                                                                        >  7:146179/1‑135 (MQ=255)                                   gAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCttt                                                                                      <  3:131932/135‑1 (MQ=255)                                         cTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCgg                                                                                >  3:88934/1‑135 (MQ=255)                                             cctcaAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATACTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCa                                                                            <  8:114162/132‑1 (MQ=255)                                               acaAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGc                                                                          >  4:54803/1‑135 (MQ=255)                                                 aagtaagtGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACGGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCAc                                                                        >  6:266727/1‑135 (MQ=255)                                                      agtGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGccc                                                                   >  8:39819/1‑135 (MQ=255)                                                         gCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCttt                                                                                      >  7:176865/1‑113 (MQ=255)                                                         gCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCttt                                                                                      <  8:176865/113‑1 (MQ=255)                                                                aCCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCAc                                                         >  8:12280/1‑135 (MQ=255)                                                                 ccATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCAcc                                                        <  5:110653/135‑1 (MQ=255)                                                                 ccATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCAcc                                                        <  4:231378/135‑1 (MQ=255)                                                                               cTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGa                                                                                                            >  4:121549/1‑69 (MQ=255)                                                                               cTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGa                                                                                                            <  3:121549/69‑1 (MQ=255)                                                                               cTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAgg                                                                                                          >  7:299687/1‑71 (MQ=255)                                                                               cTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAgg                                                                                                          <  8:299687/71‑1 (MQ=255)                                                                                ttCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGcc                                         <  1:35446/135‑1 (MQ=255)                                                                                     cgACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGc                                     >  5:283651/1‑134 (MQ=255)                                                                                     cgACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCt                                    >  3:272613/1‑135 (MQ=255)                                                                                         aaCCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGc                                <  5:183197/135‑1 (MQ=255)                                                                                         aaCCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCAACGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGc                                <  1:44152/135‑1 (MQ=255)                                                                                                      tGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCCCCTCATCGTTCAAGGCCCAGCTGAGCTACCTGtt                        >  1:2867/1‑130 (MQ=255)                                                                                                         tCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGCTACCTGTTCCAGATc                 <  6:10316/134‑1 (MQ=255)                                                                                                         tCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGCTACCTGTTCCAGATCg                >  5:294001/1‑135 (MQ=255)                                                                                                           aaGATCGCCGCCCGTATCACAGGCAGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCAcc                                                        >  3:112188/1‑93 (MQ=255)                                                                                                           aaGATCGCCGCCCGTATCACAGGCAGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCAcc                                                        >  3:107490/1‑93 (MQ=255)                                                                                                           aaGATCGCCGCCCGTATCACAGGCAGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCAcc                                                        <  4:107490/93‑1 (MQ=255)                                                                                                           aaGATCGCCGCCCGTATCACAGGCAGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCAcc                                                        <  4:112188/93‑1 (MQ=255)                                                                                                                cgccgcCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGCTACCTGTTCCAGATCGTGCGCaa         >  7:147587/1‑135 (MQ=255)                                                                                                                  ccgccCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGCTACCTGTTCCAGATCGTGCGCAAcc       <  8:147587/135‑1 (MQ=255)                                                                                                                       cggaTCACAGGCTGCCGTTCACGCGCTGAGGACGTGGGACAGGATGCCTTTTTCCGGCTGAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGCTACCTGTTCCAGATCGAGCGCAACCTGGCg  <  4:54732/132‑1 (MQ=255)                                                                                                                                   |                                                                                                                           CTCTCCCCTCTGTGACAAGGTGATTTCCATGGCGGAACAACTATCCACAAGTAAGTGCGATTCACCATTACTGCAGGCCTTCGTCGACAACCGCAGCATCCTGGTCAAGATCGCCGCCCGTATCACAGGCTGCCGTTCACGCGCTGAGGACGTGGTACAGGATGCCTTTTTCCGGCTCAGCACCGCCCCGCAGATCACCTCATCGTTCAAGGCCCAGCTGAGCTACCTGTTCCAGATCGTGCGCAACCTGGCG  >  KT2440/4825290‑4825542

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