BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,201,057	T→G	14.9%	L70L (CTT→CTG)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,201,057	0	T	G	14.9%	255.4 / 48.2	128	L70L (CTT→CTG)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base T (49/60); new base G (7/12); total (56/72)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 130 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)

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

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

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

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

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

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

        tATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTa                                                              <  1:2150604/70‑1 (MQ=255)

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

             ttCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                        <  1:481076/71‑1 (MQ=255)

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

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

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

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

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

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

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

                   gTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGc                                                                   <  1:3867789/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)

                        tGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg                                             >  1:2724003/1‑71 (MQ=255)

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

                         gACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTg                                             >  1:2265245/1‑70 (MQ=255)

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

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

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

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

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

                               tAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACg                                                           <  1:2645510/50‑1 (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)

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

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

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

                                  aGGCCACCACACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt                                   <  1:1721752/71‑3 (MQ=255)

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

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

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

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

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

                                     ccaccaCACTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATCGatt                                   <  1:5359276/68‑3 (MQ=255)

                                     ccaccaCACCGATATCGCCATTTTTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTa                                >  1:5507829/1‑71 (MQ=255)

                                       accacACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGa                                                          <  1:520764/43‑1 (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)

                                         cacaCTGATATCGCCATTATTATGGGCCTGGCGGGTAACCTGCCGGATACCGTGGATATc                                       <  1:4232881/60‑1 (MQ=255)

                                         cacaCCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAc                                                        >  1:4101675/1‑43 (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)

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

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

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

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

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

                                                atatCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGttt                      >  1:4485513/1‑70 (MQ=255)

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

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

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

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

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

                                                   tCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCg                                               >  1:181865/1‑42 (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:5637959/70‑1 (MQ=255)

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

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

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

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

                                                      ccATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc               <  1:438116/71‑1 (MQ=255)

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

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

                                                       cATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAg                                  <  1:3015882/51‑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)

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

                                                                 ggTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTagaaga    <  1:5659484/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: `‑`