BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,200,936	T→C	10.1%	V30A (GTC→GCC)	sdaA →	L‑serine deaminase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,200,936	0	T	C	10.1%	262.9 / 18.9	99	V30A (GTC→GCC)	sdaA	L‑serine deaminase I
Reads supporting (aligned to +/- strand): ref base T (39/50); new base C (7/3); total (46/53)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

Only 100 of 104 total aligned reads displayed.

TTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTG  >  minE/1200869‑1201006

tttAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTtcgtcg                                                                      <  1:147661/70‑1 (MQ=255)

  tAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCgat                                                                    >  1:2879570/1‑70 (MQ=255)

  tAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCgat                                                                    >  1:3013/1‑70 (MQ=255)

      gTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTtcgt                                                                        >  1:221722/1‑62 (MQ=255)

      gTGGGGATTGGTCCCTCATCTTCCAATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCt                                                               >  1:2071507/1‑71 (MQ=255)

       tGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTg                                                              <  1:2983885/71‑1 (MQ=255)

         gggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTgg                                                             <  1:2036541/70‑1 (MQ=255)

         gggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTgg                                                             <  1:2937608/70‑1 (MQ=255)

          ggATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTc                                                           >  1:2651147/1‑71 (MQ=255)

            aTTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTtcgtcg                                                                      >  1:1631828/1‑58 (MQ=255)

            aTTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGa                                                         <  1:2281252/71‑1 (MQ=255)

            aTTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGa                                                         <  1:1988567/71‑1 (MQ=255)

               ggTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCgatga                                                                  >  1:746460/1‑59 (MQ=255)

               ggTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATc                                                                >  1:2916292/1‑61 (MQ=255)

               ggTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaaa                                                       <  1:1937476/70‑1 (MQ=255)

               ggTACCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATc                                                                >  1:2407160/1‑61 (MQ=255)

                gTCCCTCATCTTCCCATACCTTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaaaa                                                      <  1:1065833/70‑1 (MQ=255)

                gTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaaaaa                                                     <  1:3726171/71‑1 (MQ=255)

                gTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaaaaa                                                     <  1:2151284/71‑1 (MQ=255)

                gTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaaaa                                                      <  1:636074/70‑1 (MQ=255)

                 tCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTCGAAAAAg                                                    <  1:2478714/71‑1 (MQ=255)

                 gcccTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAg                                                    <  1:1346123/70‑1 (MQ=255)

                  cccTCATCTTCCCATACCGTAGGGCGTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAg                                                    >  1:1115575/1‑70 (MQ=255)

                  cccTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAgg                                                   >  1:3469448/1‑71 (MQ=255)

                  cccTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAgg                                                   >  1:136334/1‑71 (MQ=255)

                    cTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCt                                                 <  1:3450625/71‑1 (MQ=255)

                    cTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCt                                                 <  1:2134321/71‑1 (MQ=255)

                       aTCTTCCCATACCGTAGGGCCTATGAATGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTa                                               <  1:3382132/70‑1 (MQ=255)

                       aTCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCtt                                                <  1:1176940/69‑1 (MQ=255)

                       aTCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTa                                               <  1:3549101/70‑1 (MQ=255)

                          ttCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGt                                                            >  1:3657344/1‑54 (MQ=255)

                          ttCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGt                                                            >  1:1569778/1‑54 (MQ=255)

                          ttCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTAGAAAAAGGCTTACTg                                            >  1:3207812/1‑70 (MQ=255)

                          ttCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATCATCTGGTCGAAAAAGGCTTACTg                                            >  1:938474/1‑70 (MQ=255)

                          ttCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTCGAAAAAg                                                    >  1:238145/1‑62 (MQ=255)

                          ttCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTCGAAAAAGGCTTa                                               >  1:257246/1‑67 (MQ=255)

                          ttCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTCGAAAAAGGCTTa                                               >  1:3396172/1‑67 (MQ=255)

                          ttCCCATAACGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTg                                            >  1:1565021/1‑70 (MQ=255)

                              cATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTg                                            <  1:3722115/66‑1 (MQ=255)

                              cATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTg                                            <  1:2034837/66‑1 (MQ=255)

                               aTACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGc                                      <  1:2991438/71‑1 (MQ=255)

                               aTACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGc                                      <  1:1345620/71‑1 (MQ=255)

                               aTACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGc                                      <  1:270551/71‑1 (MQ=255)

                               aTACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGc                                      <  1:1471343/71‑1 (MQ=255)

                                tACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaaa                                                       <  1:772388/53‑1 (MQ=255)

                                tACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCg                                     >  1:3105318/1‑71 (MQ=255)

                                tACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCg                                     >  1:184574/1‑71 (MQ=255)

                                tACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCg                                     >  1:3521977/1‑71 (MQ=255)

                                  ccGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGt                                    <  1:2748630/70‑1 (MQ=255)

                                  ccGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGt                                    <  1:2429369/70‑1 (MQ=255)

                                  ccGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGt                                    <  1:1198724/70‑1 (MQ=255)

                                   cGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTa                                  >  1:2134213/1‑71 (MQ=255)

                                     tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACt                                >  1:532484/1‑71 (MQ=255)

                                     tAGGGCCTATGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACt                                <  1:1285065/71‑1 (MQ=255)

                                     tAGGGCCTATGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACt                                <  1:144128/71‑1 (MQ=255)

                                      aGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaaa                                                       <  1:1693188/47‑1 (MQ=255)

                                       gggCCTATGAAGGCAGGTAAACAGTTCGTCgatga                                                                  >  1:344686/1‑35 (MQ=255)

                                       gggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTc                               <  1:3251662/70‑1 (MQ=255)

                                       gggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTc                               >  1:2401307/1‑70 (MQ=255)

                                       gggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTc                               >  1:2599743/1‑70 (MQ=255)

                                       gggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTc                               >  1:706866/1‑70 (MQ=255)

                                       gggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTc                               >  1:1297974/1‑70 (MQ=255)

                                        ggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTcgc                             <  1:483900/71‑1 (MQ=255)

                                        ggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTcgc                             <  1:2710937/71‑1 (MQ=255)

                                        ggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTc                               <  1:2959611/69‑1 (MQ=255)

                                         gCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGt                                    <  1:1461376/63‑1 (MQ=255)

                                          ccTATTAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAATAGGCTTACTGGATAGCGTTACTCGCGt                           <  1:1116225/71‑1 (MQ=255)

                                          ccTATGAAGGCAGGTAAACAGTTCGTCTATGATCTGGTCGa                                                         <  1:695841/41‑1 (MQ=255)

                                          ccTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGt                           <  1:2301903/71‑1 (MQ=255)

                                            tATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaa                                                        <  1:1400259/40‑1 (MQ=255)

                                            tATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGaa                                                        <  1:2128179/40‑1 (MQ=255)

                                            tATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTcgc                             >  1:3401294/1‑67 (MQ=255)

                                            tATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTcgc                             >  1:1240212/1‑67 (MQ=255)

                                             aTGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGt                                    <  1:944734/59‑1 (MQ=255)

                                              tGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGc                                      <  1:3047019/56‑1 (MQ=255)

                                              tGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGcc                       >  1:279724/1‑71 (MQ=255)

                                              tGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGcc                       >  1:446894/1‑71 (MQ=255)

                                              tGAAGGCAGGTAAACAGTTCGCCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGcc                       >  1:2154364/1‑71 (MQ=255)

                                                aaGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAg                                                    >  1:1001878/1‑40 (MQ=255)

                                                aaGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAg                                                    >  1:683819/1‑40 (MQ=255)

                                                  ggCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTgg                   <  1:2298826/71‑1 (MQ=255)

                                                     aGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGTATAGCGTTACTCGCGTTGCCGTGGACg                >  1:2228922/1‑71 (MQ=255)

                                                     aGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACg                >  1:3535608/1‑71 (MQ=255)

                                                      ggTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGAc                 >  1:2757954/1‑69 (MQ=255)

                                                      ggTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGt               <  1:1882991/71‑1 (MQ=255)

                                                      ggTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGt               <  1:2687233/71‑1 (MQ=255)

                                                         aaaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGt                                    <  1:1530858/47‑1 (MQ=255)

                                                         aaaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGt                                    <  1:2733706/47‑1 (MQ=255)

                                                         aaaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTa            >  1:2272869/1‑71 (MQ=255)

                                                         aaaCAGTTCGCCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTa            >  1:2353490/1‑71 (MQ=255)

                                                          aaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTc                               >  1:3652139/1‑51 (MQ=255)

                                                          aaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTc                               >  1:2648491/1‑51 (MQ=255)

                                                          aaCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTAt           >  1:1955002/1‑71 (MQ=255)

                                                           aCAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATg          >  1:1495433/1‑71 (MQ=255)

                                                            cAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATgg         >  1:3173541/1‑71 (MQ=255)

                                                              gTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGtt       >  1:2413687/1‑71 (MQ=255)

                                                              gTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGt        <  1:1330169/70‑1 (MQ=255)

                                                              gTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGt        <  1:1652188/70‑1 (MQ=255)

                                                                 cgtcgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGt                           >  1:634934/1‑48 (MQ=255)

                                                                   tcgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTg  <  1:2341192/71‑1 (MQ=255)

TTTAAGGTGGGGATTGGTCCCTCATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTG  >  minE/1200869‑1201006

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