BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,557,808	T→A	5.9%	D154E (GAT→GAA) ‡	secB →	protein export chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,557,808	0	T	A	5.9%	331.2 / 6.8	113	D154E (GAT→GAA) ‡	secB	protein export chaperone
Reads supporting (aligned to +/- strand): ref base T (52/54); new base A (5/2); total (57/56)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 118 total aligned reads displayed.

TTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTGCCA  >  Exported/2557738‑2557877

ttAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGAt                                                                       >  1:2006458/1‑71 (MQ=255)

ttAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGAt                                                                       >  1:3183030/1‑71 (MQ=255)

ttAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGAt                                                                       >  1:3476876/1‑71 (MQ=255)

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

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

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

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

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

       cGATGCGCTGTTCATGAACTATTTGCAACAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:2481717/1‑71 (MQ=255)

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

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

         aTGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgatg                                                               <  1:2223845/70‑1 (MQ=255)

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

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

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

             gcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:3713169/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)

                  tCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGTACCAACGt                                                      >  1:3071454/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:1647213/63‑1 (MQ=255)

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

                         cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTc                                              <  1:2681608/71‑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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                     aCTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGa                                   <  1:2458824/54‑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:515175/1‑56 (MQ=255)

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

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

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

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

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

                                                          agaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCAc             >  1:1846527/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:1371687/1‑52 (MQ=38)

                                                              cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            >  1:2125016/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)

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

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

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

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

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

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

TTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTGCCA  >  Exported/2557738‑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: `‑`