BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,991	T→G	7.4%	V421V (GTA→GTC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,991	0	T	G	7.4%	310.4 / 21.9	121	V421V (GTA→GTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (66/46); new base G (5/4); total (71/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 122 total aligned reads displayed.

GCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACT  >  minE/608924‑609058

gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                                                  <  1:2648335/71‑1 (MQ=255)

gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                                                  <  1:43335/71‑1 (MQ=255)

gCATGGTTTTCGCCAGGTTTGCACCCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                                                                  >  1:3004311/1‑71 (MQ=255)

   tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCGTACg                                                                  >  1:3760010/1‑68 (MQ=255)

   gcgTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                                               >  1:1398399/3‑71 (MQ=255)

     gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagca                                                             <  1:3498025/71‑1 (MQ=255)

      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgca                                                                <  1:2918905/67‑1 (MQ=255)

          cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                                         <  1:11636/70‑1 (MQ=255)

          cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                                         <  1:1670833/70‑1 (MQ=255)

          cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggca                                                         <  1:1900468/70‑1 (MQ=255)

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

               ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                                     >  1:4078426/1‑69 (MQ=255)

               ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                                     >  1:2729084/1‑69 (MQ=255)

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

               ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                                   >  1:1221440/1‑71 (MQ=255)

                gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGt                                                  <  1:1555201/71‑1 (MQ=255)

                 tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTa                                                 <  1:232203/71‑1 (MQ=255)

                  ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                        >  1:1117979/1‑63 (MQ=255)

                  ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                                      >  1:2022989/1‑65 (MQ=255)

                  ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                                                >  1:593652/1‑71 (MQ=255)

                  ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                                                >  1:40416/1‑71 (MQ=255)

                   tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                                                >  1:488209/1‑70 (MQ=255)

                   tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                               >  1:509491/1‑71 (MQ=255)

                   tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                               >  1:3324975/1‑71 (MQ=255)

                   tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                               >  1:2525102/1‑71 (MQ=255)

                   tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                               >  1:2004154/1‑71 (MQ=255)

                     cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                                                                     <  1:1916882/47‑1 (MQ=255)

                      aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                          <  1:3363694/57‑1 (MQ=255)

                      aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCAt                                            <  1:3381657/71‑1 (MQ=255)

                        gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                        >  1:282008/1‑57 (MQ=255)

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

                         tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa                                                 <  1:456201/62‑1 (MQ=255)

                         tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATAGCGTa                                                 <  1:802870/62‑1 (MQ=255)

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

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

                          gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                                        >  1:1385780/1‑69 (MQ=255)

                          gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                                        >  1:808507/1‑69 (MQ=255)

                             ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                          <  1:2944900/50‑1 (MQ=255)

                             ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                      >  1:123820/1‑70 (MQ=255)

                              cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAAt                                                       <  1:1006285/52‑1 (MQ=255)

                              cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                     >  1:2056076/1‑70 (MQ=255)

                              cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                     >  1:2061036/1‑70 (MQ=255)

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

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

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

                               gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                     >  1:1840749/1‑69 (MQ=255)

                               gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                       >  1:528020/1‑67 (MQ=255)

                               gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                    >  1:3036046/1‑70 (MQ=255)

                                aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                    <  1:618985/69‑1 (MQ=255)

                                aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                    >  1:2744738/1‑69 (MQ=255)

                                aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGaa                                  <  1:1272267/71‑1 (MQ=255)

                                 agaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAg                                 >  1:3157389/1‑71 (MQ=255)

                                  gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                                >  1:1475473/1‑71 (MQ=255)

                                    aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTcc                              <  1:1934522/71‑1 (MQ=255)

                                     cTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                          <  1:4002379/42‑1 (MQ=255)

                                     cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                        >  1:1850420/1‑44 (MQ=255)

                                     cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                             <  1:3254358/71‑1 (MQ=255)

                                     cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                             <  1:421580/71‑1 (MQ=255)

                                       gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                      >  1:1060756/1‑60 (MQ=255)

                                       gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                      >  1:765660/1‑60 (MQ=255)

                                        cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                        <  1:3988942/41‑1 (MQ=255)

                                        cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                        <  1:244663/41‑1 (MQ=255)

                                        cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                    >  1:838121/1‑45 (MQ=255)

                                        cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGAc                          >  1:2272008/1‑71 (MQ=255)

                                          tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAg                                 <  1:643970/62‑1 (MQ=255)

                                            tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCTTAGTCATCGTTGTTGaa                                  >  1:1225262/1‑59 (MQ=255)

                                            tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                  >  1:1600160/1‑59 (MQ=255)

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

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

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

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

                                             gTTTACCAACGATCTTCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                      <  1:2443336/70‑1 (MQ=255)

                                                tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa                                   >  1:1092611/1‑54 (MQ=255)

                                                tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                  >  1:3129577/1‑71 (MQ=255)

                                                    aaCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                             <  1:2976428/40‑1 (MQ=255)

                                                     aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                             <  1:2004947/55‑1 (MQ=255)

                                                      cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                           >  1:1723181/1‑40 (MQ=255)

                                                      cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                           >  1:410887/1‑56 (MQ=255)

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

                                                         tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc         >  1:3593713/1‑71 (MQ=255)

                                                          catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                       >  1:708533/1‑40 (MQ=255)

                                                           atcGGTCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                             <  1:2747710/49‑1 (MQ=255)

                                                           atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGtt          <  1:1624796/68‑1 (MQ=255)

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

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

                                                            tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                    <  1:2833879/41‑1 (MQ=255)

                                                            tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                  >  1:1983774/1‑43 (MQ=255)

                                                            tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                  >  1:1508039/1‑43 (MQ=255)

                                                            tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                  >  1:1333188/1‑43 (MQ=255)

                                                            tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa      <  1:3865728/71‑1 (MQ=255)

                                                             cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa      >  1:3344593/1‑70 (MQ=255)

                                                             cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa      >  1:49787/1‑70 (MQ=255)

                                                             cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa      >  1:3647067/1‑70 (MQ=255)

                                                             cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa      >  1:3877880/1‑70 (MQ=255)

                                                             cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa      >  1:2351553/1‑70 (MQ=255)

                                                             cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa      >  1:3761891/1‑70 (MQ=255)

                                                              gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                  >  1:3452161/1‑41 (MQ=255)

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

                                                                gCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc                     <  1:2058693/52‑1 (MQ=255)

                                                                gCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACt  >  1:2839531/1‑71 (MQ=255)

GCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACT  >  minE/608924‑609058

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