BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,013	G→A	32.1%	S24S (AGC→AGT)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,013	0	G	A	32.1%	177.5 / 112.5	140	S24S (AGC→AGT)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base G (41/54); new base A (19/26); total (60/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

GGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAT  >  minE/454944‑455080

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

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

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

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

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

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

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

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

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

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

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

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

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

                    ttAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTg                                                 <  1:414974/70‑1 (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)

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

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

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

                          gATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcgc                                           >  1:1619090/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:3263095/69‑1 (MQ=255)

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

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

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

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

                                 gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAgta                                                             >  1:507627/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)

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

                                  agaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGatat                                   <  1:3771237/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)

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

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

                                     tAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCTCAACGACATAGTa                               <  1:1553546/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)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                               gcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             >  1:494033/1‑63 (MQ=255)

                                                taaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATcgcgcgcg                    >  1:4027102/4‑71 (MQ=255)

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

                                                cagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                    >  1:925877/1‑71 (MQ=255)

                                                 agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACa                                    <  1:283032/54‑1 (MQ=255)

                                                 agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                   >  1:2968215/1‑71 (MQ=255)

                                                 agGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCgg                   >  1:1929373/1‑71 (MQ=255)

                                                  aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGAttt                                                  >  1:1893600/2‑39 (MQ=21)

                                                  aggCACGCATGCCGTAAAAACTGAAGTACTCCCAGAttt                                                  >  1:501213/2‑39 (MQ=21)

                                                     cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt               <  1:2786885/71‑1 (MQ=255)

                                                     cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt               <  1:4509645/71‑1 (MQ=255)

                                                     cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt               <  1:4503315/71‑1 (MQ=255)

                                                     cACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt               <  1:3778408/71‑1 (MQ=255)

                                                      aCGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                              >  1:5562898/1‑39 (MQ=255)

                                                       cGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa             <  1:70855/71‑1 (MQ=255)

                                                         cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAAc                                        >  1:3208739/1‑42 (MQ=255)

                                                         cATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGTTGTGATg           >  1:3128724/1‑71 (MQ=255)

                                                         cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAg                                 <  1:3572438/49‑1 (MQ=255)

                                                         cATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg           <  1:121246/71‑3 (MQ=255)

                                                           tGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                               <  1:4810159/49‑1 (MQ=255)

                                                           tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGagg           <  1:1998365/69‑3 (MQ=255)

                                                             ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCa                                          >  1:4173295/1‑36 (MQ=255)

                                                             ccGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtg       >  1:442478/1‑71 (MQ=255)

                                                                tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                <  1:1091625/43‑1 (MQ=255)

                                                                tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                <  1:941405/43‑1 (MQ=255)

                                                                tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                             <  1:4982778/46‑1 (MQ=255)

                                                                tAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCtgt               <  1:5459242/60‑1 (MQ=255)

                                                                  aaaaCTGAAGTACTCCCAGATTTGTATCGCCACGATATAGTATATCGCGCGCGGCTGTGag            >  1:3834926/1‑60 (MQ=255)

                                                                  aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAg                                 <  1:236756/40‑1 (MQ=255)

                                                                  aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc  <  1:4056179/70‑3 (MQ=255)

                                                                  aaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGA‑GGGTGttttc  <  1:1798406/70‑3 (MQ=255)

                                                                  aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACAt                                   >  1:776803/1‑38 (MQ=255)

                                                                  aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAt  >  1:3473302/1‑71 (MQ=255)

                                                                     gCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGc                  <  1:2532398/52‑1 (MQ=255)

GGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAT  >  minE/454944‑455080

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