Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	1,845,590	C→T	10.2%	A540A (GCC→GCT)	relA →	GTP pyrophosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	1,845,590	0	C	T	10.2%	114.6 / 3.8	39	A540A (GCC→GCT)	relA	GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base C (18/17);  new base T (2/2);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGCGGGCCAAGATCGTCCACTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGCCTGAGCGTATCGAGCAGATCGAGCTGGT  >  KT2440/1845460‑1845721                                                                                                                                   |                                                                                                                                    gggcgcgggcCAAGATCGTCCACTGGTTCAAATTGCAGGCCAGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGGTGATTTCGAGCGCCTGGCCga                                                                                                                                   >  6:107670/1‑133 (MQ=255)  ggcgcgggcCAAGATCGTCCACTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCAGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGagaag                                                                                                                               <  8:164340/136‑1 (MQ=255)          ccAAGATCGTCCACTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGt                                                                                                                       <  7:1303/136‑1 (MQ=255)                 cGTCCACTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGTCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGAcc                                                                                                                <  7:186111/136‑1 (MQ=255)                      aCTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGCAGACCAACGTCAAGAccgc                                                                                                              >  1:240430/1‑133 (MQ=255)                      aCTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCCACGTCAAGAccgc                                                                                                              >  1:228245/1‑133 (MQ=255)                      aCTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGAccgc                                                                                                              >  6:159507/1‑133 (MQ=255)                      aCTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGAccgc                                                                                                              >  5:92969/1‑133 (MQ=255)                       cTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGa                                                                                                           <  2:34606/135‑1 (MQ=255)                           ttCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGACAAGACCGCCGAGGAc                                                                                                       >  6:65485/1‑135 (MQ=255)                             cAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACAt                                                                                                     >  2:97531/1‑135 (MQ=255)                                aTTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCTGAGAAGACCAACGTCAAGACCGCCGAGGACa                                                                                                      >  3:199926/1‑131 (MQ=255)                                aTTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCTGAGAAGACCAACGTCAAGACCGCCGAGGACa                                                                                                      >  3:203439/1‑131 (MQ=255)                                aTTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCTGAGAAGACCAACGTCAAGACCGCCGAGGACa                                                                                                      <  4:199926/131‑1 (MQ=255)                                aTTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCTGAGAAGACCAACGTCAAGACCGCCGAGGACa                                                                                                      <  4:203439/131‑1 (MQ=255)                                     aGGGCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGGCGGCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCctc                                                                                            <  8:63539/136‑1 (MQ=255)                                       gCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCctcg                                                                                           >  6:138327/1‑135 (MQ=255)                                               ccAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGc                                                                                   <  5:159507/135‑1 (MQ=255)                                                            gcCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGcc                                                                     <  7:73521/136‑1 (MQ=255)                                                                gCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTgg                                                                  >  5:256176/1‑135 (MQ=255)                                                                gCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGGTGATTTCGCGCGCCTGGCCGAGAAGACCAACCTCAAGACCGCCGAGGTCATGTTTGCCTCGCTCGGGGCTGGCGACCTGCGCCggg                                                                  >  5:157962/1‑135 (MQ=255)                                                                  aaGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCt                                                                <  6:92969/135‑1 (MQ=255)                                                                            ttgAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGTTTTCGAGCGCCTGGCCGAGAAGACCAACGTCATGAGCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCa                                                      <  3:6214/135‑1 (MQ=255)                                                                            ttgAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCa                                                      <  1:97531/135‑1 (MQ=255)                                                                                  gcgAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGc                                               >  7:201884/1‑136 (MQ=255)                                                                                    gAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGcc                                                                     <  4:153779/112‑1 (MQ=255)                                                                                    gAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGcc                                                                     >  3:153779/1‑112 (MQ=255)                                                                                    gAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTAAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGAACTTCGCCTGGCTCATCTGTTCAACGCTGcc                                              >  6:39830/1‑135 (MQ=255)                                                                                             cgtcTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTg                                     <  2:179499/135‑1 (MQ=255)                                                                                                cTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTg                                  <  6:126216/135‑1 (MQ=255)                                                                                                    gCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGcc                             <  8:205577/136‑1 (MQ=255)                                                                                                     ccTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCa                                                      >  8:78526/1‑110 (MQ=255)                                                                                                     ccTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCa                                                      <  7:78526/110‑1 (MQ=255)                                                                                                         ccGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGCCTGAg                         >  5:54092/1‑135 (MQ=255)                                                                                                               gTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGCCTGAGCGTATc                   <  8:201884/135‑1 (MQ=255)                                                                                                               gTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGCCTGAGCGTATc                   <  3:98578/135‑1 (MQ=255)                                                                                                                  gATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGCCTGAGCGTATCGAg                >  3:68178/1‑135 (MQ=255)                                                                                                                            gcCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATCTTTGCCTCGCTCTGTTCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGCCTGAGCGTATCGAGCAGATCGAGc      >  6:208974/1‑135 (MQ=255)                                                                                                                                ggCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGCCTGAGCGTATCGAGCAGATCGAGCTGGt  >  5:60966/1‑135 (MQ=255)                                                                                                                                   |                                                                                                                                    GGGCGCGGGCCAAGATCGTCCACTGGTTCAAATTGCAGGCCCGCGACCAGAACGTTGCTGCCGGCAAGACCTTGCTTGAGCGCGAGCTCAGTCGTCTGGGCCTGCCGCAGGTTGATTTCGAGCGCCTGGCCGAGAAGACCAACGTCAAGACCGCCGAGGACATGTTTGCCTCGCTCGGTGCTGGCGACCTGCGCCTGGCTCATCTGGTCAACGCTGCCCAGCAGTTGCTGGAGCCTGAGCGTATCGAGCAGATCGAGCTGGT  >  KT2440/1845460‑1845721

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