BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	454,992	C→T	33.3%	L31L (CTG→CTA) ‡	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	454,992	0	C	T	33.3%	112.4 / 78.6	96	L31L (CTG→CTA) ‡	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (23/41); new base T (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 = 9.90e-01

Only 100 of 109 total aligned reads displayed.

GTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCG  >  Exported/454923‑455062

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

gTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagca                                                                       <  1:2973944/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:4791161/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)

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

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

                                 gTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTcc                                       >  1:342150/1‑70 (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)

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

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

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

                                               gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAatat                                          <  1:782649/53‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                 tgggTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCAc                      <  1:3862299/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:5342075/1‑45 (MQ=255)

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

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

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

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

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

                                                       agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACt                                                 <  1:447934/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:3829380/51‑1 (MQ=255)

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

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

                                                         aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGAtata                >  1:3039509/1‑69 (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)

                                                          tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa             <  1:5403146/71‑1 (MQ=255)

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

                                                           agagAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAg                                    >  1:4100467/1‑47 (MQ=255)

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

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

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

                                                            gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGata                  >  1:2970162/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:700044/70‑2 (MQ=255)

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

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

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

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

                                                            gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa           <  1:2087903/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)

                                                                   agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg    >  1:626790/1‑71 (MQ=255)

                                                                   agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg    >  1:1768082/1‑71 (MQ=255)

                                                                     cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                               <  1:2589670/42‑1 (MQ=255)

GTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCG  >  Exported/454923‑455062

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