BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,544	A→G	29.7%	N570N (AAT→AAC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,544	0	A	G	29.7%	106.9 / 54.1	114	N570N (AAT→AAC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (52/28); new base G (12/22); total (64/50)
Fisher's exact test for biased strand distribution p-value = 4.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 114 total aligned reads displayed.

CAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTC  >  minE/608478‑608606

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:2014258/69‑1 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:2459678/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:2508690/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:1455261/69‑1 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:1417187/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:2715336/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:1667174/69‑1 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:167103/69‑1 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:2306635/69‑1 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:3110821/69‑1 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:1009625/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:1276067/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:474430/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:478428/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              >  1:111156/1‑69 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:1965662/69‑1 (MQ=255)

cAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:1063826/69‑1 (MQ=255)

 aGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:1105992/68‑1 (MQ=255)

         aaTTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCAtt                                                              <  1:2851716/60‑1 (MQ=255)

           ttttCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCaa                                                      <  1:2268814/66‑1 (MQ=255)

           ttttCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCaa                                                      <  1:3384846/66‑1 (MQ=255)

                 tcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCggg                                              <  1:2246568/68‑1 (MQ=255)

                 tcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCggg                                              <  1:1350150/68‑1 (MQ=255)

                 tcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCggg                                              <  1:2734489/68‑1 (MQ=255)

                 tcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCggg                                              <  1:2183986/68‑1 (MQ=255)

                 tcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCggg                                              <  1:2447876/68‑1 (MQ=255)

                 tcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATAATTGTTACCAAACTGCggg                                              <  1:3112824/68‑1 (MQ=255)

                     gAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTc                                         <  1:210667/69‑1 (MQ=255)

                      aaaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCg                                        <  1:2838270/69‑1 (MQ=255)

                      aaaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCg                                        <  1:3279252/69‑1 (MQ=255)

                      aaaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCg                                        <  1:1204878/69‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCGAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:1293408/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:264832/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:2345610/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:2882597/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:2968811/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:602135/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:946365/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:946635/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:2301049/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:1333704/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:1099815/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:1338779/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:1553543/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:1665578/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:1745741/58‑1 (MQ=255)

                       aaCGTTCTACCAGGTCAACAGCCAGGTCATCTAAACGCGGATCGTTGTTACCAAACTg                                                  <  1:3045952/58‑1 (MQ=255)

                                ccAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTGCGGGTATTCGCCTTCGAtt                               >  1:595690/1‑68 (MQ=255)

                                ccAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTGCGGGTATTCGCCTTCGAtt                               >  1:1180807/1‑68 (MQ=255)

                                ccAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTGCGGGTATTCGCCTTCGAtt                               >  1:1511191/1‑68 (MQ=255)

                                ccAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTGCGGGTATTCGCCTTCGAtt                               >  1:2378103/1‑68 (MQ=255)

                                ccAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTGCGGGTATTCGCCTTCGAtt                               >  1:2725140/1‑68 (MQ=255)

                                ccAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTGCGGGTATTCGCCTTCGAtt                               >  1:1347381/1‑68 (MQ=255)

                                ccAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTGCGGGTATTCGCCTTCGAtt                               >  1:71791/1‑68 (MQ=255)

                                ccAGGTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTGCGGGTATTCGCCTTCGAtt                               >  1:186214/1‑68 (MQ=255)

                                   ggTCAAGAGCCAGGTCATCTACACGCGGATCATTGTTACCaa                                                      <  1:665456/42‑1 (MQ=255)

                                   ggTCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:610730/46‑1 (MQ=255)

                                     tCAACAGCCAGGTCATCTACACGCGGATCGTTGTTACCAAACTg                                                  <  1:1110268/44‑1 (MQ=255)

                                                gTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTct                                        <  1:2864401/43‑2 (MQ=255)

                                                 tCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAg              >  1:2093414/1‑68 (MQ=255)

                                                 tCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGa             >  1:1346818/1‑69 (MQ=255)

                                                 tCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGa             >  1:2190554/1‑69 (MQ=255)

                                                 tCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGa             >  1:3024079/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:2615682/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:2635401/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:2700367/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:3184173/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:1749389/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:129650/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:1778549/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:2878347/1‑69 (MQ=255)

                                                  cATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc            >  1:2254911/1‑69 (MQ=255)

                                                      tACCCGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:10318/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCCAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:83951/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGAt                    <  1:2945173/57‑1 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGAt                    <  1:1135254/57‑1 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGAt                    <  1:2213239/57‑1 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGt         >  1:1782883/1‑68 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGt         >  1:1406598/1‑68 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:1360244/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:3214488/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:1332712/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:331618/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:3344154/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:857367/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:640715/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:432584/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:14523/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:502438/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:1503088/1‑69 (MQ=255)

                                                      tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc        >  1:2501698/1‑69 (MQ=255)

                                                        cacGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt      >  1:2119154/1‑69 (MQ=255)

                                                        cacGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt      >  1:2032682/1‑69 (MQ=255)

                                                        cacGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt      >  1:3385305/1‑69 (MQ=255)

                                                        cacGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt      >  1:1186610/1‑69 (MQ=255)

                                                        cacGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt      >  1:2881538/1‑69 (MQ=255)

                                                        cacGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt      >  1:2885512/1‑69 (MQ=255)

                                                        cacGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt      >  1:2920988/1‑69 (MQ=255)

                                                        aacGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCtt      >  1:1257355/2‑69 (MQ=255)

                                                            cgGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTc  >  1:2214093/1‑69 (MQ=255)

CAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTC  >  minE/608478‑608606

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