BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,534,716	C→T	30.2%	Y103Y (TAC→TAT)	rfaG →	glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,534,716	0	C	T	30.2%	51.3 / 24.0	56	Y103Y (TAC→TAT)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base C (23/16); new base T (11/6); total (34/22)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TCGCGTTGTTGGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTAT  >  minE/2534653‑2534779

tCGCGTTGTTGGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCg                                                             >  1:1982670/1‑68 (MQ=255)

tCGCGTTGTTGGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCg                                                             >  1:3215219/1‑68 (MQ=255)

tCGCGTAGTTGGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCg                                                             >  1:2905692/1‑68 (MQ=255)

     ttgttgGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAg                                                       >  1:912334/1‑69 (MQ=255)

     ttgttgGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAg                                                       >  1:2845871/1‑69 (MQ=255)

     ttgttgGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAg                                                       >  1:2039402/1‑69 (MQ=255)

     ttgttgGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAg                                                       >  1:3032812/1‑69 (MQ=255)

     ttgttgGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAg                                                       >  1:3079026/1‑69 (MQ=255)

     ttgttgGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAg                                                       >  1:867725/1‑69 (MQ=255)

               taataaGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAg                                              >  1:3297807/1‑68 (MQ=255)

               taataaGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAg                                              >  1:2395635/1‑68 (MQ=255)

               taataaGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAg                                              >  1:2506974/1‑68 (MQ=255)

               taataaGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAg                                              >  1:2617548/1‑68 (MQ=255)

                   aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                          <  1:1605723/65‑1 (MQ=255)

                   aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                          <  1:1910814/65‑1 (MQ=255)

                   aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                          <  1:1460004/65‑1 (MQ=255)

                   aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                          <  1:2979534/65‑1 (MQ=255)

                   aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                          <  1:2885633/65‑1 (MQ=255)

                   aaaaTGCCGGGGCTGGACGTTTATTATGCCGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGagaa                                          <  1:2703397/65‑1 (MQ=255)

                    aGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAg                                        >  1:3335647/1‑69 (MQ=255)

                    aGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAg                                        >  1:479151/1‑69 (MQ=255)

                    aGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAg                                        >  1:881626/1‑69 (MQ=255)

                    aGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAg                                        >  1:917974/1‑69 (MQ=255)

                     gATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAgg                                       >  1:3159161/1‑69 (MQ=255)

                        gcctggcctgGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGtt                                     >  1:3049928/1‑68 (MQ=255)

                            ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:1081220/69‑1 (MQ=255)

                            ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:2435760/69‑1 (MQ=255)

                            ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:23520/69‑1 (MQ=255)

                            ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:2899366/69‑1 (MQ=255)

                            ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:424304/69‑1 (MQ=255)

                            ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:3380674/69‑1 (MQ=255)

                            ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:1971976/69‑1 (MQ=255)

                            ggcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGGTTTTta                                <  1:91605/69‑1 (MQ=255)

                             gcctggATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:1407054/68‑1 (MQ=255)

                                  gATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCg                           >  1:255853/1‑68 (MQ=255)

                                  gATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCg                           >  1:1186593/1‑68 (MQ=255)

                                    tGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGttt                        <  1:2011566/69‑1 (MQ=255)

                                      tttATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTaa                      <  1:1366827/69‑1 (MQ=255)

                                      tttATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTaa                      <  1:3414289/69‑1 (MQ=255)

                                      tttATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTaa                      <  1:493144/69‑1 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCgcc                         >  1:2505501/1‑54 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:2559078/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:2691857/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:1157282/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:84776/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:816921/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:2729634/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:370144/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:1741560/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:3286031/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTATGCCGAGAAAGTAGCGCAGGAGAAAGGCTTTTTCTATCGCCTGACGTCACGtta            >  1:1792444/1‑69 (MQ=255)

                                                cGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGat             <  1:2003054/68‑1 (MQ=255)

                                                    gATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGc        <  1:1419719/69‑1 (MQ=255)

                                                          tgttACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTta                                <  1:1956919/39‑1 (MQ=255)

                                                          tgttACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGc        >  1:696207/1‑63 (MQ=255)

                                                          tgttACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTAt  >  1:847416/1‑69 (MQ=255)

TCGCGTTGTTGGGTTTAATAAGATGCCTGGCCTGGATGTTTATTTTGCCGCTGATGTTTGTTACGCCGAGAAAGTTGCGCAAGAAAAAGGTTTTTTATATCGTTTAACATCACGATATCGCCATTAT  >  minE/2534653‑2534779

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