BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,254	C→G	14.3%	T14T (ACG→ACC)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,254	0	C	G	14.3%	335.6 / 41.1	133	T14T (ACG→ACC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base C (59/55); new base G (4/15); total (63/70)
Fisher's exact test for biased strand distribution p-value = 1.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 134 total aligned reads displayed.

AGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGA  >  minE/2595185‑2595322

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

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

  tttCCAGATAATAGATAATGGTCGGTTCTGTACCGCTGTTGCGGATCATCTCCAGCGTATTACGTGACGt                                                                    <  1:1088262/70‑1 (MQ=255)

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

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

      gaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGccac                                                               <  1:5354453/70‑3 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                           ttCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                          >  1:5562081/1‑71 (MQ=255)

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

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

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

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

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

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

                                tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCTGGTTGTGATAAATGGTAATg                                     <  1:68517/71‑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)

                                 gCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGt                                    >  1:3184184/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                              atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                          <  1:409249/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:2039534/1‑71 (MQ=255)

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

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

                                                catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                           >  1:4305419/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)

                                                    tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                          <  1:446750/46‑1 (MQ=255)

                                                    tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAt                     >  1:4740944/1‑67 (MQ=255)

                                                    tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAt                     >  1:2641067/1‑67 (MQ=255)

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

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

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

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

                                                     ccAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                        >  1:5326773/1‑63 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                           tATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa            <  1:877463/69‑1 (MQ=255)

                                                           tATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa            <  1:1366060/69‑1 (MQ=255)

                                                           tATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa            <  1:1251121/69‑1 (MQ=255)

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

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

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

                                                             ttACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGt        >  1:5286726/1‑71 (MQ=255)

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

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

                                                                cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa     <  1:2640712/71‑1 (MQ=255)

                                                                cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa     <  1:4988950/71‑1 (MQ=255)

                                                                 gcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa    <  1:1537887/71‑1 (MQ=255)

                                                                  cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa      >  1:3237259/1‑68 (MQ=255)

                                                                  cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa      >  1:4694044/1‑68 (MQ=255)

                                                                   gACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga  <  1:1384827/71‑1 (MQ=255)

                                                                    aCGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa    >  1:1645505/1‑68 (MQ=255)

AGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGA  >  minE/2595185‑2595322

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