BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,404,986	T→G	25.2%	T825T (ACA→ACC)	cyaA ←	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,404,986	0	T	G	25.2%	137.8 / 54.8	103	T825T (ACA→ACC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base T (34/43); new base G (13/13); total (47/56)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

TGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGA  >  minE/2404930‑2405053

tGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          <  1:1668857/68‑1 (MQ=255)

tGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          <  1:2706517/68‑1 (MQ=255)

      aGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCAc                                                    <  1:3068113/68‑1 (MQ=255)

      aGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCAc                                                    <  1:1946164/68‑1 (MQ=255)

              aTCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACg                                           <  1:1198368/69‑1 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcgga                                                              >  1:2608606/1‑49 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:476195/1‑64 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:238339/1‑64 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:2624485/1‑64 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:1571035/1‑64 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:32044/1‑64 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:1505585/1‑64 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:651183/1‑64 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGtt                                               >  1:931097/1‑64 (MQ=255)

               tCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCa                                             >  1:647273/1‑66 (MQ=255)

                 gTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGAcc                                        >  1:2751557/1‑69 (MQ=255)

                    aTCCTGATTGGCAGTCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATc                                     <  1:2532139/69‑1 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:2131739/1‑48 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:3132327/1‑48 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:3176585/1‑48 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:2790597/1‑48 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:1547821/1‑48 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:1524797/1‑48 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:604976/1‑48 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:1221195/1‑48 (MQ=255)

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:1163715/1‑48 (MQ=255)

                    aTCCTGATTGG‑AGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                          >  1:567355/1‑47 (MQ=255)

                      gcTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:243939/67‑1 (MQ=255)

                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:2727753/68‑1 (MQ=255)

                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:2482668/68‑1 (MQ=255)

                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:2038116/68‑1 (MQ=255)

                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:1868656/68‑1 (MQ=255)

                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:1856178/68‑1 (MQ=255)

                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:1707244/68‑1 (MQ=255)

                      ccTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCa                                    <  1:321248/68‑1 (MQ=255)

                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:2123791/57‑1 (MQ=255)

                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:303901/57‑1 (MQ=255)

                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:3043317/57‑1 (MQ=255)

                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:2745012/57‑1 (MQ=255)

                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:3189284/57‑1 (MQ=255)

                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:3223547/57‑1 (MQ=255)

                         gATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCAc                                            <  1:849830/57‑1 (MQ=255)

                           ttGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCtt                               >  1:103609/1‑68 (MQ=255)

                           ttGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCtt                               >  1:2707007/1‑68 (MQ=255)

                            tGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTc                              >  1:3269990/1‑68 (MQ=255)

                            tGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTc                              >  1:3149616/1‑68 (MQ=255)

                               cAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTcaca                           <  1:787507/68‑1 (MQ=255)

                               cAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTcaca                           <  1:531540/68‑1 (MQ=255)

                               cAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTcaca                           <  1:2194728/68‑1 (MQ=255)

                                  gcggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:3151045/69‑1 (MQ=255)

                                  gcggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:1513856/69‑1 (MQ=255)

                                  gcggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:2140016/69‑1 (MQ=255)

                                  gcggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:859666/69‑1 (MQ=255)

                                   cggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       <  1:1438289/68‑1 (MQ=255)

                                    ggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       >  1:832882/1‑67 (MQ=255)

                                    ggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       >  1:336464/1‑67 (MQ=255)

                                    ggcATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCt                       >  1:2373654/1‑67 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:3084391/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:2409409/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:3389026/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:42716/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:3127075/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:3080485/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:1311517/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:2983003/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGAc                                         <  1:2925019/48‑1 (MQ=255)

                                     gcATGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTGCACGAc                                         <  1:8072/48‑1 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCTAAACGGAATCACCTGTTCACGAc                                         <  1:2711608/45‑1 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGt                                                >  1:1211266/1‑38 (MQ=38)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:270163/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:640615/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:374941/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:70361/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:2504626/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1551112/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:2202101/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:2877556/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1948967/1‑68 (MQ=255)

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1768686/1‑68 (MQ=255)

                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:1400339/66‑1 (MQ=255)

                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:3124877/66‑1 (MQ=255)

                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:34046/66‑1 (MQ=255)

                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:210785/66‑1 (MQ=255)

                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:1253834/66‑1 (MQ=255)

                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:650246/66‑1 (MQ=255)

                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:80633/66‑1 (MQ=255)

                                          ttACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  <  1:785718/66‑1 (MQ=255)

                                              cGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGc            <  1:1310753/68‑1 (MQ=255)

                                              cGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGc            <  1:736466/68‑1 (MQ=255)

                                                 tAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTg                      >  1:1147242/1‑55 (MQ=255)

                                                 tAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1061324/1‑59 (MQ=255)

                                                 tAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:85542/1‑59 (MQ=255)

                                                 tAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAg                  >  1:1629389/1‑59 (MQ=255)

                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  <  1:800696/69‑1 (MQ=255)

                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:1939650/1‑69 (MQ=255)

                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  <  1:2183908/69‑1 (MQ=255)

                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:2270563/1‑69 (MQ=255)

                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:2287194/1‑69 (MQ=255)

                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:2453555/1‑69 (MQ=255)

                                                       ttGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGa  >  1:2471830/1‑69 (MQ=255)

TGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGA  >  minE/2404930‑2405053

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