BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,118	G→A	37.5%	N712N (AAC→AAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,118	0	G	A	37.5%	155.5 / 122.9	128	N712N (AAC→AAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (45/35); new base A (31/17); total (76/52)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 131 total aligned reads displayed.

AGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCAT  >  minE/608052‑608186

agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtc                                                                   <  1:3476824/70‑1 (MQ=9)

aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                                                  >  1:4707475/1‑71 (MQ=255)

 gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                                 <  1:3047823/68‑1 (MQ=9)

 gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                                 <  1:1208217/68‑1 (MQ=9)

 gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                                 <  1:1890052/68‑1 (MQ=9)

   tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                                <  1:3986411/70‑1 (MQ=255)

    acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                              >  1:378263/1‑71 (MQ=255)

     cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                                                             <  1:3994210/71‑1 (MQ=12)

     cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                                                             <  1:2878476/71‑1 (MQ=12)

     cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                                                             <  1:2610403/71‑1 (MQ=12)

     cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                                                             <  1:5515929/71‑1 (MQ=12)

     cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACg                                                             >  1:1900190/1‑71 (MQ=255)

       cGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                           <  1:724340/71‑1 (MQ=255)

        ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                           >  1:2256425/1‑70 (MQ=255)

         gATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa                                                         >  1:3215964/1‑71 (MQ=255)

             gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa                                                         <  1:388756/67‑1 (MQ=255)

               cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                                >  1:3268764/1‑58 (MQ=255)

               cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                     >  1:672858/1‑69 (MQ=255)

                aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa                                                         >  1:5346504/1‑64 (MQ=255)

                aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                  >  1:4238926/1‑71 (MQ=255)

                   ttaggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTg                                               <  1:430491/68‑1 (MQ=255)

                     cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                               <  1:3561424/69‑1 (MQ=255)

                     cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                             >  1:4815554/1‑71 (MQ=255)

                      ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                   <  1:684597/64‑1 (MQ=255)

                      ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTg                                            <  1:4759502/71‑1 (MQ=255)

                      ggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgtt                                             <  1:1529062/70‑3 (MQ=255)

                       gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                                           >  1:352683/1‑71 (MQ=12)

                       gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                                           >  1:5338729/1‑71 (MQ=12)

                        atatTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTg                                               <  1:5669114/66‑1 (MQ=255)

                         tatTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGAc                                                        >  1:28388/1‑56 (MQ=255)

                         tatTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGacc                                         >  1:3784295/1‑71 (MQ=12)

                          atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                             >  1:223045/1‑66 (MQ=255)

                          atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                        >  1:762954/1‑71 (MQ=255)

                          atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccg                                        >  1:4047357/1‑70 (MQ=255)

                          atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccg                                        >  1:2994233/1‑70 (MQ=255)

                          atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccg                                        >  1:3157662/1‑70 (MQ=255)

                          atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGacc                                         >  1:4518422/1‑70 (MQ=255)

                           tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGtt                                                                    >  1:4697980/1‑42 (MQ=255)

                           tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                   >  1:1148539/1‑59 (MQ=255)

                           tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                                       >  1:4866991/1‑71 (MQ=255)

                           tttttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGacc                                         <  1:1020302/69‑1 (MQ=255)

                            ttttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                      >  1:2690368/1‑71 (MQ=255)

                             tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTg                                            >  1:1028818/1‑64 (MQ=255)

                             tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                      <  1:41402/70‑1 (MQ=255)

                             tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                     >  1:4776960/1‑71 (MQ=255)

                             tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                     >  1:4862636/1‑71 (MQ=255)

                             tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                     >  1:5486463/1‑71 (MQ=255)

                             tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                     >  1:3596185/1‑71 (MQ=255)

                               tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTc                                   >  1:894824/1‑71 (MQ=12)

                               tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTc                                   >  1:4702231/1‑71 (MQ=12)

                               tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTc                                   >  1:2933824/1‑71 (MQ=12)

                                ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                           <  1:3789738/46‑1 (MQ=255)

                                ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                   <  1:1375712/70‑1 (MQ=255)

                                   ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                           >  1:4418540/1‑59 (MQ=255)

                                       ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                 >  1:1784602/1‑65 (MQ=255)

                                        tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat                                >  1:1564142/1‑65 (MQ=255)

                                         ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                         >  1:2694238/1‑71 (MQ=255)

                                         ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                         >  1:1309634/1‑71 (MQ=255)

                                         ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                         >  1:1071379/1‑71 (MQ=255)

                                            tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                    <  1:574831/41‑1 (MQ=255)

                                              gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                    >  1:5125089/1‑71 (MQ=255)

                                              gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgt                     <  1:4659684/70‑1 (MQ=255)

                                                gTCGAGAAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                  <  1:2142944/71‑1 (MQ=255)

                                                 tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                <  1:4177272/40‑1 (MQ=255)

                                                 tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                   <  1:5622409/53‑1 (MQ=255)

                                                  cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                 >  1:229330/1‑70 (MQ=255)

                                                  cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                 >  1:1158104/1‑70 (MQ=255)

                                                  cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg                 >  1:2851578/1‑70 (MQ=255)

                                                   gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                <  1:1019380/70‑1 (MQ=255)

                                                   gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                <  1:4259562/70‑1 (MQ=255)

                                                   gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                <  1:2140849/70‑1 (MQ=255)

                                                   gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                <  1:3282078/70‑1 (MQ=255)

                                                   gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                <  1:2369404/70‑1 (MQ=255)

                                                   gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                <  1:2385207/70‑1 (MQ=255)

                                                   cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                     >  1:5378634/2‑49 (MQ=12)

                                                   cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCga                                 >  1:449414/2‑53 (MQ=12)

                                                   cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                <  1:4463507/69‑1 (MQ=255)

                                                   cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                >  1:5280240/2‑70 (MQ=255)

                                                   cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa               <  1:2301901/70‑1 (MQ=12)

                                                     gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa             >  1:2978923/1‑70 (MQ=12)

                                                     gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa             >  1:2948280/1‑70 (MQ=12)

                                                     gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                   <  1:5141500/49‑1 (MQ=255)

                                                     gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                   <  1:2723896/49‑1 (MQ=255)

                                                     gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg             >  1:4247507/1‑71 (MQ=255)

                                                      cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                      <  1:1564587/61‑1 (MQ=255)

                                                      cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                      <  1:3371798/61‑1 (MQ=255)

                                                        gcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATac          >  1:5666752/1‑70 (MQ=12)

                                                         caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc         >  1:1334170/1‑71 (MQ=255)

                                                         caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc         <  1:3379660/71‑1 (MQ=255)

                                                           tttCACGGTTCATCACGTTAACTTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca       >  1:4915494/1‑71 (MQ=255)

                                                           tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca       >  1:2504342/1‑71 (MQ=255)

                                                           tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc        >  1:4281342/1‑70 (MQ=255)

                                                            ttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa             >  1:2831214/1‑63 (MQ=255)

                                                             tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                <  1:2085760/60‑1 (MQ=255)

                                                             tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat      <  1:700995/70‑1 (MQ=255)

                                                               aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc    <  1:2767935/70‑1 (MQ=255)

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

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

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

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

AGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCAT  >  minE/608052‑608186

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