Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	3,986,622	T→G	31.2%	D31A (GAT→GCT)	CJ019_03763 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	3,986,622	0	T	G	31.2%	45.3 / 5.3	58	D31A (GAT→GCT)	CJ019_03763	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base G (7/12);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGT  >  KT2440/3986500‑3986752                                                                                                                           |                                                                                                                                   ccATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCGGCGCggcgc                                                                                                                        >  5:389907/1‑132 (MQ=255)  cATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCAc                                                                                                                       >  2:335852/1‑135 (MQ=255)  cATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCAc                                                                                                                       >  4:46145/1‑135 (MQ=255)   aTCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACg                                                                                                                      >  6:276344/1‑135 (MQ=255)   aTCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGGCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGCGCGTGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACg                                                                                                                      <  2:60562/135‑1 (MQ=255)       gCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCCGCTCCGGCCCCGGCCCGTcag                                                                                                                  >  3:27738/1‑133 (MQ=255)       gCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGAGGAAGGGGCCGGCTCCGGCGCGGGCACGTcat                                                                                                                  >  3:32642/1‑133 (MQ=255)               gCCTGCAACTGTGCGGGGTACATGGCATTGGGGGCCAGAATCTGGTCATGGTGCTCCTGGGATAACCCGGGGTCGTTGGCTACGTGCGGGGCAGCGGAGGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTg                                                                                                           <  2:124919/134‑1 (MQ=255)               gCCCGCAACTGATCAGGCGACACGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGACCACACCGGGGCCTGGGCTACGAGCGTGGCAGCGGACGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGc                                                                                                          <  8:256203/135‑1 (MQ=255)               gCCAGAAACTGATCTGGGTCCAAGGCAAGGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGGCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTg                                                                                                           <  6:98085/134‑1 (MQ=255)                 cAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGGGCGGCGCGGGCACGGCTGCGGccgggac                                                                                                        >  6:359260/1‑130 (MQ=255)                 cAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCATGGGAGAACACCGGGTCCTTGGCTACGAGCGAGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCCAGTCTGCCGCCAGccc                                                                                                        >  2:103091/1‑133 (MQ=255)                 cAGCAACTGATCAGCGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGCAGAACACCGGGTCCTTGGCTACGAGCGTGTCAGCGGAAGGCCCCCGCACCGCCCCGGGCCCCTCTTCGCCCAGCAc                                                                                                        >  2:237908/1‑135 (MQ=255)                     aaCTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGCGCGTGGAAGCGGAAGGCGCCCGATCCCGCGCGGGCCCGTCTGCCGCCGGCCccagg                                                                                                    >  6:57250/1‑135 (MQ=255)                       cTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGACCACCCCGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCCGCCCCCGCCCCACAGCGGCCTGGAcccggcc                                                                                                  >  4:350402/1‑135 (MQ=255)                        tGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATGCGGCGCGGGCCCGTCTGCGGCCGGCAcccggccc                                                                                                 >  4:104182/1‑134 (MQ=255)                               ggaaCAAGCCACGGGGGGCCCCCATCTGGTCCAGTGGCTCCTGGCACAACACCGGGCCCTTGCCTACGAGCGGGGCAGCGGACGGCGCCGGACCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAg                                                                                          <  1:270893/132‑1 (MQ=255)                                  aCAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCCGGCCCGTCTGCGGCCTGCCCCCCGCCCGGGGCAgcgc                                                                                       >  2:10223/1‑135 (MQ=255)                                       gCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCGGCGCGGGCACGACTGCGGCCGGCACCCCGCCCCGGGCAGCGCGCAgg                                                                                  >  4:170670/1‑135 (MQ=255)                                            gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGCGCCGGCGCGGGCGCGGCGGCGGCCGGCCGCCGGCCCGGGGCGGCGCCCCGGcgcag                                                                             >  1:248241/1‑131 (MQ=255)                                             ggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGGCGGCGGCCGGCGCCCGGCCCGGGGCAGCGCGCCGGCaaa                                                                              >  5:315678/1‑133 (MQ=255)                                              gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGTCAGCGGAAGGCGCCCGAACCGGCGCGGGCACGGCTGCGGCCGGCACCCCGCCCGCGGCCGCGCTCCGGCAAAGCa                                                                           >  2:281200/1‑135 (MQ=255)                                              gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGGCGGCGCGGGCCCGGAGGCGGCCAGCGCCACCCCAGGCGCCGCGcgcacgcacacc                                                                            >  5:129171/1‑130 (MQ=255)                                              gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACCTCTTCGGCCGTCCACAGGCCACGGGCCTCGCTCAGGCAAAGCa                                                                           >  2:292677/1‑135 (MQ=255)                                              gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCCCGTCTGCGGCCTGCCCCCGGCCAGGGGCAGCGCTCCGGCCAAGCa                                                                           >  8:368163/1‑135 (MQ=255)                                                            tCCAGTTGCTCGGGGCACAACCCCGGGTCCTTGGCTACGCGCGGGGCAGCGGACGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGa                                                               <  8:198695/133‑1 (MQ=255)                                                             ccAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCCGCTCCGGCGCGGGCCCGTCTGCGGCCTGCACCCGGCCCGGGGCCGCGCGCAGGCGAACCAGGGGCGCCAGGaaa                                                             >  2:65073/1‑134 (MQ=255)                                                             ccAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCCGATCCGGCGCGGGCACCTCTTAGGGCTGCCCCAGGCCAGGGGCAGCGCGCAGGCAAAGCAGTGCAGACAGGaaa                                                             >  4:280518/1‑134 (MQ=255)                                                                gTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGATGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCaa                                                         <  6:391030/135‑1 (MQ=255)                                                                     tCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAAt                                                    <  4:63474/135‑1 (MQ=255)                                                                        ggggACAACCCGGGGCCCTTGGCTCCGAGCGGCGCAGCGGAGGGCGCCGTAGCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGc                                                 <  1:229825/134‑1 (MQ=255)                                                                          cgagaACCCCGGGTCCTTGGCTCCGAGGGTGGCTGCGCACGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCAt                                               <  3:170670/134‑1 (MQ=255)                                                                            cgcaccccggggccggggCTCCGCGGGGGGCCGGGGACGGCGCGGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATgc                                             <  2:190253/129‑1 (MQ=255)                                                                            ataaCACCGGGGCCTTGGCTAGGAGGGGGGCCGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATgc                                             <  6:194972/133‑1 (MQ=255)                                                                            agaACACCGGGGCCTTGGCTACGCGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATgc                                             <  8:293130/135‑1 (MQ=255)                                                                            acaccgccGGGGCCTGGGCCCCGAGCGTGGCCGCGGAAGGCGCGGGATCCGCCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATgc                                             <  7:307006/129‑1 (MQ=255)                                                                                 acCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCCGCGCGGGGCCGGCTGGGGCCGGCCCCCCGCCCGGGGGCGCGCGCCGGCCACCCACGGCCGCCACGAAAGCAAACAAAGGCCTGCGCCTc                                        >  4:376704/1‑135 (MQ=255)                                                                                     gggccTTGGCTCCGCGCGGGGCCGCGGACGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                    <  5:276344/132‑1 (MQ=255)                                                                                     ggTCCTGGGCTACGCGCGTGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                    <  8:172039/135‑1 (MQ=255)                                                                                          ttGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATccc                               <  4:32642/135‑1 (MQ=255)                                                                                          ttGGCTACGAGCGGGGCAGCGTAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGTGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATccc                               <  6:313178/135‑1 (MQ=255)                                                                                          tgggCTACGCGCGGGGCAGCGGCCGGCGCCGGCGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATccc                               <  7:66463/133‑1 (MQ=255)                                                                                          tgggCTACGAGCGTGGCAGCGGAAGGCGCCGGTTCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATccc                               <  4:27738/133‑1 (MQ=255)                                                                                             gCTCCGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGTCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATcccc                              <  5:163119/133‑1 (MQ=255)                                                                                              cTACGAGCGTGGCAGCGGAAGGCCCCGGCGCCCCCGCCGGCACGTCAGCCCCCAGCACCACGCCACGCGCCCCGCGCATGCAAACCACGGCCGACAACAAAAAAAAAAAAAGCATGCCCCTCGCATATACCCaag                           >  6:197638/1‑133 (MQ=255)                                                                                                    gCGTGGCAGCGGAAGGCGCCGGCTACGGCGCGGGCACGGCTGCGGCCAGCACCAGGCCCCGGGCAGCGCGCAGGCAAACCAGTGCCGACAGGCAAGCAAACAAATGCATGCGCAGAGCGTATACCCCAGAgggtg                     >  6:130134/1‑132 (MQ=255)                                                                                                        ggctggggACGGGGCCGGTTGGGGGGCGGGCCCGTGTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                 <  2:60001/129‑1 (MQ=255)                                                                                                         gctgCGGACCGCGCCGGAGCGGGGGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACcg                <  7:233114/132‑1 (MQ=255)                                                                                                          cAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGTCTGCGGCCTGCACCAGGCCAGGGGCCGCGCTCCCGGAAAGCAGGGCCGACCGGGAAGCAAACAAAGGCATGCGCCTCGCGTATCCCCATGAGGTGAGCCcgc               >  8:10665/1‑135 (MQ=255)                                                                                                            gCGGAAGGCGCCGGCTCCGGCGCGGGCCCGTCTGCGGCCGGCCCCCGGCCCGGGGGCGCGCTCCGGCCAAGCCGTGCCGCCAGGAACGCCAACAATGGCATGCGCATCGCGTATCCCCATGAGGTGAGCCcgcgg             >  2:228118/1‑134 (MQ=255)                                                                                                              ggcgggCGCCGGAGCCGCCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg           <  3:376704/131‑1 (MQ=255)                                                                                                              ggcaGGCGCCGGAGCCGGCGGGGGCACGGCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTAACGCCATGAGGTCAGACGGCGTgg           <  8:342190/132‑1 (MQ=255)                                                                                                              ggATGGCCCGGGAACGGGGGCTGGCACGTCTGCGGCCTGAACCAGGCCTGTGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCTCATAGCTTATCACCATGAGGTCAGACCGCGTgg           <  7:312966/135‑1 (MQ=255)                                                                                                              ggAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg           <  4:133938/135‑1 (MQ=255)                                                                                                                    cgcgGGAGCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAg     <  2:186876/135‑1 (MQ=255)                                                                                                                    cgcCGGATCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAg     <  2:72271/135‑1 (MQ=255)                                                                                                                       cGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGt  <  3:46145/135‑1 (MQ=255)                                                                                                                       cGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGt  <  2:278779/135‑1 (MQ=255)                                                                                                                           |                                                                                                                                   CCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGT  >  KT2440/3986500‑3986752

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