BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,201,054	T→C	17.2%	G69G (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,201,054	0	T	C	17.2%	279.3 / 46.2	129	G69G (GGT→GGC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base T (48/58); new base C (8/14); total (57/72)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 133 total aligned reads displayed.

TGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGC  >  minE/1200987‑1201124

tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACCCCTATATCGCCATTATTATGGGTCtt                                                                     >  1:4107399/1‑71 (MQ=255)

tGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCtt                                                                     <  1:4504614/71‑1 (MQ=255)

  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                                                                    >  1:5487791/1‑70 (MQ=255)

  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTg                                                                    <  1:5210681/70‑1 (MQ=255)

  gACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                                                   <  1:3776526/71‑1 (MQ=255)

  gACGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCtgg                                                                    <  1:2125531/70‑3 (MQ=255)

   aCGTGTACGGCTCGCTCTCTCTGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGc                                                                   <  1:2226196/70‑1 (MQ=255)

    cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                                                                 <  1:4697721/71‑1 (MQ=255)

    cGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                                                                 <  1:1742101/71‑1 (MQ=255)

      tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg                                                                >  1:3126734/1‑70 (MQ=255)

      tttATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAgg                                                                >  1:3736191/1‑70 (MQ=255)

        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGt                                                                        >  1:5135927/1‑60 (MQ=255)

        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                                             >  1:5435198/1‑71 (MQ=255)

            gTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                          <  1:912467/70‑1 (MQ=255)

               cACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                                      >  1:1876963/1‑71 (MQ=255)

                aCTGTCGCTGACGGGTAAAGTCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                       <  1:3951882/69‑1 (MQ=255)

                aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                       <  1:1358070/69‑1 (MQ=255)

                aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                                      <  1:106621/70‑1 (MQ=255)

                aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCt                                                      <  1:3041287/70‑1 (MQ=255)

                aCTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTg                                                     >  1:4793423/1‑71 (MQ=255)

                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                    >  1:2059994/1‑70 (MQ=255)

                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                    >  1:5175328/1‑70 (MQ=255)

                  tGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                    >  1:5539730/1‑70 (MQ=255)

                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                                                   <  1:2980099/54‑1 (MQ=255)

                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                                                   <  1:2941518/54‑1 (MQ=255)

                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                          <  1:2673851/63‑1 (MQ=255)

                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                        >  1:2793122/1‑65 (MQ=255)

                    tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAg                                                                 <  1:1518438/55‑1 (MQ=255)

                    tCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                 >  1:1145933/1‑71 (MQ=255)

                        tGACCGGTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAAc                                                            >  1:24917/1‑56 (MQ=255)

                          aCGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCAACCGTg                                             <  1:2093796/69‑1 (MQ=255)

                             ggTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                        >  1:5640380/1‑71 (MQ=255)

                              gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGatat                                        >  1:16107/1‑70 (MQ=255)

                              gTAAAGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                       <  1:3545152/71‑1 (MQ=255)

                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                                             >  1:1116551/1‑48 (MQ=255)

                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                           <  1:568434/50‑1 (MQ=255)

                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                           <  1:3137000/50‑1 (MQ=255)

                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                           <  1:4425788/50‑1 (MQ=255)

                                aaaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGata                                         >  1:2191146/1‑67 (MQ=255)

                                 aaGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGa                                     >  1:797924/1‑70 (MQ=255)

                                 aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                     <  1:5313696/70‑1 (MQ=255)

                                 aaGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                     <  1:3802516/70‑1 (MQ=255)

                                   ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                  <  1:2724485/71‑4 (MQ=255)

                                   ggCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                  <  1:822517/71‑4 (MQ=255)

                                   ggCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                               >  1:3985719/1‑58 (MQ=255)

                                    gCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                  >  1:3360804/1‑67 (MQ=255)

                                    gCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGa                                           >  1:1748441/1‑61 (MQ=255)

                                       accacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                          <  1:2920708/43‑1 (MQ=255)

                                        ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                       <  1:1903575/61‑1 (MQ=255)

                                        ccacACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGattc                                  <  1:5044630/66‑4 (MQ=255)

                                        ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                       <  1:5312026/45‑1 (MQ=255)

                                        ccacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                       <  1:4402897/45‑1 (MQ=255)

                                         cacaCCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                        >  1:1883555/1‑43 (MQ=255)

                                          acacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGa                                           >  1:1579627/1‑55 (MQ=255)

                                          acacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTccc                           >  1:2749234/1‑71 (MQ=255)

                                           cacCGATATCGCCATTATTATGGTTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                          >  1:688421/1‑71 (MQ=255)

                                           cacCGATATCGCCATTATTATGGGTCTTGCAGGTaa                                                             >  1:523184/1‑36 (MQ=255)

                                           cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                          >  1:1683737/1‑71 (MQ=255)

                                           cacCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCg                          >  1:4176217/1‑71 (MQ=255)

                                            acCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCgg                         >  1:5515268/1‑71 (MQ=255)

                                            acCGATATCGCCATTATTATGGGACTTCCACGTaa                                                             >  1:4293903/1‑35 (MQ=21)

                                             ccGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGt                        >  1:1125878/1‑71 (MQ=255)

                                                atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                    <  1:5570109/56‑1 (MQ=255)

                                                atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGAc                                    <  1:2376697/56‑1 (MQ=255)

                                                 tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                    <  1:2652223/71‑1 (MQ=255)

                                                 tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                    <  1:5695205/71‑1 (MQ=255)

                                                 tatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTa                    <  1:14276/71‑1 (MQ=255)

                                                  atCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                 >  1:5268274/1‑41 (MQ=255)

                                                  atCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                   >  1:2619947/1‑71 (MQ=255)

                                                   tCGCCATTATTATGGGTCTTGCAGGTAAGGAACCTGCCACCGTGGatat                                        >  1:8664/1‑49 (MQ=255)

                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                               >  1:799848/1‑42 (MQ=255)

                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                        >  1:3658793/1‑49 (MQ=255)

                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAtt                  >  1:4179570/1‑71 (MQ=255)

                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                   <  1:349867/70‑1 (MQ=255)

                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                   <  1:1733972/70‑1 (MQ=255)

                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTAt                   <  1:5637959/70‑1 (MQ=255)

                                                   tCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAt                   >  1:2231524/1‑70 (MQ=255)

                                                    cGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                >  1:42484/1‑56 (MQ=255)

                                                      ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGt                                              <  1:1404212/40‑1 (MQ=255)

                                                      ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGt                                              <  1:334685/40‑1 (MQ=255)

                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:4998195/51‑1 (MQ=255)

                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:3015882/51‑1 (MQ=255)

                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:166450/51‑1 (MQ=255)

                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:3332616/51‑1 (MQ=255)

                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgcg              <  1:4439293/71‑1 (MQ=255)

                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc               >  1:458109/1‑70 (MQ=255)

                                                        attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:5573431/50‑1 (MQ=255)

                                                        attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:4008774/50‑1 (MQ=255)

                                                        attattATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGa             >  1:3070030/1‑71 (MQ=255)

                                                          tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg           >  1:5206861/1‑71 (MQ=255)

                                                          tattatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACg           >  1:5009563/1‑71 (MQ=255)

                                                          tattatGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTAt                   <  1:918092/63‑1 (MQ=255)

                                                           attatGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGATTCCATCCCCAGTTTTATTc                 >  1:2963194/1‑64 (MQ=255)

                                                            ttatGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTa         <  1:4607497/71‑1 (MQ=255)

                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag        <  1:1330065/70‑1 (MQ=255)

                                                              atGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTag        <  1:5295182/70‑1 (MQ=255)

                                                               tGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaa      <  1:2686352/71‑1 (MQ=255)

                                                                  gTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAagag   <  1:4348460/71‑1 (MQ=255)

                                                                   tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc  <  1:5406685/71‑1 (MQ=255)

                                                                   tCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAgc  >  1:4706057/1‑71 (MQ=255)

TGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGC  >  minE/1200987‑1201124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`