BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,055	T→G	15.4%	V733V (GTA→GTC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,055	0	T	G	15.4%	140.4 / 24.8	65	V733V (GTA→GTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (31/24); new base G (5/5); total (36/29)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCAT  >  minE/607987‑608123

aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATac                                                                    <  1:2826648/71‑1 (MQ=255)

 aGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATac                                                                    <  1:3332902/70‑1 (MQ=255)

   gTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                 <  1:3737148/71‑1 (MQ=255)

   gTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                 <  1:2535925/71‑1 (MQ=255)

   gTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                 >  1:2028385/1‑71 (MQ=255)

    tACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                                <  1:3117734/71‑1 (MQ=255)

    tACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                                                <  1:283522/71‑1 (MQ=255)

      cGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                                               <  1:827193/70‑1 (MQ=255)

      cGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                                               <  1:2488753/70‑1 (MQ=255)

       gAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                                              >  1:591593/1‑70 (MQ=255)

         gaaataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGt                                                           >  1:1494929/3‑71 (MQ=255)

         gTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                           >  1:580131/1‑71 (MQ=255)

         gTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                           >  1:964003/1‑71 (MQ=255)

          taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                           >  1:2108558/1‑70 (MQ=255)

           aataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATaca                                                                   >  1:2941795/1‑61 (MQ=255)

            ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                        >  1:1002818/1‑71 (MQ=255)

            ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                        >  1:1282906/1‑71 (MQ=255)

               aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                      <  1:2502981/70‑1 (MQ=255)

               aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                      <  1:2694977/70‑1 (MQ=255)

               aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                      <  1:264757/70‑1 (MQ=255)

                 gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATaca                                                                   >  1:2011968/1‑55 (MQ=255)

                   ccTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                                                 <  1:705950/55‑1 (MQ=255)

                   ccTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                 >  1:145123/1‑71 (MQ=255)

                    ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                 <  1:910165/70‑1 (MQ=255)

                        tgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                                            >  1:3792193/1‑71 (MQ=255)

                        tgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                                            >  1:1300684/1‑71 (MQ=255)

                          ctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATg                                                             >  1:2330520/1‑52 (MQ=255)

                             ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                        <  1:1404730/70‑1 (MQ=255)

                             ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                        <  1:2531698/70‑1 (MQ=255)

                             ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                        <  1:2198588/70‑1 (MQ=255)

                                  tAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                           >  1:3532970/1‑46 (MQ=255)

                                  tAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGttt                                  >  1:941510/1‑71 (MQ=255)

                                     tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                        <  1:2401492/62‑1 (MQ=255)

                                       aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                                              >  1:1210512/1‑38 (MQ=255)

                                       aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                                              >  1:2091158/1‑38 (MQ=255)

                                       aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                             >  1:1828542/1‑71 (MQ=255)

                                         cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTg                                                     <  1:2570363/45‑1 (MQ=255)

                                            gatGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTACCGGGtttt                                 >  1:1958951/3‑62 (MQ=255)

                                             ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                   >  1:2265872/1‑43 (MQ=255)

                                               gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                          >  1:1997915/1‑50 (MQ=255)

                                               gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTc                       >  1:2861031/1‑69 (MQ=255)

                                               gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                      <  1:2574956/70‑1 (MQ=255)

                                                aaaCGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                      <  1:248729/37‑1 (MQ=255)

                                                aaaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCCTCGag                    <  1:589931/71‑1 (MQ=255)

                                                 aaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGat                                                <  1:1339895/42‑1 (MQ=255)

                                                 aaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                      >  1:1936877/1‑68 (MQ=255)

                                                  aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                    >  1:529629/1‑69 (MQ=255)

                                                  aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                    >  1:566372/1‑69 (MQ=255)

                                                    gTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                  >  1:3640/1‑69 (MQ=255)

                                                      aCTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAt              >  1:3378394/1‑71 (MQ=255)

                                                       cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt             <  1:3645971/69‑1 (MQ=12)

                                                       cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt             <  1:1265563/69‑1 (MQ=12)

                                                         gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt            >  1:2362056/1‑70 (MQ=255)

                                                          cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa          >  1:2963865/1‑71 (MQ=255)

                                                              gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt      <  1:3744660/71‑1 (MQ=255)

                                                              gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg                                      >  1:931519/1‑39 (MQ=9)

                                                              gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg                                      >  1:75985/1‑39 (MQ=9)

                                                              gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGtgc                                  >  1:1021577/1‑41 (MQ=9)

                                                                cAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga       <  1:2729242/68‑2 (MQ=9)

                                                                 agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCg        <  1:263141/66‑1 (MQ=9)

                                                                 aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                                   <  1:871208/39‑1 (MQ=255)

                                                                 aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca   >  1:396873/1‑71 (MQ=255)

                                                                 aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca   <  1:385533/71‑1 (MQ=255)

                                                                  taTACACGGATGGTCAGCTGCGGATATTTTTCCGGGttt                                  >  1:454934/2‑39 (MQ=255)

                                                                  gATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt            <  1:303918/61‑1 (MQ=255)

                                                                   aTACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat  >  1:3380607/1‑70 (MQ=255)

AAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCAT  >  minE/607987‑608123

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