BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	349,124	C→T	77.6%	A474A (GCC→GCT)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	349,124	0	C	T	77.6%	194.0 / 55.1	107	A474A (GCC→GCT)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand): ref base C (16/8); new base T (33/50); total (49/58)
Fisher's exact test for biased strand distribution p-value = 3.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 108 total aligned reads displayed.

GATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAG  >  minE/349058‑349193

gATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATc                                                                   >  1:5427156/1‑71 (MQ=255)

 aTTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCt                                                                  >  1:4302353/1‑71 (MQ=255)

  ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTg                                                                 <  1:3015247/71‑1 (MQ=255)

  ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTg                                                                 >  1:384995/1‑71 (MQ=255)

  ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTg                                                                 >  1:4274469/1‑71 (MQ=255)

  ttGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTg                                                                 >  1:488134/1‑71 (MQ=255)

       gagTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCt                                                            <  1:814349/71‑1 (MQ=255)

        agTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTAt                                                                    <  1:3576298/62‑1 (MQ=255)

            aGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTc                                                       >  1:4159849/1‑71 (MQ=255)

              cTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                     <  1:5021158/71‑1 (MQ=255)

              cTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCgt                                                     <  1:2902299/71‑1 (MQ=255)

              cTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCgt                                                     <  1:853784/71‑1 (MQ=255)

              cTGTTGGCGCGCAGTTCAACATTACGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCg                                                      <  1:3130808/70‑1 (MQ=255)

                 ttGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCt                                                            <  1:537883/61‑1 (MQ=255)

                 ttGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCt                                                            <  1:4005341/61‑1 (MQ=255)

                 ttGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtca                                                   <  1:803158/70‑1 (MQ=255)

                 ttGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtca                                                   <  1:3593146/70‑1 (MQ=255)

                   ggCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGa                                                          >  1:80165/1‑61 (MQ=255)

                   ggCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGa                                                          >  1:3504064/1‑61 (MQ=255)

                   ggCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTc                                                       >  1:890933/1‑64 (MQ=255)

                   ggCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtcagt                                                 <  1:1479917/70‑1 (MQ=255)

                   ggCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtcagt                                                 <  1:1681290/70‑1 (MQ=255)

                   ggCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtcagt                                                 <  1:1654593/70‑1 (MQ=255)

                   ggCGCGCAGTTCAACATTCCGAACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGa                                                          >  1:1567707/1‑61 (MQ=255)

                     cgcgcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTcac                                              <  1:4585342/71‑1 (MQ=255)

                     cgcgcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTcac                                              <  1:5284972/71‑1 (MQ=255)

                      gcgcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCt                                                            <  1:86889/56‑1 (MQ=255)

                        gcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGa                                                          >  1:2883855/1‑56 (MQ=255)

                        gcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGa                                                          >  1:3948617/1‑56 (MQ=255)

                        gcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGa                                                          >  1:5396901/1‑56 (MQ=255)

                        gcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGa                                                          >  1:4324073/1‑56 (MQ=255)

                        gcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAAc                                           >  1:4271276/1‑71 (MQ=255)

                          aGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAAcgc                                         <  1:542768/71‑1 (MQ=255)

                            ttCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAAcg                                          <  1:660983/68‑1 (MQ=255)

                             tCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTgg                                                                >  1:4593427/1‑45 (MQ=255)

                              cAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtc                                                    >  1:3948229/1‑56 (MQ=255)

                              cAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtc                                                    >  1:2556930/1‑56 (MQ=255)

                              cAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtc                                                    >  1:5189976/1‑56 (MQ=255)

                              cAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtc                                                    >  1:5196300/1‑56 (MQ=255)

                              cAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtt                                     >  1:2696464/1‑71 (MQ=255)

                               aaCATTCCGTACATCCATGTGCTGGTCAACAAGGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtt                                     <  1:2315175/70‑1 (MQ=255)

                                aCATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtttt                                   >  1:2902938/1‑71 (MQ=255)

                                aCATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtttt                                   <  1:3293100/71‑1 (MQ=255)

                                  aTTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGa                                 <  1:4844394/71‑1 (MQ=255)

                                   ttCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCtt                                     <  1:4736032/66‑1 (MQ=255)

                                   ttCCGTACATCCATGTGCTGGTCAACAACGCCTATCTggg                                                               >  1:573620/1‑40 (MQ=255)

                                     ccGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACa                               <  1:2472825/70‑1 (MQ=255)

                                      cGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATg                             <  1:2866424/71‑1 (MQ=255)

                                      cGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATg                             <  1:4130819/71‑1 (MQ=255)

                                      cGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATg                             <  1:4168380/71‑1 (MQ=255)

                                      cGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATg                             <  1:3555836/71‑1 (MQ=255)

                                       gTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                     <  1:3803340/46‑1 (MQ=255)

                                       gTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                     <  1:4585160/46‑1 (MQ=255)

                                       gTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                     <  1:2357131/46‑1 (MQ=255)

                                       gTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATgg                            <  1:2797179/71‑1 (MQ=255)

                                       gTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATgg                            <  1:47964/71‑1 (MQ=255)

                                        tACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGa                           <  1:1627715/71‑1 (MQ=255)

                                         aCATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACaa                                            >  1:215063/1‑53 (MQ=255)

                                         aCATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACaa                                            >  1:2269602/1‑53 (MQ=255)

                                         aCATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACaa                                            >  1:2664973/1‑53 (MQ=255)

                                          cATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGact                         <  1:1882733/71‑1 (MQ=255)

                                          cATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGac                          <  1:3517106/70‑1 (MQ=255)

                                          cATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGac                          >  1:4401648/1‑70 (MQ=255)

                                          cATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGac                          >  1:5544318/1‑70 (MQ=255)

                                            tCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                     <  1:5540483/41‑1 (MQ=255)

                                             ccATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGac                          <  1:3089526/67‑1 (MQ=255)

                                             ccATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGac                          <  1:2799972/67‑1 (MQ=255)

                                              cATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTg                     >  1:5359615/1‑71 (MQ=255)

                                              cATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTg                     >  1:4593348/1‑71 (MQ=255)

                                              cATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTg                     >  1:4612773/1‑71 (MQ=255)

                                               aTGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGc                    <  1:5456635/71‑1 (MQ=255)

                                                tgtgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTg                                  >  1:4991993/1‑56 (MQ=255)

                                                tgtgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTg                                  >  1:5677024/1‑56 (MQ=255)

                                                tgtgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTg                                  >  1:1490777/1‑56 (MQ=255)

                                                 gtgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCg                   >  1:3036659/1‑70 (MQ=255)

                                                  tgCTGGTCAACAACGCTTATCTGGGGCTGATTCgtca                                                   >  1:3124895/1‑37 (MQ=255)

                                                  tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTcac                                              <  1:5382114/42‑1 (MQ=255)

                                                  tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACaa                                            >  1:301230/1‑44 (MQ=255)

                                                  tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAAcgcg                                        >  1:2582397/1‑48 (MQ=255)

                                                  tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGa                           >  1:34781/1‑61 (MQ=255)

                                                  tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGt                  >  1:1027475/1‑70 (MQ=255)

                                                  tgCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTg                 >  1:3573024/1‑71 (MQ=255)

                                                  tgCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTg                 >  1:4824203/1‑71 (MQ=255)

                                                  tgCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTg                 >  1:582295/1‑71 (MQ=255)

                                                  tgCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTg                 >  1:5297173/1‑71 (MQ=255)

                                                    cTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCTTGCa               <  1:1382902/71‑1 (MQ=255)

                                                       gTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAAc                                           <  1:3927588/40‑1 (MQ=255)

                                                       gTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAAc                                           <  1:3650355/40‑1 (MQ=255)

                                                       gTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAAc                                           <  1:3435317/40‑1 (MQ=255)

                                                        tCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTc           <  1:1431972/71‑1 (MQ=255)

                                                        tCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTc           <  1:1406987/71‑1 (MQ=255)

                                                           acaacGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACa                               >  1:3483272/1‑48 (MQ=255)

                                                           acaacGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCt        >  1:4692327/1‑71 (MQ=255)

                                                           acaacGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCt        >  1:2489668/1‑71 (MQ=255)

                                                            caacGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGactact                      <  1:2011550/56‑1 (MQ=255)

                                                             aacGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCttt      <  1:2596249/71‑1 (MQ=255)

                                                              acGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGc                <  1:3990786/60‑1 (MQ=255)

                                                              acGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGc                <  1:2088721/60‑1 (MQ=255)

                                                              acGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGc                <  1:1852781/60‑1 (MQ=255)

                                                                 ccTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCgag  <  1:4699981/71‑1 (MQ=255)

GATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAG  >  minE/349058‑349193

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