BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,205	C→T	6.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
*	minE	608,205	0	C	T	6.1%	310.3 / 6.5	99	K683K (AAG→AAA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (43/49); new base T (2/4); total (46/53)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 104 total aligned reads displayed.

TCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTA  >  minE/608135‑608274

tCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc                                                                       >  1:1994647/1‑71 (MQ=255)

 cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc                                                                       >  1:1792913/1‑70 (MQ=255)

 cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc                                                                       <  1:3159160/70‑1 (MQ=255)

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

 cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGT‑GAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                      >  1:2402537/1‑70 (MQ=255)

  aGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtttt                                                                     <  1:2231230/71‑1 (MQ=255)

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

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

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

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

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

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

  aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTGGTCt                                                                      <  1:3119832/70‑1 (MQ=255)

   ggTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTa                                                                    <  1:3409193/71‑1 (MQ=255)

   ggTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTa                                                                    <  1:3824925/71‑1 (MQ=255)

   ggTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTa                                                                    <  1:3880655/71‑1 (MQ=255)

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

    gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                                   >  1:923689/1‑71 (MQ=255)

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

          accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAAc                                                               <  1:3631156/69‑1 (MQ=255)

          accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                              >  1:1538295/1‑70 (MQ=255)

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

           ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                                   >  1:2649502/1‑64 (MQ=255)

           ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                                   >  1:1479166/1‑64 (MQ=255)

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

            caccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcg                                                           >  1:2952013/1‑71 (MQ=255)

              ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc                                                         <  1:2284422/71‑1 (MQ=255)

              ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc                                                         <  1:924198/71‑1 (MQ=255)

                 tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                                      >  1:3731260/1‑71 (MQ=255)

                  cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAGGAACTTcgtcgt                                                       <  1:2293151/69‑1 (MQ=255)

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

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

                  cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                                      <  1:3724881/70‑1 (MQ=255)

                  cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                                      <  1:1797071/70‑1 (MQ=255)

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

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

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

                   gatggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                                    >  1:3933899/1‑71 (MQ=255)

                     tggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                                      <  1:2676588/51‑1 (MQ=255)

                       gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                                                   >  1:2492549/1‑52 (MQ=255)

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

                           cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                                    <  1:151775/47‑1 (MQ=255)

                           cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                                    <  1:334105/47‑1 (MQ=255)

                           cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                                       <  1:463558/60‑1 (MQ=255)

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

                           cTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgc                                            <  1:504485/71‑2 (MQ=255)

                             tCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACtgg                                                            >  1:2273165/1‑51 (MQ=255)

                             tCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACtgg                                                            >  1:3358305/1‑51 (MQ=255)

                              cGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGaa                                                                >  1:2059806/1‑48 (MQ=255)

                               gtggtgGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                                        >  1:1672782/1‑71 (MQ=255)

                                ttgtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                              <  1:824844/46‑1 (MQ=255)

                                  gtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTc                                      <  1:546230/70‑1 (MQ=255)

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

                                    ggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                                    <  1:2843820/54‑1 (MQ=255)

                                    ggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                                    >  1:677342/1‑70 (MQ=255)

                                     gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                                               >  1:637500/1‑58 (MQ=255)

                                      aaGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                                 >  1:571959/1‑71 (MQ=255)

                                       aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                                 >  1:2924643/1‑70 (MQ=255)

                                         tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                            >  1:3342509/1‑41 (MQ=255)

                                         tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                                            >  1:3751140/1‑41 (MQ=255)

                                         tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                                       >  1:622982/1‑46 (MQ=255)

                                         tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                                 >  1:1036054/1‑52 (MQ=255)

                                          aaCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                                      >  1:856909/1‑46 (MQ=255)

                                          aaCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                                     >  1:2703262/1‑47 (MQ=255)

                                            ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                                         <  1:2832041/57‑1 (MQ=255)

                                            ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                                         <  1:1975745/57‑1 (MQ=255)

                                            ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa                             >  1:1611263/1‑69 (MQ=255)

                                            ccatccatCAGACCACCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                                         <  1:1263972/57‑1 (MQ=255)

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

                                               tccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag                          >  1:203586/1‑69 (MQ=255)

                                                 catcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa                             >  1:3905764/1‑64 (MQ=255)

                                                  atcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTg                                           <  1:126335/49‑1 (MQ=255)

                                                   tcAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                                        <  1:2829526/35‑1 (MQ=255)

                                                   tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAGCGATAGAGAAGgtg                    >  1:1607805/1‑71 (MQ=255)

                                                   tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa                             <  1:1007353/62‑1 (MQ=255)

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

                                                    cAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTc                                      >  1:2140048/1‑52 (MQ=255)

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

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

                                                     aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg                    >  1:151327/1‑69 (MQ=255)

                                                       aCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAgg                <  1:2461678/71‑1 (MQ=255)

                                                       aCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAgg                <  1:435437/71‑1 (MQ=255)

                                                        ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGtt                                       >  1:1256042/1‑47 (MQ=255)

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

                                                           gccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGAt            <  1:1203182/71‑1 (MQ=255)

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

                                                               ggTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGgaga             >  1:331444/1‑66 (MQ=255)

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

                                                                gTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCa        >  1:3724385/1‑70 (MQ=255)

                                                                gTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCAt       >  1:2121433/1‑71 (MQ=255)

                                                                gTTGGTATTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt                   >  1:689505/1‑59 (MQ=255)

                                                                  tGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGaa                        >  1:2476657/1‑52 (MQ=255)

                                                                   ggTTTTACGTACTGGATCTTCTTTGCCCAGCGCCGCAGGAACGATTGAGAAGGTGTacg                <  1:2164387/59‑3 (MQ=255)

                                                                   ggTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCtt    >  1:2248047/1‑71 (MQ=255)

                                                                    gTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCttt   <  1:3630105/71‑1 (MQ=255)

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

                                                                     tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgt                     >  1:2908105/1‑52 (MQ=255)

                                                                     tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTa                  <  1:183149/55‑1 (MQ=255)

                                                                     tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTa  <  1:3187426/71‑1 (MQ=255)

                                                                     tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTa  <  1:3325699/71‑1 (MQ=255)

TCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTA  >  minE/608135‑608274

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