BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	316,051	T→G	38.9%	L349L (CTT→CTG)	kefA →	fused mechanosensitive channel proteins

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	316,051	0	T	G	38.9%	60.4 / 121.9	126	L349L (CTT→CTG)	kefA	fused mechanosensitive channel proteins
Reads supporting (aligned to +/- strand): ref base T (57/20); new base G (25/24); total (82/44)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 126 total aligned reads displayed.

GGAACGCAATATTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGAT  >  minE/315988‑316117

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

           ttAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACaaa                                                    >  1:1594545/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:3164471/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                           ctgctgTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATg                    <  1:1163584/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:1215895/1‑60 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                          aTCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATcgcg     <  1:2139028/69‑1 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGAcc                 >  1:2403208/1‑54 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGAcc                 >  1:729252/1‑54 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGAcc                 >  1:1320916/1‑54 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGAcc                 >  1:3358346/1‑54 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGAcc                 >  1:235076/1‑54 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGa   <  1:3207616/68‑1 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGa   <  1:947831/68‑1 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGa   <  1:2452521/68‑1 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGa   <  1:2746980/68‑1 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGAt  >  1:1356804/1‑69 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGAt  >  1:184278/1‑69 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGAt  >  1:2707913/1‑69 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGAt  >  1:251252/1‑69 (MQ=255)

                                                             cTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGAt  >  1:908136/1‑69 (MQ=255)

                                                             cTTTACAAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGAcc                 >  1:1444275/1‑54 (MQ=255)

GGAACGCAATATTAAAGAGCAGATTGCCGTCCTGAAGGGCAGCCTGCTGTTGTCTCGTATCCTTTACCAGCAACAACAAACGCTGCCCTCGGCGGATGAACTGGAAAACATGACCAACCGCATCGCGGAT  >  minE/315988‑316117

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