BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,202,140	C→A	21.8%	T431T (ACC→ACA)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,202,140	0	C	A	21.8%	289.3 / 46.1	110	T431T (ACC→ACA)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base C (37/49); new base A (9/15); total (46/64)
Fisher's exact test for biased strand distribution p-value = 8.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 110 total aligned reads displayed.

GATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGT  >  minE/1202071‑1202206

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

   ggCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGt                                                                 >  1:1703545/1‑70 (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)

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

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

       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: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)

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

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

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

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

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

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

                         ccGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGa                                                       >  1:2673068/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:4482551/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                 ctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCa                                                 <  1:3526690/56‑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:3343172/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

GATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGT  >  minE/1202071‑1202206

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