Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	4,677,569	A→C	17.7%	K187N (AAA→AAC)	CJ019_04418 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	4,677,569	0	A	C	17.7%	66.1 / 4.0	42	K187N (AAA→AAC)	CJ019_04418	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/20);  new base C (2/6);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CCATCACCTCGGCTGACGTGGCCTTCTCGTACCGCACGCTGCTCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGCGGCCCCTACCG  >  KT2440/4677436‑4677694                                                                                                                                      |                                                                                                                              ccATCACCTCGGCTGACGTGGCCTTCTCGTACCGCACGCTGCTCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAAc                                                                                                                              >  1:84994/1‑135 (MQ=255)            gCTGACGTGGCCTTGTCGTACCGCACGCTGCTCAAGGATGGCCACCCGATATACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCaa                                                                                                                    <  4:8455/134‑1 (MQ=255)            gCTGACGTGGCCTGCTCGTACCGCACGCTGCTCACGGATGGCCACCCGATCTACCGCACCAACCGGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTGCGTGTTCAAACGCGCCGGCa                                                                                                                     <  5:40755/133‑1 (MQ=255)                aCGTGGCCTTCTCGTACCGCACGCTGCTCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCgct                                                                                                              <  8:64715/136‑1 (MQ=255)                 cggggCCGTCTCGTCCCGCACGCTGCTCAAGGAGGCCCACCCGAGCTACCGCACCAACCTGCAGGAAGTCCAGCGGGGAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCgct                                                                                                              <  3:213421/132‑1 (MQ=255)                 cGTGGCCTTCTCGTACCGCACGCTGCTCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGGGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTCCAACCGCGCCGGCAACCCgct                                                                                                              <  1:182304/135‑1 (MQ=255)                                  gCACGCTGCTCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTggg                                                                                            <  8:189740/136‑1 (MQ=255)                                   cccgcgGCTCAAGGATGGCCACCCGATCTACCGCACCAACCTGCCGGAAGGGCAGCGGGTAGACATCCTCGGCCCGCTACGCAGCGGCTTCGTGTTAAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTggg                                                                                            <  6:31615/133‑1 (MQ=255)                                   cACGCTGCGCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTggg                                                                                            <  2:147269/135‑1 (MQ=255)                                      gctgctCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAACCGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCg                                                                                          >  4:241804/1‑134 (MQ=255)                                      gctgctCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCg                                                                                          >  2:57776/1‑134 (MQ=255)                                             aaGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCg                                                                                  >  2:177801/1‑135 (MQ=255)                                                  tGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCTGCGGGTAGACATCCTCGGCCCGCTACGCATCGGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAATTGCCGGTACTg                                                                            <  7:159566/136‑1 (MQ=255)                                                   ggCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTg                                                                            <  3:272713/135‑1 (MQ=255)                                                     ccACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAACCGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGc                                                                           >  2:111912/1‑134 (MQ=255)                                                      cACCCTATCTCCCGCACCAACCTGCAGGACGGGCAGCGGGTAGACATCCGCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGccc                                                                         <  5:20245/135‑1 (MQ=255)                                                      cACCCGATCTACCGCACCCACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGccc                                                                         <  3:191797/135‑1 (MQ=255)                                                       aCCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGccc                                                                         >  4:54267/1‑134 (MQ=255)                                                       aCCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGccc                                                                         >  4:141880/1‑134 (MQ=255)                                                       aCCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGccc                                                                         >  2:217640/1‑134 (MQ=255)                                                       aCCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCa                                                                        >  8:129023/1‑135 (MQ=255)                                                       aCCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCa                                                                        >  7:62150/1‑135 (MQ=255)                                                       aCCCGATCTACCGCACCAACCTGCAGGAAGTGAAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACTCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCa                                                                        >  7:244476/1‑135 (MQ=255)                                                                        aaCCTGCCGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAg                                                       <  5:55479/135‑1 (MQ=255)                                                                         aCCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCgc                                                                                                               >  8:145098/1‑78 (MQ=255)                                                                         aCCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCgc                                                                                                               <  7:145098/78‑1 (MQ=255)                                                                            tGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAACCACTACTGGCAGAAg                                                       >  2:227824/1‑131 (MQ=255)                                                                            tGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAg                                                       <  1:227824/131‑1 (MQ=255)                                                                                                 aCATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCg                              <  3:54267/135‑1 (MQ=255)                                                                                                 aCATCCGCGGCCCGCTACGCATCCGCTGCGGGCTCACCCGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCg                              <  3:60971/135‑1 (MQ=255)                                                                                                  cATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCg                              <  1:57776/134‑1 (MQ=255)                                                                                                         ggCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTGCGAGcccacc                      >  2:230179/1‑135 (MQ=255)                                                                                                          gCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGACACCACCTTCGAGCCCCCCt                     >  4:190594/1‑135 (MQ=255)                                                                                                            ccGCTACCCATCCGCTTCGTGTTCAACCGCTCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCATCGAGCCCCCCTTgg                  <  8:244476/136‑1 (MQ=255)                                                                                                                  cGCATCCGCTTCGTGTTCAACCGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGc             <  1:111912/135‑1 (MQ=255)                                                                                                                  cGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACCCTACTGGCAGAAGCGCGACTTCAAGGCCACCCCCTTCGAGCCCCCCTTGGGCAGc             >  4:242521/1‑135 (MQ=255)                                                                                                                      gcggCTGCGTGTCCAACCGCGCCGGCAACCCGCGGCTGACCCTGCGCCTGGGCGAACTGCCGGAACTGCCCACACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGCGGcc         <  6:151949/132‑1 (MQ=255)                                                                                                                       ccGCTTCTGGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGCGGccc        <  5:258880/135‑1 (MQ=255)                                                                                                                       ccGCTTCGTGTCCAAACGCGCCGGCAACCCGCGGCTGATCCTGCGCCGGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACATCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGCGGccc        <  4:104893/135‑1 (MQ=255)                                                                                                                             cGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGCGGCCCCTACCg  <  7:129023/135‑1 (MQ=255)                                                                                                                             cGTGTGCCACCGCGCAGGCAACCGGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCCGAAGCGCGACTTCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGCGGCCCCTACCg  <  2:263409/135‑1 (MQ=255)                                                                                                                              ggtttCAAACGCGCCGGCAAGCCGCTGCTGATCGTGCGCCTGGGCGAAATCCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGCGGCCCCTACCg  <  1:63217/131‑1 (MQ=255)                                                                                                                                      |                                                                                                                              CCATCACCTCGGCTGACGTGGCCTTCTCGTACCGCACGCTGCTCAAGGATGGCCACCCGATCTACCGCACCAACCTGCAGGAAGTGCAGCGGGTAGACATCCTCGGCCCGCTACGCATCCGCTTCGTGTTCAAACGCGCCGGCAACCCGCTGCTGATCCTGCGCCTGGGCGAAATGCCGGTACTGCCCAAACACTACTGGCAGAAGCGCGACTTCAAGGCCACCACCTTCGAGCCCCCCTTGGGCAGCGGCCCCTACCG  >  KT2440/4677436‑4677694

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