BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,929,501	G→T	8.3%	V424F (GTT→TTT)	sdaB →	L‑serine deaminase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,929,501	0	G	T	8.3%	424.3 / 12.2	131	V424F (GTT→TTT)	sdaB	L‑serine deaminase II
Reads supporting (aligned to +/- strand): ref base G (70/50); new base T (3/8); total (73/58)
Fisher's exact test for biased strand distribution p-value = 6.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 131 total aligned reads displayed.

TAAAAGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGG  >  W3110S.gb/2929433‑2929568

tAAAAGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGtt                                                                   <  1:4857065/71‑1 (MQ=255)

  aaaGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTa                                                                  <  1:718699/70‑1 (MQ=255)

   aaGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATc                                                                <  1:7027248/71‑1 (MQ=255)

    aGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATc                                                                >  1:5528805/1‑70 (MQ=255)

    aGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATc                                                                >  1:5839882/1‑70 (MQ=255)

     gCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGa                                                              <  1:1203676/71‑1 (MQ=255)

      cGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATc                                                                >  1:586335/1‑68 (MQ=255)

      cGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCATCGATAAAGTTATCGaa                                                             <  1:4587840/71‑1 (MQ=255)

          gAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCa                                                         >  1:699722/1‑71 (MQ=255)

           aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCa                                                         >  1:5034058/1‑70 (MQ=255)

            aCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATg                                                       >  1:46447/1‑71 (MQ=255)

             cgccgcACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATg                                                       >  1:1738498/1‑70 (MQ=255)

               ccgcACGTATGGCGCTGCGCCGTACCAGCGCGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTAc                                                    >  1:6694455/1‑71 (MQ=255)

               ccgcACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCAt                                                        <  1:197884/67‑1 (MQ=255)

                 gcACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACg                                                   >  1:6558379/1‑70 (MQ=255)

                    cGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                >  1:5446248/1‑70 (MQ=255)

                    cGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                >  1:3948175/1‑70 (MQ=255)

                     gTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACa                                              <  1:712133/71‑1 (MQ=255)

                      tATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAg                                             <  1:835677/71‑1 (MQ=255)

                      tATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAg                                             <  1:3538357/71‑1 (MQ=255)

                        tGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAgg                                            >  1:7099247/1‑70 (MQ=255)

                         ggCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGt                                                                    <  1:6773974/45‑1 (MQ=255)

                          gcgcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTa                                          >  1:4266160/1‑70 (MQ=255)

                           cgcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaa                                        <  1:2727544/71‑1 (MQ=255)

                           cgcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTa                                          >  1:940751/1‑69 (MQ=255)

                            gcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       <  1:1392129/71‑1 (MQ=255)

                             cTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       <  1:2791887/70‑1 (MQ=255)

                               gcgcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGtt                                                                   <  1:5296372/40‑1 (MQ=255)

                                cgcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatat                                   >  1:3833386/1‑71 (MQ=255)

                                 gcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAc                                                           <  1:6614253/46‑1 (MQ=255)

                                 gcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatat                                   >  1:7733265/1‑70 (MQ=255)

                                 gcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata                                    >  1:677522/1‑69 (MQ=255)

                                  ccGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGa                                 <  1:1233290/71‑1 (MQ=255)

                                  ccGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGa                                 <  1:7384173/71‑1 (MQ=255)

                                   cGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCa                                                         >  1:3277740/1‑46 (MQ=255)

                                   cGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCa                                                         >  1:3351825/1‑46 (MQ=255)

                                   cGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGa                                                  <  1:2046179/53‑1 (MQ=255)

                                   cGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGa                                 <  1:5545298/70‑1 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                >  1:1541778/1‑53 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                >  1:1532428/1‑53 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGa                                      <  1:6028287/63‑1 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACg                              >  1:1761350/1‑71 (MQ=255)

                                         aGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAAc                               >  1:3000381/1‑66 (MQ=255)

                                          gCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGATACAGGTAAAGATATGAACGcc                            >  1:3390138/1‑68 (MQ=255)

                                          gCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       >  1:6468716/1‑57 (MQ=255)

                                          gCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAg                         >  1:2913585/1‑71 (MQ=255)

                                          gCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATAAACGCCAAg                         >  1:679844/1‑71 (MQ=255)

                                           cGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAAc                                               <  1:3830034/48‑1 (MQ=255)

                                              gCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGc                             <  1:1068568/63‑1 (MQ=255)

                                              gCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGc                             <  1:871746/63‑1 (MQ=255)

                                               ccGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatg                                    >  1:6996491/1‑54 (MQ=255)

                                                cgcgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATg                                  <  1:2528305/56‑1 (MQ=255)

                                                cgcgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   <  1:2721011/71‑1 (MQ=255)

                                                 gcgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:4570855/1‑70 (MQ=255)

                                                  cgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGa                 <  1:7384204/71‑1 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGCAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAAc              >  1:3255967/1‑71 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGat                                     <  1:4997287/48‑1 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGat                                     <  1:5728343/48‑1 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGat                                     <  1:3477806/48‑1 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAAc              >  1:4961947/1‑71 (MQ=255)

                                                       cTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACct            >  1:7637019/1‑71 (MQ=255)

                                                         gCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTAcc                     >  1:4820783/1‑60 (MQ=255)

                                                         gCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTAcc                     >  1:4811241/1‑60 (MQ=255)

                                                         gCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:4580813/1‑62 (MQ=255)

                                                         gCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:2982016/1‑62 (MQ=255)

                                                         gCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:3580518/1‑62 (MQ=255)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:426866/65‑2 (MQ=25)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:7398388/65‑2 (MQ=25)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:2601710/65‑2 (MQ=25)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:1069425/65‑2 (MQ=25)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:6437223/65‑2 (MQ=25)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:6524508/65‑2 (MQ=25)

                                                          gcTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       >  1:4662449/2‑41 (MQ=255)

                                                          ccTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGa                                 <  1:2006492/47‑1 (MQ=255)

                                                          ccTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACctc           <  1:4130382/69‑1 (MQ=255)

                                                          ccTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACct            <  1:2203227/68‑1 (MQ=255)

                                                           ttCGATAACTTTATCGAGACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcg        >  1:7591951/2‑71 (MQ=21)

                                                           cTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcg        <  1:4448795/71‑1 (MQ=255)

                                                            tCGATAACTTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:4956487/1‑59 (MQ=21)

                                                            tCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTAc                      >  1:5261191/1‑56 (MQ=255)

                                                            tCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc       >  1:4734270/1‑71 (MQ=255)

                                                            tCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc       >  1:1726656/1‑71 (MQ=255)

                                                             cGATAACTTTATCGAAACCATGTACGAAACAGGTAAAGatat                                   >  1:2014821/1‑42 (MQ=38)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatat                                   >  1:7041894/1‑42 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatat                                   >  1:6031749/1‑42 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata                                    >  1:3450165/1‑41 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata                                    >  1:7862862/1‑41 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata                                    >  1:5890227/1‑41 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc       >  1:6215451/1‑70 (MQ=255)

                                                              gATAACTTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgcg      <  1:5247900/70‑1 (MQ=255)

                                                              gATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACct            <  1:5956545/64‑1 (MQ=255)

                                                              gATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc       <  1:7578641/69‑1 (MQ=255)

                                                              gATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCgcgg     <  1:3310168/71‑1 (MQ=255)

                                                               aTAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGcc                            <  1:1392488/47‑1 (MQ=255)

                                                                tAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       <  1:1398492/35‑1 (MQ=255)

                                                                tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCgcggcg   >  1:6041998/1‑71 (MQ=255)

                                                                tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCgcggcg   <  1:1631415/71‑1 (MQ=255)

                                                                 aaaGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTAcc                     >  1:5731181/1‑52 (MQ=255)

                                                                 aaaGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTAcc                     >  1:7816489/1‑52 (MQ=255)

                                                                 aaaGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGcggcgg  >  1:3462198/1‑71 (MQ=255)

TAAAAGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGG  >  W3110S.gb/2929433‑2929568

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