BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	1,928,608	T→G	8.9%	D649A (GAC→GCC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	1,928,608	0	T	G	8.9%	247.9 / 21.5	101	D649A (GAC→GCC)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base T (47/45); new base G (3/6); total (50/51)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

Only 100 of 108 total aligned reads displayed.

CAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTG  >  AP009048/1928538‑1928677

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

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

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

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

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

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

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

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

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

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

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

   aaaaGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGt                                                                       >  1:5661749/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:1555241/71‑1 (MQ=255)

      aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                                 <  1:2443703/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:2968366/1‑70 (MQ=255)

          tatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                              >  1:2708979/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)

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

                 tCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                                 <  1:6150082/60‑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)

                         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)

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

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

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

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

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

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

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

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

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

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

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

                             cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                          >  1:1316849/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:6940503/1‑43 (MQ=255)

                                   aTTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                >  1:835325/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)

                                       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: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)

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

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

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

                                                          acATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGttt                                <  1:323805/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: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)

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

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

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

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

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

                                                                    aGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAtttt   >  1:3665055/1‑71 (MQ=255)

                                                                     gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTa               >  1:7852679/1‑58 (MQ=255)

                                                                     gTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTg  <  1:2474147/71‑1 (MQ=255)

CAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTG  >  AP009048/1928538‑1928677

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