Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,573,449	A→T	53.3%	T94T (ACT→ACA)	gltX ←	glutamyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,573,449	0	A	T	53.3%	6.4 / 68.9	60	T94T (ACT→ACA)	gltX	glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (11/17);  new base T (12/20);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACG  >  minE/1573387‑1573516                                                               |                                                                    ttCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGccc                                                               >  1:2176627/1‑69 (MQ=255) ttCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGccc                                                               >  1:3338483/1‑69 (MQ=255) ttCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGccc                                                               <  1:3189730/69‑1 (MQ=255) ttCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGccc                                                               <  1:3017674/69‑1 (MQ=255) ttCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGccc                                                               <  1:2368495/69‑1 (MQ=255) ttCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGcc                                                                >  1:794498/1‑68 (MQ=255)          tGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGc                                                      >  1:2423972/1‑69 (MQ=255)          tGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGc                                                      >  1:1989247/1‑69 (MQ=255)               ttCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTg                                                 >  1:3388404/1‑69 (MQ=255)                 cgcgCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGccc                                                               <  1:2916442/53‑1 (MQ=255)                 cgcgCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTga                                                >  1:2924493/1‑68 (MQ=255)                 cgcgCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGACCTCTTCCAGCATCTga                                                >  1:1947463/1‑68 (MQ=255)                           ttCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACcgc                                     <  1:87098/69‑1 (MQ=255)                           ttCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACcgc                                     <  1:1702803/69‑1 (MQ=255)                           ttCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACcgc                                     <  1:2446819/69‑1 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAg                              >  1:2917352/1‑68 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:2395971/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:251745/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:221319/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:3313515/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:600124/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:745027/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:1693374/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:749961/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:1623322/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:1581296/1‑69 (MQ=255)                                   gcTCTTTAGAGCAATAGCATTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGc                             >  1:1493552/1‑69 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:3199270/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:3266097/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:99293/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:426777/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:3169925/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:675765/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:797642/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:2855880/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:2731/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:2581605/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:1522818/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:1669136/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:1853576/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:2068180/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:2138542/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:216566/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:2242387/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:2263173/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:242082/68‑1 (MQ=255)                                                     aTTTATAAGCTGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTc            <  1:2670440/68‑1 (MQ=255)                                                       ttATAAGCAGTGCCCTCTTCCAGCATCTGATCTATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTgg         <  1:117812/69‑1 (MQ=255)                                                       ttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTgg         <  1:2167098/69‑1 (MQ=255)                                                       ttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTgg         <  1:207584/69‑1 (MQ=255)                                                       ttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTgg         <  1:1966731/69‑1 (MQ=255)                                                       ttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTgg         <  1:1914519/69‑1 (MQ=255)                                                       ttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTgg         <  1:1416007/69‑1 (MQ=255)                                                       ttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTgg         <  1:853216/69‑1 (MQ=255)                                                       ttATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTgg         <  1:1192860/69‑1 (MQ=255)                                                        tataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGgt        >  1:2263609/1‑69 (MQ=255)                                                        tataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGgt        >  1:2150918/1‑69 (MQ=255)                                                        tataAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGgt        >  1:2920457/1‑69 (MQ=255)                                                              cAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACg  <  1:2499345/69‑1 (MQ=255)                                                              cAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACg  <  1:2568650/69‑1 (MQ=255)                                                               |                                                                    TTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTGGTCTGGTAGTACG  >  minE/1573387‑1573516

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑