BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,922,600	C→T	38.2%	E874E (GAG→GAA)	valS ←	valyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,922,600	0	C	T	38.2%	23.4 / 27.3	47	E874E (GAG→GAA)	valS	valyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base C (26/3); new base T (3/15); total (29/18)
Fisher's exact test for biased strand distribution p-value = 6.66e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

TTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACG  >  minE/2922537‑2922660

ttcttTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGc                                                         <  1:2576004/69‑1 (MQ=255)

 tcttTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCa                                                        >  1:3408433/1‑69 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:971165/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:2486873/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:236087/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:2360572/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:3063977/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:3168138/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:1988646/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:1540404/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:14071/68‑1 (MQ=255)

    ttCGCCAGACGTGCCAGCTCATCTTCTTTGTTGATGAGGCCTGCCATCGGGATCAGCAGTTCTGCGcc                                                      <  1:643375/68‑1 (MQ=255)

                      tCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAAc                                    >  1:350411/1‑68 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCg                                        >  1:3058676/1‑59 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACgg                                  >  1:1818107/1‑65 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:505767/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:3394954/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:666703/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:685872/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:859560/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:914168/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:2942721/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:924784/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:2563691/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:1051516/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:1069004/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:1119377/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:2052751/1‑69 (MQ=255)

                           ttcttTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAAc                              >  1:1029694/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:1289411/66‑1 (MQ=255)

                                  ttgATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGa                          <  1:3035970/66‑1 (MQ=255)

                                            ctgccaTCGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTatcatc               >  1:757306/3‑67 (MQ=255)

                                            ctgccaTCGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTatcatc               >  1:828199/3‑67 (MQ=255)

                                            ctgccaTCGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTatcatc               >  1:2423807/3‑67 (MQ=255)

                                                   cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag      <  1:2036910/69‑1 (MQ=255)

                                                   cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag      <  1:307729/69‑1 (MQ=255)

                                                   cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag      <  1:1150533/69‑1 (MQ=255)

                                                   cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag      <  1:2357994/69‑1 (MQ=255)

                                                   cGGGATCAGCAGTTCTGCGCCGTCGACGATCTTTGTAACGGAAACCGGACCTTTATCATCAgcaggcag      <  1:970384/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)

TTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACG  >  minE/2922537‑2922660

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