BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	1,928,614	0	T	G	7.3%	248.0 / 14.5	96	D647A (GAC→GCC)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base T (48/41); new base G (3/4); total (51/45)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 102 total aligned reads displayed.

AGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTT  >  AP009048/1928544‑1928683

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

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

 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)

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

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

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

      ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa                                                                  <  1:7132902/70‑1 (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)

            cGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                                <  1:4888769/66‑1 (MQ=255)

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

               tACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAAt                                                         >  1:2789412/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:5632970/71‑1 (MQ=255)

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

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

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

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

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

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

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

                       cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                >  1:1004967/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)

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

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

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

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

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

                                                     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:4876954/1‑66 (MQ=255)

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

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

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

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

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

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

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

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

                                                                 ccATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                   >  1:7372419/1‑42 (MQ=255)

                                                                 ccATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCa                                   >  1:3354966/1‑42 (MQ=255)

                                                                  cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAg                                  >  1:56593/1‑42 (MQ=255)

                                                                  cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAAttt                      <  1:2612137/54‑1 (MQ=255)

                                                                  cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc       <  1:4676448/69‑1 (MQ=255)

                                                                  cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCg     <  1:2837661/71‑1 (MQ=255)

                                                                  cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCg     >  1:3729379/1‑71 (MQ=255)

                                                                   aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGa                 >  1:7866021/1‑58 (MQ=255)

                                                                   aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAtttt         >  1:3633281/1‑66 (MQ=255)

                                                                   aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAtttt         >  1:1404964/1‑66 (MQ=255)

                                                                   aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg    <  1:7046640/71‑1 (MQ=255)

                                                                   aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg    <  1:7373488/71‑1 (MQ=255)

                                                                   aTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg    <  1:5361023/71‑1 (MQ=255)

                                                                    tGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGt   <  1:3906738/71‑1 (MQ=255)

                                                                    tGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGt   <  1:1430142/71‑1 (MQ=255)

                                                                    tGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGt   <  1:5585130/71‑1 (MQ=255)

                                                                     gTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAtttt         <  1:2906590/64‑1 (MQ=255)

                                                                     gTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGtt  >  1:4627986/1‑71 (MQ=255)

                                                                      tCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGtt  <  1:6772256/70‑1 (MQ=255)

AGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTT  >  AP009048/1928544‑1928683

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