BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,202,149	C→T	16.6%	D434D (GAC→GAT)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,202,149	0	C	T	16.6%	247.0 / 31.3	96	D434D (GAC→GAT)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base C (38/42); new base T (5/11); total (43/53)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 102 total aligned reads displayed.

CGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCT  >  Exported/1202081‑1202219

cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                                                      <  1:5233609/71‑1 (MQ=255)

 gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                                                     <  1:2945933/71‑1 (MQ=255)

 gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                                                     <  1:1053086/71‑1 (MQ=255)

 gccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                                                     <  1:3386637/71‑1 (MQ=255)

  ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                                    <  1:485031/71‑1 (MQ=255)

  ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                                    <  1:3133687/71‑1 (MQ=255)

  ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                                    <  1:2237289/71‑1 (MQ=255)

  ccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                                    <  1:740483/71‑1 (MQ=255)

   cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATg                                                                     <  1:122387/69‑1 (MQ=255)

   cgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa                                                                   >  1:2482802/1‑71 (MQ=255)

     cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg                                                                 >  1:3152718/1‑71 (MQ=255)

     cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg                                                                 >  1:2128450/1‑71 (MQ=255)

       ccAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc                                                               <  1:2282120/71‑1 (MQ=255)

       ccAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc                                                               <  1:1873247/71‑1 (MQ=255)

        cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc                                                               >  1:4645873/1‑70 (MQ=255)

          gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                            <  1:5386880/71‑1 (MQ=255)

          gTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                            <  1:2388045/71‑1 (MQ=255)

           tGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAAt                                                           <  1:3202003/71‑1 (MQ=255)

            gCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                          >  1:1085806/1‑71 (MQ=255)

             cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                          <  1:588897/70‑1 (MQ=255)

              gccGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATg                                                                     <  1:3728018/57‑1 (MQ=2)

              gccGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATg                                                                     <  1:4211468/57‑1 (MQ=2)

              gccGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATg                                                                     <  1:5078341/57‑1 (MQ=2)

               ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                                    >  1:2673068/1‑58 (MQ=255)

               ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                                    >  1:5550167/1‑58 (MQ=255)

               ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                              <  1:2307687/64‑1 (MQ=255)

                  cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcgcg                                                     >  1:2217578/1‑70 (MQ=255)

                  cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                                    <  1:3744168/71‑1 (MQ=255)

                  cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                                    <  1:4482551/71‑1 (MQ=255)

                  cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGa                                                    <  1:4411502/71‑1 (MQ=255)

                      tctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcgc                                                      <  1:5520954/65‑1 (MQ=255)

                       ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa                                                                   >  1:222754/1‑51 (MQ=255)

                       ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa                                                                   >  1:3029786/1‑51 (MQ=255)

                       ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                              <  1:453443/56‑1 (MQ=255)

                       ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                              <  1:3526690/56‑1 (MQ=255)

                       ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                              <  1:584898/56‑1 (MQ=255)

                        tcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                            >  1:3928558/1‑57 (MQ=255)

                        tcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaa                                                             >  1:2460716/1‑56 (MQ=255)

                         cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCa                                             >  1:4829874/1‑71 (MQ=255)

                           cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCaaa                                                            >  1:2878054/1‑54 (MQ=255)

                           cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           <  1:4143883/71‑1 (MQ=255)

                           cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           >  1:5574699/1‑71 (MQ=255)

                           cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           <  1:466880/71‑1 (MQ=255)

                           cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           <  1:3343172/71‑1 (MQ=255)

                           cTGGATAAGGGCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           <  1:2006494/71‑1 (MQ=255)

                            tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCa                                                              <  1:4834383/49‑1 (MQ=2)

                            tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                                          >  1:1473167/3‑71 (MQ=9)

                            tcgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc                                          >  1:4221071/3‑71 (MQ=9)

                             ggATAAGGTCATCGAAACGATGTACGAAACCGTTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgc                                          >  1:5586504/1‑70 (MQ=255)

                             ggATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                              <  1:2158246/50‑1 (MQ=255)

                             cgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatat                                                                      >  1:1029108/2‑42 (MQ=16)

                             cgATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatat                                                                      >  1:612330/2‑42 (MQ=16)

                                tAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtg                                      <  1:2025406/71‑1 (MQ=255)

                                  aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg                                         <  1:3738303/66‑1 (MQ=255)

                                  aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtgg                                     <  1:5193472/70‑1 (MQ=255)

                                    gttaTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCt                                  <  1:4544821/68‑1 (MQ=12)

                                     tCATCGAAACGATGTAAGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                           <  1:5067371/61‑1 (MQ=255)

                                      cATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTgg                                >  1:1888669/1‑71 (MQ=255)

                                       aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGc                               <  1:3744487/71‑1 (MQ=255)

                                        tCGAAACCATGTACGAAACAGGTAAAGATATGAAc                                                                  <  1:1564200/35‑1 (MQ=9)

                                         cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaa                                                   <  1:5697285/49‑1 (MQ=255)

                                         cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                             >  1:2265199/1‑71 (MQ=255)

                                         cGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                             >  1:759792/1‑71 (MQ=255)

                                         cGAAACGATGTACGAAACCGGTAAGGACATGAACGACAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                             >  1:5452808/1‑71 (MQ=255)

                                          gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATa                                                          >  1:4248574/1‑41 (MQ=255)

                                          gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCaa                             >  1:160257/1‑70 (MQ=255)

                                          gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                            >  1:5074494/1‑71 (MQ=255)

                                          gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                            >  1:3688775/1‑71 (MQ=255)

                                          gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                            >  1:3277154/1‑71 (MQ=255)

                                          gAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAt                            <  1:4977722/71‑1 (MQ=17)

                                           aaaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCgg                                        >  1:2875460/1‑58 (MQ=255)

                                            aaCCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtga                          <  1:3170059/71‑3 (MQ=12)

                                              ccaTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                        <  1:2610426/69‑2 (MQ=12)

                                               caTGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaaga                       <  1:761970/70‑3 (MQ=12)

                                                 tGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGaa                                                   <  1:1006210/41‑1 (MQ=255)

                                                 tGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTc                     <  1:3174331/71‑1 (MQ=255)

                                                   tACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCgg                                        >  1:499273/1‑50 (MQ=17)

                                                   tACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCAAtg                           <  1:1496302/63‑2 (MQ=17)

                                                     cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaaa                        >  1:1115203/1‑64 (MQ=255)

                                                     cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCaaa                        >  1:2100008/1‑64 (MQ=255)

                                                     cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                 <  1:5020958/71‑1 (MQ=255)

                                                     cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                 >  1:1861590/1‑71 (MQ=255)

                                                     cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                 >  1:48996/1‑71 (MQ=255)

                                                     cGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCaa                             >  1:240513/1‑59 (MQ=17)

                                                     cGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGGCCTGGCaa                             >  1:5320277/1‑59 (MQ=17)

                                                       aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt               >  1:2959377/1‑71 (MQ=255)

                                                       aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt               <  1:2330294/71‑1 (MQ=255)

                                                       aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt               <  1:2329887/71‑1 (MQ=255)

                                                       aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt               <  1:3417817/71‑1 (MQ=255)

                                                             gTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa         <  1:1662673/71‑1 (MQ=255)

                                                              tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa         >  1:2895223/1‑70 (MQ=255)

                                                              tAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa         >  1:4780250/1‑70 (MQ=255)

                                                               aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGAc            >  1:1363211/1‑66 (MQ=255)

                                                               aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa       >  1:4117496/1‑71 (MQ=255)

                                                               aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa       >  1:721408/1‑71 (MQ=255)

                                                                  gACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActt    <  1:5143428/71‑1 (MQ=255)

                                                                    cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGTCAATCAAAGTCCAGTGTGACTaa         <  1:2364631/64‑1 (MQ=255)

                                                                    cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct  <  1:3857742/71‑1 (MQ=255)

                                                                    cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtct  <  1:4912416/71‑1 (MQ=255)

                                                                    cATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtc   <  1:1245795/70‑1 (MQ=255)

CGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCT  >  Exported/1202081‑1202219

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