BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,955	G→A	24.0%	F433F (TTC→TTT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,955	0	G	A	24.0%	268.7 / 80.2	146	F433F (TTC→TTT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (42/69); new base A (20/15); total (62/84)
Fisher's exact test for biased strand distribution p-value = 5.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 168 total aligned reads displayed.

TTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGA  >  minE/608887‑609025

  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                                     <  1:5081114/70‑1 (MQ=255)

  cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                                     <  1:3217052/70‑1 (MQ=255)

   aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                                   <  1:2485728/71‑1 (MQ=255)

     tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACt                                                                 <  1:3820837/71‑1 (MQ=255)

      ttctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                                                                >  1:5586913/4‑71 (MQ=255)

       tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                                >  1:2813492/1‑70 (MQ=255)

       tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                                >  1:3031028/1‑70 (MQ=255)

        ctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                               <  1:5301396/69‑1 (MQ=255)

         tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                             >  1:2057723/1‑71 (MQ=255)

         tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                             >  1:1163239/1‑71 (MQ=255)

         tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                                   >  1:3078734/1‑65 (MQ=255)

            tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg                                                          <  1:1694310/71‑1 (MQ=255)

            tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg                                                          <  1:5268998/71‑1 (MQ=255)

              aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                        <  1:898099/71‑1 (MQ=255)

              aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                                        >  1:2402722/1‑71 (MQ=255)

               aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACTCGCACCGAAGAACTGCATTTGtt                                                        <  1:2118390/70‑1 (MQ=255)

                ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                                              <  1:1155144/61‑1 (MQ=255)

                 cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                                     >  1:3756992/2‑71 (MQ=255)

                   cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                                     <  1:4665827/53‑1 (MQ=255)

                   cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                                     <  1:1526179/53‑1 (MQ=255)

                   cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCa                                                   <  1:2065457/71‑1 (MQ=255)

                   cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                                                        >  1:1081659/1‑64 (MQ=255)

                   cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                                                        >  1:3426410/1‑64 (MQ=255)

                   cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                    >  1:4221099/1‑70 (MQ=255)

                   cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                                    >  1:5186119/1‑70 (MQ=255)

                     ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                                      >  1:5516307/1‑66 (MQ=255)

                     ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaat                                                 <  1:2375314/71‑2 (MQ=255)

                      cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg                                                <  1:2771617/71‑1 (MQ=255)

                        ttGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                              >  1:4718237/1‑71 (MQ=255)

                         tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaa                                                                      >  1:3720903/1‑46 (MQ=255)

                          gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                             <  1:4320014/70‑1 (MQ=255)

                          aaTTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                                             >  1:4177436/2‑54 (MQ=255)

                           aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                            <  1:261160/70‑1 (MQ=255)

                            ttGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                                     <  1:3908913/44‑1 (MQ=12)

                               cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                        <  1:493001/70‑1 (MQ=255)

                                 tACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                     <  1:4018357/71‑1 (MQ=255)

                                   cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                                  >  1:5458941/1‑40 (MQ=255)

                                   cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                                  >  1:3498844/1‑40 (MQ=255)

                                   cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                   >  1:5572364/1‑71 (MQ=255)

                                    aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                                  <  1:196137/71‑1 (MQ=255)

                                    aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                                  <  1:1128710/71‑1 (MQ=255)

                                      cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                   <  1:5202882/68‑1 (MQ=255)

                                      cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                 >  1:2583007/1‑70 (MQ=255)

                                       aTGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                                   >  1:3132667/1‑67 (MQ=255)

                                        tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag                              >  1:814231/1‑71 (MQ=255)

                                          gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                             >  1:1569553/1‑70 (MQ=255)

                                          gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc                             <  1:1461504/70‑1 (MQ=255)

                                           ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                            <  1:576990/54‑1 (MQ=255)

                                           ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                            <  1:3517010/54‑1 (MQ=255)

                                           ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                            <  1:399597/54‑1 (MQ=255)

                                           ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                                    <  1:1943784/62‑1 (MQ=255)

                                           ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa                           <  1:5109958/71‑1 (MQ=255)

                                           ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa                           <  1:4812571/71‑1 (MQ=255)

                                            tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                          >  1:3055222/1‑71 (MQ=255)

                                            tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                          >  1:5117703/1‑71 (MQ=255)

                                             ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                          >  1:2663018/1‑70 (MQ=255)

                                             ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                          >  1:2391923/1‑70 (MQ=255)

                                             ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                          >  1:5332640/1‑68 (MQ=255)

                                             ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg                          >  1:32789/1‑68 (MQ=255)

                                              tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                     <  1:5587630/58‑1 (MQ=255)

                                              tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                         <  1:5547664/70‑1 (MQ=255)

                                              tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                         <  1:2032018/70‑1 (MQ=255)

                                              tCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggc                         <  1:4334380/70‑1 (MQ=255)

                                                gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                       <  1:4224293/70‑1 (MQ=255)

                                                gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                       <  1:4653418/70‑1 (MQ=255)

                                                  cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaatc                                                <  1:1912273/43‑3 (MQ=19)

                                                  cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaatc                                                <  1:1905524/43‑3 (MQ=19)

                                                   aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                   >  1:4947648/1‑71 (MQ=255)

                                                   aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc                   <  1:4923163/71‑1 (MQ=255)

                                                    ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                             >  1:1341037/1‑44 (MQ=255)

                                                    ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                         >  1:4642752/1‑64 (MQ=255)

                                                    ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                  <  1:224252/71‑1 (MQ=255)

                                                     gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                     <  1:4772410/51‑1 (MQ=255)

                                                     gTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                 <  1:2925223/71‑1 (MQ=255)

                                                      tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                                 <  1:1089070/38‑1 (MQ=255)

                                                         gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                             <  1:5612882/39‑1 (MQ=255)

                                                         gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                             <  1:2935523/39‑1 (MQ=255)

                                                         gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                        <  1:4919366/44‑1 (MQ=255)

                                                         gCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                        <  1:2511788/44‑1 (MQ=255)

                                                           aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                           >  1:4914727/1‑39 (MQ=255)

                                                           aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                           >  1:1657053/1‑39 (MQ=255)

                                                           aCGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                            <  1:252889/38‑1 (MQ=255)

                                                             gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg         <  1:3806412/71‑1 (MQ=255)

                                                              cgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgt        <  1:1980355/71‑1 (MQ=255)

                                                               gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt       >  1:4528867/1‑71 (MQ=255)

                                                               gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt       >  1:3664825/1‑71 (MQ=255)

                                                               gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgt        >  1:140475/1‑70 (MQ=255)

                                                               gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAgcaggcaa                       <  1:4796127/55‑1 (MQ=255)

                                                                cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                   <  1:4355125/58‑1 (MQ=255)

                                                                cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg      <  1:2154030/71‑1 (MQ=255)

                                                                cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc                     <  1:4774095/56‑2 (MQ=255)

                                                                cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg      >  1:3185884/1‑71 (MQ=255)

                                                                 aCCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                             <  1:5172893/47‑1 (MQ=16)

                                                                  ccGAAGAACTGCATTTGTTTACCAACGATCATCGGCCTTACgcag                              <  1:771054/45‑1 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgca                               <  1:3427250/43‑1 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg   <  1:2209975/71‑1 (MQ=255)

                                                                   cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg   <  1:4459885/71‑1 (MQ=255)

                                                                    gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                            >  1:2836770/1‑45 (MQ=255)

                                                                    gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                            >  1:4493100/1‑45 (MQ=255)

                                                                    gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa  <  1:217997/71‑1 (MQ=255)

TTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGA  >  minE/608887‑609025

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