BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	982,617	C→T	24.7%	D21D (GAC→GAT)	ydfA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	982,617	0	C	T	24.7%	284.4 / 93.6	162	D21D (GAC→GAT)	ydfA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (54/68); new base T (15/25); total (69/93)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 168 total aligned reads displayed.

AGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGG  >  Exported/982547‑982686

aGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                                                       >  1:560113/1‑71 (MQ=255)

 ggTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACg                                                                      <  1:254718/71‑1 (MQ=255)

     aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                                                  >  1:4293902/1‑71 (MQ=255)

      gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACg                                                                 <  1:1306891/71‑1 (MQ=255)

       cATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGt                                                                >  1:1875911/1‑71 (MQ=255)

         tGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTc                                                              <  1:3036222/71‑1 (MQ=255)

           gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                                                            >  1:5263158/1‑71 (MQ=255)

           gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                                                            <  1:2339448/71‑1 (MQ=255)

           gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                                                            <  1:2359212/71‑1 (MQ=255)

            aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                                                       >  1:767283/1‑59 (MQ=255)

            aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAcc                                                           >  1:93754/1‑71 (MQ=255)

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

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

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

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

                tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCa                                                             >  1:5086101/1‑65 (MQ=255)

                aaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGTAACGTTCACTGcc                                                        >  1:2488450/2‑69 (MQ=255)

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

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

                  tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGtt                                                               >  1:103855/1‑61 (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:2927463/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

                          ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAGGTTCACCGCAATGACGTATAc                                             >  1:5165757/1‑71 (MQ=255)

                          ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  <  1:5432214/66‑1 (MQ=255)

                          ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  <  1:1106693/66‑1 (MQ=255)

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

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

                            caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                         <  1:342520/57‑1 (MQ=255)

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

                            caacaaCGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTc                                                              >  1:1368259/1‑52 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                      tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGc                                                                    <  1:672586/36‑1 (MQ=255)

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

                                      tccctGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCaa                                    <  1:3567720/65‑1 (MQ=255)

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

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

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

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

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

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

                                            gtgtGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                           <  1:5040189/71‑1 (MQ=255)

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

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

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

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

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

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

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

                                             tgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaa                          <  1:5295608/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                    cGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaa                   <  1:1957919/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                  aGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                    <  1:4531709/56‑1 (MQ=255)

                                                                  aGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                    <  1:2312492/56‑1 (MQ=255)

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

                                                                     aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctgg  <  1:3627447/71‑1 (MQ=255)

                                                                      cGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                  >  1:4481924/1‑54 (MQ=255)

AGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGG  >  Exported/982547‑982686

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