BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,030	C→A	14.3%	P408P (CCG→CCT) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,030	0	C	A	14.3%	365.9 / 29.7	125	P408P (CCG→CCT) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (45/62); new base A (9/9); total (54/71)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 128 total aligned reads displayed.

ATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTC  >  Exported/608965‑609099

  ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                                                      <  1:659474/65‑1 (MQ=255)

  ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                                                      <  1:5176578/65‑1 (MQ=255)

   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                                               >  1:268613/1‑71 (MQ=255)

       tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                                            <  1:1626106/70‑1 (MQ=255)

        aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                                                           <  1:3253912/70‑1 (MQ=255)

         ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt                                                          <  1:1329495/70‑1 (MQ=255)

             cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                     >  1:822877/1‑71 (MQ=255)

             cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                     <  1:4051844/71‑1 (MQ=255)

               accattGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt                                                   <  1:2958484/69‑1 (MQ=255)

                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                         >  1:947131/1‑64 (MQ=255)

                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                         >  1:1007848/1‑64 (MQ=255)

                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                         >  1:1557215/1‑64 (MQ=255)

                tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                                  <  1:2113837/71‑1 (MQ=255)

                 catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                 >  1:1482523/1‑71 (MQ=255)

                 catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                 >  1:5542797/1‑71 (MQ=255)

                  atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                                            <  1:2042020/59‑1 (MQ=255)

                  atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                 <  1:439156/70‑1 (MQ=255)

                   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                               >  1:3290511/3‑71 (MQ=255)

                   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                               >  1:4556021/3‑71 (MQ=255)

                   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                               >  1:5684195/3‑71 (MQ=255)

                     gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                              <  1:4952891/70‑1 (MQ=255)

                     gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                                             >  1:230515/1‑71 (MQ=255)

                           aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                          >  1:4535093/1‑68 (MQ=255)

                            cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                      <  1:2738835/71‑1 (MQ=255)

                             gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                                     <  1:3895082/71‑4 (MQ=255)

                             gcagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                         <  1:373776/51‑1 (MQ=255)

                              cagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga                                    <  1:2118396/71‑1 (MQ=255)

                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                                    >  1:1982035/1‑38 (MQ=255)

                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                                    >  1:3094628/1‑38 (MQ=255)

                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGATCTCATACTGCAGagaa                                   <  1:1570117/71‑1 (MQ=255)

                                gcaggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAggaag                                  <  1:2366258/71‑1 (MQ=255)

                                gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                               <  1:3215043/58‑1 (MQ=255)

                                gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                  <  1:2007039/71‑1 (MQ=255)

                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                                    >  1:3979835/1‑35 (MQ=255)

                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                                    >  1:4328136/1‑35 (MQ=255)

                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                                    >  1:3000644/1‑35 (MQ=255)

                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                   >  1:515667/1‑68 (MQ=255)

                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                   >  1:1216229/1‑68 (MQ=255)

                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                <  1:316184/71‑1 (MQ=255)

                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                >  1:3078264/1‑71 (MQ=255)

                                   agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                <  1:3330683/70‑1 (MQ=255)

                                    gcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgt                              <  1:4835971/71‑1 (MQ=255)

                                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                              <  1:419717/71‑1 (MQ=255)

                                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                              <  1:1435471/71‑1 (MQ=255)

                                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                              >  1:4337394/1‑71 (MQ=255)

                                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                              <  1:1461097/71‑1 (MQ=255)

                                     caaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATc                                                         >  1:41501/1‑43 (MQ=12)

                                      aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga                            <  1:4618897/71‑2 (MQ=255)

                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCa                                                             <  1:2824134/37‑1 (MQ=255)

                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCa                                                        <  1:4323708/42‑1 (MQ=255)

                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                           <  1:5588921/71‑1 (MQ=255)

                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                           <  1:1551045/71‑1 (MQ=255)

                                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                           <  1:1082422/71‑1 (MQ=255)

                                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                               <  1:1385924/50‑1 (MQ=255)

                                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                      <  1:1803017/59‑1 (MQ=255)

                                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                      <  1:1820379/59‑1 (MQ=255)

                                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg                          <  1:5383566/71‑1 (MQ=255)

                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                                            >  1:3662086/1‑35 (MQ=255)

                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                       >  1:5055751/1‑40 (MQ=255)

                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                       >  1:1150438/1‑40 (MQ=255)

                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGtt                                                   >  1:938943/1‑44 (MQ=255)

                                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt                        <  1:2751295/71‑1 (MQ=255)

                                          gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa                         <  1:940828/70‑1 (MQ=255)

                                          gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGAt                        >  1:2137280/1‑71 (MQ=255)

                                             tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                                    <  1:1850512/40‑1 (MQ=255)

                                             tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGa                     <  1:3301044/71‑1 (MQ=255)

                                              aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                                            <  1:5126632/47‑1 (MQ=255)

                                              aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                                            <  1:2223316/47‑1 (MQ=255)

                                               gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                              <  1:2293534/44‑1 (MQ=255)

                                                tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                  <  1:1115463/68‑1 (MQ=255)

                                                tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                     >  1:1970741/4‑68 (MQ=255)

                                                tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                     >  1:3663707/4‑68 (MQ=255)

                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                     <  1:5290414/36‑1 (MQ=255)

                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                     <  1:4903047/36‑1 (MQ=255)

                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                 <  1:235236/40‑1 (MQ=255)

                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                 <  1:3550239/40‑1 (MQ=255)

                                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                     <  1:5667610/52‑1 (MQ=255)

                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                                             >  1:770698/1‑42 (MQ=255)

                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa                                        <  1:4137198/47‑1 (MQ=255)

                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                >  1:3379225/1‑71 (MQ=255)

                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                >  1:445539/1‑71 (MQ=255)

                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                >  1:3179307/1‑71 (MQ=255)

                                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                >  1:705335/1‑71 (MQ=255)

                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                  <  1:2456467/52‑1 (MQ=255)

                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                  <  1:1334793/52‑1 (MQ=255)

                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                  <  1:5640415/52‑1 (MQ=255)

                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                >  1:3116851/1‑70 (MQ=255)

                                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                >  1:3110437/1‑70 (MQ=255)

                                                    cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                  <  1:2280545/51‑1 (MQ=255)

                                                    cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                >  1:2781174/1‑69 (MQ=255)

                                                     gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa              <  1:845802/70‑1 (MQ=255)

                                                     gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg             >  1:5548385/1‑71 (MQ=255)

                                                       tgttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg           >  1:456444/1‑71 (MQ=255)

                                                         ttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaaga                                 <  1:5177378/47‑1 (MQ=255)

                                                          tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                     <  1:2901027/42‑1 (MQ=255)

                                                          tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                     <  1:4411910/42‑1 (MQ=255)

                                                           gAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAgg                               >  1:2996215/1‑47 (MQ=255)

                                                            aaGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                   >  1:2942945/1‑42 (MQ=255)

                                                                 cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc  <  1:5337889/70‑1 (MQ=255)

                                                                 cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc  <  1:2198931/70‑1 (MQ=255)

ATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTC  >  Exported/608965‑609099

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