BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,202,128	G→C	18.3%	T427T (ACG→ACC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,202,128	0	G	C	18.3%	276.7 / 44.1	109	T427T (ACG→ACC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base G (39/50); new base C (6/14); total (45/64)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 109 total aligned reads displayed.

TAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCA  >  Exported/1202059‑1202195

tAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                                                    >  1:573100/1‑71 (MQ=255)

 aaCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGa                                                                    <  1:2635090/70‑1 (MQ=255)

  aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                                                   <  1:1416027/70‑1 (MQ=255)

  aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                                                   <  1:5636342/70‑1 (MQ=255)

  aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                                                  <  1:214850/71‑1 (MQ=255)

  aCGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                                                  >  1:1591409/1‑71 (MQ=255)

      cgcacgTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACg                                                              <  1:2094673/71‑1 (MQ=12)

       gcgcgGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGa                                                             <  1:3272937/71‑1 (MQ=255)

          cgtaTGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAAc                                                          <  1:3611663/68‑1 (MQ=12)

            gATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCg                                                        >  1:3287383/1‑71 (MQ=255)

                gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                                                     <  1:374468/70‑1 (MQ=255)

                gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                                                     <  1:602788/70‑1 (MQ=255)

                gCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTa                                                     <  1:5028452/70‑1 (MQ=255)

                  gctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaag                                                  <  1:3752269/70‑3 (MQ=12)

                  gctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaag                                                  <  1:155195/70‑3 (MQ=12)

                   ctGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaa                                                    >  1:1338565/1‑68 (MQ=16)

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

                    tGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                                                               <  1:279462/56‑1 (MQ=255)

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

                      cgccgcACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGa                                                             >  1:4681059/1‑56 (MQ=255)

                      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:740483/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:2128450/1‑71 (MQ=255)

                           cACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACg                                         >  1:3152718/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)

                               aGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTa                                                                >  1:3994514/1‑44 (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)

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

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

                                        cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGaaa                                                           >  1:5416344/1‑40 (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)

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

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

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

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

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

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

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

                                              tcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAgg                                                  <  1:5378160/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:466880/71‑1 (MQ=255)

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

                                                 cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                   >  1:5574699/1‑71 (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)

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

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

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

                                                    gATAAAGTTATCGAAACCATGTACGAAACAGGTaaa                                                   <  1:69973/36‑2 (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: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)

TAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCA  >  Exported/1202059‑1202195

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