BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,136	C→G	33.8%	L706L (CTG→CTC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,136	0	C	G	33.8%	189.0 / 121.4	136	L706L (CTG→CTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (55/35); new base G (36/10); total (91/45)
Fisher's exact test for biased strand distribution p-value = 5.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 144 total aligned reads displayed.

AGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCT  >  Exported/608068‑608206

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

   ttaggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTg                                                                   <  1:430491/68‑1 (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)

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

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

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

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

          atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGacc                                                             >  1:4518422/1‑70 (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)

             tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                                                  >  1:4919348/1‑61 (MQ=255)

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

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

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

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

             tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                                         >  1:4862636/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)

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

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

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

                     ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                                 >  1:4004942/1‑71 (MQ=255)

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

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

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

                         ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                             >  1:2694238/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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                     gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg                                 >  1:5649424/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)

                                         caTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                             <  1:3871174/71‑3 (MQ=9)

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

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

                                          aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                             >  1:1072629/1‑68 (MQ=9)

                                           tttCACGGTTCATCACGTTAACGTTCAGGTGCTGac                                                              >  1:3132951/1‑36 (MQ=255)

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

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

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

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

                                             tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc                         >  1:701406/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                     catcaCGTTAACGTTCAGGTGCTGACCACCTTCga                                                     >  1:2808846/1‑35 (MQ=255)

                                                     catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                             <  1:1334738/59‑1 (MQ=255)

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

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

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

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

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

                                                         aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                       >  1:1947077/1‑61 (MQ=255)

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

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

                                                          cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                    <  1:4196643/47‑1 (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)

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

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

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

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

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

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

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

AGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCT  >  Exported/608068‑608206

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