BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	454,989	C→G	33.3%	L32L (CTG→CTC)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	454,989	0	C	G	33.3%	111.6 / 82.2	96	L32L (CTG→CTC)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (23/41); new base G (13/19); total (36/60)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

ACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGC  >  minE/454920‑455059

aCAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATca                                                                       <  1:5640953/71‑1 (MQ=255)

aCAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATc                                                                        <  1:1336033/70‑1 (MQ=255)

   gTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                    <  1:2973944/71‑1 (MQ=255)

   gTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                    <  1:4067698/71‑1 (MQ=255)

           tAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                    <  1:4240205/63‑1 (MQ=255)

           tAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagc                                                                     <  1:4193495/62‑1 (MQ=255)

             aaTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagcag                                                                   >  1:579540/1‑62 (MQ=255)

             aaTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                          >  1:4850534/1‑71 (MQ=255)

               tAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCGGCAGGGCACGCATGc                                                        >  1:5521668/1‑71 (MQ=255)

                  cTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGt                                                     <  1:2548878/71‑1 (MQ=255)

                        ggCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGc                                               >  1:1243234/1‑71 (MQ=255)

                           gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAATAGCTaa                                            <  1:2381/71‑1 (MQ=255)

                           gTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                            <  1:3644714/71‑1 (MQ=255)

                            tGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagc                                                                     <  1:2320098/45‑1 (MQ=255)

                            tGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaa                                           <  1:1362548/71‑1 (MQ=255)

                              aGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAat                                         <  1:3626119/71‑1 (MQ=255)

                                tATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                       <  1:3961095/71‑1 (MQ=255)

                                 aTCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                            <  1:3515809/65‑1 (MQ=255)

                                  tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATca                                                                       >  1:143323/1‑37 (MQ=255)

                                  tCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                            <  1:46847/64‑1 (MQ=255)

                                   cGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                            <  1:4189223/63‑1 (MQ=255)

                                    gTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                    >  1:2946076/1‑70 (MQ=255)

                                    gTTGTATTTTAGATGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATAtt                                      >  1:5247165/1‑68 (MQ=255)

                                      tGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                    >  1:776356/1‑68 (MQ=255)

                                      tGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCa                                  <  1:5350036/70‑1 (MQ=255)

                                           tttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                    <  1:1029401/63‑1 (MQ=255)

                                            ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                            <  1:414974/70‑1 (MQ=255)

                                               gTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                         <  1:5079031/70‑1 (MQ=255)

                                                ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                            <  1:2818561/50‑1 (MQ=255)

                                                ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTa                                             <  1:5027261/49‑1 (MQ=255)

                                                ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                       >  1:2372739/1‑71 (MQ=255)

                                                ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                         <  1:1774744/69‑1 (MQ=255)

                                                ttGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                         <  1:2972474/69‑1 (MQ=255)

                                                 tGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                       >  1:1587338/1‑70 (MQ=255)

                                                  gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                          >  1:1941796/1‑50 (MQ=255)

                                                  gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                          >  1:5095262/1‑50 (MQ=255)

                                                  gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                          >  1:2646619/1‑50 (MQ=255)

                                                  gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                      >  1:596092/1‑70 (MQ=255)

                                                  gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                      >  1:1619090/1‑70 (MQ=255)

                                                  gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                      >  1:3082149/1‑70 (MQ=255)

                                                    ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                          <  1:686172/48‑1 (MQ=255)

                                                    ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                          <  1:1316045/48‑1 (MQ=255)

                                                    ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                          <  1:4013472/48‑1 (MQ=255)

                                                    ttGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                   >  1:2081188/1‑71 (MQ=255)

                                                    tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                   <  1:3862299/69‑1 (MQ=255)

                                                    tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                   <  1:3353582/69‑1 (MQ=255)

                                                    tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                   <  1:3263095/69‑1 (MQ=255)

                                                     gggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACtg                                             <  1:1240955/43‑2 (MQ=25)

                                                     gggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTc                                     >  1:5201728/2‑52 (MQ=255)

                                                       gTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAt                               >  1:4363715/1‑56 (MQ=255)

                                                        tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGAttt                             <  1:3099311/57‑1 (MQ=255)

                                                        tGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa               <  1:1790594/71‑1 (MQ=255)

                                                         gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                        >  1:5553697/1‑45 (MQ=255)

                                                         gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                        >  1:5420728/1‑45 (MQ=255)

                                                         gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                        >  1:507627/1‑45 (MQ=255)

                                                         gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                        >  1:5342075/1‑45 (MQ=255)

                                                         gagaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaa                                           <  1:1035054/42‑1 (MQ=255)

                                                          agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACt                                              <  1:447934/38‑1 (MQ=21)

                                                          agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACt                                              <  1:1822730/38‑1 (MQ=21)

                                                          agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACt                                              <  1:2422257/38‑1 (MQ=21)

                                                          agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat              <  1:3771237/70‑1 (MQ=255)

                                                          agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat              <  1:3291639/70‑1 (MQ=255)

                                                           gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg            <  1:265496/71‑1 (MQ=255)

                                                           gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg            >  1:3082007/1‑71 (MQ=255)

                                                            aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTgaa                                           >  1:4833002/1‑39 (MQ=16)

                                                            aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGAt                               <  1:2419711/51‑1 (MQ=255)

                                                            aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGAt                               <  1:2896687/51‑1 (MQ=255)

                                                            aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGAt                               <  1:3829380/51‑1 (MQ=255)

                                                            aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata             >  1:3039509/1‑69 (MQ=255)

                                                            aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata             <  1:2225637/69‑1 (MQ=255)

                                                            aTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg            >  1:2010975/1‑70 (MQ=255)

                                                             tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata             <  1:5071689/68‑1 (MQ=255)

                                                             tAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt           <  1:4606702/70‑1 (MQ=255)

                                                             tAGAGAATCAGCAGTGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa          <  1:1399617/71‑1 (MQ=255)

                                                             tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa                                            <  1:2885980/37‑1 (MQ=255)

                                                             tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCTCAACGACATAGTa          <  1:1553546/71‑1 (MQ=255)

                                                             tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTACCAGATTTGCAGCGCAACGACATAGTa          <  1:4030605/71‑1 (MQ=255)

                                                              agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTa          >  1:136720/1‑70 (MQ=255)

                                                              agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa          <  1:4914839/70‑1 (MQ=255)

                                                               gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata               >  1:3129678/1‑62 (MQ=255)

                                                               gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt           >  1:1350989/1‑68 (MQ=255)

                                                               gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGt           >  1:4002479/1‑68 (MQ=255)

                                                               gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTata        <  1:700401/71‑3 (MQ=255)

                                                               gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTat         <  1:5468723/70‑2 (MQ=255)

                                                               gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTat         <  1:700044/70‑2 (MQ=255)

                                                               gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa        <  1:762769/71‑1 (MQ=255)

                                                               gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa        <  1:2087903/71‑1 (MQ=255)

                                                               gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa        <  1:4845549/71‑1 (MQ=255)

                                                               gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa        <  1:4334378/71‑1 (MQ=255)

                                                               gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa        <  1:3796277/71‑1 (MQ=255)

                                                               gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa        <  1:3352644/71‑1 (MQ=255)

                                                               gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa        >  1:328063/1‑71 (MQ=255)

                                                                agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatat       <  1:386759/71‑1 (MQ=255)

                                                                agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa        >  1:4696071/1‑70 (MQ=255)

                                                                 gaATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa                                          <  1:1031592/35‑1 (MQ=255)

                                                                  aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa             >  1:1265823/1‑63 (MQ=255)

                                                                  aaTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa             >  1:4607506/1‑63 (MQ=255)

                                                                   aTCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg     <  1:1708427/70‑1 (MQ=255)

                                                                     cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc  >  1:4318729/1‑71 (MQ=255)

                                                                     cagcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgc  >  1:2903474/1‑71 (MQ=255)

ACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGC  >  minE/454920‑455059

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