BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,202,185	T→C	11.9%	G446G (GGT→GGC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,202,185	0	T	C	11.9%	244.0 / 18.8	92	G446G (GGT→GGC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base T (44/37); new base C (2/9); total (46/46)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 101 total aligned reads displayed.

GTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTA  >  minE/1202117‑1202252

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

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

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

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

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

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

        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)

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

                   aaaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAgtgt                                                <  1:2330294/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:4780250/1‑70 (MQ=255)

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

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

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

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

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

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

                                 aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGa                                              >  1:850636/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:1091376/66‑1 (MQ=255)

                                               ccGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCa                         <  1:2010905/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:2064428/40‑1 (MQ=255)

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

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

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

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

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

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

                                                            cgcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                                        >  1:874680/1‑38 (MQ=255)

                                                            cgcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATa                                        >  1:3197837/1‑38 (MQ=255)

                                                            cgcgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCg        <  1:4435621/70‑1 (MQ=255)

                                                              cgGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAAt     <  1:4796274/71‑1 (MQ=255)

                                                                gtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:4335264/61‑1 (MQ=255)

                                                                gtggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:1051568/61‑1 (MQ=255)

                                                                 tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGAt             <  1:5036189/60‑1 (MQ=255)

                                                                 tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTta  <  1:5513497/71‑1 (MQ=255)

                                                                 tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTta  <  1:1874762/71‑1 (MQ=255)

                                                                  ggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTta  >  1:2159459/1‑70 (MQ=255)

                                                                   gtCTGGCAATCAAAGTCCAGTGTGACTAATActtc                                    <  1:5605383/35‑1 (MQ=255)

                                                                   gtCTGGCAATCAAAGTCCAGTGTGACTAATActtc                                    <  1:1465285/35‑1 (MQ=255)

                                                                   gtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTa                                 >  1:3221902/1‑38 (MQ=255)

GTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTA  >  minE/1202117‑1202252

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