BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,404,992	G→A	26.9%	F823F (TTC→TTT)	cyaA ←	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,404,992	0	G	A	26.9%	110.9 / 37.3	96	F823F (TTC→TTT)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base G (28/42); new base A (13/13); total (41/55)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 119 total aligned reads displayed.

TGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTTG  >  minE/2404930‑2405058

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

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

tGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                               <  1:785880/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:32044/1‑64 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                    aTCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGcggaacgg                                                               >  1:1524797/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:321248/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                        tGTTACCGATAGATTTGGTGCGAAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATa                        >  1:2780453/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                              cGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGc                 <  1:56367/68‑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:2287194/1‑69 (MQ=255)

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

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

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

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

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

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

                                                            cggaacggaaTCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTtg  >  1:1160725/1‑69 (MQ=255)

                                                            cggaacggaaTCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTtg  >  1:838998/1‑69 (MQ=255)

                                                            cggaacggaaTCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTtg  >  1:1908820/1‑69 (MQ=255)

                                                            cggaacggaaTCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTtg  >  1:1278876/1‑69 (MQ=255)

                                                            cggaacggaaTCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTtg  >  1:577903/1‑69 (MQ=255)

                                                            cggaacggaaTCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTtg  >  1:3320899/1‑69 (MQ=255)

                                                            cggaacggaaTCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTtg  >  1:2918657/1‑69 (MQ=255)

                                                            cggaacggaaTCACCTGTTAACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTtg  >  1:3070263/1‑69 (MQ=255)

TGTAATAGCGGCGTATCGTGATCCTGATTGGCAGGCGGCATGTTACCGATAGATTTTGTGCGGAACGGAATCACCTGTTCACGACCATCAACCTTCACAATCTGATAGAACTGCGGCAGGTTGAAGTTG  >  minE/2404930‑2405058

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