BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,036,471	T→C	35.1%	I201I (ATT→ATC)	ygiC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,036,471	0	T	C	35.1%	71.4 / 64.5	108	I201I (ATT→ATC)	ygiC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (9/61); new base C (22/16); total (31/77)
Fisher's exact test for biased strand distribution p-value = 1.51e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 108 total aligned reads displayed.

GCGACACGGTGGAAGATCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAGGTGAAAA  >  minE/2036407‑2036515

gcgACACGGTGGAAGATCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:1924334/68‑1 (MQ=255)

gcgACACGGTGGAAGATCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:3168771/68‑1 (MQ=255)

gcgACACGGTGGAAGATCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:1184274/68‑1 (MQ=255)

gcgACACGGTGGAAGATCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:2501844/68‑1 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:1526844/52‑1 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:1835970/52‑1 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:45924/52‑1 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:2040740/52‑1 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:156880/52‑1 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:2519567/52‑1 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2424617/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:1022197/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2329393/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2144031/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:3185465/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2131061/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:415585/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2718984/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2005623/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2792417/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:282495/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2893847/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:925888/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:852893/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:1556492/1‑68 (MQ=255)

                tCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:1062685/1‑68 (MQ=255)

                tCGCGGAACCATTAAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATCGCTACCGAGTTCCTCTACa                           >  1:2019962/1‑68 (MQ=255)

                 cgcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCt                                           <  1:1365547/51‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTTAAATTGCTACTGAGTTCCTCTACATCg                        <  1:1529786/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:3038624/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:3335341/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:1093432/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:2268678/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:351848/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:2832573/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:2076411/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:1606436/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:3050430/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:2044815/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:1633886/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:2972645/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:3016324/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:1659035/69‑1 (MQ=255)

                  gcgGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:672878/69‑1 (MQ=255)

                  gcagAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCg                        <  1:3420950/66‑1 (MQ=255)

                    ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                      <  1:3091003/69‑1 (MQ=255)

                    ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                      <  1:1065673/69‑1 (MQ=255)

                    ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                      <  1:1734481/69‑1 (MQ=255)

                    ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                      <  1:3310052/69‑1 (MQ=255)

                    ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                      <  1:2290628/69‑1 (MQ=255)

                    ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                      <  1:79056/69‑1 (MQ=255)

                          aTTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATc                <  1:1415372/69‑1 (MQ=255)

                          aTTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATc                <  1:2096735/69‑1 (MQ=255)

                            tCAGTTTTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:959009/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:3325002/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:3413827/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:316250/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:347117/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:826034/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:433943/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:3039850/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:1054349/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:816666/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:770255/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:280680/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:2791057/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:1901502/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:1782722/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:1649324/69‑1 (MQ=255)

                            tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTCCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:2887573/69‑1 (MQ=255)

                             cAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGc                                            >  1:639281/1‑38 (MQ=255)

                             cAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGc                                            >  1:2656507/1‑38 (MQ=255)

                             cAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGc                                            >  1:384515/1‑38 (MQ=255)

                             cAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGc                                            >  1:1647866/1‑38 (MQ=255)

                             cAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              >  1:682187/1‑68 (MQ=255)

                                  tttGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                      >  1:457685/1‑55 (MQ=255)

                                    tGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAgg        >  1:2881199/1‑67 (MQ=255)

                                    tGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAgg        >  1:1337629/1‑67 (MQ=255)

                                    tGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAGGt       >  1:2733605/1‑68 (MQ=255)

                                        ggACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAGGTGaaaa  <  1:1577372/69‑1 (MQ=255)

                                        ggACTGCGCAACGGAAGCTGAAATTGCTACTGACTTCCTCTACATCGATGATATCGGGTTAGGTGaaaa  <  1:2072818/69‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTat         <  1:375185/52‑2 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:799926/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:695014/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:1470767/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:1197088/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:506638/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:475340/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:2047522/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:2186501/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:2253314/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:2274175/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:237013/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:2541532/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:2588552/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:3294844/52‑1 (MQ=255)

                                                  aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg         <  1:2699877/52‑1 (MQ=255)

                                                   cGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:1980001/46‑1 (MQ=255)

                                                   cGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:791378/46‑1 (MQ=255)

                                                         cTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg              <  1:3263245/40‑1 (MQ=255)

GCGACACGGTGGAAGATCGCGGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAGGTGAAAA  >  minE/2036407‑2036515

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