BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,928,605	G→C	8.5%	S650* (TCA→TGA)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,928,605	0	G	C	8.5%	311.6 / 13.4	107	S650* (TCA→TGA)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base G (48/48); new base C (3/6); total (52/55)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 107 total aligned reads displayed.

CCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAT  >  W3110S.gb/1928537‑1928673

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

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

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

 cAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                   >  1:3178725/1‑71 (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:1670060/1‑70 (MQ=255)

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

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

    aaaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                   >  1:6682251/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:4360636/1‑70 (MQ=255)

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

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

             ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTcc                                                                 >  1:702698/1‑61 (MQ=255)

             ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa                                                                >  1:2915954/1‑62 (MQ=255)

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

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

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

                   cGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                      <  1:4888769/66‑1 (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)

                         gCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACagag                                                   >  1:615600/1‑63 (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:2146654/71‑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:2134147/49‑1 (MQ=255)

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

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

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

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

                             tCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                       >  1:6701188/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:4700849/71‑1 (MQ=255)

                                  cGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAt                                  <  1:2892934/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)

                                       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:323805/51‑1 (MQ=21)

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

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

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

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

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

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

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

                                                             aTGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                          >  1:2031024/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)

                                                                  cTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGa       >  1:2197952/1‑66 (MQ=255)

                                                                  cTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAt  >  1:1858847/1‑71 (MQ=255)

                                                                  cTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAt  >  1:1569986/1‑71 (MQ=255)

                                                                  cTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAt  <  1:2620031/71‑1 (MQ=255)

CCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAT  >  W3110S.gb/1928537‑1928673

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