BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	316,042	T→A	44.5%	S346S (TCT→TCA)	kefA →	fused mechanosensitive channel proteins

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	316,042	0	T	A	44.5%	37.5 / 117.2	110	S346S (TCT→TCA)	kefA	fused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand): ref base T (45/16); new base A (25/24); total (70/40)
Fisher's exact test for biased strand distribution p-value = 1.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GGCGCTGCAATCGGAACGCAATATTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATG  >  minE/315976‑316099

            ggAACGCAATATTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCg                                                         >  1:834285/1‑57 (MQ=255)

            ggAACGCAATATTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTAcc                                              >  1:1310432/1‑68 (MQ=255)

            ggAACGCAATATTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTAc                                               >  1:76699/1‑67 (MQ=255)

                      atTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCGTGCTGTTGTCTCGTATCCTTTACCAGcaacaacaa                                   >  1:1369128/1‑69 (MQ=255)

                      atTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGcaacaacaa                                   >  1:3121673/1‑69 (MQ=255)

                      atTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGcaacaacaa                                   >  1:961571/1‑69 (MQ=255)

                      atTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGcaacaacaa                                   >  1:1929554/1‑69 (MQ=255)

                      atTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGcaacaacaa                                   >  1:1524687/1‑69 (MQ=255)

                      atTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCCGCTGTTGTCTCGTATCCTTTACCAGcaacaacaa                                   >  1:877283/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGcaacaacaa                                   >  1:768853/1‑68 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:2427784/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:1621204/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:421133/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:1985398/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:2157236/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:2209371/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:228125/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:2365534/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:1032008/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:2604993/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:2650860/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:2687803/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:3164471/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:2951849/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:1032268/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:1594545/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:156441/1‑69 (MQ=255)

                       ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                  >  1:1546029/1‑69 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCg                                                         >  1:2359585/1‑44 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAAc                                 >  1:3246014/1‑68 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACg                                >  1:3088300/1‑69 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACg                                >  1:349776/1‑69 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACg                                >  1:1440022/1‑69 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACg                                >  1:1268574/1‑69 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACg                                >  1:2996390/1‑69 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACg                                >  1:193941/1‑69 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACg                                >  1:1726888/1‑69 (MQ=255)

                         aaaGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACg                                >  1:1607598/1‑69 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:2646247/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:1183341/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:2542005/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:1145748/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:1085865/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:2708759/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:3150912/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:339824/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:2994846/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:365077/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:1474330/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:211595/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:418723/67‑3 (MQ=255)

                           ggagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:1888165/67‑3 (MQ=255)

                            gagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:748295/67‑3 (MQ=255)

                            gagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:2786511/67‑3 (MQ=255)

                            gagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:1083779/67‑3 (MQ=255)

                            gagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:461258/67‑3 (MQ=255)

                            gagCAGATTGCCGTTCTGAAGGGCAGTCTGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:3285046/67‑3 (MQ=255)

                            gagCAGATTGCCGTTCTGAAGGGCAGTCCGCTGTTGTCACGTATCCTGTACCAGCAACAACAAAcac                               <  1:2309280/67‑3 (MQ=255)

                                         tCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATg                <  1:1813787/69‑1 (MQ=255)

                                         tCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATg                <  1:3090097/69‑1 (MQ=255)

                                         tCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATg                <  1:3235003/69‑1 (MQ=255)

                                                  gCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGaa        >  1:1596826/1‑68 (MQ=255)

                                                  gCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGaa        >  1:2336515/1‑68 (MQ=255)

                                                   cAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAGACGCTGCCCTCGGCGGATGAACTGGaaaa      >  1:3195869/1‑69 (MQ=255)

                                                   cAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGaaaa      <  1:272972/69‑1 (MQ=255)

                                                   cAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGaaaa      <  1:305868/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:975545/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:560466/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:367116/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:1163584/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:2879662/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:1388203/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:1459865/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:1931101/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:193410/69‑1 (MQ=255)

                                                       ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg  <  1:2189493/69‑1 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGGAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:498141/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATTAACTGGaa        >  1:110501/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTTGaa        >  1:10087/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:3396181/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:3029657/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:1040817/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:931767/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:2974368/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:2893353/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:809880/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:1215895/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:1226882/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:2812717/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:2746342/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:1656394/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:247195/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:2061044/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:2335357/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:3412188/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGaa        >  1:2391205/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGa         >  1:2710240/1‑59 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCCCCTCGGCAGATGAACTGGa         >  1:2421099/1‑59 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCACCTCGGCAGATGAACTGGaa        >  1:3238085/1‑60 (MQ=255)

                                                          ctgTTGTCACGTATCCTGTACCAGCAACAACAAACACTCACCTCGGCAGATGAACTGGaa        >  1:23363/1‑60 (MQ=255)

GGCGCTGCAATCGGAACGCAATATTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATG  >  minE/315976‑316099

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