Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	5,853,869	T→G	27.9%	Q561H (CAA→CAC)	CJ019_05481 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	5,853,869	0	T	G	27.9%	73.9 / 5.8	55	Q561H (CAA→CAC)	CJ019_05481	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/36);  new base G (13/2);  total (17/38)
Fisher's exact test for biased strand distribution p-value = 1.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACCGTCAGTTGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCGGCATCGTCCATGTTCAGCACGGCAAT  >  KT2440/5853736‑5854003                                                                                                                                      |                                                                                                                                       aGACCGTCAGTTGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTTGCGCCAGGCtt                                                                                                                                       >  4:374283/1‑135 (MQ=255) aGACCGTCAGTTGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGGTGGCGGCCGGCGTg                                                                                                                                      >  7:12705/1‑136 (MQ=255)  gACCGTCAGTTGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTgg                                                                                                                                     <  8:39529/136‑1 (MQ=255)   cccGTCAGTTGGGTCTCGGCGGCCGACGCCTGAAGGCTGAACACGGCGTAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGCCGGTCGCGGCGCGTTGGCGCCAGGCTTgg                                                                                                                                     <  5:409324/134‑1 (MQ=255)   aCCGTCAGTTGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTgg                                                                                                                                     <  3:347674/135‑1 (MQ=255)      gTCAGTTGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAg                                                                                                                                  <  4:117653/135‑1 (MQ=255)        cAGTTGGGTATCGGCGGCCGAAGCCTGCAGGCGGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAgcg                                                                                                                                <  2:79840/135‑1 (MQ=255)           ttGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGCCGGTCGCGGCGCGGTGGGGCCCCGCCTGGGGgccgcg                                                                                                                             >  6:30025/1‑135 (MQ=255)           ttGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGGGCGGTGGCGCCCGGCGTGggggcggggc                                                                                                                            >  7:148327/1‑136 (MQ=255)           ttGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTTGCGCCAGGCGTTggcgcggcgc                                                                                                                            >  8:133966/1‑136 (MQ=255)           ttGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCGTGGGCgcggcg                                                                                                                             >  3:141300/1‑135 (MQ=255)           ttGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGGTGGCGCCACGCGTGGGGgccgcg                                                                                                                             >  6:444650/1‑135 (MQ=255)           ttGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTAGCGGCGGGTTTGGGCCCGGCGTGGGGGGGgcgc                                                                                                                            >  8:14168/1‑136 (MQ=255)                 aTCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCCGGCGTGGGCGCGGcgccgccc                                                                                                                       >  4:165373/1‑133 (MQ=255)                      cggcCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCcgc                                                                                                                  <  6:339662/135‑1 (MQ=255)                         ccGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTGCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCGTGGGAGCGGCGCAGCAGGCcgcgcg                                                                                                               <  1:174140/135‑1 (MQ=255)                              gCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGGCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCt                                                                                                          <  5:340153/135‑1 (MQ=255)                                    aGGCTGAACACGGCGGAACCGGCAACGGCAAGGGCAAGGGGCTTGCGCATGGGTAGGTGCCTTACAGGTCAGACGGGCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGcagcc                                                                                                    <  3:232104/135‑1 (MQ=255)                                      gCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGGGCGTTTGGGCCCGGCGTGGGGGCGGCGCAGCCCGCCGCGCGCCGCGGcagcccgc                                                                                                 >  8:137893/1‑136 (MQ=255)                                           acacGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTTGCGCCACGCTTTGGCGCGGCGCCGCCGGCCCCGCGCGGCTGCAGCCCGCACCa                                                                                             >  4:350852/1‑135 (MQ=255)                                           acacGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGGCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAg                                                                                            <  8:287140/136‑1 (MQ=255)                                            cacGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAg                                                                                            <  4:337720/135‑1 (MQ=255)                                            cacGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAg                                                                                            <  4:394877/135‑1 (MQ=255)                                            cacGGCGGAAACGGCAACGGCAAGTGAAAGGTGCTTGTGCATGGGTCGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAg                                                                                            <  6:44564/135‑1 (MQ=255)                                             acGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGGGGTGGCGCCAGGCGTTGGGGCGGGGCCGCCGGCCGCGCGAGGCTGGCGCCAGCAACact                                                                                           >  4:72539/1‑133 (MQ=255)                                                gcggAATCGGCAAGTGCAAGTGCAAGGTGGTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGTCCGCGCGAGGCTGCAGCCAGCAACAGTGAgg                                                                                       <  8:89472/136‑1 (MQ=255)                                                 cggAAACGTCAACGGCAAGTGCAAGGTGCTTGTGGATGGTTAGGTTCCTTACAGGTCAGACGGTCGCGTCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAgg                                                                                       <  2:108516/135‑1 (MQ=255)                                                 cggAAACGGCAACGGCACGTGCAAGGGGCTTGTGCATGGGTAGGGTACTTCCAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAgg                                                                                       <  3:345129/135‑1 (MQ=255)                                                    aatcggcaagggcaaGTGCAAGGTGCTTGTGCAGGGGTAGGTGCCTTACAGGGCAGCCGGGCGCGGCGGGGTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGgcgc                                                                                    <  5:121608/132‑1 (MQ=255)                                                     aacggcaacggcAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCg                                                                                  <  8:117382/136‑1 (MQ=255)                                                      acggcaacggcaAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCg                                                                                  <  4:433548/135‑1 (MQ=255)                                                      acggcaacggcaAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCg                                                                                  <  6:419187/135‑1 (MQ=255)                                                      acgccaacggcaAGTGCAAGGTGCTGGTGCAGGGGTTGGTTCCTTACAGGTCCGCCGGTCGCGGCGCGTGGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCg                                                                                  <  6:254970/135‑1 (MQ=255)                                                              ggcaaGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGGTTGCGCCCGGGGTGGGCGCGGCGCCGCCCGCCGCGCCCGGCGGCCGCCCGCCACACGGGGGCGCCGGCCCGGGt                                                                          >  3:308837/1‑135 (MQ=255)                                                              ggcaaGTGCAAGGGGCTTGTGCATGGGTAGGTTCCTTACAGGGCCGAGGGGCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGt                                                                          <  2:390293/135‑1 (MQ=255)                                                                  aGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGGCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTcagca                                                                     <  7:14168/136‑1 (MQ=255)                                                                                gtgCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGa                                                        <  8:65568/135‑1 (MQ=255)                                                                                 tgCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGACGTGGCGCCAGCGCTGGGCCCGCGCCCGCACGCCCCGCCCGGCGGGCGCCACCCACACAGAGGCGGCGGCACAGGGCACCCGGAGCAAgcgtgc                                                        >  2:350860/1‑129 (MQ=255)                                                                                        gtcggTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATggcgg                                                <  3:298598/132‑1 (MQ=255)                                                                                        gTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCCGTGGCGCCAGGCGAGGGAGCCGCCCCCCCGGCCGCGCGAGGCTTCAGCCAGCAACATTGAGGGGCCGGCCCCCGTGAGCACGATCAAGGGGGCCAAagcggc                                               >  8:360421/1‑136 (MQ=255)                                                                                           ggTTCCTTACAGGTCAGACGGTCGCGGCGCGGTGGCGCCCGGCTTGGGGGCGGCGCCGCAGGCCGCGCGAGGCTGCCGCCCGCCACACTGAGGCGCCGGCCCTGGTGAGCACGCTCAATGTGGACAAggcggcgg                                             >  3:202349/1‑135 (MQ=255)                                                                                                    cAGGTCCGACGGGCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGGGGACATGGCGGCGGCACCGcc                                      <  1:322980/133‑1 (MQ=255)                                                                                                          aGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCgcagc                                                                                                                         <  7:142778/44‑1 (MQ=255)                                                                                                          aGACGGTCGCGGCGCGTTGGCGCCAGGCGTGGGGGCGGCGCCGCAGGCCGCCCGCGGCTGCAGCCAGCCACAGTGAGGCCCCGGCCCTGGTGAGCAGGCTCAAGGTGGACAAGGCCGCGGCCCCCCCCAccgtgc                              >  2:312733/1‑130 (MQ=255)                                                                                                          aGACGGTCGCGGCGCGTTGGCGCCAGGCGTGGGAGCGGCgcagc                                                                                                                         >  8:142778/1‑44 (MQ=255)                                                                                                          aGACGGGCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGcc                             <  8:148327/136‑1 (MQ=255)                                                                                                           gccGGGCGCGGCGCGTGGGCGCCTGGCTGGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGGGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGcc                             <  5:349330/133‑1 (MQ=255)                                                                                                           gACGTTCGCGGCGCGGTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCGGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGcc                             <  5:105095/135‑1 (MQ=255)                                                                                                            aCGGTCGCGGCGCGTTGGCGCCAGGCTGGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCgg                           <  7:360421/136‑1 (MQ=255)                                                                                                             gggTCGGGGCGGGTGGGCGCGAGGGGTGGTCGCGGCGCAGCTGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCACGGTGGACATGGCGGCGGCAGCGCCAACGTTGCCgg                           <  6:372906/134‑1 (MQ=255)                                                                                                             cGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCgg                           <  2:109113/135‑1 (MQ=255)                                                                                                             cGGGCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCgg                           <  4:202349/135‑1 (MQ=255)                                                                                                                       gcgTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCGGCATCGTCCAt                 <  3:39887/135‑1 (MQ=255)                                                                                                                                   ggCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCGGCATCGTCCATGTTCAGCACGGCa    <  8:269993/136‑1 (MQ=255)                                                                                                                                    gCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCGGCATCGTCCATGTTCAGCACGGCa    <  4:436529/135‑1 (MQ=255)                                                                                                                                     cTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCGGCATCGTCCATGTTCAGCACGGCaa   <  4:308837/135‑1 (MQ=255)                                                                                                                                      ttGGGAGCGGCGCAGCAGGCCGCGCGAGGCAGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCGGCATCGTCCATGTTCAGCACGGCAAt  <  6:428839/135‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       AGACCGTCAGTTGGGTATCGGCGGCCGAAGCCTGCAGGCTGAACACGGCGGAAACGGCAACGGCAAGTGCAAGGTGCTTGTGCATGGGTAGGTTCCTTACAGGTCAGACGGTCGCGGCGCGTTGGCGCCAGGCTTGGGAGCGGCGCAGCAGGCCGCGCGAGGCTGCAGCCAGCAACAGTGAGGCGCCGGCACTGGTCAGCAGGATCAAGGTGGACATGGCGGCGGCACCGCCAACGTTGCCGGCATCGTCCATGTTCAGCACGGCAAT  >  KT2440/5853736‑5854003

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