Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	252,598	G→A	6.4%	G362D (GGC→GAC)	CJ019_00227 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	252,598	0	G	A	6.4%	202.6 / 3.4	62	G362D (GGC→GAC)	CJ019_00227	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

GTGGCCTGGACACCGGCAAGTTCATCGACGACCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCAGGAGGGCAAGGACATCCGC  >  KT2440/252466‑252731                                                                                                                                     |                                                                                                                                      gTGGCCTGGACACCGGCAAGTTCATCGACGACCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGc                                                                                                                                      >  4:311237/1‑134 (MQ=255) gTGGCCGGGACACCGGCAAGTTCATCGACGACCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGcc                                                                                                                                     <  6:187873/135‑1 (MQ=255)            acCGGCAAGTTCATCGACGACCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCAcc                                                                                                                          >  5:84378/1‑135 (MQ=255)                      tCATCGACGACCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAgcc                                                                                                                <  1:312060/135‑1 (MQ=255)                          cgacgaCCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCt                                                                                                            <  1:102016/135‑1 (MQ=255)                            acgacCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTg                                                                                                           <  8:51945/134‑1 (MQ=255)                            acgacCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTg                                                                                                           <  7:187494/134‑1 (MQ=255)                                         ccGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGt                                                                                             <  6:60521/135‑1 (MQ=255)                                         ccGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGt                                                                                             <  2:66305/135‑1 (MQ=255)                                         ccGGGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGt                                                                                             <  1:328554/135‑1 (MQ=255)                                           gTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACACGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTgc                                                                                           <  2:122475/135‑1 (MQ=255)                                            tGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTgcg                                                                                          >  2:188555/1‑135 (MQ=255)                                             gccgccTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTgcgc                                                                                         <  7:120553/135‑1 (MQ=255)                                             gccgccTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTgcgc                                                                                         <  7:269205/135‑1 (MQ=255)                                               cgccTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCgg                                                                                       >  4:42330/1‑135 (MQ=255)                                               cgccTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCgg                                                                                       >  7:7631/1‑135 (MQ=255)                                                 ccTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGa                                                                                                    <  4:72637/120‑1 (MQ=255)                                                 ccTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGa                                                                                                    >  3:72637/1‑120 (MQ=255)                                                 ccTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGt                                                                                      >  4:324375/1‑134 (MQ=255)                                                 ccTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGt                                                                                      >  3:323717/1‑134 (MQ=255)                                                 ccTTCAAAGAGGCAGGCCATGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGt                                                                                      >  3:27962/1‑134 (MQ=255)                                                    tCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGACCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAg                                                                                  >  4:19774/1‑135 (MQ=255)                                                         gagGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCAGTCAAGGCGGATGACGCGCTGACGGGCCAACATATCACCGTCTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             <  7:55485/135‑1 (MQ=255)                                                         gagGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             <  3:324375/135‑1 (MQ=255)                                                         gagGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             <  5:334622/135‑1 (MQ=255)                                                         gagGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGACCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             <  3:19774/135‑1 (MQ=255)                                                         gagGCAGGCCGTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             <  5:37577/135‑1 (MQ=255)                                                            gcaggcCTTGACTACGAGAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCTGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTc                                                                          <  2:323480/135‑1 (MQ=255)                                                               ggcCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCg                                                                        >  6:301905/1‑134 (MQ=255)                                                                gcCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGcc                                                                      >  6:62221/1‑135 (MQ=255)                                                                gcCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGACCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             <  7:291720/128‑1 (MQ=255)                                                                gcCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGACCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             >  8:291720/1‑128 (MQ=255)                                                                 ccTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCa                                                                     >  6:225407/1‑135 (MQ=255)                                                                        tACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACt                                                              >  4:124660/1‑135 (MQ=255)                                                                             cAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAgcc                                                                                                                >  8:166505/1‑80 (MQ=255)                                                                             cAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAgcc                                                                                                                <  7:166505/80‑1 (MQ=255)                                                                             cAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             >  4:255591/1‑115 (MQ=255)                                                                             cAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAg                                                                             <  3:255591/115‑1 (MQ=255)                                                                             cAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGcc                                                         >  8:40238/1‑135 (MQ=255)                                                                               ctgcgcTGAAAGCCTCCGACGTACTGCCACTCACGGCCGATGACGCGCTGCCGGGCCAACCGAGCCGCGACTTCAGCCGCCGGGCGCAGAGCTGGGGGCGCGGGGAAGACAAGGGCCGCCACTACACTTCGCCgg                                                       <  6:217953/133‑1 (MQ=255)                                                                               aGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCgg                                                       <  4:131892/135‑1 (MQ=255)                                                                                ggCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGa                                                      <  7:40238/135‑1 (MQ=255)                                                                                 gcgcTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAg                                                     <  4:27962/135‑1 (MQ=255)                                                                                   gcTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGgc                                                   >  3:267405/1‑135 (MQ=255)                                                                                   gcTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGCAGgc                                                   >  8:51623/1‑135 (MQ=255)                                                                                        aaGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGccc                                              >  4:96080/1‑135 (MQ=255)                                                                                              acgacgTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTg                                        >  2:132854/1‑135 (MQ=255)                                                                                               cgacgTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGc                                       >  1:97465/1‑135 (MQ=255)                                                                                                  cgTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCtt                                    >  8:221759/1‑135 (MQ=255)                                                                                                        gCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGt                                                                                             >  3:307640/1‑72 (MQ=255)                                                                                                        gCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGt                                                                                             <  4:307640/72‑1 (MQ=255)                                                                                                        gCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCa                              >  3:197487/1‑135 (MQ=255)                                                                                                              cAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGa                        >  4:286744/1‑135 (MQ=255)                                                                                                               aaGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTggcg                       >  1:170409/1‑133 (MQ=255)                                                                                                                aGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGc                      <  8:246627/135‑1 (MQ=255)                                                                                                                 ggCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCa                     >  6:3514/1‑135 (MQ=255)                                                                                                                    cGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCagga                  >  7:17074/1‑135 (MQ=255)                                                                                                                         aCGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCAGGAggg               >  1:163996/1‑133 (MQ=255)                                                                                                                         aCGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCAGGAggg               >  3:261428/1‑133 (MQ=255)                                                                                                                                gACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCAGGAGGGCAAGGACAt      >  4:299044/1‑135 (MQ=255)                                                                                                                                  ccGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCAGGAGGGCAAGGACATcc    <  4:75065/135‑1 (MQ=255)                                                                                                                                   cGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCAGGAGGGCAAGGACATCcg   <  5:3514/135‑1 (MQ=255)                                                                                                                                    ggCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGACGACAAGGTCCGCCACTCCACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCAGGAGGGCAAGGACCTCcgc  >  1:66510/1‑135 (MQ=255)                                                                                                                                     |                                                                                                                                      GTGGCCTGGACACCGGCAAGTTCATCGACGACCAGTACATCCGTGCCGCCTTCAAAGAGGCAGGCCTTGACTACGACAAGGCGCTGAAAGACTACGACGTACTGCCACTCAAGGCCGATGACGCGCTGACCGGCCAACCGATCACCGACTTCAGCCGCCTGGCGCAGATCTGGGTGCGCGGTGAAGACAAGGTCCGCCACTACACTTCGCCGGAGGCCGCCCTTGGTGCCTTGGCCCAGCTGGAGCAGGAGGGCAAGGACATCCGC  >  KT2440/252466‑252731

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