BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,557,809	G→T	5.9%	A155S (GCC→TCC)	secB →	protein export chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,557,809	0	G	T	5.9%	388.7 / 4.6	113	A155S (GCC→TCC)	secB	protein export chaperone
Reads supporting (aligned to +/- strand): ref base G (50/56); new base T (5/2); total (55/58)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

TAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTGCCA  >  minE/2557739‑2557877

tAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATg                                                                      >  1:1936576/1‑71 (MQ=255)

tAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATg                                                                      >  1:127235/1‑71 (MQ=255)

 aaCTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGc                                                                     >  1:1963380/1‑71 (MQ=255)

    ttCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCt                                                                   >  1:1791758/1‑70 (MQ=255)

     tCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtga                                                                 >  1:1010840/1‑71 (MQ=255)

      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:629047/1‑71 (MQ=255)

      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:1631171/1‑71 (MQ=255)

      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:3119868/1‑71 (MQ=255)

      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:1500639/1‑71 (MQ=255)

      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtg                                                                  <  1:1247075/69‑1 (MQ=255)

       gATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGCATCCtgat                                                                <  1:3737989/70‑1 (MQ=255)

       gATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                <  1:3464701/70‑1 (MQ=255)

         tGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGaa                                                             >  1:2181313/1‑71 (MQ=255)

           cgcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAAcc                                                           <  1:1043723/71‑1 (MQ=255)

            gcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:3713169/71‑1 (MQ=255)

            gcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:3401574/71‑1 (MQ=255)

            gcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:3269516/71‑1 (MQ=255)

             cTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCaa                                                         <  1:1368688/71‑1 (MQ=255)

             cTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCaa                                                         <  1:480809/71‑1 (MQ=255)

             cTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:1105298/70‑1 (MQ=255)

                ttCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACg                                                       >  1:1123821/1‑70 (MQ=255)

                   aTGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTaa                                                    <  1:3021589/70‑1 (MQ=255)

                    tGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGaa                                                             >  1:1423354/1‑60 (MQ=255)

                     gAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGc                                                 >  1:3674928/1‑71 (MQ=255)

                        cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGt                                                      <  1:2575221/63‑1 (MQ=255)

                        cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGt                                                      <  1:2579604/63‑1 (MQ=255)

                        cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGt                                                      <  1:746977/63‑1 (MQ=255)

                        cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTc                                              <  1:916661/71‑1 (MQ=255)

                          aTTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACTTAATGCTTCaa                                            >  1:644038/1‑71 (MQ=255)

                          aTTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTaa                                                    >  1:624539/1‑63 (MQ=255)

                          aTTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTaa                                                    >  1:794151/1‑63 (MQ=255)

                             tGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTT‑AATGAc                                        >  1:674477/1‑71 (MQ=255)

                               cagcagcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGAc                                        <  1:2143029/70‑1 (MQ=255)

                                  cagcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgtg                                    <  1:143893/71‑1 (MQ=255)

                                  cagcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                     >  1:3489215/1‑70 (MQ=255)

                                  cagcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                     >  1:1571835/1‑70 (MQ=255)

                                    gcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGa                                   >  1:872594/1‑70 (MQ=255)

                                     cagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTTATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:1744750/71‑1 (MQ=255)

                                      agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                                >  1:2181527/1‑71 (MQ=255)

                                      agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                                >  1:1217168/1‑71 (MQ=255)

                                      agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                                >  1:1559478/1‑71 (MQ=255)

                                      agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                                >  1:507549/1‑71 (MQ=255)

                                        gCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGt                              <  1:1545981/71‑1 (MQ=255)

                                          tGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtg                                      >  1:3169690/1‑61 (MQ=255)

                                          tGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            >  1:1976305/1‑71 (MQ=255)

                                            gCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            <  1:1276473/69‑1 (MQ=255)

                                            gCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCg                          >  1:2141495/1‑71 (MQ=255)

                                              gAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAAAGACTGTGATCGGTGCCgg                         >  1:573549/1‑70 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTc                                              <  1:288534/47‑1 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtg                                      >  1:3634398/1‑55 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtg                                      >  1:3629965/1‑55 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:1569690/1‑71 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:2989645/1‑71 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:2989161/1‑71 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:2829414/1‑71 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:2720194/1‑71 (MQ=255)

                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:827047/1‑71 (MQ=255)

                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:2711177/71‑1 (MQ=255)

                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:3406379/71‑1 (MQ=255)

                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:1180309/71‑1 (MQ=255)

                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:3156958/71‑1 (MQ=255)

                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     >  1:714420/1‑71 (MQ=255)

                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:2941791/71‑1 (MQ=255)

                                                  gTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                     <  1:2430258/54‑1 (MQ=255)

                                                  gTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                     <  1:1241235/54‑1 (MQ=255)

                                                  gTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 >  1:3048263/1‑58 (MQ=255)

                                                    aCTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGa                                   <  1:2458824/54‑1 (MQ=255)

                                                    aCTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGa                                   >  1:1879032/1‑54 (MQ=255)

                                                     cTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTAc                  <  1:1771923/70‑1 (MQ=255)

                                                      tGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCa                                             <  1:2996987/42‑1 (MQ=255)

                                                      tGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCa                                             <  1:1668620/42‑1 (MQ=255)

                                                       gaagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCa                                             <  1:2274385/41‑1 (MQ=255)

                                                       gaagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCa                                             <  1:1745229/41‑1 (MQ=255)

                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTg                             >  1:1040093/1‑56 (MQ=255)

                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTg                             >  1:515175/1‑56 (MQ=255)

                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACgg                <  1:1398186/69‑1 (MQ=255)

                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACgg                <  1:2495991/69‑1 (MQ=255)

                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACgg                <  1:855288/69‑1 (MQ=255)

                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCa              >  1:1288020/1‑71 (MQ=255)

                                                         agaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCAc             >  1:1846527/1‑71 (MQ=255)

                                                         agaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCAc             >  1:2494856/1‑71 (MQ=255)

                                                          gaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTAc                  >  1:3344252/1‑65 (MQ=255)

                                                            aCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGc          >  1:2592756/1‑71 (MQ=255)

                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGcc                           >  1:3057337/1‑53 (MQ=38)

                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            >  1:2125016/1‑52 (MQ=38)

                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            >  1:1371687/1‑52 (MQ=38)

                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGt                    >  1:1664242/1‑60 (MQ=255)

                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTa                   >  1:2123739/1‑61 (MQ=255)

                                                                cAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCtt      >  1:2032265/1‑71 (MQ=255)

                                                                 aGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTg     <  1:1352153/71‑1 (MQ=255)

                                                                   gATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      <  1:2909443/52‑1 (MQ=255)

                                                                   gATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      <  1:2317709/52‑1 (MQ=255)

                                                                    aggCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:3663940/38‑1 (MQ=255)

                                                                    aTGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:518848/40‑1 (MQ=255)

                                                                    aTGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:2434030/40‑1 (MQ=255)

                                                                    aTGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:3153104/40‑1 (MQ=255)

                                                                    aTGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:3445742/40‑1 (MQ=255)

                                                                     tGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:3715663/39‑1 (MQ=255)

                                                                     tGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTGCCa  <  1:2263256/70‑1 (MQ=255)

                                                                      gCCTGATGAACCAACGTAATGCTTCAATGACtgtg                                    <  1:3215041/35‑1 (MQ=255)

TAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTGCCA  >  minE/2557739‑2557877

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`