BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,205	C→T	5.1%	K683K (AAG→AAA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,205	0	C	T	5.1%	328.1 / 5.6	98	K683K (AAG→AAA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (38/55); new base T (2/3); total (40/58)
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.99e-01

Only 100 of 102 total aligned reads displayed.

AGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTT  >  Exported/608137‑608272

aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                    <  1:2370936/70‑1 (MQ=255)

aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                    <  1:1380264/70‑1 (MQ=255)

aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTGGTCtg                                                                   <  1:828901/71‑2 (MQ=39)

    gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCTATCAGACCAGCCAGGTTGGTCTTACg                                                                >  1:3258761/1‑70 (MQ=255)

     cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGa                                                               <  1:434264/70‑1 (MQ=255)

      tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAAc                                                             >  1:1377271/1‑71 (MQ=255)

        accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACtt                                                           <  1:3205139/71‑1 (MQ=255)

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

         ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                          >  1:831423/1‑71 (MQ=255)

         ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                          <  1:599943/71‑1 (MQ=255)

         ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                          <  1:1914503/71‑1 (MQ=255)

         ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                          <  1:1910382/71‑1 (MQ=255)

          caccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcg                                                         <  1:2452448/71‑1 (MQ=255)

             cTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                                      >  1:1584842/1‑71 (MQ=255)

              ttCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCtt                                                                   <  1:769596/57‑1 (MQ=255)

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

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

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

                 gatggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                                  <  1:2778301/71‑1 (MQ=255)

                  atggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                          <  1:547807/62‑1 (MQ=255)

                   tggatgCTTCGTGGTGTAAGTAACCATCCATCAGACCAGCCAGCTTGGTCTTACGa                                                               >  1:1817731/1‑56 (MQ=39)

                    ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                               >  1:978465/1‑71 (MQ=255)

                    ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                               >  1:1084660/1‑71 (MQ=255)

                    ggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                               >  1:1056648/1‑71 (MQ=255)

                     gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTAGTCTTTAcc                                              <  1:1333530/71‑1 (MQ=255)

                      atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                    <  1:2632856/48‑1 (MQ=255)

                      atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                    <  1:2407937/48‑1 (MQ=255)

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

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

                        gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                                  >  1:151684/1‑48 (MQ=255)

                        gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                                  >  1:3438648/1‑48 (MQ=255)

                        gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                                  >  1:1416806/1‑48 (MQ=255)

                        gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                                  >  1:1191422/1‑48 (MQ=255)

                        gCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                                  >  1:621257/1‑48 (MQ=255)

                         cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                                          <  1:2006578/71‑1 (MQ=255)

                          ttCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                                         >  1:3317690/1‑69 (MQ=255)

                          ttCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgcg                                         >  1:555346/1‑69 (MQ=255)

                              tggtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCa                                            >  1:1810574/1‑64 (MQ=255)

                              tggtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCa                                            >  1:2123364/1‑64 (MQ=255)

                              ctgtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                                      <  1:1079716/68‑1 (MQ=255)

                                 tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                    <  1:242005/37‑1 (MQ=255)

                                 tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                                      <  1:1129711/67‑1 (MQ=255)

                                 tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                                      <  1:1732766/67‑1 (MQ=255)

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

                                  ggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCtt                                                                   <  1:651479/37‑1 (MQ=255)

                                   gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                    >  1:2901128/1‑35 (MQ=255)

                                   gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                    >  1:271989/1‑35 (MQ=255)

                                   gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAAc                                                             >  1:3705305/1‑42 (MQ=255)

                                   gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                          <  1:2564473/45‑1 (MQ=255)

                                   gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                  <  1:81576/69‑1 (MQ=255)

                                     aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACGCAGTGCGTTCGGAAc                               >  1:3145403/1‑70 (MQ=255)

                                     aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                  <  1:3243649/67‑1 (MQ=255)

                                     aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                  <  1:2212706/67‑1 (MQ=255)

                                         aCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                                    <  1:2868811/45‑1 (MQ=255)

                                           catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                               <  1:2684017/64‑1 (MQ=255)

                                            atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                        >  1:2396999/1‑70 (MQ=255)

                                            atccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                        >  1:2197010/1‑70 (MQ=255)

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

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

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

                                              ccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAg                     >  1:2451130/1‑71 (MQ=255)

                                               catcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAgg                    >  1:3194138/1‑71 (MQ=255)

                                                  cAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg                  >  1:2079794/1‑70 (MQ=255)

                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt                 <  1:1355717/70‑1 (MQ=255)

                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt                 <  1:2219934/70‑1 (MQ=255)

                                                   aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt                 <  1:541844/70‑1 (MQ=255)

                                                    gACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg               <  1:189502/71‑1 (MQ=255)

                                                      ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                                   >  1:2560959/1‑49 (MQ=255)

                                                      ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                                   >  1:179306/1‑49 (MQ=255)

                                                      ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCg                                   >  1:2847790/1‑49 (MQ=255)

                                                      ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg               >  1:470641/1‑69 (MQ=255)

                                                       cagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGGACGATAGAGAAGGTGTAg               >  1:341976/1‑68 (MQ=255)

                                                        agccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAg               >  1:3269415/1‑67 (MQ=255)

                                                        agccagGTTGGTCTTACGAACTTCGTCGTCTATACCCAGTGCGTTCGGAACGATAGAGAAGGTGTa                >  1:2142812/1‑66 (MQ=255)

                                                         gccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGgaga           >  1:2956301/1‑70 (MQ=255)

                                                          ccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                        <  1:2029097/56‑1 (MQ=255)

                                                          ccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATa         <  1:3022282/71‑1 (MQ=255)

                                                           cagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAc        >  1:2464540/1‑71 (MQ=255)

                                                            agGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGgaga           <  1:1493229/67‑1 (MQ=255)

                                                            agGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc       <  1:1577235/71‑1 (MQ=255)

                                                            agGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATAcc       <  1:601451/71‑1 (MQ=255)

                                                             ggTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCa      <  1:2778570/71‑1 (MQ=255)

                                                             ggTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCa      >  1:1274257/1‑71 (MQ=255)

                                                              gTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCAt     <  1:2998289/71‑1 (MQ=255)

                                                              gTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCAt     <  1:632466/71‑1 (MQ=255)

                                                               ttGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATTGAGAAGGTGTACGAGATCCCATc    <  1:2991555/71‑1 (MQ=255)

                                                               ttGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATTGAGAAGGTGTACGAGATCCCATc    <  1:672893/71‑1 (MQ=255)

                                                               ttGGTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATTGAGAAGGTGTACGAGATCCCATc    <  1:3163479/71‑1 (MQ=255)

                                                               ttGGTCTTACGAACTTGGTCGTCTTTACCCAGTGCTTTCGGAACGATAGAGAACGTGTAGGAGATACCa      <  1:3379116/69‑1 (MQ=255)

                                                               ttGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATc    >  1:3214630/1‑71 (MQ=255)

                                                               ttGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATc    >  1:1002129/1‑71 (MQ=255)

                                                                tGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCt   >  1:3175971/1‑71 (MQ=255)

                                                                 ggTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCt   <  1:3048909/70‑1 (MQ=255)

                                                                  gTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt  <  1:3229215/70‑1 (MQ=255)

                                                                   tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATaga                         <  1:225981/46‑1 (MQ=255)

                                                                   tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt  <  1:2306891/69‑1 (MQ=255)

                                                                   tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCt   <  1:343722/68‑1 (MQ=255)

                                                                    cTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATag                          <  1:2567041/44‑1 (MQ=255)

                                                                    cTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATag                          <  1:3274462/44‑1 (MQ=255)

                                                                    cTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATag                          <  1:1673309/44‑1 (MQ=255)

AGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTT  >  Exported/608137‑608272

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