BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,727	G→A	34.2%	A509A (GCC→GCT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,727	0	G	A	34.2%	166.7 / 130.6	152	A509A (GCC→GCT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (53/47); new base A (19/33); total (72/80)
Fisher's exact test for biased strand distribution p-value = 6.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 167 total aligned reads displayed.

ACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAG  >  minE/608659‑608795

aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGc                                                                    >  1:225641/1‑71 (MQ=255)

  ggACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                                                  <  1:203083/71‑1 (MQ=255)

          cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                                                  >  1:4797998/1‑63 (MQ=255)

          cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                                                  >  1:789345/1‑63 (MQ=255)

          cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                                                  >  1:4804471/1‑63 (MQ=255)

          cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                                          <  1:4899071/71‑1 (MQ=255)

            gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                                        >  1:1757715/1‑71 (MQ=255)

              gATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTAc                                                       >  1:4218048/1‑70 (MQ=255)

              gATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACt                                                      >  1:4981778/1‑71 (MQ=255)

               aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTc                                                                 >  1:2933757/1‑59 (MQ=255)

                  ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                                          <  1:4109552/63‑1 (MQ=255)

                    acacGCCATGGTGCTGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGt                                                <  1:5588762/71‑1 (MQ=255)

                     catgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg                                               <  1:3102444/68‑1 (MQ=255)

                     catgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg                                               <  1:3557553/68‑1 (MQ=255)

                     catgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg                                               <  1:4083297/68‑1 (MQ=255)

                     cacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg                                               >  1:3314966/1‑71 (MQ=255)

                       tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa                                             >  1:568748/2‑71 (MQ=255)

                       tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa                                             >  1:3838314/2‑71 (MQ=255)

                         ccATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATg                                           <  1:4580196/71‑1 (MQ=255)

                          cATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                          <  1:4603465/71‑1 (MQ=255)

                           aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg                                               <  1:5161873/65‑1 (MQ=255)

                                     tA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                                                >  1:3158117/1‑38 (MQ=255)

                                     tA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                                                >  1:1718547/1‑38 (MQ=255)

                                     tA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                               >  1:2485028/1‑71 (MQ=255)

                                      a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGTCTTc                                                                 >  1:1685677/1‑36 (MQ=255)

                                      a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                              >  1:3993825/1‑71 (MQ=255)

                                        acatCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                             <  1:4782496/68‑1 (MQ=255)

                                        aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                             <  1:392572/71‑1 (MQ=255)

                                         catCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                             >  1:3770710/3‑70 (MQ=255)

                                         cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                             >  1:3245952/1‑70 (MQ=255)

                                          atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata                               >  1:1286047/2‑66 (MQ=255)

                                          atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata                               >  1:3689863/2‑66 (MQ=255)

                                          atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTca                           <  1:3606150/70‑1 (MQ=255)

                                          atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTca                           <  1:4062699/70‑1 (MQ=255)

                                          atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTca                           <  1:2054779/70‑1 (MQ=255)

                                           tCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTa                                         <  1:4340230/56‑1 (MQ=255)

                                           tCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                             <  1:2958806/68‑1 (MQ=255)

                                            cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                          >  1:4652470/1‑54 (MQ=255)

                                            cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                    >  1:5182136/1‑60 (MQ=255)

                                            cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                    >  1:3317384/1‑60 (MQ=255)

                                            cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagt                         >  1:3647712/1‑71 (MQ=255)

                                            cACGATCGTGCAGCTCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTg                                      <  1:4873007/58‑1 (MQ=255)

                                            cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgat                                >  1:5504193/1‑64 (MQ=255)

                                             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                        <  1:2905936/71‑1 (MQ=255)

                                             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                        >  1:2192726/1‑71 (MQ=255)

                                             aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                        >  1:3443817/1‑69 (MQ=255)

                                              cgatCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcag                          >  1:1028952/4‑68 (MQ=255)

                                                gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag                     >  1:2970207/1‑71 (MQ=255)

                                                 tCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                    <  1:3292929/71‑4 (MQ=255)

                                                 tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCAt                                            <  1:2779826/47‑1 (MQ=255)

                                                 tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCAt                                            <  1:963286/47‑1 (MQ=255)

                                                 tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                    <  1:3191561/71‑1 (MQ=255)

                                                  cGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                    <  1:989415/70‑4 (MQ=255)

                                                   gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                   <  1:1358362/70‑1 (MQ=255)

                                                   gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                  >  1:1885736/1‑71 (MQ=255)

                                                     gCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                >  1:5679924/1‑71 (MQ=255)

                                                     gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                >  1:5244114/1‑71 (MQ=255)

                                                     gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                >  1:755507/1‑71 (MQ=255)

                                                     gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                >  1:2277605/1‑71 (MQ=255)

                                                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                <  1:3837388/70‑1 (MQ=255)

                                                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                <  1:4888514/70‑1 (MQ=255)

                                                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                <  1:1965890/70‑1 (MQ=255)

                                                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGt               <  1:470009/71‑1 (MQ=255)

                                                      cAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                <  1:2351756/70‑1 (MQ=255)

                                                      cAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                <  1:1620676/70‑1 (MQ=255)

                                                        gcgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                             <  1:1676327/55‑1 (MQ=255)

                                                         cgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt            <  1:1975400/71‑1 (MQ=255)

                                                         cgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt            <  1:2579433/71‑1 (MQ=255)

                                                         cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                 >  1:714355/1‑50 (MQ=255)

                                                         cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                 >  1:5423308/1‑50 (MQ=255)

                                                         cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                 >  1:1561869/1‑50 (MQ=255)

                                                         cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                 >  1:5415543/1‑50 (MQ=255)

                                                          gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg           <  1:3739247/71‑1 (MQ=255)

                                                          gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg           <  1:2673486/71‑1 (MQ=255)

                                                          gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg           <  1:3488505/71‑1 (MQ=255)

                                                            cATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgca         <  1:1111821/71‑3 (MQ=255)

                                                            cATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgca         <  1:2729436/71‑3 (MQ=255)

                                                            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                              <  1:317435/50‑1 (MQ=255)

                                                            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg           >  1:4474950/1‑69 (MQ=255)

                                                            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg           >  1:2101143/1‑69 (MQ=255)

                                                            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt         <  1:609636/71‑1 (MQ=255)

                                                             aTCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                    <  1:1812992/59‑1 (MQ=255)

                                                              tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcag                          >  1:3070682/1‑52 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                          <  1:20060/36‑1 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                          <  1:902813/36‑1 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                          <  1:4623851/36‑1 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatg                                  <  1:5555646/44‑1 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                      <  1:2762348/56‑1 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                      <  1:5468068/56‑1 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                      <  1:1877917/56‑1 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                      <  1:217065/56‑1 (MQ=255)

                                                               cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTag      >  1:1534630/1‑71 (MQ=255)

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                           >  1:2807347/1‑49 (MQ=255)

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc     >  1:4505278/1‑71 (MQ=255)

                                                                 gagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTg                                      <  1:2529322/37‑1 (MQ=255)

                                                                   gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                    >  1:2416177/1‑37 (MQ=255)

                                                                   gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                    >  1:1248271/1‑37 (MQ=255)

                                                                    agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                        <  1:3539066/48‑1 (MQ=255)

                                                                    agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                        <  1:928533/48‑1 (MQ=255)

                                                                     ggCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                              <  1:4176974/41‑1 (MQ=255)

ACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAG  >  minE/608659‑608795

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