BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	455,045	C→T	23.1%	V14I (GTC→ATC)	ybgH ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	455,045	0	C	T	23.1%	199.8 / 56.0	108	V14I (GTC→ATC)	ybgH	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (36/47); new base T (9/16); total (45/63)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 114 total aligned reads displayed.

TGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGT  >  minE/454976‑455115

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

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

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

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

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

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

      agagAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTa                                                                  <  1:4914839/70‑1 (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)

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

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

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

       gagaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                >  1:328063/1‑71 (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:2087903/71‑1 (MQ=255)

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

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

        agaATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                                                >  1:4696071/1‑70 (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: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)

                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)

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

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

                    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)

                        gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGag                                               >  1:2119683/1‑70 (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)

                           tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTatat                                                               <  1:4975880/52‑1 (MQ=255)

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

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

                                tAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGt                                                                   <  1:1091625/43‑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)

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

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

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

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

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

                                       tAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCtgt                                                  <  1:515501/53‑1 (MQ=255)

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

                                            tatTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt                                        <  1:1532654/58‑1 (MQ=255)

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

                                                 ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                        >  1:4394748/1‑69 (MQ=255)

                                                 ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac                      <  1:5599699/71‑1 (MQ=255)

                                                 ccAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca                       >  1:3020669/1‑70 (MQ=255)

                                                 ccAGATTTGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCca                       >  1:1175136/1‑70 (MQ=255)

                                                  cAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcaca                     <  1:1510223/71‑1 (MQ=255)

                                                   aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                >  1:2925206/1‑43 (MQ=255)

                                                   aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGa                                                >  1:3796198/1‑43 (MQ=255)

                                                   aGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACaa                    >  1:4977867/1‑71 (MQ=255)

                                                    gATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAt                   >  1:5553338/1‑71 (MQ=255)

                                                     aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                  >  1:1924153/1‑71 (MQ=255)

                                                     aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                  <  1:3458865/71‑1 (MQ=255)

                                                     aTTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                  <  1:5153035/71‑1 (MQ=255)

                                                       ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                >  1:5031569/1‑71 (MQ=255)

                                                       ttGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                >  1:2460373/1‑71 (MQ=255)

                                                       ttGCAGAGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                >  1:737611/1‑71 (MQ=255)

                                                         gCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGa                <  1:5020296/69‑1 (MQ=255)

                                                          cAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACgg             >  1:3548632/1‑71 (MQ=255)

                                                           aGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc            <  1:1001850/71‑1 (MQ=255)

                                                            gcgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGc            <  1:4602369/70‑1 (MQ=255)

                                                             cgcAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCaa          >  1:4018329/1‑71 (MQ=255)

                                                                 aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGtt                                        <  1:940477/37‑1 (MQ=255)

                                                                 aCGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTcc                        <  1:2095013/53‑1 (MQ=255)

                                                                  cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                  <  1:5311807/58‑1 (MQ=255)

                                                                  cGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAtt                  <  1:3388961/58‑1 (MQ=255)

                                                                   gACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCcac                      >  1:1080846/1‑53 (MQ=255)

                                                                   gACATAATAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCAc    <  1:1200868/71‑1 (MQ=255)

                                                                    aCATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTc      <  1:4935546/68‑1 (MQ=255)

                                                                     cATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGt  <  1:4355098/71‑1 (MQ=255)

TGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGT  >  minE/454976‑455115

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