BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,979	G→C	8.2%	I425M (ATC→ATG) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,979	0	G	C	8.2%	329.2 / 12.6	109	I425M (ATC→ATG) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (56/44); new base C (5/4); total (61/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

CGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCG  >  minE/608909‑609048

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

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

    gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                                                                   >  1:587614/1‑71 (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)

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

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

                         cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                               <  1:11636/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:2729084/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                         gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                              >  1:1211478/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)

                                             cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                         <  1:2124285/56‑1 (MQ=255)

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

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

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

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

                                             cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                          <  1:2539150/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:618985/69‑1 (MQ=255)

                                               aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                          >  1:3905432/1‑69 (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: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:244663/41‑1 (MQ=255)

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

                                                       cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                              <  1:3988942/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:4081085/1‑37 (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:2830485/1‑71 (MQ=255)

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

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

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

CGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCG  >  minE/608909‑609048

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