BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,064	G→T	27.5%	T730T (ACC→ACA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,064	0	G	T	27.5%	152.6 / 37.5	69	T730T (ACC→ACA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (27/23); new base T (8/11); total (35/34)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

AGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAA  >  Exported/607995‑608131

aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg                                                                    <  1:87811/71‑1 (MQ=255)

aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg                                                                    <  1:3814349/71‑1 (MQ=255)

   aataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGCATTGTCa                                                                 <  1:569736/71‑1 (MQ=255)

    ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                                <  1:472346/71‑1 (MQ=255)

    ataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                                >  1:1064665/1‑71 (MQ=255)

    ataaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAg                                                                >  1:4040379/1‑71 (MQ=255)

     taaCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                                                 <  1:3450180/69‑1 (MQ=255)

       aCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                                                 >  1:2340301/1‑67 (MQ=255)

         gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                           >  1:3563115/1‑71 (MQ=255)

         gTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                                                           >  1:5527442/1‑71 (MQ=255)

          tCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                                   <  1:5225139/62‑1 (MQ=255)

            ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                                                                 >  1:753132/1‑62 (MQ=255)

            ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                         >  1:2966344/1‑70 (MQ=255)

            ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                         >  1:3375138/1‑70 (MQ=255)

            ctgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                         >  1:1178530/1‑70 (MQ=255)

            ctgctgcagTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                                                         >  1:562899/1‑70 (MQ=255)

             tgctgctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                                                              <  1:1232043/64‑1 (MQ=255)

                  ctgTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                                  <  1:1119205/71‑1 (MQ=255)

                     ttCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                                <  1:1488244/70‑1 (MQ=255)

                      tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                              >  1:2144533/1‑71 (MQ=255)

                      tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                              >  1:5385737/1‑71 (MQ=255)

                      tCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg                                              >  1:3291109/1‑71 (MQ=255)

                       cTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                                                                >  1:4205704/1‑52 (MQ=255)

                           aGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                                                  <  1:5332446/62‑1 (MQ=255)

                            gTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCgg                                              >  1:3904211/1‑65 (MQ=9)

                            gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc                                        >  1:4713787/1‑71 (MQ=255)

                            gTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc                                        <  1:2381112/71‑1 (MQ=255)

                             tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc                                       <  1:4341197/71‑1 (MQ=255)

                             tCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc                                       <  1:1067665/71‑1 (MQ=255)

                               aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                                           <  1:862411/65‑1 (MQ=255)

                               aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt                                           <  1:5478999/65‑1 (MQ=255)

                               aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                     >  1:4465790/1‑71 (MQ=255)

                                gCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                                                   <  1:2426693/40‑1 (MQ=255)

                                 cGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATcg                                   >  1:2599153/1‑71 (MQ=255)

                                   cgTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTcc                                                <  1:5500462/55‑1 (MQ=12)

                                     ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                                                 <  1:1855915/53‑1 (MQ=255)

                                     ttGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc                                                 <  1:1169838/53‑1 (MQ=255)

                                       gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCa                                      >  1:3090896/1‑62 (MQ=255)

                                       gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                              <  1:1841906/70‑1 (MQ=255)

                                         aaCGTACTGCGTAGCCAGATACAGGGATGGTCAGCTGCGGATATTTTTcc                                                <  1:4784597/50‑1 (MQ=255)

                                          aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                          <  1:2574704/71‑1 (MQ=255)

                                          aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                            >  1:4723853/1‑69 (MQ=255)

                                             tACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                                                    <  1:1550897/42‑1 (MQ=255)

                                               cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                     >  1:2496566/3‑71 (MQ=12)

                                               cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAt                      >  1:2145640/3‑70 (MQ=12)

                                                 gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc                                                <  1:4827450/42‑1 (MQ=255)

                                                 gCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTc                   >  1:3894100/1‑71 (MQ=16)

                                                  cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                     >  1:4987806/1‑52 (MQ=255)

                                                  cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                    <  1:2534990/69‑1 (MQ=255)

                                                  cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                  >  1:686846/1‑71 (MQ=255)

                                                     aGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga               >  1:611247/1‑70 (MQ=12)

                                                      gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCgat              >  1:204410/1‑69 (MQ=12)

                                                         agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtc            <  1:3476824/70‑1 (MQ=9)

                                                         aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca           >  1:4707475/1‑71 (MQ=255)

                                                          gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat          <  1:1890052/68‑1 (MQ=9)

                                                          gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat          <  1:1208217/68‑1 (MQ=9)

                                                          gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat          <  1:3047823/68‑1 (MQ=9)

                                                            tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc         <  1:4459557/70‑1 (MQ=255)

                                                            tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc         <  1:3986411/70‑1 (MQ=255)

                                                            gacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcatta        >  1:2583400/2‑69 (MQ=9)

                                                             acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtc            >  1:2752349/1‑66 (MQ=255)

                                                             acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc       >  1:125761/1‑71 (MQ=255)

                                                             acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc       >  1:378263/1‑71 (MQ=255)

                                                              cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg      <  1:5515929/71‑1 (MQ=12)

                                                              cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg      <  1:3994210/71‑1 (MQ=12)

                                                              cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg      <  1:2610403/71‑1 (MQ=12)

                                                              cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg      <  1:2878476/71‑1 (MQ=12)

                                                              cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACg      >  1:1900190/1‑71 (MQ=255)

                                                                cGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt    <  1:724340/71‑1 (MQ=255)

                                                                 ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt    >  1:2256425/1‑70 (MQ=255)

                                                                  gATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa  >  1:3215964/1‑71 (MQ=255)

AGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAA  >  Exported/607995‑608131

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