BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,981	T→C	9.2%	I425V (ATC→GTC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,981	0	T	C	9.2%	277.2 / 16.7	98	I425V (ATC→GTC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (49/40); new base C (5/4); total (54/44)
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.98e-01

Only 100 of 110 total aligned reads displayed.

TGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTC  >  minE/608912‑609051

tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat                                                                        >  1:657558/1‑70 (MQ=255)

 gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                      >  1:1852293/1‑71 (MQ=255)

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

 gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTCGTTTACCAACGatca                                                                      >  1:685199/1‑71 (MQ=255)

   ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                                                    <  1:910646/71‑1 (MQ=255)

   ttGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcatc                                                                    <  1:1823578/71‑1 (MQ=255)

         aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                                                                     >  1:111947/1‑64 (MQ=255)

         aCAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAAAGatcat                                                                     >  1:1450022/1‑64 (MQ=255)

            gCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                                                               >  1:1288974/1‑67 (MQ=255)

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

               tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                                  <  1:1040688/61‑1 (MQ=255)

               tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                                  <  1:4123172/61‑1 (MQ=255)

               tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                                                                  <  1:2636176/61‑1 (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)

                    ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                                       <  1:1990208/51‑1 (MQ=255)

                    ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                                                                       >  1:3624840/1‑51 (MQ=255)

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

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

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

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

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

                         cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                              <  1:1472220/71‑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)

                             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:2525102/1‑71 (MQ=255)

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

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

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

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

                                    gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                 >  1:2363300/1‑57 (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:808507/1‑69 (MQ=255)

                                      gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                                 >  1:1385780/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)

                                          cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                            <  1:2124285/56‑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)

                                           gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGAAGCAAGCAATAGCGTAGTCATCgttg                                >  1:2618494/1‑67 (MQ=255)

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

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

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

                                            aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                           <  1:2852718/71‑1 (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:421580/71‑1 (MQ=255)

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

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

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

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

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

                                                    cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                 <  1:1817824/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)

                                                      tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGc                 <  1:3672883/71‑1 (MQ=255)

                                                        tGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                 >  1:2006482/1‑37 (MQ=255)

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

                                                        tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTc                        >  1:56076/1‑62 (MQ=255)

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

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

                                                        tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt               >  1:3799987/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)

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

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

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

TGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTC  >  minE/608912‑609051

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