BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,090	T→G	16.4%	N722H (AAC→CAC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,090	0	T	G	16.4%	113.5 / 22.2	55	N722H (AAC→CAC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (26/20); new base G (2/7); total (28/27)
Fisher's exact test for biased strand distribution p-value = 7.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGA  >  minE/608021‑608159

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

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

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

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

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

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

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

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

                aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                        >  1:529629/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)

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

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

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

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

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

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

                                gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                      <  1:803100/68‑1 (MQ=9)

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

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

                                   acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                   >  1:1149639/1‑71 (MQ=255)

                                    cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                   <  1:251882/70‑1 (MQ=255)

                                     acGGATGGTCAGCTGCGGAAATTTTTCCGGGTTTTCCATCGCGTCGagcagc                                                    >  1:25087/1‑52 (MQ=255)

                                       ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg                               <  1:1812077/71‑2 (MQ=255)

                                       ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                >  1:2612834/1‑70 (MQ=255)

                                            gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                <  1:731714/49‑1 (MQ=255)

                                             tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                             <  1:160098/51‑1 (MQ=255)

                                             tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                         <  1:928150/71‑1 (MQ=255)

                                              cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                          >  1:254837/1‑69 (MQ=255)

                                               aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                       >  1:1010417/1‑71 (MQ=255)

                                               aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                       >  1:3155460/1‑71 (MQ=255)

                                               aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                       <  1:3241491/71‑1 (MQ=255)

                                                gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                     >  1:2916952/1‑56 (MQ=255)

                                                gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                        <  1:1912623/69‑1 (MQ=255)

                                                 cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                     <  1:2538351/71‑1 (MQ=255)

                                                 cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                     <  1:3583132/71‑1 (MQ=255)

                                                 cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                     >  1:40621/1‑71 (MQ=255)

                                                   taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                   >  1:1934025/3‑70 (MQ=255)

                                                   gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                   >  1:3690568/1‑71 (MQ=255)

                                                     ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                             >  1:2967272/1‑43 (MQ=255)

                                                      gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                <  1:1929460/71‑1 (MQ=255)

                                                      gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                >  1:616967/1‑71 (MQ=12)

                                                         atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGt                                 <  1:1328372/51‑1 (MQ=255)

                                                         atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca             >  1:564393/1‑71 (MQ=255)

                                                          tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac            <  1:2149654/71‑1 (MQ=255)

                                                           ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc           >  1:855116/1‑71 (MQ=255)

                                                           ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc           >  1:534188/1‑71 (MQ=255)

                                                           ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc           >  1:3251582/1‑71 (MQ=255)

                                                            tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc           <  1:3278536/70‑1 (MQ=255)

                                                                cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga      <  1:868636/71‑1 (MQ=255)

                                                                cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga      <  1:1278827/71‑1 (MQ=255)

                                                                 gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                      <  1:558239/38‑1 (MQ=255)

                                                                 gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat     <  1:1627687/71‑1 (MQ=255)

                                                                   gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg    >  1:1553947/1‑70 (MQ=255)

                                                                    ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                       <  1:1953432/50‑1 (MQ=255)

                                                                    ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                       <  1:237535/50‑1 (MQ=255)

                                                                    ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga  >  1:1659397/1‑71 (MQ=255)

TAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGA  >  minE/608021‑608159

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