BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,248	A→G	14.8%	R16R (CGT→CGC)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,248	0	A	G	14.8%	303.4 / 38.5	128	R16R (CGT→CGC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base A (51/58); new base G (4/15); total (55/73)
Fisher's exact test for biased strand distribution p-value = 4.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 133 total aligned reads displayed.

GGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTG  >  minE/2595180‑2595317

ggCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAc                                                                      <  1:1590286/70‑1 (MQ=255)

 gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                                                    <  1:2444658/71‑1 (MQ=255)

 gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcg                                                                     <  1:772138/70‑1 (MQ=255)

   ggAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                                                  <  1:3161115/71‑1 (MQ=255)

   ggAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                                                  <  1:26413/71‑1 (MQ=255)

     aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                                                >  1:214205/1‑71 (MQ=255)

      gTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGt                                                               <  1:3291164/71‑1 (MQ=255)

        ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGAc                                                                 <  1:32204/67‑1 (MQ=255)

         tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                                                             >  1:3430818/4‑70 (MQ=255)

         tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                                            >  1:532932/1‑71 (MQ=255)

          ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                                           >  1:2515271/1‑71 (MQ=255)

           gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGcca                                                           <  1:4622529/69‑2 (MQ=255)

           gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGcca                                                           <  1:3918693/69‑2 (MQ=255)

           gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGcca                                                           <  1:1903951/69‑2 (MQ=255)

                 aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCa                                                         >  1:4934723/1‑66 (MQ=255)

                 aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCa                                                         >  1:3825895/1‑66 (MQ=255)

                    gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGt                                                               >  1:446952/1‑57 (MQ=255)

                    gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGt                                                 >  1:2619882/1‑71 (MQ=255)

                     ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                >  1:2615485/1‑71 (MQ=255)

                      taataGTCGGTTCTGTGCC‑CTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                              <  1:552968/71‑1 (MQ=255)

                        attgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                              >  1:765989/4‑70 (MQ=255)

                        attgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                              <  1:2990478/67‑1 (MQ=255)

                        ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                              <  1:1234973/70‑1 (MQ=255)

                         taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                                       >  1:2270696/1‑44 (MQ=255)

                         taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                                            <  1:5624812/55‑1 (MQ=255)

                         taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                                            <  1:5451187/55‑1 (MQ=255)

                           gTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGa                                            >  1:1025140/1‑69 (MQ=255)

                            tCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                         <  1:930599/71‑1 (MQ=255)

                            tCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                         <  1:1132343/71‑1 (MQ=255)

                            tCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATa                                          <  1:3370188/70‑1 (MQ=255)

                             cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa                                        >  1:5471374/1‑71 (MQ=255)

                              ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGt                                                 >  1:4083819/1‑61 (MQ=255)

                              ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAAt                                       <  1:1609557/71‑1 (MQ=255)

                              ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAAt                                       >  1:5296498/1‑71 (MQ=255)

                               gTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                                           <  1:3248690/50‑1 (MQ=255)

                                ttCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                     <  1:5046617/71‑1 (MQ=255)

                                  cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                    >  1:2132299/1‑70 (MQ=255)

                                  cTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                    <  1:1929783/70‑1 (MQ=255)

                                  cTGTGCCGCTGTTGCAGATCATCTCCAGCGTATTACGCGACGTGCCGCa                                                         >  1:3609444/1‑49 (MQ=255)

                                   tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                                                  <  1:4425046/39‑1 (MQ=255)

                                   tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCg                                                    <  1:2059901/53‑1 (MQ=255)

                                   gttgCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTaa                                  <  1:3594924/69‑1 (MQ=255)

                                     tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGAt                                           <  1:166696/60‑1 (MQ=255)

                                     tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTa                                   <  1:1049420/68‑1 (MQ=255)

                                     tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTa                                   <  1:3335105/68‑1 (MQ=255)

                                     tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                <  1:4748671/71‑1 (MQ=255)

                                     tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                <  1:5096028/71‑1 (MQ=255)

                                       ccGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGtt                              >  1:5296390/1‑71 (MQ=255)

                                       ccGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGtt                              >  1:2829626/1‑71 (MQ=255)

                                        cGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTg                             <  1:4973189/71‑1 (MQ=255)

                                           tGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                              <  1:5333956/51‑1 (MQ=255)

                                           tGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                              <  1:2260799/51‑1 (MQ=255)

                                           tGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt                           >  1:2473706/1‑70 (MQ=255)

                                             ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                         <  1:5669419/54‑1 (MQ=255)

                                             ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                         <  1:4345273/54‑1 (MQ=255)

                                              tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata                       >  1:3934983/1‑71 (MQ=255)

                                               gCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGc                            >  1:5166536/1‑65 (MQ=255)

                                               gCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGc                            >  1:613852/1‑65 (MQ=255)

                                                cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCatat                      <  1:2095784/70‑1 (MQ=255)

                                                cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCatat                      <  1:928998/70‑1 (MQ=255)

                                                cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATc                     <  1:4791176/71‑1 (MQ=255)

                                                 ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCa                    <  1:2723350/71‑1 (MQ=255)

                                                   atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                     <  1:136605/52‑1 (MQ=255)

                                                   atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                     <  1:409249/52‑1 (MQ=255)

                                                   atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                     <  1:1932403/52‑1 (MQ=255)

                                                   atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                   >  1:5677860/1‑70 (MQ=255)

                                                   agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                              <  1:4023910/41‑1 (MQ=255)

                                                   agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                              <  1:2933133/41‑1 (MQ=255)

                                                    tcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa                 >  1:1044466/1‑71 (MQ=255)

                                                    tcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa                 >  1:2039534/1‑71 (MQ=255)

                                                    gcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATag                                         <  1:1504319/46‑2 (MQ=255)

                                                     catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                      >  1:4305419/1‑49 (MQ=255)

                                                     catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                      >  1:1088043/1‑49 (MQ=255)

                                                     catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                      >  1:3673613/1‑49 (MQ=255)

                                                      atcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtc               >  1:932610/1‑71 (MQ=255)

                                                       tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt                  <  1:3404029/67‑1 (MQ=255)

                                                         tcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                              <  1:132499/37‑1 (MQ=16)

                                                         tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa            >  1:454615/1‑71 (MQ=255)

                                                         tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa            >  1:5332946/1‑71 (MQ=255)

                                                           cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtct              <  1:2660411/67‑1 (MQ=255)

                                                           cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt           >  1:2048545/1‑70 (MQ=255)

                                                           cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt           >  1:1323531/1‑70 (MQ=255)

                                                           cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt           >  1:557149/1‑70 (MQ=255)

                                                           cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt           >  1:2904947/1‑70 (MQ=255)

                                                             gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTa         >  1:512151/1‑70 (MQ=255)

                                                             gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa            >  1:5641697/1‑67 (MQ=255)

                                                             gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa            <  1:3064936/67‑1 (MQ=255)

                                                              cGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATgg                                     <  1:1165287/41‑1 (MQ=255)

                                                              cGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGTTAATGTTGCTCATATCAGTATCTCATTAc        >  1:2334989/1‑70 (MQ=255)

                                                              cGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTAc        >  1:219959/1‑70 (MQ=255)

                                                              cGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTAc        >  1:3241682/1‑70 (MQ=255)

                                                              cGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTAc        >  1:2941594/1‑70 (MQ=255)

                                                               gTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa      <  1:1465342/71‑1 (MQ=255)

                                                               gTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa      <  1:72747/71‑1 (MQ=255)

                                                                  ttACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa                 >  1:4813361/1‑57 (MQ=255)

                                                                  ttACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAtt          >  1:5307268/1‑64 (MQ=255)

                                                                  ttACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAtt          >  1:1536322/1‑64 (MQ=255)

                                                                    aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                    <  1:316472/36‑1 (MQ=255)

                                                                    aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa            >  1:2854054/1‑60 (MQ=255)

                                                                    aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg  >  1:3929362/1‑70 (MQ=255)

GGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTG  >  minE/2595180‑2595317

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