BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,160	T→C	22.3%	A698A (GCA→GCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,160	0	T	C	22.3%	244.8 / 62.3	121	A698A (GCA→GCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (60/34); new base C (16/11); total (76/45)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 136 total aligned reads displayed.

CCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCC  >  minE/608093‑608229

ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                                                    >  1:1309634/1‑71 (MQ=255)

ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                                                    >  1:2694238/1‑71 (MQ=255)

     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                                                               >  1:5125089/1‑71 (MQ=255)

     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgt                                                                <  1:4659684/70‑1 (MQ=255)

       gTCGAGAAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                                             <  1:2142944/71‑1 (MQ=255)

         cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                                                            >  1:1158104/1‑70 (MQ=255)

         cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                                                            >  1:229330/1‑70 (MQ=255)

          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:1019380/70‑1 (MQ=255)

          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2140849/70‑1 (MQ=255)

          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2369404/70‑1 (MQ=255)

          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:2385207/70‑1 (MQ=255)

          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:3282078/70‑1 (MQ=255)

          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:4259562/70‑1 (MQ=255)

          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                                                           >  1:5280240/2‑70 (MQ=255)

          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                                                           <  1:4463507/69‑1 (MQ=255)

          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                                                          <  1:2301901/70‑1 (MQ=12)

            gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2948280/1‑70 (MQ=12)

            gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2978923/1‑70 (MQ=12)

             cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                                                                 <  1:3371798/61‑1 (MQ=255)

              agcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATa                                                      >  1:4622590/1‑71 (MQ=12)

               gcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATac                                                     >  1:5666752/1‑70 (MQ=12)

                caTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                                                    <  1:3871174/71‑3 (MQ=9)

                caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                                                    >  1:1334170/1‑71 (MQ=255)

                 aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                                                    >  1:1072629/1‑68 (MQ=9)

                  tttCACGGTTCATCACGTTAACTTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                  >  1:4915494/1‑71 (MQ=255)

                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                                   >  1:4274676/1‑70 (MQ=255)

                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                                   >  1:4281342/1‑70 (MQ=255)

                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAaca                                                  >  1:3366388/1‑71 (MQ=255)

                   ttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                                        >  1:2831214/1‑63 (MQ=255)

                    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc                                                >  1:701406/1‑71 (MQ=255)

                    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat                                                 <  1:700995/70‑1 (MQ=255)

                      aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc                                               <  1:2767935/70‑1 (MQ=255)

                       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:3110125/71‑1 (MQ=255)

                       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:3082681/71‑1 (MQ=255)

                       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:4829242/71‑1 (MQ=255)

                        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:3961175/1‑70 (MQ=255)

                        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:5401046/1‑70 (MQ=255)

                        ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             >  1:5186034/1‑70 (MQ=255)

                            cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc                                        >  1:1950331/1‑69 (MQ=255)

                            cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc                                        >  1:5406626/1‑69 (MQ=255)

                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                        >  1:5157747/1‑71 (MQ=255)

                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                        >  1:4958019/1‑71 (MQ=255)

                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           <  1:5129726/64‑1 (MQ=9)

                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           <  1:2189330/64‑1 (MQ=9)

                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                                      <  1:5048742/69‑2 (MQ=255)

                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                                      <  1:3477946/69‑2 (MQ=255)

                              tcaCGTTAACGTTCAGGTCCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAg                                          <  1:244548/67‑1 (MQ=255)

                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                                                         >  1:3474597/1‑50 (MQ=255)

                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                              >  1:3737946/1‑61 (MQ=255)

                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                              >  1:1697490/1‑61 (MQ=255)

                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                                              >  1:1947077/1‑61 (MQ=255)

                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                                     <  1:1231279/70‑1 (MQ=255)

                                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                                           <  1:4196643/47‑1 (MQ=255)

                                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                                    >  1:1497237/1‑70 (MQ=255)

                                     gaCGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           >  1:847201/2‑59 (MQ=9)

                                     gaCGTTGAGGTGTTGACAGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                           >  1:2651620/2‑59 (MQ=255)

                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAcc                                       >  1:1661064/1‑63 (MQ=255)

                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAgg                                <  1:2450950/70‑1 (MQ=255)

                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:3673025/1‑71 (MQ=255)

                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:3487062/1‑71 (MQ=255)

                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:2440528/1‑71 (MQ=255)

                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt                               >  1:4374389/1‑71 (MQ=255)

                                       cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGATCAGCCAGGTTg                             >  1:1063503/1‑71 (MQ=255)

                                        gTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTgg                            >  1:1924003/1‑71 (MQ=255)

                                          tgaGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtt                          >  1:3347066/3‑70 (MQ=255)

                                           cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:4177110/1‑71 (MQ=255)

                                              gtgtTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatccag                                             <  1:5696337/48‑2 (MQ=255)

                                              gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                      <  1:1212158/71‑1 (MQ=255)

                                              gTGATGACCACCTTCGATGGATGCTTCGTGGTGGAAGt                                                       >  1:4922195/1‑38 (MQ=255)

                                               tGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                                  >  1:4328242/1‑42 (MQ=255)

                                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                      <  1:4439662/69‑1 (MQ=255)

                                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:2530585/1‑70 (MQ=255)

                                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:5183963/1‑70 (MQ=255)

                                                gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:1493020/1‑70 (MQ=255)

                                                  tGACCGCCTTCGACATACGCTTCGTGGTGGAAATACACATCCAGCAGGCCGACAAGGTTGGTTTTACGTAc                  >  1:1881386/1‑71 (MQ=255)

                                                   gACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     <  1:4104449/67‑1 (MQ=255)

                                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:935576/1‑66 (MQ=255)

                                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                     >  1:3098408/1‑66 (MQ=255)

                                                    accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                 <  1:1572301/70‑1 (MQ=255)

                                                    accaccTACGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACtt                >  1:4563970/1‑71 (MQ=255)

                                                        ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGa                    >  1:1552807/1‑63 (MQ=255)

                                                        ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCAAGGTTGGTCTTACGa                    >  1:3175673/1‑63 (MQ=255)

                                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:5352547/1‑50 (MQ=255)

                                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:2331445/1‑50 (MQ=255)

                                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:5594752/1‑50 (MQ=255)

                                                           tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                              >  1:4045519/1‑50 (MQ=255)

                                                            cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt          <  1:4803447/69‑1 (MQ=255)

                                                            cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt        >  1:4551955/1‑71 (MQ=255)

                                                              atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt      <  1:1041323/71‑1 (MQ=255)

                                                              atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt       >  1:1830649/1‑70 (MQ=255)

                                                               tggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt      <  1:5088166/70‑1 (MQ=255)

                                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGt                           <  1:571495/47‑1 (MQ=255)

                                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc   <  1:3715514/71‑1 (MQ=255)

                                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc   >  1:4544237/1‑71 (MQ=255)

                                                                 gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc   >  1:756104/1‑71 (MQ=255)

                                                                  atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc  <  1:3626318/71‑1 (MQ=255)

                                                                  atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc  <  1:260989/71‑1 (MQ=255)

                                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:3246150/1‑47 (MQ=255)

                                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:3721636/1‑47 (MQ=255)

                                                                   tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         >  1:4405462/1‑47 (MQ=255)

CCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCC  >  minE/608093‑608229

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