BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,928,602	C→A	8.4%	G651V (GGC→GTC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,928,602	0	C	A	8.4%	333.6 / 8.8	104	G651V (GGC→GTC)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base C (46/49); new base A (3/6); total (49/55)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

CGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACC  >  W3110S.gb/1928532‑1928670

cgcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGc                                                                      >  1:2351159/1‑71 (MQ=255)

    aCCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                                                                    <  1:3521143/69‑1 (MQ=255)

     ccAGAAAAGCATATTCCATCGCTATGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAg                                                                 <  1:3247438/71‑1 (MQ=255)

      cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAgg                                                                <  1:663710/71‑2 (MQ=255)

      cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                >  1:3178725/1‑71 (MQ=255)

      cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                <  1:545332/71‑1 (MQ=255)

      cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                <  1:4225787/71‑1 (MQ=255)

      cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                <  1:4025628/71‑1 (MQ=255)

       aGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                               <  1:4657678/71‑1 (MQ=255)

       aGAAAAACATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                               >  1:8313/1‑71 (MQ=255)

        gAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                               <  1:4954123/70‑1 (MQ=255)

        gAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                               >  1:2157472/1‑70 (MQ=255)

        gAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                                                               >  1:1670060/1‑70 (MQ=255)

         aaaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                >  1:6682251/1‑68 (MQ=255)

         aaaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                >  1:5661749/1‑68 (MQ=255)

         aaaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                >  1:740579/1‑68 (MQ=255)

         aaaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa                                                             >  1:7262287/1‑71 (MQ=255)

          aaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTTAGTCCAt                                                            <  1:1960510/71‑1 (MQ=255)

           aaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATg                                                           >  1:7911079/1‑71 (MQ=255)

            aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                          <  1:2443703/71‑1 (MQ=255)

            aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                          <  1:1555241/71‑1 (MQ=255)

            aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATgg                                                          <  1:7877596/71‑2 (MQ=37)

             gCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                         <  1:689055/71‑1 (MQ=255)

             gCATATTCCATCGCTACGCCTTCGTACAATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                         >  1:2265705/1‑71 (MQ=255)

               atatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                          <  1:5935737/68‑1 (MQ=255)

               atatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                       <  1:454492/71‑1 (MQ=255)

               atatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATgg                                                          <  1:592074/68‑2 (MQ=37)

                tatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                       >  1:2708979/1‑70 (MQ=255)

                tatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                       >  1:4360636/1‑70 (MQ=255)

                tatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                       >  1:2968366/1‑70 (MQ=255)

                  ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                                                                    <  1:7161686/55‑1 (MQ=255)

                  ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                     <  1:7132902/70‑1 (MQ=255)

                    ccATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGcc                                                                     >  1:2916499/1‑52 (MQ=255)

                    ccATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                                                                    >  1:1245169/1‑53 (MQ=255)

                     cATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACag                                                  <  1:2885223/70‑1 (MQ=255)

                       tCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                          <  1:4298361/60‑1 (MQ=255)

                       tCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                          <  1:6150082/60‑1 (MQ=255)

                       tCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACagag                                                >  1:1723800/1‑70 (MQ=255)

                           tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTcc                                                              <  1:4807404/52‑1 (MQ=255)

                           tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTcc                                                              <  1:6804676/52‑1 (MQ=255)

                           tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAAt                                            >  1:2789412/1‑70 (MQ=255)

                           tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAAt                                            >  1:5735107/1‑70 (MQ=255)

                            aCGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCaa                                             >  1:1952725/1‑68 (MQ=255)

                              gCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                        <  1:3526555/71‑1 (MQ=255)

                               ccTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaa                                          >  1:3698953/1‑68 (MQ=255)

                               ccTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAg                                       <  1:2538803/71‑1 (MQ=255)

                                cTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCg                                                        <  1:1538986/53‑1 (MQ=255)

                                cTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGa                                      <  1:6782834/71‑1 (MQ=255)

                                cTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGa                                      <  1:5632970/71‑1 (MQ=255)

                                 ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATg                                                           <  1:4062577/49‑1 (MQ=255)

                                 ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATg                                                           <  1:6655226/49‑1 (MQ=255)

                                 ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATg                                                           <  1:3531160/49‑1 (MQ=255)

                                 ttCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATg                                                           <  1:2134147/49‑1 (MQ=255)

                                  tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGc                                                                      <  1:2337544/37‑1 (MQ=255)

                                  tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                    >  1:4709657/1‑71 (MQ=255)

                                  tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                    >  1:6701188/1‑71 (MQ=255)

                                  tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                    >  1:1017646/1‑71 (MQ=255)

                                   cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGAt                                     <  1:884588/69‑1 (MQ=255)

                                   cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                   >  1:1004967/1‑71 (MQ=255)

                                   cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                   >  1:1316849/1‑71 (MQ=255)

                                   cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                   >  1:2722004/1‑71 (MQ=255)

                                       cGATTTAAATCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCa                                <  1:4782895/70‑1 (MQ=255)

                                       cGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAt                               <  1:443345/71‑1 (MQ=255)

                                       cGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAt                               <  1:2892934/71‑1 (MQ=255)

                                       cGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAt                               <  1:4700849/71‑1 (MQ=255)

                                         aTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                         >  1:835325/1‑43 (MQ=255)

                                         aTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                         >  1:6940503/1‑43 (MQ=255)

                                          tttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                       <  1:5253432/44‑1 (MQ=255)

                                          tttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAt                               >  1:5627417/1‑68 (MQ=255)

                                          tttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCgg                            >  1:4934504/1‑71 (MQ=255)

                                           ttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTc                                 <  1:6981021/65‑1 (MQ=255)

                                            tAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGtt                          >  1:2573765/1‑71 (MQ=255)

                                             aaaGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGtt                          <  1:7522801/70‑1 (MQ=255)

                                              aaGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGtttt                        <  1:7304424/71‑1 (MQ=255)

                                              aaGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGttt                         <  1:175657/70‑1 (MQ=255)

                                               aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCg                             >  1:1439454/1‑65 (MQ=255)

                                               aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGtttt                        >  1:6356176/1‑70 (MQ=255)

                                               aGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGtttt                        >  1:3382133/1‑70 (MQ=255)

                                                gCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCg                             >  1:6791550/1‑64 (MQ=255)

                                                     cAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcg                 >  1:897288/1‑71 (MQ=255)

                                                     cAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTc                  >  1:6235794/1‑70 (MQ=255)

                                                            aCCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAtt          >  1:2400774/1‑71 (MQ=255)

                                                             ccGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAttt         <  1:3711200/71‑1 (MQ=255)

                                                              cGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAtt          >  1:1972397/1‑69 (MQ=255)

                                                                ccATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAtt          <  1:722221/67‑1 (MQ=255)

                                                                acATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGttt                         <  1:252898/51‑1 (MQ=21)

                                                                acATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGttt                         <  1:323805/51‑1 (MQ=21)

                                                                acATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGttt                         <  1:3182955/51‑1 (MQ=21)

                                                                 cATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGc      >  1:4107156/1‑70 (MQ=255)

                                                                 cATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg     <  1:3499132/71‑1 (MQ=255)

                                                                 cATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCg     <  1:3639058/71‑1 (MQ=255)

                                                                  aTGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                       >  1:6829061/1‑52 (MQ=25)

                                                                  aTGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                       >  1:2031024/1‑52 (MQ=25)

                                                                  aTGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                       >  1:5260896/1‑52 (MQ=25)

                                                                   tGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGt                           <  1:3997724/47‑1 (MQ=255)

                                                                   tGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTc                  <  1:5804886/56‑1 (MQ=255)

                                                                   tGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGAc   >  1:7780178/1‑71 (MQ=255)

                                                                    ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTa        >  1:6068105/1‑65 (MQ=255)

                                                                    ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGATCGCGCAATTTAGCGa    >  1:4717049/1‑69 (MQ=255)

                                                                     gCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGAcc  <  1:6521574/70‑1 (MQ=255)

CGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACC  >  W3110S.gb/1928532‑1928670

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