BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,724	A→C	34.2%	S510S (TCT→TCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,724	0	A	C	34.2%	104.0 / 159.2	152	S510S (TCT→TCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (53/47); new base C (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.93e-01

Only 100 of 169 total aligned reads displayed.

ACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGC  >  minE/608659‑608792

 ccgACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCt                                                                >  1:2536424/3‑71 (MQ=255)

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

          ccgCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTg                                                       >  1:2590246/3‑71 (MQ=255)

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

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

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

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

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

                  ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagag                                                                     >  1:1411553/1‑49 (MQ=255)

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

                  ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                                                <  1:4063157/70‑1 (MQ=255)

                  ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTc                                               <  1:4524482/71‑1 (MQ=255)

                   cacaCGCCATCGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCaga                                                                      >  1:4626461/1‑47 (MQ=255)

                    acacGCCATGGTGCGGATG‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGt                                             <  1:5425087/71‑1 (MQ=255)

                    acacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagag                                                                     >  1:5594075/1‑47 (MQ=255)

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

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

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

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

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

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

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

                              gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTg                                   <  1:1283842/71‑1 (MQ=255)

                                  ggATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAg                                     >  1:697820/1‑65 (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:5158606/1‑71 (MQ=255)

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

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

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

                                          atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                          <  1:1976578/68‑1 (MQ=255)

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

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

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

                                           tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                          <  1:5486016/68‑1 (MQ=255)

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

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

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

                                            cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTc                         <  1:378611/68‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                     gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                              >  1:5197420/1‑54 (MQ=255)

                                                     gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                              >  1:4880499/1‑54 (MQ=255)

                                                     gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                              >  1:600218/1‑54 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                          gcCATCAGAGAGTCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg        <  1:3744199/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                             aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgcac     <  1:4290162/71‑4 (MQ=255)

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

                                                              tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatga                              <  1:4764929/45‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAg                                     >  1:4959340/1‑36 (MQ=255)

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

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt      >  1:2912605/1‑67 (MQ=255)

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

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

ACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGC  >  minE/608659‑608792

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