BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,022	T→C	34.7%	E21E (GAA→GAG)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,022	0	T	C	34.7%	115.8 / 106.3	118	E21E (GAA→GAG)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base T (38/39); new base C (15/26); total (53/65)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.97e-01

Only 100 of 123 total aligned reads displayed.

TATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTC  >  minE/454952‑455092

    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)

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

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

                 tGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                                                        >  1:1587338/1‑70 (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)

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

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

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

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

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

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

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

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

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

                           gaTAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAg                                             <  1:265496/71‑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)

                             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:2970162/1‑62 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                      agcagGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                                  >  1:626790/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)

                                         aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                               >  1:1049272/3‑71 (MQ=255)

                                         aaggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCgg                               >  1:7948/3‑71 (MQ=255)

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

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

                                            gCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                            <  1:3150876/71‑1 (MQ=255)

                                             cACGCATGCCGTAAAAGCTAAAATATTTCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                           <  1:1217570/71‑1 (MQ=255)

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

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

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

                                             cACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtg                            <  1:1933962/70‑1 (MQ=255)

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

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

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

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

                                               cGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGa                         <  1:294826/71‑1 (MQ=255)

                                                gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                        >  1:2119683/1‑70 (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:4221607/1‑36 (MQ=255)

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

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

                                                        tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATa                                              <  1:4168569/41‑1 (MQ=255)

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

                                                        tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                            <  1:1091625/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:3652356/70‑3 (MQ=255)

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

                                                          aaaGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                      >  1:1839905/1‑47 (MQ=255)

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

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

                                                                aaaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTa               >  1:4483838/1‑64 (MQ=255)

                                                                 aaaTATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat       <  1:4807160/71‑1 (MQ=255)

                                                                      ttCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTc  <  1:4033631/71‑1 (MQ=255)

TATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTC  >  minE/454952‑455092

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