BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,378,355	T→C	66.7%	K188E (AAA→GAA)	yjeJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,378,355	0	T	C	66.7%	84.2 / 62.6	69	K188E (AAA→GAA)	yjeJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (15/8); new base C (23/23); total (38/31)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTAT  >  W3110S.gb/4378292‑4378419

cTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTcc                                                           >  1:1899495/1‑71 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:3648793/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:3559927/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:3387058/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:3181598/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:259399/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:2570884/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:2366443/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:2320179/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:163695/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:1454114/71‑1 (MQ=255)

       aCTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTACGGTTggg                                                    <  1:1169511/71‑1 (MQ=255)

            cATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATa                                                <  1:2357674/70‑1 (MQ=255)

                  tAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCa                                         >  1:3688843/1‑71 (MQ=255)

                  tAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCa                                         >  1:1881045/1‑71 (MQ=255)

                  tAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTAGGAATAGGTGTCa                                         >  1:1225044/1‑71 (MQ=255)

                    aaGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCata                                       <  1:1512223/71‑1 (MQ=255)

                    aaGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCata                                       <  1:925599/71‑1 (MQ=255)

                    aaGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCata                                       <  1:596753/71‑1 (MQ=255)

                    aaGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCata                                       <  1:3415467/71‑1 (MQ=255)

                          tAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTc                                    >  1:1851793/1‑68 (MQ=255)

                          tAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTc                                    >  1:1212700/1‑68 (MQ=255)

                          tAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTc                                    >  1:2063350/1‑68 (MQ=255)

                          tAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCa                                  <  1:2902967/70‑1 (MQ=255)

                          tAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCa                                  <  1:1240515/70‑1 (MQ=255)

                          tAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCa                                  <  1:3267770/70‑1 (MQ=255)

                           aaaGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTc                                    >  1:3279199/1‑67 (MQ=255)

                           aaaGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTc                                    >  1:2712722/1‑67 (MQ=255)

                           aaaGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTc                                    >  1:2355073/1‑67 (MQ=255)

                           aaaGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTc                                    >  1:2103838/1‑67 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCAtt                                                             >  1:2932105/1‑41 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGa                                >  1:1837200/1‑70 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGAt                               >  1:3588025/1‑71 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGAt                               >  1:1207575/1‑71 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGa                                >  1:3498639/1‑70 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGa                                >  1:586758/1‑70 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGa                                >  1:1302560/1‑70 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGAt                               <  1:1265060/71‑1 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGAt                               <  1:2674074/71‑1 (MQ=255)

                            aaGCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGAt                               <  1:2784528/71‑1 (MQ=255)

                            aaGCACCGCCAGATAATGATCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGa                                >  1:927988/1‑70 (MQ=255)

                              gCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:2776829/48‑1 (MQ=255)

                              gCACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGa                                >  1:3145116/1‑68 (MQ=255)

                               cACCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTggg                                                    <  1:2166323/47‑1 (MQ=255)

                                aCCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGaa                                                  >  1:2708320/1‑48 (MQ=255)

                                aCCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGaa                                                  >  1:2371733/1‑48 (MQ=255)

                                aCCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGaa                                                  >  1:2969981/1‑48 (MQ=255)

                                aCCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGaa                                                  >  1:3697724/1‑48 (MQ=255)

                                aCCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGaa                                                  >  1:1941968/1‑48 (MQ=255)

                                aCCGCCAGATAATGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGaa                                                  >  1:863597/1‑48 (MQ=255)

                                           aTGGTCGAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGc                <  1:3559669/71‑1 (MQ=255)

                                               tCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGt             >  1:2629470/1‑70 (MQ=255)

                                                 gAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCaa          >  1:2720552/1‑71 (MQ=255)

                                                 gAACAAGTTATGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCaa          >  1:2617107/1‑71 (MQ=255)

                                                  aaCAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCa                                  <  1:1259412/46‑1 (MQ=255)

                                                  aaCAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCa                                  <  1:849038/46‑1 (MQ=255)

                                                   aCAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAAcg        <  1:1775931/71‑1 (MQ=255)

                                                    cAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAAcgt       >  1:1210033/1‑71 (MQ=255)

                                                    cAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAAcgt       >  1:3179634/1‑71 (MQ=255)

                                                    cAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAAcgt       >  1:485810/1‑71 (MQ=255)

                                                    cAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAAcgt       >  1:1724619/1‑71 (MQ=255)

                                                    cAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAAcgt       >  1:1554086/1‑71 (MQ=255)

                                                          aTGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCAt                         >  1:3107304/1‑47 (MQ=255)

                                                          aTGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCAt                         >  1:2205868/1‑47 (MQ=255)

                                                          aTGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCAt                         >  1:2195499/1‑47 (MQ=255)

                                                          aTGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCAt                         >  1:1350387/1‑47 (MQ=255)

                                                          aTGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGtat  <  1:3104654/70‑1 (MQ=255)

                                                          aTGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGtat  <  1:1141900/70‑1 (MQ=255)

                                                          aTGTTCCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGtat  <  1:2079112/70‑1 (MQ=255)

CTTTGCCACTTTCATCTTTAAAGCGGTAAAGCACCGCCAGATAATGGTCGAACAAGTTATGTTTCCATTCCGGTTGGGAATAGGTGTCATATTCCAGATTGCCATTGTCCTGGCAGTCAACGTCGTAT  >  W3110S.gb/4378292‑4378419

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