BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,284,104	C→T	42.2%	V1961V (GTC→GTT)	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,284,104	0	C	T	42.2%	56.2 / 105.6	142	V1961V (GTC→GTT)	yeeJ	adhesin
Reads supporting (aligned to +/- strand): ref base C (52/30); new base T (14/46); total (66/76)
Fisher's exact test for biased strand distribution p-value = 3.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 151 total aligned reads displayed.

GGTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATCGC  >  minE/1284036‑1284163

ggTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTc                                                             <  1:1645082/69‑1 (MQ=255)

ggTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTc                                                             <  1:2153290/69‑1 (MQ=255)

ggTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTc                                                             <  1:737677/69‑1 (MQ=255)

ggTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTCAAAGTc                                                             <  1:435699/69‑1 (MQ=255)

      tataCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGc                                                       >  1:2612205/1‑69 (MQ=255)

      tataCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGc                                                       >  1:1810959/1‑69 (MQ=255)

      tataCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGc                                                       >  1:2712082/1‑69 (MQ=255)

       ataCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCa                                                      <  1:3237250/69‑1 (MQ=255)

        taCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCaa                                                     >  1:2251937/1‑69 (MQ=255)

        taCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCaa                                                     >  1:2678014/1‑69 (MQ=255)

        taCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCaa                                                     >  1:853672/1‑69 (MQ=255)

        taCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCaa                                                     >  1:687200/1‑69 (MQ=255)

        taCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCa                                                      <  1:743578/68‑1 (MQ=255)

            gTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCa                                                 <  1:1168627/69‑1 (MQ=255)

            gTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCa                                                 <  1:1571072/69‑1 (MQ=255)

                       tttCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      >  1:1241583/1‑69 (MQ=255)

                           cATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:3397578/65‑1 (MQ=255)

                           cATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:2792111/65‑1 (MQ=255)

                           cATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:2397105/65‑1 (MQ=255)

                               aCGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCt                               >  1:1815947/1‑68 (MQ=255)

                               aCGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTa                              >  1:1706493/1‑69 (MQ=255)

                                 ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:1889943/1‑69 (MQ=255)

                                 ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:281650/1‑69 (MQ=255)

                                 ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:1705178/1‑69 (MQ=255)

                                 ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:1090379/1‑69 (MQ=255)

                                 ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:841766/1‑69 (MQ=255)

                                 ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:426715/1‑69 (MQ=255)

                                 ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            <  1:1263562/69‑1 (MQ=255)

                                    gtcaCAATACAGACACAGATAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:2204535/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:2788779/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:848103/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:2903242/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:601545/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:3004057/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1690874/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:3157958/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:3195167/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:523499/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:3271160/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1193529/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1867221/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1834730/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:244441/66‑1 (MQ=255)

                                    gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1615067/66‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:95081/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:3066497/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:81608/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:522026/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:45563/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:3308477/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:980217/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:2587025/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1515444/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1575628/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1652571/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1942884/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1983894/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1993847/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1291701/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:2451561/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:2261118/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1196447/69‑1 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:771484/1‑69 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:786642/1‑69 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:1924680/1‑69 (MQ=255)

                                       aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:2296253/1‑69 (MQ=255)

                                            acagacacagacAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGa                 >  1:5151/1‑69 (MQ=255)

                                             cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGc                                         >  1:1444211/1‑44 (MQ=255)

                                             cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:3099334/60‑1 (MQ=255)

                                             cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:3383530/60‑1 (MQ=255)

                                             cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:380671/60‑1 (MQ=255)

                                             cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:3119245/1‑67 (MQ=255)

                                             cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:2410764/1‑67 (MQ=255)

                                             cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:82770/1‑67 (MQ=255)

                                             cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:502280/1‑67 (MQ=255)

                                               gacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATcc              <  1:797289/69‑1 (MQ=255)

                                                   cagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:2203479/54‑1 (MQ=255)

                                                     gacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGAc        >  1:1734703/1‑69 (MQ=255)

                                                     gacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGAc        >  1:2070052/1‑69 (MQ=255)

                                                     gacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGAc        >  1:1432185/1‑69 (MQ=255)

                                                       caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:2238324/1‑69 (MQ=255)

                                                       caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:1398825/1‑69 (MQ=255)

                                                       caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:2735281/1‑69 (MQ=255)

                                                       caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:317396/1‑69 (MQ=255)

                                                         aCCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:1026425/69‑1 (MQ=255)

                                                         aCCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:2695576/69‑1 (MQ=255)

                                                         aCCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:2485235/69‑1 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1447427/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1595098/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1386683/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:3096995/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:790081/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:3088546/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1326864/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:201990/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1301504/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:3004898/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1171693/1‑58 (MQ=255)

                                                          ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATc               >  1:3370054/1‑57 (MQ=255)

                                                           cGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATcgc  >  1:1545416/1‑69 (MQ=255)

GGTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATCGC  >  minE/1284036‑1284163

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