BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,316	T→C	8.3%	V646V (GTA→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,316	0	T	C	8.3%	168.5 / 6.6	60	V646V (GTA→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (30/25); new base C (2/3); total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACG  >  minE/608250‑608381

gaAGGAGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                                  >  1:116959/1‑68 (MQ=255)

 aaGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGc                                                              <  1:2397246/71‑1 (MQ=255)

  aGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCa                                                             >  1:66239/1‑71 (MQ=255)

       tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCttt                                                        >  1:1331584/1‑71 (MQ=255)

              tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                                   <  1:433874/69‑1 (MQ=255)

               aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                >  1:3711575/1‑71 (MQ=255)

               aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                >  1:2276328/1‑71 (MQ=255)

               aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                >  1:1997276/1‑71 (MQ=255)

                   tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCt                                                                   <  1:612356/48‑1 (MQ=255)

                   tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                                      <  1:931244/61‑1 (MQ=255)

                   tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                                      <  1:57827/61‑1 (MQ=255)

                   tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                                    <  1:2276715/63‑1 (MQ=255)

                   tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                            >  1:2824860/1‑71 (MQ=255)

                     tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                                          >  1:86378/1‑71 (MQ=255)

                       tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                         >  1:2871084/1‑70 (MQ=255)

                       tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                         >  1:162429/1‑70 (MQ=255)

                       tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                         >  1:2274216/1‑70 (MQ=255)

                       tAGCGTAAGCAAACCGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                         >  1:3170732/1‑70 (MQ=255)

                         gCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                                               <  1:2989302/46‑1 (MQ=255)

                           gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                                    <  1:3083121/71‑1 (MQ=255)

                           gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                                    >  1:3206440/1‑71 (MQ=255)

                            tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                                    <  1:1532937/54‑1 (MQ=255)

                              aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                 >  1:1218867/1‑71 (MQ=255)

                              aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                 >  1:1485182/1‑71 (MQ=255)

                              aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                 >  1:3070467/1‑71 (MQ=255)

                                  aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                                <  1:51290/68‑1 (MQ=255)

                                   aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                                  >  1:1588358/1‑49 (MQ=255)

                                    cGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAc                           >  1:2351881/1‑71 (MQ=255)

                                     ggCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc                          <  1:1757453/71‑1 (MQ=255)

                                       cAGTTTCGCCACCGAGGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCg                        >  1:244307/1‑71 (MQ=255)

                                         gTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCggc                      >  1:442481/1‑70 (MQ=255)

                                          tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                              <  1:3423367/62‑1 (MQ=255)

                                          tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg                              <  1:3418642/62‑1 (MQ=255)

                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                                   >  1:2472862/1‑55 (MQ=255)

                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                    >  1:2777947/1‑70 (MQ=255)

                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                    >  1:2686875/1‑70 (MQ=255)

                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                    >  1:1308062/1‑70 (MQ=255)

                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                    >  1:249141/1‑70 (MQ=255)

                                            tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc                    >  1:816383/1‑70 (MQ=255)

                                             aGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt                                <  1:3357585/57‑1 (MQ=255)

                                              gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGa                  <  1:2898274/70‑1 (MQ=255)

                                               cAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCGCGACCGTGCATCgg                                   <  1:2989599/52‑1 (MQ=255)

                                                aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                                   >  1:1076369/1‑51 (MQ=255)

                                                aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg               >  1:3648859/1‑71 (MQ=255)

                                                aaCGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACg               <  1:1243177/71‑1 (MQ=255)

                                                  cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACTAcc                                            >  1:3876436/1‑40 (MQ=255)

                                                    gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc             <  1:18356/69‑1 (MQ=255)

                                                     aaGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcg          >  1:3561611/1‑71 (MQ=255)

                                                       gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                 <  1:2751199/46‑1 (MQ=255)

                                                       gTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGcc        <  1:3416744/71‑1 (MQ=255)

                                                        tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                         <  1:2555676/37‑1 (MQ=255)

                                                        tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                         <  1:735118/37‑1 (MQ=255)

                                                        tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAg                              <  1:2509218/48‑1 (MQ=255)

                                                         cAGAGAGGCTACTGCACCTTTCTGGTCACGAACGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg      >  1:867294/1‑71 (MQ=255)

                                                          agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                 >  1:3127725/1‑59 (MQ=255)

                                                          agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg      <  1:2643037/70‑1 (MQ=255)

                                                          agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg      <  1:3095710/70‑1 (MQ=255)

                                                          agagagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGc     <  1:1417372/71‑1 (MQ=255)

                                                             cagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacg  <  1:2164731/70‑1 (MQ=255)

                                                              agGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacg  <  1:1962260/70‑1 (MQ=255)

                                                                gCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAg                              <  1:813800/40‑1 (MQ=255)

GAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACG  >  minE/608250‑608381

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`