BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,202,176	A→T	13.7%	S443S (TCA→TCT)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,202,176	0	A	T	13.7%	265.7 / 29.3	102	S443S (TCA→TCT)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base A (45/43); new base T (5/9); total (50/52)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 103 total aligned reads displayed.

CGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATG  >  Exported/1202106‑1202242

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

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

  cTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcg                                                                  <  1:466880/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)

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

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

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

       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)

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

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

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

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

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

                 gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAAt                                                   >  1:5074494/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)

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

                            cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgt                                        <  1:5020958/71‑1 (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:3417817/71‑1 (MQ=255)

                              aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                      <  1:2330294/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: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)

                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTgaga                                    >  1:1917519/1‑59 (MQ=255)

                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                                    >  1:850636/1‑59 (MQ=255)

                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                                    >  1:1910741/1‑59 (MQ=255)

                                            aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                >  1:600639/1‑63 (MQ=255)

                                             tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtctt                        <  1:131133/70‑1 (MQ=255)

                                              gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTa                       >  1:2953457/1‑70 (MQ=255)

                                              gAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTAc                      >  1:2256805/1‑71 (MQ=255)

                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATActtctt                        <  1:134453/68‑1 (MQ=255)

                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTAc                      <  1:1087350/70‑1 (MQ=255)

                                               aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACt                     >  1:5120370/1‑71 (MQ=255)

                                                  gCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTcc                                           <  1:3840604/46‑1 (MQ=255)

                                                   ccAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                                          <  1:3666635/46‑1 (MQ=255)

                                                    cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                                         >  1:2886856/1‑46 (MQ=255)

                                                    cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAg                                         >  1:4025890/1‑46 (MQ=255)

                                                    cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGccc                <  1:1561643/71‑1 (MQ=255)

                                                     aaaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa               >  1:201241/1‑71 (MQ=255)

                                                      aaTACCGCGAAACCTCACGCGGTTGTCTGGCAATCAAAGTCCAgt                                        <  1:1155662/45‑1 (MQ=255)

                                                      aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa               >  1:1569990/1‑70 (MQ=255)

                                                      aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCAt              >  1:514951/1‑71 (MQ=255)

                                                       aTACCGCGAAACCTCACGCGGTGGTCTGGCAATCa                                                 <  1:463645/35‑1 (MQ=255)

                                                       aTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                              >  1:2129067/1‑54 (MQ=255)

                                                       aTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                              >  1:2004569/1‑54 (MQ=255)

                                                        tACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                                               <  1:1328273/36‑2 (MQ=12)

                                                        tACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                                               <  1:5181447/36‑2 (MQ=12)

                                                        tACCGCGAAACCTCTCGCGGCGGCCTGGCAATGaag                                               <  1:4040918/36‑2 (MQ=12)

                                                        tACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCt            >  1:2090788/1‑71 (MQ=255)

                                                         aCCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                >  1:3560348/1‑50 (MQ=255)

                                                         aCCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTaa                                >  1:4849932/1‑50 (MQ=255)

                                                          ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCTATCTGc          <  1:739819/71‑1 (MQ=255)

                                                          ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa               <  1:2010905/66‑1 (MQ=255)

                                                          ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa               <  1:1091376/66‑1 (MQ=255)

                                                          ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGc          <  1:2620639/71‑1 (MQ=255)

                                                            gcgAAACCTCACGCGGTGTTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCaa        <  1:5013659/71‑1 (MQ=255)

                                                            gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                              >  1:5338672/1‑49 (MQ=255)

                                                            gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCaa        <  1:2207420/71‑1 (MQ=255)

                                                             cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCa                                          >  1:2993534/1‑36 (MQ=255)

                                                             cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                      <  1:5617285/40‑1 (MQ=255)

                                                             cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                      <  1:2064428/40‑1 (MQ=255)

                                                             cgAAACCTAACGCGGTGGTCTGGCAATCAAAGTCCa                                          >  1:2164211/1‑36 (MQ=255)

                                                               aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGccc                <  1:2483240/60‑1 (MQ=255)

                                                               aaaCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACgg     <  1:1692655/71‑1 (MQ=255)

                                                                 aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt   <  1:2737675/71‑1 (MQ=255)

                                                                 aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt   <  1:4266285/71‑1 (MQ=255)

                                                                 aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt   <  1:3426651/71‑1 (MQ=255)

                                                                  ccTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATg  <  1:2294783/71‑1 (MQ=255)

                                                                  ccTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATg  <  1:719396/71‑1 (MQ=255)

                                                                    tCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                                    >  1:5195265/1‑35 (MQ=255)

                                                                    tCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                                    >  1:2376557/1‑35 (MQ=255)

CGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATG  >  Exported/1202106‑1202242

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