BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,913	G→A	14.7%	I447I (ATC→ATT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,913	0	G	A	14.7%	259.9 / 30.1	102	I447I (ATC→ATT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (50/37); new base A (6/9); total (56/46)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 104 total aligned reads displayed.

CATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGAT  >  minE/608846‑608981

cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                                   >  1:266883/1‑71 (MQ=255)

cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                                   <  1:3893848/71‑1 (MQ=255)

cATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                                   <  1:3944075/71‑1 (MQ=255)

    gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                                >  1:3283948/1‑70 (MQ=255)

    gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                                >  1:3545287/1‑70 (MQ=255)

    gCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                                >  1:685629/1‑70 (MQ=255)

      cTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTAc                                                             <  1:2993798/71‑1 (MQ=255)

         tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                           >  1:3720852/1‑70 (MQ=255)

         tGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                           >  1:401538/1‑70 (MQ=255)

          gATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGc                                                          <  1:3535708/70‑1 (MQ=255)

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

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

              ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGt                                                     >  1:2842195/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:743499/1‑71 (MQ=255)

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

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

                    gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                                    <  1:2513308/50‑1 (MQ=255)

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

                      cTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                                   >  1:3858894/1‑49 (MQ=255)

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

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

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

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

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

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

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

                          ggACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                          <  1:1811441/70‑1 (MQ=255)

                           gACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAtt                                                                   >  1:3177881/1‑44 (MQ=255)

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

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

                               aaCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                                      <  1:3071286/37‑1 (MQ=255)

                                aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                                    >  1:2564407/1‑38 (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)

                                 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)

                                         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)

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

                                           cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                                            >  1:2858225/1‑35 (MQ=255)

                                           cAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACa                                                            >  1:659043/1‑35 (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:2912248/71‑1 (MQ=255)

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

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

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

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

                                                       aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGt             <  1:2327653/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)

                                                              ccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc     <  1:3946145/71‑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)

CATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGAT  >  minE/608846‑608981

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