BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,136	0	C	G	13.5%	290.9 / 30.1	111	L706L (CTG→CTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (49/47); new base G (2/13); total (51/60)
Fisher's exact test for biased strand distribution p-value = 1.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 118 total aligned reads displayed.

TCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCT  >  minE/608066‑608206

tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                                        <  1:928150/71‑1 (MQ=255)

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

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

  aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                                      <  1:3241491/71‑1 (MQ=255)

   gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                                       <  1:1912623/69‑1 (MQ=255)

    cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                                    <  1:3583132/71‑1 (MQ=255)

    cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                                    >  1:40621/1‑71 (MQ=255)

    cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                                                    <  1:2538351/71‑1 (MQ=255)

      taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                                                  >  1:1934025/3‑70 (MQ=255)

      gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                                                                  >  1:3690568/1‑71 (MQ=255)

         gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                               <  1:1929460/71‑1 (MQ=255)

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

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

             tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                                                           <  1:2149654/71‑1 (MQ=255)

              ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                                          >  1:3251582/1‑71 (MQ=255)

              ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                                          >  1:534188/1‑71 (MQ=255)

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

                   cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                                     <  1:1278827/71‑1 (MQ=255)

                   cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                                     <  1:868636/71‑1 (MQ=255)

                      gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                                   >  1:1553947/1‑70 (MQ=255)

                       ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                                      <  1:237535/50‑1 (MQ=255)

                       ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                                      <  1:1953432/50‑1 (MQ=255)

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

                          tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                                              >  1:3689004/1‑71 (MQ=255)

                           ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTGAGGTGCTGACCACCTTCGATGGATGc                                              <  1:3419300/70‑1 (MQ=255)

                           ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                                              <  1:1213134/70‑1 (MQ=255)

                            cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAgg                                                                     <  1:904983/46‑1 (MQ=255)

                            cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCtt                                            >  1:3171065/1‑71 (MQ=255)

                             aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                                        >  1:1568133/1‑58 (MQ=255)

                             aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                                        >  1:1577473/1‑58 (MQ=255)

                             aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                                           >  1:350886/1‑71 (MQ=255)

                              tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                                          >  1:136118/1‑71 (MQ=255)

                               cgcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                                      >  1:1356558/1‑42 (MQ=255)

                                gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                                   >  1:1872647/1‑44 (MQ=255)

                                gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                                   >  1:2517574/1‑44 (MQ=255)

                                gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                                                   >  1:3904503/1‑44 (MQ=255)

                                  gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                      >  1:2158204/1‑71 (MQ=255)

                                     gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                                                 >  1:3094586/1‑41 (MQ=255)

                                     gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                                           <  1:37952/63‑1 (MQ=255)

                                     gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                    >  1:1551415/1‑70 (MQ=255)

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

                                     gAGCAGCAATTCACGGTTCATCACGTTAACGTTCAg                                                                      >  1:1227772/1‑36 (MQ=255)

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

                                     cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                                   <  1:1099635/70‑1 (MQ=255)

                                     cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                                   <  1:1307701/70‑1 (MQ=255)

                                      agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                                          <  1:1502184/47‑1 (MQ=255)

                                      agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaa                                  <  1:2845961/71‑1 (MQ=12)

                                      agcagcATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                               >  1:2861022/1‑42 (MQ=255)

                                       gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                 <  1:2209198/71‑2 (MQ=255)

                                          gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                              >  1:690240/1‑71 (MQ=255)

                                          gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                              >  1:686513/1‑71 (MQ=255)

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

                                              ttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat                          >  1:1077683/1‑71 (MQ=255)

                                               tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                                                <  1:2052940/48‑1 (MQ=255)

                                                cACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                      >  1:923323/1‑57 (MQ=255)

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

                                                  cGGTTCATCACGTTAACGTTCAGGTGCTTACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                      <  1:3467065/71‑1 (MQ=255)

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

                                                    gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                      >  1:2550594/1‑69 (MQ=255)

                                                    gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                    <  1:743175/71‑1 (MQ=255)

                                                    gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                    <  1:1093349/71‑1 (MQ=255)

                                                    gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                    <  1:3007658/71‑1 (MQ=255)

                                                     ttCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                    <  1:400357/70‑1 (MQ=255)

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

                                                       catcaCGTTAACGATCAGGTGCTGACCACCTACgatgga                                                 >  1:2430576/1‑39 (MQ=37)

                                                        atcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                 <  1:1054824/70‑1 (MQ=255)

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

                                                         tcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAcca               <  1:3408674/71‑1 (MQ=255)

                                                          tacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaa                                  >  1:556602/2‑51 (MQ=255)

                                                          caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg                                 >  1:3902437/1‑52 (MQ=255)

                                                           aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc             >  1:2987300/1‑71 (MQ=255)

                                                           aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc             >  1:1065758/1‑71 (MQ=255)

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

                                                           aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag              >  1:462972/1‑70 (MQ=255)

                                                            cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                                  >  1:2506238/1‑49 (MQ=255)

                                                             gTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatccag                      <  1:2263703/60‑2 (MQ=255)

                                                             gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca           <  1:866714/71‑1 (MQ=255)

                                                             gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca           <  1:10364/71‑1 (MQ=255)

                                                             gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca           <  1:1102502/71‑1 (MQ=255)

                                                              ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGt                                <  1:2996603/49‑1 (MQ=255)

                                                              ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                              >  1:1356578/1‑51 (MQ=255)

                                                               tAACGTTCAGGTGCTGACCACCTTCGTTGGATGCtt                                            <  1:2838288/36‑1 (MQ=255)

                                                               tAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                             >  1:448236/1‑51 (MQ=255)

                                                               tAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                             >  1:3099265/1‑51 (MQ=255)

                                                                aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc                        >  1:1893065/1‑55 (MQ=255)

                                                                aaCGTTCAGGGGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTa                               <  1:3781202/48‑1 (MQ=255)

                                                                 aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt       >  1:1799271/1‑71 (MQ=255)

                                                                 aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt       <  1:3703610/71‑1 (MQ=255)

                                                                  cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg      <  1:3223304/71‑1 (MQ=255)

                                                                  cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg      <  1:2821012/71‑1 (MQ=255)

                                                                  cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg      <  1:668188/71‑1 (MQ=255)

                                                                   gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt       >  1:3239501/1‑69 (MQ=255)

                                                                    ttCAGGTGCTGACCACCTTCGATGGATGCTTCgtgg                                       <  1:2320480/36‑1 (MQ=255)

                                                                    ttCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                      <  1:1175710/53‑1 (MQ=255)

                                                                    ttCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                 <  1:3874000/58‑1 (MQ=255)

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

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

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

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

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

TCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCT  >  minE/608066‑608206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`