BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,284,101	A→G	38.2%	K1960K (AAA→AAG)	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,284,101	0	A	G	38.2%	83.0 / 107.5	157	K1960K (AAA→AAG)	yeeJ	adhesin
Reads supporting (aligned to +/- strand): ref base A (52/45); new base G (14/46); total (66/91)
Fisher's exact test for biased strand distribution p-value = 2.29e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

CGGTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATCGC  >  minE/1284035‑1284163

cGGTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGt                                                              <  1:613902/69‑1 (MQ=255)

cGGTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGt                                                              <  1:2537351/69‑1 (MQ=255)

cGGTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGt                                                              <  1:3333901/69‑1 (MQ=255)

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

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

       tataCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTgg                                                        >  1:2227035/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

         taCGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCa                                                      >  1:1439573/1‑68 (MQ=255)

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

           cGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTc                                                             <  1:347305/59‑1 (MQ=255)

           cGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAAc                                                    <  1:1477524/68‑1 (MQ=255)

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

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

                  tGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCAc                                            >  1:1871924/1‑69 (MQ=255)

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

                        gttCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      >  1:3119208/2‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                        aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:786642/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:380671/60‑1 (MQ=255)

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

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

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

                                              cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:1556537/1‑67 (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:280502/1‑67 (MQ=255)

                                              cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCCAGCttt                         <  1:435191/60‑1 (MQ=255)

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

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

                                                      gacaACCGTCAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGa         >  1:2198153/1‑68 (MQ=255)

                                                        caACCGTGAAAGTCACTGGCAACTCAAGCACCGccc                                       >  1:1993294/1‑36 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                           ccGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:549640/68‑1 (MQ=255)

                                                           ccGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:2623362/68‑1 (MQ=255)

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

                                                              tGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:2257091/47‑1 (MQ=255)

CGGTACATATACGGTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATCGC  >  minE/1284035‑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: `‑`