BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,637	T→G	9.2%	A539A (GCA→GCC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,637	0	T	G	9.2%	266.6 / 19.3	108	A539A (GCA→GCC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (47/51); new base G (2/8); total (49/59)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGT  >  Exported/608567‑608707

cGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGAtt                                                                        >  1:4668290/1‑71 (MQ=255)

 gCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGAtt                                                                        <  1:4925840/70‑1 (MQ=255)

    ttCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTg                                                                       >  1:2852850/1‑68 (MQ=255)

     tCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGa                                                                   >  1:3383365/1‑71 (MQ=255)

     tCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGa                                                                   >  1:2561910/1‑71 (MQ=255)

     tCGATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCTGGCATATTTGATGGCAGa                                                                   <  1:2195776/71‑1 (MQ=255)

      cGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGc                                                                      >  1:2313723/1‑67 (MQ=255)

      cGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGa                                                                   <  1:457391/70‑1 (MQ=255)

      cGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGa                                                                   <  1:5508069/70‑1 (MQ=255)

      cGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGa                                                                   >  1:3989754/1‑70 (MQ=255)

      cGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGa                                                                   <  1:3957706/70‑1 (MQ=255)

       gATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACa                                                                 <  1:1533572/71‑1 (MQ=255)

       gATTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACa                                                                 <  1:3924634/71‑1 (MQ=255)

        aTTTCAAAGTCCACCGCCAGGCCGTTTTCGTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAg                                                                <  1:3500287/71‑1 (MQ=255)

         tttCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAgg                                                               <  1:1425660/71‑1 (MQ=255)

          ttCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAgg                                                               >  1:3310988/1‑70 (MQ=255)

          ttCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAgg                                                               >  1:3805755/1‑70 (MQ=255)

            cGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGa                                                             <  1:3324675/70‑1 (MQ=255)

            cGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGa                                                             <  1:3331302/70‑1 (MQ=255)

             gAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGt                                                           >  1:4827992/1‑71 (MQ=255)

              aaGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTc                                                          <  1:3012058/71‑1 (MQ=255)

              aaGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTc                                                          <  1:3694775/71‑1 (MQ=255)

                gTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcag                                                        >  1:3606078/1‑71 (MQ=255)

                   gATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGa                                                             >  1:453125/1‑63 (MQ=255)

                   gATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagcag                                                     >  1:2594567/1‑71 (MQ=255)

                   gATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagca                                                      <  1:1947354/70‑1 (MQ=255)

                    aTAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCagcagc                                                    <  1:4490458/71‑1 (MQ=255)

                     tAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagcag                                                     >  1:2958687/1‑69 (MQ=255)

                     tAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCagcagc                                                    >  1:3295102/1‑70 (MQ=255)

                        ccAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACg                                                <  1:827431/71‑1 (MQ=255)

                         cAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCaa                                                  >  1:4241703/1‑68 (MQ=255)

                          aGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                              >  1:4964313/1‑71 (MQ=255)

                           gACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                              >  1:2606141/1‑70 (MQ=255)

                              cGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagca                                                      >  1:798317/1‑59 (MQ=255)

                              cGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                            >  1:5132053/1‑69 (MQ=255)

                                tCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTg                                                                       <  1:5400373/40‑1 (MQ=255)

                                tCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTg                                                                       <  1:4662849/40‑1 (MQ=255)

                                tCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTg                                                                       <  1:1049976/40‑1 (MQ=255)

                                tCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTg                                                                       <  1:1520065/40‑1 (MQ=255)

                                tCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTg                                                                       <  1:1200173/40‑1 (MQ=255)

                                tCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCagcagc                                                    <  1:2299178/59‑1 (MQ=255)

                                  ttCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAg                                      >  1:2221473/1‑71 (MQ=255)

                                  ttCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAg                                      >  1:5559560/1‑71 (MQ=255)

                                  ttCGTCACGAATCGGTTTAACTTTAGCATATTTGATTGCAGa                                                                   <  1:4900488/42‑1 (MQ=255)

                                    cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                     <  1:525254/70‑1 (MQ=255)

                                    cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                     <  1:829012/70‑1 (MQ=255)

                                    cGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                     <  1:5559538/70‑1 (MQ=255)

                                     gTCACGGATTGGTTTCACGCGGGCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGc                                     <  1:5654942/69‑1 (MQ=255)

                                      tCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGAt                                  <  1:4089282/71‑1 (MQ=255)

                                        aCGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagc                                                       >  1:1254233/1‑48 (MQ=255)

                                        aCGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATaa                                >  1:4233397/1‑70 (MQ=255)

                                          gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACg                                                <  1:914416/53‑1 (MQ=255)

                                          gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACca                              <  1:463864/71‑1 (MQ=255)

                                          gAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACca                              <  1:2480866/71‑1 (MQ=255)

                                          gAATCGATTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGAc                                             <  1:4183711/56‑1 (MQ=255)

                                            aTCGGTTTAACTTTCGCATATTTTATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACcaca                            <  1:3384186/71‑1 (MQ=255)

                                            aTCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACcaca                            <  1:2769197/71‑1 (MQ=255)

                                             tCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCacac                           >  1:5277265/1‑71 (MQ=255)

                                             tCGGTTTAACTTTCCCATATTTGATTGCAGACAGGGa                                                             <  1:1368351/37‑1 (MQ=255)

                                              cGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGt                                                           <  1:4138769/38‑1 (MQ=255)

                                               ggTTTAACTTTCGCATATTTGATTGCAGACAGGGAGt                                                           >  1:2592090/1‑37 (MQ=255)

                                               ggTTTAACTTTCGCATATTTGATTGCAGACAGGGAGt                                                           >  1:2722996/1‑37 (MQ=255)

                                                gTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGcc                        <  1:1521978/71‑1 (MQ=255)

                                                gTTAAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGc                         <  1:3387994/70‑1 (MQ=255)

                                                  ttAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCAt                      >  1:5109093/1‑71 (MQ=255)

                                                   tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                     >  1:3085328/1‑71 (MQ=255)

                                                   tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                     >  1:115620/1‑71 (MQ=255)

                                                   tAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                     >  1:5075723/1‑71 (MQ=255)

                                                    aaCTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATg                     >  1:2842431/1‑70 (MQ=255)

                                                      cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                       <  1:4423978/50‑1 (MQ=255)

                                                      cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                       <  1:4223321/50‑1 (MQ=255)

                                                      cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCa                                       <  1:758341/50‑1 (MQ=255)

                                                      cTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTg                  <  1:117589/71‑1 (MQ=255)

                                                       tttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                 >  1:2713393/1‑71 (MQ=255)

                                                        ttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTg                  >  1:5359604/1‑69 (MQ=255)

                                                        ttCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCg                >  1:4107728/1‑71 (MQ=255)

                                                         tCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                 <  1:2585078/69‑1 (MQ=255)

                                                         gggCATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGATGCCGCATGCCATAGTGCga               <  1:373464/69‑2 (MQ=255)

                                                          cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                    >  1:1824049/1‑65 (MQ=255)

                                                          cGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATgg                    >  1:1881742/1‑65 (MQ=255)

                                                           gCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGAt             <  1:937516/71‑1 (MQ=255)

                                                            cttatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa            <  1:2431491/69‑1 (MQ=255)

                                                            cATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATa            <  1:920176/71‑1 (MQ=255)

                                                            cATATTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGCGa                                   >  1:1089330/1‑48 (MQ=255)

                                                              tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                            <  1:4210037/37‑1 (MQ=255)

                                                              tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACa                                            <  1:3466546/37‑1 (MQ=255)

                                                              tatTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGAt                                  <  1:1127564/47‑1 (MQ=255)

                                                              tatTTGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGcccg                                      <  1:1613128/43‑3 (MQ=12)

                                                               atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                     <  1:574690/43‑1 (MQ=255)

                                                               atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACg         >  1:1149481/1‑71 (MQ=255)

                                                                tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAg                                           >  1:3417255/1‑36 (MQ=255)

                                                                tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAg                                           >  1:1767876/1‑36 (MQ=255)

                                                                tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGa                                   <  1:218253/44‑1 (MQ=255)

                                                                tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt        >  1:4624953/1‑71 (MQ=255)

                                                                tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt        >  1:4120457/1‑71 (MQ=255)

                                                                tttGATGGCAGACAGGGAGTCCGTCGCCACCGACAGGCCCGc                                     >  1:1975401/1‑42 (MQ=255)

                                                                 ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt        >  1:3549713/1‑70 (MQ=255)

                                                                 ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt        >  1:4482926/1‑70 (MQ=255)

                                                                     ttGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACgg   <  1:709857/71‑1 (MQ=255)

                                                                      tGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGt  >  1:3881570/1‑71 (MQ=255)

CGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGT  >  Exported/608567‑608707

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