Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	429,936	A→G	38.0%	S181S (TCA→TCG)	glnS →	glutamyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	429,936	0	A	G	38.0%	33.2 / 45.3	71	S181S (TCA→TCG)	glnS	glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base G (25/2);  total (47/24)
Fisher's exact test for biased strand distribution p-value = 2.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

AAAATGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGCA  >  minE/429871‑430003                                                                  |                                                                    aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2994747/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2266410/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2261406/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2299093/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2403011/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2518431/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2544054/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:258885/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2701220/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2775794/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2783343/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:294460/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2961367/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:3108354/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:1734548/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:1623218/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:33587/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:1562050/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:3397771/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:3405021/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:389397/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:773522/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:795863/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:866878/1‑69 (MQ=255) aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:1087211/1‑69 (MQ=255)     tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:1109641/1‑69 (MQ=255)     tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:3143217/1‑69 (MQ=255)     tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:2612619/1‑69 (MQ=255)     tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:702240/1‑69 (MQ=255)     tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:2438195/1‑69 (MQ=255)     tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:1981752/1‑69 (MQ=255)     tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:1565534/1‑69 (MQ=255)     tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:1786528/1‑69 (MQ=255)         tGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGt                                                          <  1:1959853/69‑1 (MQ=255)                ggtTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATgcgc                                                   >  1:642476/1‑69 (MQ=255)                            ggTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCGTTTATCGTGATGCGCGATCCGGTGCTg                                       <  1:768106/69‑1 (MQ=255)                            ggTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCGTTTATCGTGATGCGCGATCCGGTGCTg                                       <  1:3114090/69‑1 (MQ=255)                              tAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTa                                     >  1:2369693/1‑69 (MQ=255)                              tAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTa                                     >  1:383954/1‑69 (MQ=255)                              tAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTa                                     >  1:165257/1‑69 (MQ=255)                              tAAAGCATGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTa                                     >  1:52923/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:2259646/41‑1 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:872574/41‑1 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:1431230/41‑1 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:760422/41‑1 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:1528355/41‑1 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:2389804/41‑1 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:2053594/41‑1 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:2732282/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:1480220/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:386252/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:2792997/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:380953/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:2206456/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:1248436/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:946653/1‑69 (MQ=255)                               aaaGCCTGCCTGCGTGCGAAAAGCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:3142454/1‑69 (MQ=255)                                aaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:579818/40‑1 (MQ=255)                                             tgcgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTg                       <  1:1467800/68‑1 (MQ=255)                                             tgcgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTg                       <  1:1763067/68‑1 (MQ=255)                                             tgcgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTg                       <  1:195312/68‑1 (MQ=255)                                             tgcgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGAGATCCGGTGCTGTACCGTATTAAATTTg                       <  1:2134135/68‑1 (MQ=255)                                               cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:261337/69‑1 (MQ=255)                                               cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1332111/69‑1 (MQ=255)                                               cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1553901/69‑1 (MQ=255)                                               cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1883385/69‑1 (MQ=255)                                               cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1870957/69‑1 (MQ=255)                                               cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:3093585/69‑1 (MQ=255)                                               cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1825980/69‑1 (MQ=255)                                                gAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:634950/68‑1 (MQ=255)                                                         cATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACAccacca          <  1:570933/69‑1 (MQ=255)                                                                  aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:2672017/1‑68 (MQ=255)                                                                  aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:3193641/1‑68 (MQ=255)                                                                  aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:1786949/1‑68 (MQ=255)                                                                  aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:1806644/1‑68 (MQ=255)                                                                  aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:227283/1‑68 (MQ=255)                                                                  aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:2917844/1‑68 (MQ=255)                                                                  |                                                                    AAAATGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGCA  >  minE/429871‑430003

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