BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,985	C→T	9.2%	P423P (CCG→CCA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,985	0	C	T	9.2%	299.5 / 15.6	98	P423P (CCG→CCA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (49/40); new base T (5/4); total (54/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 122 total aligned reads displayed.

TTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCA  >  minE/608915‑609054

ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                                                       <  1:1823578/71‑1 (MQ=255)

ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                                                       <  1:910646/71‑1 (MQ=255)

         gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                                  >  1:1288974/1‑67 (MQ=255)

         gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                                              <  1:2648335/71‑1 (MQ=255)

         gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                                              <  1:43335/71‑1 (MQ=255)

         gCATGGTTTTCGCCAGGTTTGCACCCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                                              >  1:3004311/1‑71 (MQ=255)

            tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                                     <  1:1040688/61‑1 (MQ=255)

            tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                                     <  1:2636176/61‑1 (MQ=255)

            tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                                     <  1:4123172/61‑1 (MQ=255)

            gcgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                                           >  1:1398399/3‑71 (MQ=255)

              gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                                         <  1:3498025/71‑1 (MQ=255)

               ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgca                                                            <  1:2918905/67‑1 (MQ=255)

                   cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                                     <  1:1199145/70‑1 (MQ=255)

                   cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                                     <  1:1670833/70‑1 (MQ=255)

                   cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggca                                                     <  1:1900468/70‑1 (MQ=255)

                      cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                                 <  1:1472220/71‑1 (MQ=255)

                       aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                >  1:4119712/1‑71 (MQ=255)

                        ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                                 >  1:2729084/1‑69 (MQ=255)

                        ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                                 >  1:4078426/1‑69 (MQ=255)

                        ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                               <  1:2076321/71‑1 (MQ=255)

                        ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                               >  1:1221440/1‑71 (MQ=255)

                          tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTa                                             <  1:232203/71‑1 (MQ=255)

                           ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                    >  1:1117979/1‑63 (MQ=255)

                           ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                                  >  1:2022989/1‑65 (MQ=255)

                           ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                                            >  1:40416/1‑71 (MQ=255)

                            tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                           >  1:509491/1‑71 (MQ=255)

                            tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                           >  1:2525102/1‑71 (MQ=255)

                            tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                           >  1:3324975/1‑71 (MQ=255)

                            tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                           >  1:2004154/1‑71 (MQ=255)

                               aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                      <  1:3363694/57‑1 (MQ=255)

                               aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCAt                                        <  1:3381657/71‑1 (MQ=255)

                                 gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                    >  1:282008/1‑57 (MQ=255)

                                 gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg                                      <  1:2744558/71‑1 (MQ=255)

                                  tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa                                             <  1:456201/62‑1 (MQ=255)

                                   gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                    >  1:3364338/1‑71 (MQ=255)

                                   gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                    >  1:1211478/1‑71 (MQ=255)

                                   gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                                    >  1:1385780/1‑69 (MQ=255)

                                   gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                                    >  1:808507/1‑69 (MQ=255)

                                      ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                      <  1:2944900/50‑1 (MQ=255)

                                      ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                  >  1:123820/1‑70 (MQ=255)

                                       cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAAt                                                   <  1:1006285/52‑1 (MQ=255)

                                       cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                 >  1:2061036/1‑70 (MQ=255)

                                       cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                 >  1:2056076/1‑70 (MQ=255)

                                       cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                <  1:2539150/71‑1 (MQ=255)

                                       cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                <  1:151874/71‑1 (MQ=255)

                                        gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                 >  1:1840749/1‑69 (MQ=255)

                                        gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                   >  1:528020/1‑67 (MQ=255)

                                        gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                >  1:3036046/1‑70 (MQ=255)

                                         aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                <  1:618985/69‑1 (MQ=255)

                                         aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                >  1:3905432/1‑69 (MQ=255)

                                          agaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAg                             >  1:3157389/1‑71 (MQ=255)

                                           gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                            >  1:1475473/1‑71 (MQ=255)

                                             aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTcc                          <  1:1934522/71‑1 (MQ=255)

                                              cTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                      <  1:4002379/42‑1 (MQ=255)

                                              cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                    >  1:1850420/1‑44 (MQ=255)

                                              cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                         <  1:421580/71‑1 (MQ=255)

                                              cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                         <  1:3254358/71‑1 (MQ=255)

                                                gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                  >  1:765660/1‑60 (MQ=255)

                                                gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                  >  1:1060756/1‑60 (MQ=255)

                                                 cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                    <  1:1817824/41‑1 (MQ=255)

                                                 cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                    <  1:244663/41‑1 (MQ=255)

                                                 cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                >  1:838121/1‑45 (MQ=255)

                                                 cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGAc                      >  1:2272008/1‑71 (MQ=255)

                                                   tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAg                             <  1:643970/62‑1 (MQ=255)

                                                   tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGc                    <  1:3672883/71‑1 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                    >  1:2006482/1‑37 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                    >  1:4081085/1‑37 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCTTAGTCATCGTTGTTGaa                              >  1:1225262/1‑59 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                              >  1:1600160/1‑59 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTc                           >  1:56076/1‑62 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                  >  1:2190831/1‑71 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                  >  1:3799987/1‑71 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                  >  1:259203/1‑71 (MQ=255)

                                                     tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                  >  1:2830485/1‑71 (MQ=255)

                                                      gTTTACCAACGATCTTCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                  <  1:2443336/70‑1 (MQ=255)

                                                         tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa                               >  1:1092611/1‑54 (MQ=255)

                                                         tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg              >  1:3129577/1‑71 (MQ=255)

                                                             aaCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                         <  1:2976428/40‑1 (MQ=255)

                                                               cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                       >  1:2747646/1‑40 (MQ=255)

                                                               cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                       >  1:1723181/1‑40 (MQ=255)

                                                                  tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc     >  1:3593713/1‑71 (MQ=255)

                                                                   catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                   >  1:708533/1‑40 (MQ=255)

                                                                    atcGGTCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                         <  1:2747710/49‑1 (MQ=255)

                                                                    atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCa           <  1:1326218/63‑1 (MQ=255)

                                                                    atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGtt      <  1:1624796/68‑1 (MQ=255)

                                                                    atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct    <  1:938850/70‑1 (MQ=255)

                                                                    atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct    <  1:3065100/70‑1 (MQ=255)

                                                                     tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                <  1:2833879/41‑1 (MQ=255)

                                                                     tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                              >  1:1333188/1‑43 (MQ=255)

                                                                     tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                              >  1:1508039/1‑43 (MQ=255)

                                                                     tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                              >  1:1983774/1‑43 (MQ=255)

                                                                     tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  <  1:3865728/71‑1 (MQ=255)

                                                                     tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  <  1:1935341/71‑1 (MQ=255)

                                                                      cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  >  1:2351553/1‑70 (MQ=255)

                                                                      cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  >  1:49787/1‑70 (MQ=255)

                                                                      cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  >  1:789443/1‑70 (MQ=255)

                                                                      cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  >  1:3344593/1‑70 (MQ=255)

                                                                      cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  >  1:3647067/1‑70 (MQ=255)

                                                                      cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  >  1:3761891/1‑70 (MQ=255)

                                                                      cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa  >  1:3877880/1‑70 (MQ=255)

TTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCA  >  minE/608915‑609054

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