Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,364,205	T→C	18.6%	D565D (GAT→GAC)	metG →	methionyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,364,205	0	T	C	18.6%	44.3 / 3.6	32	D565D (GAT→GAC)	metG	methionyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (18/8);  new base C (2/4);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCTGCGCG  >  minE/1364139‑1364257                                                                   |                                                     ggtggAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTTGCAGATGATcc                                                    >  1:3327713/1‑69 (MQ=255) ggtggAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGCTcc                                                    >  1:1344619/1‑69 (MQ=255)        gCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCa                                              >  1:2330480/1‑68 (MQ=255)        gCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCa                                              >  1:373435/1‑68 (MQ=255)        gCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCa                                              >  1:1363116/1‑68 (MQ=255)         ccTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCa                                              >  1:1938370/1‑67 (MQ=255)         ccTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAg                                             >  1:1022159/1‑68 (MQ=255)           tctAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGaa                                          <  1:1430837/69‑1 (MQ=255)           tctAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGaa                                          <  1:3256921/69‑1 (MQ=255)             caaaGAGGAAGTAAAAGCCACAGCCGCACCGGTAACTGGCCCACTGGCTGACGACCCGATTCAGGaa                                          <  1:3374652/66‑1 (MQ=255)             caaaGAGGAAGTAAAAGCCACAGCCGCACCGGTAACTGGCCCACTGGCTGACGACCCGATTCAGGaa                                          <  1:1680202/66‑1 (MQ=255)             caaaGAGGAAGTAAAAGCCACAGCCGCACCGGTAACTGGCCCACTGGCTGACGACCCGATTCAGGaa                                          <  1:490429/66‑1 (MQ=255)             caaaGAGGAAGTAAAAGCCACAGCCGCACCGGTAACTGGCCCACTGGCTGACGACCCGATTCAGGaa                                          <  1:2959922/66‑1 (MQ=255)                   agaagTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCAc                                  <  1:1165001/69‑1 (MQ=255)                    gaagTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCAcc                                 >  1:2405138/1‑69 (MQ=255)                          aaaGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTgac                           >  1:140623/1‑69 (MQ=255)                          aaaGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTgac                           <  1:2169708/69‑1 (MQ=255)                          aaaGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTgac                           >  1:2329490/1‑69 (MQ=255)                                           tgTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGa          <  1:3291308/68‑1 (MQ=255)                                           ggTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGa          <  1:2077338/69‑1 (MQ=255)                                           ggTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGa          <  1:2081416/69‑1 (MQ=255)                                            gTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGa          <  1:1059361/68‑1 (MQ=255)                                              aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt       >  1:3035919/1‑69 (MQ=255)                                              aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt       >  1:2714475/1‑69 (MQ=255)                                              aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt       >  1:2082150/1‑69 (MQ=255)                                              aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt       >  1:1371364/1‑69 (MQ=255)                                              aaCTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCt       >  1:1093275/1‑69 (MQ=255)                                                  ggCCCGCTGGCAGATGATCCGATTCAGGAAACCATCGCCTTTGACGACTTCGCTAAAGTTGACCTgcgc   >  1:2550720/1‑69 (MQ=255)                                                  ggCCCGCTGGCAGATGATCCGATTCAGGAAACCATCGCCTTTGACGACTTCGCTAAAGTTGACCTgcgc   >  1:1229400/1‑69 (MQ=255)                                                  ggCCCGCTGGCAGATGATCCGATTCAGGAAACCATCGCCTTTGACGACTTCGCTAAAGTTGACCTgcgc   >  1:849407/1‑69 (MQ=255)                                                   gCCCACTGGCTGACGACCCGATTCAGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTg      >  1:379189/1‑65 (MQ=255)                                                   gCCCACTGGCTGACGACCCGATTCAGGAAACCATCACCTTTGACGACTTCGCCAAAGTTGACCTGcgcg  >  1:3395374/1‑69 (MQ=255)                                                                   |                                                     GGTGGAAGCCTCTAAAGAAGAAGTAAAAGCCGCTGCCGCGCCGGTAACTGGCCCGCTGGCAGATGATCCGATTCAGGAAACCATCACCTTTGACGACTTCGCTAAAGTTGACCTGCGCG  >  minE/1364139‑1364257

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