BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,186	T→G	15.4%	M690L (ATG→CTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,186	0	T	G	15.4%	197.7 / 34.4	91	M690L (ATG→CTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (49/28); new base G (7/7); total (56/35)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 112 total aligned reads displayed.

CGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTG  >  Exported/608116‑608256

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

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

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

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

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

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

  attCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                                                                      <  1:4609794/70‑1 (MQ=12)

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

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

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

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

      attacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGcc                                                                  <  1:4386120/68‑1 (MQ=9)

       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)

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

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

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

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

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

          cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccaga                                                               >  1:3046544/1‑69 (MQ=255)

            ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagccag                                                            >  1:5154923/1‑71 (MQ=255)

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

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

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

              aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATcagcccagcc                                                              >  1:1906388/1‑67 (MQ=255)

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

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

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

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

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

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

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

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

                 gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg                                                       <  1:4275958/71‑1 (MQ=255)

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

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

                      ggTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                  >  1:2654879/1‑71 (MQ=255)

                       gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                  >  1:2324261/1‑70 (MQ=255)

                       gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                 <  1:1212158/71‑1 (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)

                          cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGaa                                              <  1:12246/71‑1 (MQ=255)

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

                           tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                                                   >  1:1664053/1‑49 (MQ=255)

                           tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                  <  1:4780395/66‑1 (MQ=255)

                           tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATAAGAc                                                                   >  1:3797227/1‑49 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                     cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                   >  1:330922/1‑71 (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:756104/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

                                            tgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAcc                              >  1:64198/1‑69 (MQ=255)

                                             gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                    >  1:774782/1‑62 (MQ=255)

                                             gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                    >  1:1995157/1‑62 (MQ=255)

                                              cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc                                       >  1:251112/1‑58 (MQ=255)

                                               ttCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCttt                                 <  1:725121/63‑1 (MQ=255)

                                                tCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                         >  1:4196184/1‑68 (MQ=255)

                                                 cGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                                <  1:2983188/46‑1 (MQ=255)

                                                 cGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg                                                <  1:460381/46‑1 (MQ=255)

                                                     gtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                   >  1:5161950/1‑71 (MQ=255)

                                                      tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                  <  1:2945434/71‑1 (MQ=255)

                                                           gTAACCATCCATCAGACCAGCCAGGTTGGTCTTTCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGa             <  1:4596744/71‑1 (MQ=255)

                                                              aCCATCTATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa           <  1:5236490/70‑1 (MQ=255)

                                                              aCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGa             >  1:1746875/1‑68 (MQ=255)

                                                               ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATaga         <  1:304087/71‑1 (MQ=255)

                                                               ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATaga         <  1:2638319/71‑1 (MQ=255)

                                                                catccatcAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                      <  1:3559112/41‑1 (MQ=255)

                                                                catccatcAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                      <  1:208932/41‑1 (MQ=255)

                                                                catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag        <  1:5529521/71‑1 (MQ=255)

                                                                catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag        <  1:837849/71‑1 (MQ=255)

                                                                 atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAgaga       <  1:3290948/71‑1 (MQ=255)

                                                                  tccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                   >  1:1827767/1‑58 (MQ=255)

                                                                   ccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                  >  1:3014024/1‑42 (MQ=255)

                                                                    catcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCa                            >  1:2290821/1‑47 (MQ=255)

                                                                      tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg  <  1:151381/71‑1 (MQ=255)

CGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTG  >  Exported/608116‑608256

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