BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,381,100	A→C	34.5%	L299L (CTT→CTG)	yeiB ←	conserved inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,381,100	0	A	C	34.5%	17.7 / 20.2	29	L299L (CTT→CTG)	yeiB	conserved inner membrane protein
Reads supporting (aligned to +/- strand): ref base A (7/12); new base C (4/6); total (11/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCGTGGTACAAATCAGCGTTTGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGCGGCCAGAAGCCATAAAACAGCGACGCATAGCCA  >  minE/1381033‑1381160

gCGTGGTACAAATCAGCGTTTGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCa                                                              >  1:923349/1‑68 (MQ=255)

gCGTGGTACAAATCAGCGTTTGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCa                                                              >  1:590054/1‑68 (MQ=255)

        cAAATCAGCGTTTGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAg                                                      >  1:1341603/1‑68 (MQ=255)

        cAAATCAGCGTTTGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAg                                                      >  1:2109197/1‑68 (MQ=255)

              aGCGTTTGCAATAGATAGTTGGTCAGTGCCATCCGTCCGACACAGGCGATTGCCAGCACCAGCTTAAAg                                               <  1:1130118/69‑1 (MQ=255)

              aGCGTTTGCAATAGATAGTTGGTCAGTGCCATCCGTCCGACACAGGCGATTGCCAGCACCAGCTTAAAg                                               <  1:1080308/69‑1 (MQ=255)

              aGCGTTTGCAATAGATAGTTGGTCAGTGCCATCCGTCCGACACAGGCGATTGCCAGCACCAGCTTAAAg                                               <  1:171347/69‑1 (MQ=255)

              aGCGTTTGCAATAGATAGTTGGTCAGTGCCATCCGTCCGACACAGGCGATTGCCAGCACCAGCTTAAAg                                               <  1:742659/69‑1 (MQ=255)

              aGCGTTTGCAATAGATAGTTGGTCAGTGCCATCCGTCCGACACAGGCGATTGCCAGCACCAGCTTAAAg                                               <  1:2213161/69‑1 (MQ=255)

              aGCGTTTGCAATAGATAGTTGGTCAGTGCCATCCGTCCGACACAGGCGATTGCCAGCACCAGCTTAAAg                                               <  1:3423196/69‑1 (MQ=255)

                  tttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:987612/69‑1 (MQ=255)

                  tttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:860887/69‑1 (MQ=255)

                  tttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:707166/69‑1 (MQ=255)

                  tttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:619192/69‑1 (MQ=255)

                  tttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:593365/69‑1 (MQ=255)

                  tttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:3244820/69‑1 (MQ=255)

                  tttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:3232771/69‑1 (MQ=255)

                  tttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:2338805/69‑1 (MQ=255)

                   ttGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:1271162/68‑1 (MQ=255)

                          agatagTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGc                                   >  1:387978/1‑69 (MQ=255)

                          agatagTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGc                                   >  1:2345604/1‑69 (MQ=255)

                          agatagTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGc                                   >  1:1928512/1‑69 (MQ=255)

                          agatagGTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGc                                   >  1:440677/1‑69 (MQ=255)

                                ttGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGc                                           <  1:2702866/55‑1 (MQ=255)

                                     cAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGCGGCCAGAAGc                         <  1:364233/68‑1 (MQ=255)

                                      aGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGCGGCCAGAAGc                         <  1:778631/67‑1 (MQ=255)

                                               cgtccgACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGCGGCCAGAAGCCATAAAACAGc              >  1:378523/1‑69 (MQ=255)

                                                           gCGATTGCCAGCACCAGCTTAAAGCGGCTTAATTGCGGCCAGAAGCCATAAAACAGCGACGCATAGCCa  >  1:2981246/1‑69 (MQ=255)

                                                           gCGATTGCCAGCACCAGCTTAAAGCGGCTTAATTGCGGCCAGAAGCCATAAAACAGCGACGCATAGCCa  >  1:2047514/1‑69 (MQ=255)

                                                           gCGATTGCCAGCACCAGCTTAAAGCGGCTTAATTGCGGCCAGAAGCCATAAAACAGCGACGCATAGCCa  >  1:1988465/1‑69 (MQ=255)

                                                           gCGATTGCCAGCACCAGCTTAAAGCGGCTTAATTGCGGCCAGAAGCCATAAAACAGCGACGCATAGCCa  >  1:1480871/1‑69 (MQ=255)

GCGTGGTACAAATCAGCGTTTGCAATAGATAGTTGGTCAGCGCCATCCGTCCGACGCAGGCGATCGCAAGCACCAGCTTAAAGCGGCTCAATTGCGGCCAGAAGCCATAAAACAGCGACGCATAGCCA  >  minE/1381033‑1381160

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