BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,255,185	C→T	26.0%	G146G (GGC→GGT)	yecD →	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,255,185	0	C	T	26.0%	186.5 / 59.5	153	G146G (GGC→GGT)	yecD	predicted hydrolase
Reads supporting (aligned to +/- strand): ref base C (48/65); new base T (29/11); total (77/76)
Fisher's exact test for biased strand distribution p-value = 1.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

GTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTA  >  minE/1255120‑1255249

gtgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                               <  1:2966223/69‑1 (MQ=255)

gtgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                               <  1:17626/69‑1 (MQ=255)

gtgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                               <  1:3345164/69‑1 (MQ=255)

gtgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                               <  1:2999262/69‑1 (MQ=255)

gtgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                               <  1:1090122/69‑1 (MQ=255)

gtgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                               <  1:1663395/69‑1 (MQ=255)

gtgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                               <  1:570120/69‑1 (MQ=255)

gtgtTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                               <  1:2337439/69‑1 (MQ=255)

   ttATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAAt                                                            >  1:3130073/1‑69 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCGCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:3213854/69‑1 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:2645561/69‑1 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:3433240/69‑1 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:2136528/69‑1 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:3186430/69‑1 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:2267351/69‑1 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:57368/69‑1 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:309192/69‑1 (MQ=255)

        tgGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:3091631/69‑1 (MQ=255)

          ggATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTg                                                      <  1:606074/68‑1 (MQ=255)

            aTCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:1063687/65‑1 (MQ=255)

            aTCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:1970676/65‑1 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTgg                                                        >  1:1076989/1‑58 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTgg                                                        >  1:1757216/1‑58 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:2357782/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:2926229/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:3325402/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:1008413/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:2092598/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:1833626/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:842943/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:1032619/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:738554/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:1119282/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:446748/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:2697425/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:1202863/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:2685233/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:1431414/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:1408981/1‑65 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCGaa                                             >  1:3383814/1‑69 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTAGGAACTCGGTTTTAATCTGGTg                                                      >  1:535073/1‑60 (MQ=255)

                  aCCAATATCGGTGTTGAATCCACCGCACGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGc                                                 >  1:687843/1‑65 (MQ=255)

                   ccAATATCGGTTTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCg                                               <  1:2076598/66‑1 (MQ=255)

                   ccAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCg                                               <  1:3151535/66‑1 (MQ=255)

                   ccAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCg                                               <  1:2741559/66‑1 (MQ=255)

                   ccAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCg                                               <  1:2548746/66‑1 (MQ=255)

                   ccAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCg                                               <  1:2057629/66‑1 (MQ=255)

                   ccAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCg                                               <  1:2763446/66‑1 (MQ=255)

                   ccAATATCAGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCg                                               <  1:2430255/66‑1 (MQ=255)

                    cAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTaa                                                             >  1:1573721/1‑51 (MQ=255)

                    cAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAg                                            >  1:1646446/1‑68 (MQ=255)

                    cAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGa                                           >  1:2877486/1‑69 (MQ=255)

                    cAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGa                                           >  1:2670689/1‑69 (MQ=255)

                    cAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGa                                           >  1:2364703/1‑69 (MQ=255)

                    cAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGa                                           >  1:348135/1‑69 (MQ=255)

                    cAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGa                                           >  1:652634/1‑69 (MQ=255)

                    cAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTCGGTTTTAATCTGGTGATTGCCg                                               <  1:1088688/65‑1 (MQ=255)

                         tCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGAtt                                                   <  1:728785/56‑1 (MQ=255)

                         tCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGAtt                                                   <  1:971988/56‑1 (MQ=255)

                         tCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGAtt                                                   <  1:2605505/56‑1 (MQ=255)

                         tCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGAtt                                                   <  1:1715090/56‑1 (MQ=255)

                         tCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGAtt                                                   <  1:24967/56‑1 (MQ=255)

                         tCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGAtt                                                   <  1:1116971/56‑1 (MQ=255)

                           ggTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCTAAGACGACTGt                                    <  1:1181169/69‑1 (MQ=255)

                           ggTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGt                                    <  1:817787/69‑1 (MQ=255)

                           ggTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGt                                    <  1:191634/69‑1 (MQ=255)

                           ggTGTTGAATCCACCGCCCGCAATGCCTGGGAACATGGCTTTAATCTGGTGATTGCCGAAGACGCCTGt                                    <  1:540128/69‑1 (MQ=255)

                            gtgtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGt                                    <  1:1147353/68‑1 (MQ=255)

                             tgtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGt                                    <  1:1508192/67‑1 (MQ=255)

                                   aTCCACCGCCCGCAATGCCTGTGAACTTGGCTTTAATCTGGt                                                       <  1:1243755/42‑1 (MQ=255)

                                    tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:2146176/41‑1 (MQ=255)

                                    tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:1317096/41‑1 (MQ=255)

                                    tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:71855/41‑1 (MQ=255)

                                    tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:2644486/41‑1 (MQ=255)

                                    tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       >  1:758824/1‑41 (MQ=255)

                                    tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGt                                                       <  1:2523918/41‑1 (MQ=255)

                                                aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCGGAAGACGCCTGt                                    <  1:2049091/48‑1 (MQ=255)

                                                aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGagcag               <  1:580270/69‑1 (MQ=255)

                                                           aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGGCGAGCAGCACAATAACAg    >  1:3195856/1‑69 (MQ=255)

                                                           aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAg    >  1:379247/1‑69 (MQ=255)

                                                           aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAg    >  1:628336/1‑69 (MQ=255)

                                                           aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAg    >  1:2196763/1‑69 (MQ=255)

                                                           aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAg    >  1:3193862/1‑69 (MQ=255)

                                                           aCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCAAAATAACAg    >  1:1025969/1‑69 (MQ=255)

                                                            cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGagcag               >  1:102331/1‑57 (MQ=255)

                                                            cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGt   >  1:2350327/1‑69 (MQ=255)

                                                            cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGt   >  1:2833376/1‑69 (MQ=255)

                                                            cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGt   <  1:2975887/69‑1 (MQ=255)

                                                            cTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGt   >  1:3380180/1‑69 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:377856/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:137050/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:1782953/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:400668/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:1880097/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:2194495/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:2265634/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:2299161/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:2322555/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:2810424/1‑68 (MQ=255)

                                                              tGGCTTTAATCTGGTGATTGCCGAACACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa  >  1:2662134/1‑68 (MQ=255)

GTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTA  >  minE/1255120‑1255249

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`