BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,557,807	A→C	6.3%	D154A (GAT→GCT) ‡	secB →	protein export chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,557,807	0	A	C	6.3%	296.1 / 12.2	110	D154A (GAT→GCT) ‡	secB	protein export chaperone
Reads supporting (aligned to +/- strand): ref base A (55/48); new base C (5/2); total (60/50)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTG  >  Exported/2557737‑2557874

gTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGa                                                                     >  1:3315/1‑71 (MQ=255)

gTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGa                                                                     <  1:3331613/71‑1 (MQ=255)

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

 ttAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGAt                                                                    >  1:3476876/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:210844/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)

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

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

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

         gATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGCATCCtgat                                                             <  1:3737989/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)

               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)

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

                          cTATTTGCAGCAGCAGGCTGGCGAAGGTGCTGAAGAACATCAGGATGCCTGATGAACCAACGt                                                   <  1:614248/63‑1 (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:916661/71‑1 (MQ=255)

                            aTTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACTTAATGCTTCaa                                         >  1:644038/1‑71 (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)

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

                                      gcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGAt                               <  1:373756/71‑1 (MQ=255)

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

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

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

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

                                        agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                             >  1:1217168/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)

                                                 aaGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGc                     <  1:1295375/70‑1 (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:827047/1‑71 (MQ=255)

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

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

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

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

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

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

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

                                                    gTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                  <  1:2979666/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)

                                                          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)

                                                          aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACgg             <  1:2495991/69‑1 (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)

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

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

                                                               cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGt                 >  1:1664242/1‑60 (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)

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

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

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

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

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

GTTAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTG  >  Exported/2557737‑2557874

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