BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	165,064	T→C	37.9%	V804V (GTA→GTG)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	165,064	0	T	C	37.9%	28.2 / 35.9	58	V804V (GTA→GTG)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand): ref base T (18/18); new base C (18/4); total (36/22)
Fisher's exact test for biased strand distribution p-value = 2.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCT  >  minE/165002‑165125

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCg                                                         >  1:2388337/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:3213313/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:2814772/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:2713087/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:2706768/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:2585705/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:3369489/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:3395144/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:365495/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:1939566/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:1675091/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:498600/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:3040229/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:3274047/1‑69 (MQ=255)

ggTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:1218304/1‑69 (MQ=255)

ggTTGGCCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCg                                                         >  1:3275513/1‑69 (MQ=255)

  ttGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGt                                                       <  1:2277141/69‑1 (MQ=255)

  ttGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGt                                                       <  1:2534859/69‑1 (MQ=255)

  ttGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGt                                                       <  1:2616080/69‑1 (MQ=255)

  ttGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGt                                                       <  1:3053436/69‑1 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTTTGGGTCGGCAAAAACGTTACTTCGGt                                                       >  1:2877791/1‑68 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGt                                                       >  1:1529997/1‑68 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGt                                                       >  1:2019539/1‑68 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGt                                                       >  1:2975946/1‑68 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGt                                                       >  1:786389/1‑68 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCGGtt                                                      <  1:2817272/69‑1 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCGGtt                                                      <  1:3123337/69‑1 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCGGtt                                                      <  1:2183582/69‑1 (MQ=255)

   tGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTCACTTCGGtt                                                      <  1:47778/69‑1 (MQ=255)

                     tCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGt                                    <  1:3312468/69‑1 (MQ=255)

                     tCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGt                                    <  1:1492525/69‑1 (MQ=255)

                     tCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGt                                    <  1:1421190/69‑1 (MQ=255)

                     tCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGt                                    <  1:3258875/69‑1 (MQ=255)

                     tCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGt                                    <  1:836458/69‑1 (MQ=255)

                     tCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGt                                    <  1:1066036/69‑1 (MQ=255)

                                 ttGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGt                                    <  1:426185/57‑1 (MQ=255)

                                       tcggtATGGGTCGACAAAAACGTCACTTCGGTTTCAACAGTAAAATGGCGCAATTTAGCGGGTTGGCGa                  >  1:681427/1‑69 (MQ=255)

                                       tcggtATGGGTCGACAAAAACGTCACTTCGGTTTCAACAGTAAAATGGCGCAATTTAGCGGGTTGGCGa                  >  1:981939/1‑69 (MQ=255)

                                       tcggtATGGGTCGACAAAAACGTCACTTCGGTTTCAACAGTAAAATGGCGCAATTTAGCGGGTTGGCGa                  >  1:1081935/1‑69 (MQ=255)

                                             tGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGt                                    <  1:967203/45‑1 (MQ=255)

                                                   ggCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGt      <  1:1168714/69‑1 (MQ=255)

                                                   ggCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGt      <  1:219166/69‑1 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGc                    <  1:1236626/49‑1 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGc                    <  1:2223492/49‑1 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGc                    <  1:315483/49‑1 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGc                    <  1:155619/49‑1 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:2059430/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:870650/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:1099992/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:590604/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:1523883/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:411538/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:2520253/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:2863577/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:296321/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:3045636/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:3050832/1‑67 (MQ=255)

                                                         aaCGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCt  >  1:3116798/1‑67 (MQ=255)

GGTTGGTCGAGGGCGATCAGTTCGAGGAACGATTTGCGGTCGGTATGGGTCGGCAAAAACGTTACTTCGGTTTCAACAGTAAAATGGCGTAATTTGGCGGGTTGGCGACGGGGCTGCGGTGGCT  >  minE/165002‑165125

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