Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	4,042,017	A→C	24.4%	D308A (GAC→GCC)	quiA →	Quinate/shikimate dehydrogenase (quinone)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	4,042,017	0	A	C	24.4%	43.2 / 2.9	41	D308A (GAC→GCC)	quiA	Quinate/shikimate dehydrogenase (quinone)
Reads supporting (aligned to +/- strand):  ref base A (25/5);  new base C (5/6);  total (30/11)
Fisher's exact test for biased strand distribution p-value = 4.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTGGCATGGCCTACTTCGACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCCTACTACCAGCTGTCCTCGGCACCGCTGAT  >  KT2440/4041886‑4042152                                                                                                                                    |                                                                                                                                        ccGTGGCATGGCCTACTTCGACGCCACTGCCGCCATCGCCCCGCCAACCCAGCCGCACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACgg                                                                                                                                      >  5:197705/1‑135 (MQ=255) ccGTGGCATGGCCTACTTCGACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACgg                                                                                                                                      >  7:294113/1‑135 (MQ=255) ccGTGGCATGGCCTACTTCGACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCGGCGGCCCCACACCATAGACGgc                                                                                                                                     >  3:242456/1‑136 (MQ=255)      gCATGGCCTACTTCGACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGgccgcc                                                                                                                                 <  8:269752/135‑1 (MQ=255)         tGGCCTACTTCGACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGa                                                                                                                              >  8:171871/1‑135 (MQ=255)              tACTTCGACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGcc                                                                                                                         >  7:27436/1‑135 (MQ=255)                   cGACGCCACTGCCGCCATCGCCCAGCCACCCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTCCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGAc                                                                                                                   >  3:58290/1‑136 (MQ=255)                   cGACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGa                                                                                                                    >  2:83299/1‑135 (MQ=255)                    gACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACCGCCGCCGGATTGCCGTCGAc                                                                                                                   >  5:360497/1‑135 (MQ=255)                        ccACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACCCCATCGACGGCCGCCTGATTGCCGTCGACGCCg                                                                                                               >  5:274920/1‑135 (MQ=255)                        ccACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGa                                                                                                              >  3:285280/1‑136 (MQ=255)                             gccgccATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGCCCAACCCCATCGACGGCCGCCTGATTGCCGTCGACGCCgccccc                                                                                                          >  5:11485/1‑130 (MQ=255)                              ccgccATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGCCCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCg                                                                                                         >  8:204138/1‑135 (MQ=255)                                   aTCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGCCCCCGGCgcgg                                                                                                   >  4:198511/1‑133 (MQ=255)                                      gCCCAGCCATCCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCt                                                                                                >  4:131202/1‑136 (MQ=255)                                      gCCCAGCCACCCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGCCCAACCCCATCGCCCGCCGCCTGATTGCCGTCGACGCCGACACCGGCgcggtc                                                                                                 >  1:116918/1‑130 (MQ=255)                                      gCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGGTCt                                                                                                >  3:326881/1‑136 (MQ=255)                                      gCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCGGCTGCCCAACACCATCGACGGCCGCCGGATTGCCGTCGACGCCGACACCGGCGACTTcg                                                                                                >  3:326884/1‑135 (MQ=255)                                        ccAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCACCTGCCAGCGTCGCCTGCTGCCCAACACCATCGCCGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTg                                                                                               >  8:275231/1‑135 (MQ=255)                                              aaCCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGCCACCGGCGAGTTCTGCCAGGGc                                                                                        >  4:7691/1‑136 (MQ=255)                                                        aaCAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGCCCAACACCATCGCCGGCCGCCTGATTGCCGTCGACGCCGACCCCCGCGAGGTATGCCACGGCGTCgggcac                                                                               >  1:316979/1‑131 (MQ=255)                                                             ctctCCGATCACCGGCGTTTGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCGGCTGACCAACACCATCGACGGCCGCCGGAATGCCTTCGACGCGGACACCGGCGCGTTCTGCAAGGGCTTCCGCAACCACgg                                                                          >  1:41449/1‑135 (MQ=255)                                                                  cGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCCTCGCCTGCTGACCAACACCATCGACGGCCGCCTGAATGCCCTCGACGCCGCCCCCGGCGAGTTCTTCCAGCGCTACGGCCACCACGGCCacct                                                                    >  4:321806/1‑133 (MQ=255)                                                                     cccccGGGGTTCGCTTCGGGGCGGGGGCAACCTGCCCGCGCCGCCGGCCCCCCACCCCCACCCCCGGCCGCCTTATGCCCGGCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGa                                                                 <  4:130852/133‑1 (MQ=255)                                                                                       cggcggGTGCCACCGCCCCGCGCGCCCTGCTCCCCAACACCAGCGCGGGCCCCCTTTTTGCCGGCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTgg                                                <  8:240250/135‑1 (MQ=255)                                                                                       cgGCCGGTGCAACCTGCCAGCGGCGCCTGCTCACCAACACCCGCGACGGCCGCCTGATTGCCGACGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTgg                                                <  8:204062/135‑1 (MQ=255)                                                                                        ggCGGGTCCCACCTCCCCGCGGCGCCGGCTCACAAACCCCAGCGCCGGCCGCGTGATGGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTggg                                               <  8:22941/135‑1 (MQ=255)                                                                                         gCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGCCGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGCCCTCAAGGCAGGCCTGGGt                                              >  2:2210/1‑135 (MQ=255)                                                                                         gCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTggg                                               >  8:22986/1‑134 (MQ=255)                                                                                               gCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTc                                        >  7:300591/1‑135 (MQ=255)                                                                                                ccAACTGCCCCCGCCCCCGGCTGCCCACCACCAGCCCCGCCCGCCTTATGCCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTccc                                      <  3:200411/136‑1 (MQ=255)                                                                                                ccAACTCCCCGCGGCCCCTGCTCACCAACCCCAGCCCCGGCCGCCTGATGGCCGGCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTccc                                      <  3:245842/136‑1 (MQ=255)                                                                                                   aCTGCCGGCGGCGCCGGCTCACAACCCCAACCCACCGCCCCCTTATTGCCGGCGACGCCGACACCGGCGAGTGCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGa                                   <  3:281999/136‑1 (MQ=255)                                                                                                           cGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCcta                             >  5:183382/1‑134 (MQ=255)                                                                                                             tCGCCTGCTGACCAACACCATCGCCGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCACCAACGGCCAGGTCGACCTCAAGGAAGGCCTGGGTGATGTCCCGGACTCctacta                          >  6:59386/1‑135 (MQ=255)                                                                                                             tCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGCCTCctacta                          >  6:312461/1‑135 (MQ=255)                                                                                                              cgctggctcACCAACCCCAGCGCCGCCCGCCTGAGTGCCGCCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCCtactac                         <  1:224267/130‑1 (MQ=255)                                                                                                               gCCGGCTCACAAACCCCAGCCACGGCCGCCTGAGTCCGGGCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCCTACTACCa                       <  4:339167/136‑1 (MQ=255)                                                                                                                ccTGCTCACAACCACCATCGAGGGCCGCCTGATTGCCGGCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCCTACTACCa                       <  8:300591/135‑1 (MQ=255)                                                                                                                  tgctgACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGCCACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCCTACTACCAGc                     >  2:298293/1‑135 (MQ=255)                                                                                                                  tgctgACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCCTACTACCAGCt                    >  4:263115/1‑136 (MQ=255)                                                                                                                                    aCGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCCTACTACCAGCTGTCCTCGGCACCGCTGAt  >  3:69795/1‑136 (MQ=255)                                                                                                                                    |                                                                                                                                        CCGTGGCATGGCCTACTTCGACGCCACTGCCGCCATCGCCCAGCCAACCCAGCCGAACAGCTCGCCGATCACCGGCGTTAGCGTAGCGGCCGGTGCCAACTGCCAGCGTCGCCTGCTGACCAACACCATCGACGGCCGCCTGATTGCCGTCGACGCCGACACCGGCGAGTTCTGCCAGGGCTTCGGCAACAACGGCCAGGTCGACCTCAAGGCAGGCCTGGGTGATGTCCCGGACTCCTACTACCAGCTGTCCTCGGCACCGCTGAT  >  KT2440/4041886‑4042152

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