Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	2,553,519	T→G	19.0%	W121G (TGG→GGG)	CJ019_02406 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	2,553,519	0	T	G	19.0%	78.6 / 2.2	50	W121G (TGG→GGG)	CJ019_02406	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base G (4/7);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCAC  >  KT2440/2553388‑2553651                                                                                                                                    |                                                                                                                                      gTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCctggc                                                                                                                                    <  7:20701/135‑1 (MQ=255) gTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCctggc                                                                                                                                    <  3:367999/135‑1 (MQ=255)           tgtgtCGTTGAAGCAACTACAAGCCGAGCAACCGCCGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGGGCGGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCtt                                                                                                                          <  6:160365/135‑1 (MQ=255)                cGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGCCa                                                                                                                     >  4:275172/1‑135 (MQ=255)                cGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCATCGACa                                                                                                                     >  8:307677/1‑135 (MQ=255)                       gCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACCGGCCggg                                                                                                              >  7:310512/1‑135 (MQ=255)                       gCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCCGGCGTCGACAGGCCggg                                                                                                              >  6:312322/1‑135 (MQ=255)                            tACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGAct                                                                                                         <  6:99620/135‑1 (MQ=255)                                 gCCGAGCAACCGCAGGCGTTGGCCGAATGGGTGCAGGATCCGGCCAGTGCCGTGCATGGGGGTGAGTCGTTCGCGGCGCTATGCCAGCGCCTGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGActggctg                                                                                                    <  5:306554/135‑1 (MQ=255)                                  ccGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGActggctgg                                                                                                   >  7:241250/1‑135 (MQ=255)                                      gCAACCGCAGGCGTTGGCCGATTGGCTGAAGGAACCGGCCAGTGCCGTGCATGGTGGGGAGTCGTTCGGTGCGCTATGCCAGCGCCTGGCCGCCTGGCGGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGt                                                                                              <  4:159280/135‑1 (MQ=255)                                        aaCCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCCCGCCGCCGGGCTGGCGGCCGTCGGCCGGGCGGGGGGCTGGCGGGC‑GGTgg                                                                                            >  5:173909/1‑134 (MQ=255)                                         aCCGCTGGGGTTGGCCGATTGGCTGCAGTATCCGGCCAGTGCCGTGCTTGGGGGTGAGGCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCGGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGac                                                                                           <  3:365311/135‑1 (MQ=255)                                         aCCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATTCCAGCGCCCGGCCGCCTGGCTGGCCGCCGTCGACAGGCCGGGTGACAGGCTGGC‑GGTGac                                                                                           >  7:229858/1‑135 (MQ=255)                                              aGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGacccacc                                                                                      <  3:256092/135‑1 (MQ=255)                                                  gTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCGATGCCAGCGCCCGGCCGCCGGGCTGGCGGCCGTGGCCCGGGCGGGGGGCGGGCTGGG‑GGGGGCCCCCCCgct                                                                                  >  1:308025/1‑133 (MQ=255)                                                      gCCGATGGGCTGCAGGAGCGGGCCAGTGCGGGGCTTGGTGGGGAGTCGTTCGCTGCGCTATGCCAGGGCCAGGCCGCGTGGGTGGGGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTc                                                                              <  8:144353/135‑1 (MQ=255)                                                        cGAAGGGCTGCAGGATCCGGCCAGTGCCGGGCATGGTGGGGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCAt                                                                            <  4:155291/135‑1 (MQ=255)                                                           ctGGCGGCAGGAGCCGGCCAGGGCCGTGCAGGGGGGTGAGTCGTTCGCGGCGCTATGCCAGCGCCAGGCCGCCGGGCGGGCGGCCGTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATgcg                                                                         <  4:340828/134‑1 (MQ=255)                                                           cggggtGCAGGATCCGGCCCGGGCCGTGCTGGGCGGGTAGGCGTGCGCGGCGCTAGGCCCGCCCCTGGCCGCGGGGGGGGGGGCCCTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATgcg                                                                         <  6:37892/130‑1 (MQ=255)                                                             ggCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCCCCCg                                                                                    >  6:184777/1‑122 (MQ=255)                                                             ggCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCg                                                                                    <  5:184777/122‑1 (MQ=255)                                                                tGCAGGATCCGGCCAGTGCCGTGCATGGTGGGGAGGCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCa                                                                    <  3:91385/135‑1 (MQ=255)                                                                tGCAGGAGCCGGCCAGTGCCGTGCAGGGTGGGGAGGCGTCCGCTGCGCTAGGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCa                                                                    <  3:360430/135‑1 (MQ=255)                                                                   aGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTAGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGc                                                                 >  7:316164/1‑135 (MQ=255)                                                                      agcgggCCCGTGCCGGGCTTGGGGGTGTGGCTGGGGCTGCGCTAGGCCCGCGCCTGGCGGCCGGGCGGGGGGCGGTCGACAGGCCGGGGGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTgg                                                              <  1:178820/131‑1 (MQ=255)                                                                            ccAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCCGGCCGCCTGGCGGGCGGCCTTCGACCGGCCGGGTGACTGGCTGGC‑GGTGCCCCCCCCGATGGGCCTGCGCGCCATGCTGGTACAgg                                                        >  5:180028/1‑135 (MQ=255)                                                                             ctgTGCCGTGCTTGGTGGTGAGTCGTTCGCTGCGCTAGGCCAGCGCCAGGCGGCCGGGCGGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGt                                                       <  5:66875/133‑1 (MQ=255)                                                                                 gCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCtt                                                   <  8:273708/135‑1 (MQ=255)                                                                                  ccGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCCGGCCGCCTGGCTGGCGGCCGTGCGCAGGGCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGGCATGCGCGCCAAGCTGGTACAgggggatg                                                  >  8:134603/1‑129 (MQ=255)                                                                                   cGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCCGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTgg                                                 >  3:178948/1‑135 (MQ=255)                                                                                    gTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTggg                                                <  8:316164/135‑1 (MQ=255)                                                                                        aTGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGCCAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTc                                            >  8:260206/1‑135 (MQ=255)                                                                                        aTGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCGGGCTGGCGGCCGTCGCCAGGCCGGGTGACTGGCTGGC‑GGTGACCCCCCCGATGGTCATGCGCGCCCTGCTGGTACAGGTGCTTGGgtgtg                                            >  3:68258/1‑134 (MQ=255)                                                                                         tggtggGGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCGTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTcc                                           <  2:150284/135‑1 (MQ=255)                                                                                         tgggggTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGGTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTcc                                           <  4:336748/135‑1 (MQ=255)                                                                                                  gTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCGGGCTGGCGGCCGTCGCCAGGCCGGGGGGCTGGCGGGCGGGGGCCCCCCCCCAGGGACAGCCGCGCCATGGGGGAACCGGTGGTGGGGGGTCGGAtgtgtg                                   >  2:343941/1‑135 (MQ=255)                                                                                                   tCGTTCGCTGCGCTATGCCAGGGCCAGGCCGCCTGGGTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGcc                                 <  8:57617/135‑1 (MQ=255)                                                                                                        cggtGCGCTAGGCCAGCGCCAGGCCGCCTGGCGGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCa                            <  8:264234/132‑1 (MQ=255)                                                                                                         gCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAg                           >  7:925/1‑135 (MQ=255)                                                                                                           tGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCg                         <  8:241250/135‑1 (MQ=255)                                                                                                                    ccAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACg                 <  8:229858/134‑1 (MQ=255)                                                                                                                     cgcggcCTGGCCGCCGGGGGGGGGGCCGTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGGCCGATGAGTGCCGGGCAGCGTATCGACGtg               <  6:49963/131‑1 (MQ=255)                                                                                                                     cAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGtg               >  7:135718/1‑135 (MQ=255)                                                                                                                       gcgcCAGGCCGCCGGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGtgct             <  3:84447/135‑1 (MQ=255)                                                                                                                          ccAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCggg                                                                                                              >  3:80722/1‑36 (MQ=255)                                                                                                                          ccAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCggg                                                                                                              <  4:80722/36‑1 (MQ=255)                                                                                                                          ccAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCCCCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGCCGTGCTGcc          >  4:80119/1‑135 (MQ=255)                                                                                                                          ccAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTGATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGGTGcc          >  5:31853/1‑135 (MQ=255)                                                                                                                          ccAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGcc          >  2:20959/1‑135 (MQ=255)                                                                                                                               ccgccTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTg      >  3:93114/1‑134 (MQ=255)                                                                                                                                  ccTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCAc  <  4:178948/135‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      GTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGC‑GGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCAC  >  KT2440/2553388‑2553651

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