Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	2,696,303	T→C	7.1%	S135P (TCG→CCG)	CJ019_02551 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	2,696,303	0	T	C	7.1%	167.8 / 2.9	57	S135P (TCG→CCG)	CJ019_02551	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/21);  new base C (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCAAAACCGTCGATATCGCCTGCTTTACCTAGCACGTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGCTTGAAGGCGAGGGGG  >  KT2440/2696171‑2696432                                                                                                                                     |                                                                                                                                  cAGCCAAAACCGTCGATATCGCCTGCTTTACCTAGCACGTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCg                                                                                                                                 >  7:300408/1‑135 (MQ=255)            gTCGATATCGCCTGCTTTACCTAGCACGTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCAc                                                                                                                      >  1:404050/1‑135 (MQ=255)              cGATCTCGCCTGCTTACCCTACCACGACACCCCGCACCGCCAGCACCCCGACGCTTCTCTCGGCGGCACTGCAAACCTCCCCCGAAAAGGAAATGCAGCCTGGCGGGTGCGCGGCCGGGTCGCAGTTGGTCACTGc                                                                                                                   <  7:163141/136‑1 (MQ=255)                atatCGCCTGCTTTACCTAGCACGTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCCCTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTg                                                                                                                    >  4:218398/1‑133 (MQ=255)                atatCGCCTGCTTTACCTAGCACGTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTg                                                                                                                    >  6:157447/1‑133 (MQ=255)                           tttACCTAGCACGTCAACACGCATCGCCCCCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCAc                                                                                                         <  1:48033/133‑1 (MQ=255)                                   gCACGTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTAGCAGTTTGTCACTGCCAGCGGGCACCTGCTGCCGc                                                                                               <  8:447103/135‑1 (MQ=255)                                    cACGTCACCACCCAGCGTCAGCCCCCACACGCCTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCCGCCTGGCGGGTGCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCg                                                                                                <  3:101063/133‑1 (MQ=255)                                       gTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAg                                                                                           >  2:404655/1‑135 (MQ=255)                                       gTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAg                                                                                           >  3:384988/1‑135 (MQ=255)                                       gTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAg                                                                                           >  2:276903/1‑135 (MQ=255)                                       gTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGCCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAg                                                                                           >  4:59290/1‑135 (MQ=255)                                       gTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGCCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAg                                                                                           >  4:53557/1‑135 (MQ=255)                                          aaCCCGCATCGCAAGACCCCACACGCTTCTCTCGGCGGCACTGCAAAACTTCCCATACAAGGAAATGCAGCCTGGCTGGTTAGCGGCCGGTTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTg                                                                                         <  7:169781/134‑1 (MQ=255)                                          aaCACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTgc                                                                                        >  1:292468/1‑135 (MQ=255)                                                   cGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCttt                                                                              <  8:61660/136‑1 (MQ=255)                                                    gTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCCCCTGCTGCCGCTAAGTGCGCAAGCCttt                                                                              >  2:284753/1‑135 (MQ=255)                                                    gTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTg                                                                             >  8:348672/1‑136 (MQ=255)                                                     tCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGt                                                                            <  8:300408/136‑1 (MQ=255)                                                            ccAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCt                                                                      >  2:109124/1‑135 (MQ=255)                                                            ccAGACGCTTCTCTCGCCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGCTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTATGTCCGccgt                                                                     >  7:33062/1‑134 (MQ=255)                                                                aCGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGTAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGc                                                                  <  6:165126/135‑1 (MQ=255)                                                                 cGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCGTTGTCCGCCTTTGCa                                                                 >  4:106916/1‑135 (MQ=255)                                                                  gCTTCTCTCGGCGGCACTGCATAACTTCCACTAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCag                                                                <  5:380833/135‑1 (MQ=255)                                                                      ctctcGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCCGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCagag                                                              >  5:49366/1‑133 (MQ=255)                                                                      ctctcGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGAc                                                            >  3:118681/1‑135 (MQ=255)                                                                        ctcGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACAt                                                          >  3:176667/1‑135 (MQ=255)                                                                         tcGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCTGGTCGCAGTTTGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCGTTGCAGAGACAt                                                          >  2:249646/1‑134 (MQ=255)                                                                          cggcgtCACTCCATAACTTCAACTAAAAGTAAATGCTGTCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACAtt                                                         <  6:264149/134‑1 (MQ=255)                                                                           ggcggcACTGCATAACTTCCACGAAACGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACAtt                                                         <  8:132603/133‑1 (MQ=255)                                                                           ggcggcACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACAtt                                                         <  7:224311/133‑1 (MQ=255)                                                                              ggcACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGatat                                                    <  5:376596/135‑1 (MQ=255)                                                                              ggcACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGCCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCGGAGACATTGatat                                                    <  3:53557/135‑1 (MQ=255)                                                                              ggcACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGCCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCGGAGACATTGatat                                                    <  3:59290/135‑1 (MQ=255)                                                                                 aCTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTg                                                  >  4:48569/1‑134 (MQ=255)                                                                                   tGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCa                                               >  2:8636/1‑135 (MQ=255)                                                                                        aaCTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGt                                          >  8:418960/1‑135 (MQ=255)                                                                                                cGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAg                                  >  3:105109/1‑135 (MQ=255)                                                                                                  aaaaGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCa                                               >  4:228009/1‑120 (MQ=255)                                                                                                  aaaaGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCa                                               <  3:228009/120‑1 (MQ=255)                                                                                                  aaaaGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGcc                                >  6:297755/1‑135 (MQ=255)                                                                                                  aaaaGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGcc                                >  1:275566/1‑135 (MQ=255)                                                                                                  aaaaGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCa                               >  7:267073/1‑136 (MQ=255)                                                                                                       gAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACAtt                                                         >  7:139695/1‑105 (MQ=255)                                                                                                       gAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACAtt                                                         <  8:139695/105‑1 (MQ=255)                                                                                                        aaaTGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATAATGGCAGGCGTTATCAAAGCCATTTGc                          >  8:157271/1‑135 (MQ=255)                                                                                                          aTGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCaa                        <  3:48569/135‑1 (MQ=255)                                                                                                              aGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCagag                                                              >  6:306443/1‑93 (MQ=255)                                                                                                              aGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCagag                                                              <  5:306443/93‑1 (MQ=255)                                                                                                                 gtgggtggTTCGCGGCCGGGGCGCAGTTGGTCACTGCCAGCGGGCACGTGCTGCCGCTAAGGGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGAAGGCGTTATCAAAGCGATTTGCAATCGAAga                 <  1:233203/134‑2 (MQ=255)                                                                                                                 ctggctggTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGc                 <  2:43309/135‑1 (MQ=255)                                                                                                                 ctggctggTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGc                 <  1:404655/135‑1 (MQ=255)                                                                                                                   ggctggTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGCt                >  3:233017/1‑134 (MQ=255)                                                                                                                   ggctggTTCGCGGCCGGGTCGCAGTTGGTCACCGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGCt                >  3:343521/1‑134 (MQ=255)                                                                                                                     ctggTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGCtt               <  4:192385/133‑1 (MQ=255)                                                                                                                      gggggcgcgGCCGGGGCGCAGTTGGGCACTGCCAGCGGGCACCGGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGCTTGaa            <  1:284753/130‑1 (MQ=255)                                                                                                                                ccGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGCTTGAAGGCGAggggg  <  5:297755/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  CAGCCAAAACCGTCGATATCGCCTGCTTTACCTAGCACGTCAACACGCATCGTCAGCACCCAGACGCTTCTCTCGGCGGCACTGCATAACTTCCACGAAAAGGAAATGCAGCCTGGCTGGTTCGCGGCCGGGTCGCAGTTGGTCACTGCCAGCGGGCACCTGCTGCCGCTAAGTGCGCAAGCCTTTGTCCGCCTTTGCAGAGACATTGATATTGGCAGGCGTTATCAAAGCCATTTGCAATCGAAGCTTGAAGGCGAGGGGG  >  KT2440/2696171‑2696432

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