BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,928	G→C	16.6%	T442T (ACC→ACG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,928	0	G	C	16.6%	242.3 / 34.9	102	T442T (ACC→ACG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (41/44); new base C (8/9); total (49/53)
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.07e-01

Only 100 of 106 total aligned reads displayed.

TCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAG  >  minE/608858‑608997

tCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                                       >  1:934303/1‑71 (MQ=255)

 cGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                                                      >  1:3949552/1‑71 (MQ=255)

  ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                                                     <  1:798967/71‑1 (MQ=255)

    ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                                                   <  1:3929044/71‑1 (MQ=255)

      cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                                                  <  1:2184315/70‑1 (MQ=255)

      cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                                                 >  1:941882/1‑71 (MQ=255)

      cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                                                 >  1:743499/1‑71 (MQ=255)

      cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                                                 >  1:479043/1‑71 (MQ=255)

        gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                                                               <  1:2220937/71‑1 (MQ=255)

          cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                                                      >  1:2574848/1‑62 (MQ=255)

            tCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                                            <  1:2889250/70‑1 (MQ=255)

             cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                                          >  1:2065075/1‑71 (MQ=255)

             cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                                          >  1:1660998/1‑71 (MQ=255)

             cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                                          >  1:2574866/1‑71 (MQ=255)

             cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                                          >  1:2579454/1‑71 (MQ=255)

              ggACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCATGtt                                                         <  1:3209922/71‑1 (MQ=255)

                 ccgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGc                                                      <  1:485992/68‑1 (MQ=255)

                   gaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgc                                                    <  1:2869311/70‑3 (MQ=255)

                    aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                                           >  1:1137796/1‑63 (MQ=255)

                    aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                                          >  1:2764278/1‑64 (MQ=255)

                    aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                   >  1:1475953/1‑71 (MQ=255)

                    aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAg                                                            <  1:2463701/62‑1 (MQ=255)

                    aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                   <  1:1521300/71‑1 (MQ=255)

                     ccTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgc                                                  <  1:2996212/71‑1 (MQ=255)

                     ccTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                                          <  1:3059050/63‑1 (MQ=255)

                      cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                                 <  1:3665939/71‑1 (MQ=255)

                         cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                               <  1:2835929/70‑1 (MQ=255)

                         cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                                              >  1:1781102/1‑71 (MQ=255)

                           ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                                                  <  1:833394/49‑1 (MQ=255)

                           ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                            <  1:921355/71‑1 (MQ=255)

                            cttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGt                                                          <  1:3800219/55‑1 (MQ=255)

                             ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATgg                                                                      <  1:3159143/43‑1 (MQ=255)

                             ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                                              >  1:3646875/1‑51 (MQ=255)

                             ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                                              >  1:3053695/1‑51 (MQ=255)

                             ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                                              >  1:1858069/1‑51 (MQ=255)

                             ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCga                                           <  1:3563169/70‑1 (MQ=255)

                             ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                            >  1:3240407/1‑69 (MQ=255)

                              tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                                                   >  1:569910/1‑45 (MQ=255)

                               cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                         <  1:675343/70‑1 (MQ=255)

                               cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                         >  1:2215457/1‑70 (MQ=255)

                                aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                                                         <  1:3144932/53‑1 (MQ=255)

                                aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                       >  1:3825802/1‑71 (MQ=255)

                                aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                       >  1:922123/1‑71 (MQ=255)

                                 gTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGttt                                                        <  1:2342304/53‑1 (MQ=255)

                                  tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                       >  1:2451095/1‑53 (MQ=255)

                                  tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTg                                                       >  1:3416380/1‑53 (MQ=255)

                                  tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                    >  1:570867/1‑56 (MQ=255)

                                   ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                    >  1:3185564/1‑71 (MQ=255)

                                    tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCa                                                             >  1:1163429/1‑45 (MQ=255)

                                    tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                    >  1:3422021/1‑70 (MQ=255)

                                     ttCGTCAACGCCGCCGTTGATTGCGTACATCATGGTTTTCGcc                                                              <  1:1492905/43‑1 (MQ=255)

                                     ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCga                                           >  1:1259392/1‑62 (MQ=255)

                                     ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                  >  1:3152497/1‑71 (MQ=255)

                                     ttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                  >  1:3789051/1‑71 (MQ=255)

                                      tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                 >  1:311948/1‑71 (MQ=255)

                                       cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                 <  1:1852187/70‑1 (MQ=255)

                                        gTCCACGCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                <  1:1468722/70‑1 (MQ=255)

                                        gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                <  1:2778101/70‑1 (MQ=255)

                                         tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                               <  1:2557423/54‑1 (MQ=255)

                                          cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAc                                                    <  1:1149445/48‑1 (MQ=255)

                                          cAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                             >  1:206006/1‑71 (MQ=255)

                                           aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                            >  1:3542483/1‑71 (MQ=255)

                                           aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                            >  1:2708089/1‑71 (MQ=255)

                                           aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                            <  1:2912248/71‑1 (MQ=255)

                                           aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                             <  1:2048483/70‑1 (MQ=255)

                                           aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                             <  1:2327653/70‑1 (MQ=255)

                                           aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt                             <  1:3233016/70‑1 (MQ=255)

                                           aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGG‑TTGCACGCGCAcc                                             >  1:489912/1‑53 (MQ=255)

                                            aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGttt                           <  1:2576364/71‑1 (MQ=255)

                                             cgccgcCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaaga                                        <  1:2662450/57‑1 (MQ=255)

                                              cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGAGCGTGCACCAAAGAACTGCATATgc                             >  1:682230/2‑66 (MQ=255)

                                               ccgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                                          <  1:1694967/37‑1 (MQ=255)

                                                cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgt                                                  >  1:2889951/1‑43 (MQ=16)

                                                cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                        >  1:1643088/1‑70 (MQ=255)

                                                 gccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCaa                      >  1:2542308/1‑71 (MQ=255)

                                                  ccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                         <  1:3901684/51‑1 (MQ=255)

                                                    gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                 <  1:2242866/57‑1 (MQ=255)

                                                    gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGa                   >  1:3502839/1‑71 (MQ=255)

                                                    gTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGa                   >  1:2423550/1‑71 (MQ=255)

                                                    gTTGATTGCGTACACCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                 <  1:2564288/57‑1 (MQ=255)

                                                       gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                  >  1:1525207/1‑69 (MQ=255)

                                                        aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                            <  1:823223/42‑1 (MQ=255)

                                                        aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                            <  1:2077600/42‑1 (MQ=255)

                                                        aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat               <  1:3354044/71‑1 (MQ=255)

                                                        aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                <  1:3315875/70‑1 (MQ=255)

                                                        aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                <  1:2226859/70‑1 (MQ=255)

                                                        aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                <  1:3074143/70‑1 (MQ=255)

                                                           gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                            <  1:2344491/39‑1 (MQ=255)

                                                           gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCg             <  1:3931339/70‑1 (MQ=255)

                                                            cGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                >  1:148422/1‑50 (MQ=255)

                                                            cGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                >  1:3137387/1‑50 (MQ=255)

                                                            cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCggg           >  1:1808979/1‑71 (MQ=255)

                                                             gTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                      <  1:1013124/43‑1 (MQ=255)

                                                             gTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGc          >  1:1335388/1‑71 (MQ=255)

                                                              tACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt         <  1:1903957/71‑1 (MQ=255)

                                                                cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa       <  1:3093794/71‑1 (MQ=255)

                                                                  gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                            <  1:2632525/48‑1 (MQ=255)

                                                                  gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCggg           <  1:1882474/65‑1 (MQ=255)

                                                                  gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg            >  1:3197307/1‑64 (MQ=255)

                                                                     tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag  >  1:2136562/1‑71 (MQ=255)

TCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAG  >  minE/608858‑608997

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