BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,202,146	G→A	18.8%	K433K (AAG→AAA)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,202,146	0	G	A	18.8%	246.7 / 39.0	101	K433K (AAG→AAA)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base G (40/42); new base A (7/12); total (47/54)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 106 total aligned reads displayed.

CTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATA  >  Exported/1202078‑1202214

ctGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                    >  1:5130055/1‑71 (MQ=255)

   cgccgtACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata                                                                  <  1:1741280/70‑3 (MQ=9)

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

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

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

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

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

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

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

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

          ccAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGcc                                                          <  1:2282120/71‑1 (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:4211468/57‑1 (MQ=2)

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

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

                  ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                               >  1:852721/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)

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

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

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

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

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

                          ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                                                      <  1:4255399/43‑1 (MQ=255)

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

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

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

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

                          ctcGCTGGATAAGGTCATCGAAACGATGCACGAAACCGGTAAg                                                                      <  1:5162626/43‑1 (MQ=255)

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

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

                            cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAgg                                                                     <  1:1468633/42‑1 (MQ=255)

                            cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                    >  1:5576121/1‑43 (MQ=255)

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

                             gCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGa                                                                    <  1:438422/42‑1 (MQ=255)

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

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

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

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

                              cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                      <  1:4143883/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:612330/2‑42 (MQ=16)

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

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

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

                                             gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                       >  1:5074494/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:2100008/1‑64 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

CTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATA  >  Exported/1202078‑1202214

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