Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,922,614	T→C	19.8%	I870V (ATC→GTC)	valS ←	valyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,922,614	0	T	C	19.8%	55.6 / 6.9	40	I870V (ATC→GTC)	valS	valyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (25/7);  new base C (3/5);  total (28/12)
Fisher's exact test for biased strand distribution p-value = 3.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

TCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGA  >  minE/2922559‑2922664                                                        |                                                   tCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAAc                                        >  1:350411/1‑68 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCg                                            >  1:3058676/1‑59 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACgg                                      >  1:1818107/1‑65 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:1051516/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:924784/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:914168/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:859560/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:685872/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:666703/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:505767/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:3394954/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:2942721/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:2563691/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:2052751/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:1029694/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:1069004/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:1119377/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                                  >  1:1415255/1‑69 (MQ=255)      ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTATGCACCGTCGATGATCTTCGTAACGGaaa                                   >  1:759315/1‑68 (MQ=255)            gttgATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGAcc                            >  1:2260632/1‑69 (MQ=255)            gttgATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGAcc                            >  1:188537/1‑69 (MQ=255)            gttgATGAGGCCAGCCATAGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGAcc                            >  1:527183/1‑69 (MQ=255)             ttgATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGa                              <  1:3035970/66‑1 (MQ=255)             ttgATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGa                              <  1:1289411/66‑1 (MQ=255)                       ctgccaTCGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTatcatc                   >  1:2423807/3‑67 (MQ=255)                       ctgccaTCGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTatcatc                   >  1:828199/3‑67 (MQ=255)                       ctgccaTCGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTatcatc                   >  1:757306/3‑67 (MQ=255)                              cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag          <  1:970384/69‑1 (MQ=255)                              cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag          <  1:2036910/69‑1 (MQ=255)                              cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag          <  1:1150533/69‑1 (MQ=255)                              cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag          <  1:307729/69‑1 (MQ=255)                              cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag          <  1:2357994/69‑1 (MQ=255)                                   tCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAACCCGGACCATTGTCATCGGCAGGCAGCAc       >  1:1358558/1‑67 (MQ=255)                                   tCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCg                                >  1:1350276/1‑42 (MQ=255)                                   tCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACg      >  1:1929916/1‑68 (MQ=255)                                           tctGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGa  <  1:2838239/64‑1 (MQ=255)                                           tctGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGa  <  1:2708384/64‑1 (MQ=255)                                           tctGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGa  <  1:1991032/64‑1 (MQ=255)                                           tctGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGa  <  1:2345310/64‑1 (MQ=255)                                           tctGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGa  <  1:2200645/64‑1 (MQ=255)                                                        |                                                   TCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGA  >  minE/2922559‑2922664

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