Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	5,022,441	G→A	14.3%	A295T (GCC→ACC)	doeA →	Ectoine hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	5,022,441	0	G	A	14.3%	62.5 / 4.8	28	A295T (GCC→ACC)	doeA	Ectoine hydrolase
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base A (2/2);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

TTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCAC  >  KT2440/5022309‑5022574                                                                                                                                     |                                                                                                                                      ttcGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGCAGGTATCGCTGCCGGCCGCGAGgccgcc                                                                                                                                     >  5:160737/1‑135 (MQ=255)              gctgctACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAaccc                                                                                                                       >  8:137470/1‑134 (MQ=255)                              tCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCg                                                                                                         <  7:92833/134‑1 (MQ=255)                                       gCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGCGGCCGCCAAGCCTGGCAACACCACTGGGGCTATAGCACGTGCctt                                                                                               >  4:247877/1‑135 (MQ=255)                                       gCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCa                                                                                                                           >  8:117315/1‑107 (MQ=255)                                       gCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCa                                                                                                                           <  7:117315/107‑1 (MQ=255)                                           tCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCttcttc                                                                                           >  4:69546/1‑135 (MQ=255)                                              cGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATGGCACGTGCCGTCTTCGGc                                                                                        >  1:131350/1‑135 (MQ=255)                                               gAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCACCAAGCCTGGCAACACCACTGGTg                                                                                                               >  1:37981/1‑111 (MQ=255)                                               gAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCACCAAGCCTGGCAACACCACTGGTg                                                                                                               <  2:37981/111‑1 (MQ=255)                                               gAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCACCAAGCCTGGCAACACCACTGGTg                                                                                                               <  6:132084/111‑1 (MQ=255)                                               gAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCACCAAGCCTGGCAACACCACTGGTg                                                                                                               >  5:132084/1‑111 (MQ=255)                                                      ttACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTg                                                                                 <  1:5813/134‑1 (MQ=255)                                                        aCCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACCCCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTga                                                                                >  2:70220/1‑133 (MQ=255)                                                        aCCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTga                                                                                >  4:64031/1‑133 (MQ=255)                                                        aCCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTga                                                                                >  6:200325/1‑133 (MQ=255)                                                              gCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCg                                                                        <  1:112204/135‑1 (MQ=255)                                                                aaGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGat                                                                      >  4:158788/1‑135 (MQ=255)                                                                      ccGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCcac                                                                <  4:62065/135‑1 (MQ=255)                                                                                           gagaAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTAcc                                          <  7:135993/136‑1 (MQ=255)                                                                                                aGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATc                                     <  7:202625/136‑1 (MQ=255)                                                                                                      cgtcgAAGGTATCGCTGCCGGCCTCGAGGCCGCCAACCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCCATCGGTa                                 >  1:62022/1‑134 (MQ=255)                                                                                                       gtcgAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATc                               >  1:144227/1‑135 (MQ=255)                                                                                                        tcgAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCa                              >  4:158140/1‑135 (MQ=255)                                                                                                             ggTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTAt                         <  3:247877/135‑1 (MQ=255)                                                                                                             ggTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTAt                         <  1:42815/135‑1 (MQ=255)                                                                                                                               gAGGCCGCCTAGCCTGGCAACAGCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTTTCCACCAGAGTGGGGCGa        <  7:84942/134‑1 (MQ=255)                                                                                                                               gAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGa        <  8:48977/134‑1 (MQ=255)                                                                                                                                aGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGa        <  5:237552/133‑1 (MQ=255)                                                                                                                                  gccgccAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAgcgc    >  4:56312/1‑135 (MQ=255)                                                                                                                                    cgccAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAgcgctc  >  5:33039/1‑133 (MQ=255)                                                                                                                                     |                                                                                                                                      TTCGAAATCGCCGGCTGCTACAAGCGCTATCACTGCCCGCTGTCCCGAACTGTTTACCTGGGCAAGCCACCGCAGCATTTCATCGACGGCGAGAAAGCAGTCGTCGAAGGTATCGCTGCCGGCCTCGAGGCCGCCAAGCCTGGCAACACCACTGGTGATATCGCACGTGCCTTCTTCGGCGTGCTTGAGAAGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCAC  >  KT2440/5022309‑5022574

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