BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	982,629	C→T	29.7%	T25T (ACC→ACT)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	982,629	0	C	T	29.7%	234.4 / 122.0	168	T25T (ACC→ACT)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (47/71); new base T (21/29); total (68/100)
Fisher's exact test for biased strand distribution p-value = 8.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

Only 100 of 170 total aligned reads displayed.

ATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGA  >  Exported/982559‑982695

 tACTATCGATCTTGGCAACAACGAATCTCTGGTGTACTGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                                   <  1:1293860/71‑1 (MQ=255)

 tACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                                   <  1:1743708/71‑1 (MQ=255)

   cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                                  <  1:2441991/70‑1 (MQ=255)

   cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                                  <  1:3261743/70‑1 (MQ=255)

   cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                                  <  1:165066/70‑1 (MQ=255)

   cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                                  <  1:5578374/70‑1 (MQ=255)

   cTATCGATCTTGGCAACAACGAATCTCTGGTGTACAGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                                  <  1:4647135/70‑1 (MQ=255)

    aaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGcca                                                                >  1:5656746/2‑69 (MQ=255)

     aTCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                               >  1:4801204/1‑71 (MQ=255)

      tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                               >  1:2985504/1‑70 (MQ=255)

      tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                               >  1:2598837/1‑70 (MQ=255)

      tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                               >  1:5238324/1‑70 (MQ=255)

       cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATg                                                              >  1:2933886/1‑70 (MQ=255)

       cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                                                             >  1:3132958/1‑71 (MQ=255)

       cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                                                             >  1:2444799/1‑71 (MQ=255)

       cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                                                             >  1:1089754/1‑71 (MQ=255)

        gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                            >  1:486061/1‑71 (MQ=255)

           ctcggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACg                                                           <  1:3665795/66‑1 (MQ=255)

           cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                         >  1:4194718/1‑71 (MQ=255)

           cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                         >  1:1714743/1‑71 (MQ=255)

             cggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGa                                                             >  1:2261816/2‑65 (MQ=255)

               gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                         <  1:3138491/67‑1 (MQ=255)

               gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                         <  1:1698663/67‑1 (MQ=255)

               gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                                     <  1:4958940/71‑1 (MQ=255)

                caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                                     >  1:4368369/1‑70 (MQ=255)

                 aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                                   <  1:5425847/71‑1 (MQ=255)

                 aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                                   <  1:587474/71‑1 (MQ=255)

                 aacaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCaa                                                   <  1:5612030/71‑1 (MQ=255)

                 aacaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCaa                                                   <  1:3829569/71‑1 (MQ=255)

                  acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaaa                                                  <  1:3543836/71‑1 (MQ=255)

                    aacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                >  1:1635342/1‑55 (MQ=255)

                    aacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCa                                                                 <  1:2332322/54‑1 (MQ=255)

                       gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                             <  1:147057/71‑1 (MQ=255)

                       gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                             <  1:348569/71‑1 (MQ=255)

                        aaTCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                             >  1:5408693/1‑70 (MQ=255)

                         aTCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCaaaa                                           <  1:4875973/71‑1 (MQ=255)

                          tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGc                                               <  1:4885688/66‑1 (MQ=255)

                          tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGc                                               <  1:4505428/66‑1 (MQ=255)

                           ctctGGTGTACGTCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                >  1:693582/1‑48 (MQ=255)

                           ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                            >  1:2656630/1‑68 (MQ=255)

                           ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                            >  1:2412398/1‑68 (MQ=255)

                            tctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata                                                       >  1:1999018/1‑56 (MQ=255)

                             ctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTAAACCAAAAGCAAAACGtt                                       >  1:476508/1‑71 (MQ=255)

                              tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                                         <  1:3554955/68‑1 (MQ=255)

                              tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                                         <  1:1984123/68‑1 (MQ=255)

                                gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                    <  1:4165288/71‑1 (MQ=255)

                                gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                    <  1:5059234/71‑1 (MQ=255)

                                gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                     <  1:5164881/70‑1 (MQ=255)

                                 tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaa                                   <  1:3289837/71‑1 (MQ=255)

                                 tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaa                                   <  1:718535/71‑1 (MQ=255)

                                 tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaa                                   <  1:4568437/71‑1 (MQ=255)

                                 tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                                    >  1:157976/1‑70 (MQ=255)

                                   tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                           <  1:3186065/45‑1 (MQ=255)

                                     cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                >  1:477412/1‑38 (MQ=255)

                                     cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                >  1:3903751/1‑38 (MQ=255)

                                     cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCg                               <  1:2057759/71‑1 (MQ=255)

                                     cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCg                               <  1:4455141/71‑1 (MQ=255)

                                     cGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCg                               <  1:2175992/71‑1 (MQ=255)

                                     cGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCg                               <  1:1420245/71‑1 (MQ=255)

                                      ggCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                           <  1:2668792/42‑1 (MQ=255)

                                      ggCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                           <  1:2238668/42‑1 (MQ=255)

                                      ggCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGa                              <  1:3707637/71‑1 (MQ=255)

                                          tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGt                                        <  1:3731539/57‑1 (MQ=255)

                                          tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGt                                        <  1:3591019/57‑1 (MQ=255)

                                          tgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaaa                                  >  1:2095174/1‑63 (MQ=255)

                                            tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGttt                                      >  1:4496144/1‑57 (MQ=255)

                                             ttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                            <  1:790443/50‑1 (MQ=255)

                                              tCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                         >  1:1068631/1‑36 (MQ=255)

                                              tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaaa                                  >  1:3939251/1‑59 (MQ=255)

                                                cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                           >  1:1941308/1‑48 (MQ=255)

                                                cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                    >  1:578358/1‑71 (MQ=255)

                                                cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                    >  1:5012760/1‑71 (MQ=255)

                                                cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                    >  1:4718072/1‑71 (MQ=255)

                                                cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc                     <  1:1170512/70‑1 (MQ=255)

                                                 caaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAttt                                      >  1:5544429/2‑52 (MQ=255)

                                                 caaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAttt                                      >  1:3571703/2‑52 (MQ=255)

                                                 caaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAttt                                      >  1:1347487/2‑52 (MQ=255)

                                                 caaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAttt                                      >  1:2897144/2‑52 (MQ=255)

                                                  aaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcg                  <  1:4637559/71‑1 (MQ=255)

                                                  aaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcg                  <  1:5188924/71‑1 (MQ=255)

                                                  aaCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                                         <  1:2747890/48‑1 (MQ=255)

                                                     cAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAAca                                         <  1:1400611/45‑2 (MQ=25)

                                                      aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCga              <  1:3247907/71‑2 (MQ=255)

                                                        gACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                           >  1:2560055/1‑40 (MQ=255)

                                                        gACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                             <  1:667521/54‑1 (MQ=255)

                                                             cACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                             <  1:1716495/49‑1 (MQ=255)

                                                             cACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                             <  1:1800783/49‑1 (MQ=255)

                                                             cACGTTCACCACAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                             <  1:4415941/49‑1 (MQ=255)

                                                              aCGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                             <  1:3474791/48‑1 (MQ=255)

                                                              aCGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                             <  1:528210/48‑1 (MQ=255)

                                                              aCGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                             <  1:3086563/48‑1 (MQ=255)

                                                               cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAAATGGTGCGCGCCgatggttg        >  1:3871693/1‑68 (MQ=255)

                                                               cGTTCACCGCAATGACGTATACCAATAGCAAAACGTTTAAAACCGAAAATGGTGcc                    <  1:8980/56‑1 (MQ=255)

                                                                 ttCACTGCCATGACGTATACCAAAAGCAAAACAtt                                       <  1:3068998/35‑1 (MQ=25)

                                                                 ttCACTGCCATGACGTATACCAAAAGCAAAACAtt                                       <  1:4011734/35‑1 (MQ=25)

                                                                  tCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGaaa                            >  1:2965193/1‑45 (MQ=38)

                                                                   cACCGCAATGACGTATACCAAAAGCAAAACTTTTAAAAcc                                <  1:5110408/40‑1 (MQ=255)

                                                                   cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga  <  1:4180364/70‑1 (MQ=255)

                                                                    aCCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                    <  1:1145679/35‑1 (MQ=255)

                                                                    aCCGCAATGACGTATACCAAAAGCAAAACGTTTAAAAc                                 >  1:5433123/1‑38 (MQ=255)

ATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGA  >  Exported/982559‑982695

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