BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	997,578	G→A	37.9%	S63S (TCG→TCA)	ynfM →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	997,578	0	G	A	37.9%	73.4 / 77.3	122	S63S (TCG→TCA)	ynfM	predicted transporter
Reads supporting (aligned to +/- strand): ref base G (44/31); new base A (36/10); total (80/42)
Fisher's exact test for biased strand distribution p-value = 3.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 122 total aligned reads displayed.

CTCTGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTG  >  minE/997515‑997641

ctctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGa                                                            >  1:3338255/1‑69 (MQ=255)

ctctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGa                                                            >  1:1621888/1‑69 (MQ=255)

ctctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGa                                                            >  1:731252/1‑69 (MQ=255)

ctctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGa                                                            >  1:1967101/1‑69 (MQ=255)

ctctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGa                                                            >  1:1764741/1‑69 (MQ=255)

ctctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGa                                                            >  1:1343334/1‑69 (MQ=255)

ctctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGa                                                            >  1:3246166/1‑69 (MQ=255)

  ctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAg                                                           <  1:2015417/68‑1 (MQ=255)

  ctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAg                                                           <  1:2458138/68‑1 (MQ=255)

  ctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAg                                                           <  1:1102636/68‑1 (MQ=255)

  ctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAg                                                           <  1:2712988/68‑1 (MQ=255)

  ctGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAg                                                           <  1:1395611/68‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCCCAAGAGTTTGGCt                                                    <  1:3348274/69‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    <  1:2745609/69‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    <  1:3312241/69‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    <  1:1002121/69‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    <  1:2410798/69‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    <  1:235831/69‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    <  1:1493262/69‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    <  1:1075723/69‑1 (MQ=255)

        gACTGGCAACTTTCGCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    <  1:503350/69‑1 (MQ=255)

          cTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTa                                                  <  1:1531382/69‑1 (MQ=255)

          cTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTa                                                  >  1:1968446/1‑69 (MQ=255)

          cTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTa                                                  <  1:678123/69‑1 (MQ=255)

          cTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTa                                                  <  1:2484182/69‑1 (MQ=255)

          cTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTa                                                  >  1:526383/1‑69 (MQ=255)

          cTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTa                                                  >  1:3047093/1‑69 (MQ=255)

          cTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTa                                                  <  1:3376239/69‑1 (MQ=255)

                     tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                         >  1:2905883/1‑67 (MQ=255)

                     tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                         >  1:501031/1‑67 (MQ=255)

                     tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                         >  1:2623455/1‑67 (MQ=255)

                     tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                         >  1:1033373/1‑67 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:2996629/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:3143182/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:3204058/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:3241320/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:3319515/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:911876/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:1006319/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:2280636/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:1337833/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:1760076/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:1785849/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:1756165/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                      >  1:1968442/1‑69 (MQ=255)

                      gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGACTTTGGCTTAACGCCCGCGaa                                       >  1:2656387/1‑68 (MQ=255)

                           tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                 <  1:1672131/69‑1 (MQ=255)

                           tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                 <  1:1909582/69‑1 (MQ=255)

                           tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                 <  1:733758/69‑1 (MQ=255)

                           tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                 <  1:3194678/69‑1 (MQ=255)

                           tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                 <  1:825562/69‑1 (MQ=255)

                           tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                 <  1:1678886/69‑1 (MQ=255)

                            ctctATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                 <  1:103315/68‑1 (MQ=255)

                                 ttGTTTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:2227199/1‑69 (MQ=255)

                                 ttGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTc                           >  1:2745964/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGttt                                                        >  1:3238189/1‑40 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGTCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:2104548/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCt                                                    >  1:3051768/1‑44 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTAttt                            >  1:556764/1‑68 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTAttt                            >  1:777747/1‑68 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTAttt                            >  1:254129/1‑68 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:2692217/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:1585429/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:1675805/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:1207872/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:2993245/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:1840184/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:2562564/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:794286/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:1958118/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:2062395/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:2276344/1‑69 (MQ=255)

                                 ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                           >  1:914507/1‑69 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCtt                                                   >  1:2617879/1‑41 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACt                        >  1:2256555/1‑68 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                       >  1:665632/1‑69 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                       >  1:2159005/1‑69 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                       >  1:1548838/1‑69 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                       >  1:882835/1‑69 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                       >  1:3138947/1‑69 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                       >  1:3125851/1‑69 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                       >  1:3159921/1‑69 (MQ=255)

                                     gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                       >  1:2527035/1‑69 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                    >  1:1522378/1‑69 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                    >  1:2146625/1‑69 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                    >  1:1817517/1‑69 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                    >  1:1096765/1‑69 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                    >  1:2446567/1‑69 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                     <  1:1022141/68‑1 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                     <  1:2783062/68‑1 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                     <  1:2115582/68‑1 (MQ=255)

                                        cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                     <  1:655152/68‑1 (MQ=255)

                                        cAGCCTATCCTTACGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                     <  1:1852525/68‑1 (MQ=255)

                                                 cTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACg           >  1:2490550/1‑69 (MQ=255)

                                                 cTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACg           >  1:1900442/1‑69 (MQ=255)

                                                 cTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACg           >  1:1246199/1‑69 (MQ=255)

                                                       gTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGAt      <  1:642634/68‑1 (MQ=255)

                                                          cTTTCGCAGGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTg  >  1:2488653/1‑69 (MQ=255)

                                                          cTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTg  >  1:1390630/1‑69 (MQ=255)

                                                          cTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTg  >  1:916539/1‑69 (MQ=255)

CTCTGCCGGACTGGCAACATTTGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTG  >  minE/997515‑997641

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