BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,483,170	C→T	21.1%	L442L (CTG→CTA)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,483,170	0	C	T	21.1%	65.3 / 9.4	47	L442L (CTG→CTA)	trmE	GTPase
Reads supporting (aligned to +/- strand): ref base C (21/16); new base T (3/7); total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAATGCGGCTTCGTAAGCGCGGTTACTTACCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGCAA  >  minE/2483108‑2483229

aGAATGCGGCTTCGTAAGCGCGGTTACTTACCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGgtcg                                                       >  1:1396719/1‑69 (MQ=255)

         cTTCGTAAGCGCGGTTACTTACCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTa                                              <  1:901871/69‑1 (MQ=255)

          ttCGTAAGCGCGGTTACTTACCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTa                                              <  1:681894/68‑1 (MQ=255)

           tCGCAAGCGCGGTTACTTACCGATACAGAAACTGGAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTaaa                                            <  1:3053210/69‑1 (MQ=255)

           tCGCAAGCGCGGTTACTTACCGATACAGAAACTGGAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTaaa                                            <  1:738072/69‑1 (MQ=255)

           tCGCAAGCGCGGTTACTTACCGATACAGAAACTGGAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTaaa                                            <  1:3431337/69‑1 (MQ=255)

           tCGCAAGCGCGGTTACTTACCGATACAGAAACTGGAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTaaa                                            <  1:3026925/69‑1 (MQ=255)

           tCGCAAGCGCGGTTACTTACCGATACAGAAACTGGAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTaaa                                            <  1:260881/69‑1 (MQ=255)

           tCGCAAGCGCGGTTACTTACCGATACAGAAACTGGAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTaaa                                            <  1:2142243/69‑1 (MQ=255)

                      gTTACTTACCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGa                                 <  1:3425160/69‑1 (MQ=255)

                       ttacttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGa                                 <  1:2572142/68‑1 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAttt                              >  1:606081/1‑69 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAttt                              >  1:2434723/1‑69 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAttt                              >  1:577345/1‑69 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAttt                              >  1:3415877/1‑69 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAttt                              >  1:1852001/1‑69 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAtt                               >  1:1065912/1‑68 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAtt                               >  1:422342/1‑68 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAtt                               >  1:1013607/1‑68 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGAtt                               >  1:2859289/1‑68 (MQ=255)

                         acttacCAATACAGAAGCTGGAGAAAATCAGCCCcagc                                                             >  1:1506102/1‑38 (MQ=255)

                                  tACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGt                     <  1:69355/69‑1 (MQ=255)

                                  tACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGt                     <  1:473958/69‑1 (MQ=255)

                                  tACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGt                     <  1:1975219/69‑1 (MQ=255)

                                  tACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGt                     <  1:2730967/69‑1 (MQ=255)

                                  tACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGt                     <  1:2727024/69‑1 (MQ=255)

                                   aCAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGt                     <  1:2468594/68‑1 (MQ=255)

                                   aCAGAAACTGGAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTaaa                                            <  1:45772/45‑1 (MQ=255)

                                      gAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGt                     >  1:1575674/1‑65 (MQ=255)

                                      gAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTctg                 >  1:3381056/1‑69 (MQ=255)

                                      gAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTctg                 >  1:3133241/1‑69 (MQ=255)

                                      gAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTctg                 >  1:407361/1‑69 (MQ=255)

                                      gAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTctg                 >  1:713378/1‑69 (MQ=255)

                                      gAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTctg                 >  1:2784497/1‑69 (MQ=255)

                                            ggAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTAAATTCCCCGGTGATCTCGCTTAAGTTCTGCtgt             >  1:2769309/1‑67 (MQ=255)

                                            ggAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTAAATTCCCCGGTGATCTCGCTTAAGTTCTGCtgt             >  1:2469146/1‑67 (MQ=255)

                                            ggAGAAAATCCGCCCCAGTAGGTCGTCTGAAGTAAATTCACCGGTGATCTCGCTTAAGTTCTGCtgt             >  1:2774805/1‑67 (MQ=255)

                                               gaAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCg                                     >  1:1581577/1‑40 (MQ=255)

                                               gaAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCa         >  1:679972/1‑68 (MQ=255)

                                                   aTCCGCCCCAGCAGGTCGTCTGAAGTAAATTGCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGc    <  1:2193197/69‑1 (MQ=255)

                                                   aTCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGc    <  1:632250/69‑1 (MQ=255)

                                                   aTCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGc    <  1:2212065/69‑1 (MQ=255)

                                                   aTCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGc    <  1:540371/69‑1 (MQ=255)

                                                   aTCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGc    <  1:704766/69‑1 (MQ=255)

                                                   aTCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGc    <  1:2839027/69‑1 (MQ=255)

                                                     ccGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGtt                    >  1:409098/1‑51 (MQ=255)

                                                     ccGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGCaa  >  1:327563/1‑69 (MQ=255)

                                                     ccGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGCaa  >  1:826974/1‑69 (MQ=255)

AGAATGCGGCTTCGTAAGCGCGGTTACTTACCAATACAGAAGCTGGAGAAAATCCGCCCCAGCAGGTCGTCTGAAGTAAATTCCCCGGTGATTTCGCTTAAGTTCTGCTGTGCCAGACGCAA  >  minE/2483108‑2483229

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