BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,222	C→T	11.4%	S25N (AGC→AAC)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,222	0	C	T	11.4%	301.9 / 22.5	114	S25N (AGC→AAC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base C (44/57); new base T (3/10); total (47/67)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 114 total aligned reads displayed.

ATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTG  >  minE/2595153‑2595290

aTGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCt                                                                     <  1:1599215/71‑1 (MQ=255)

aTGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTGCTGTGCCGCt                                                                     <  1:5434576/71‑1 (MQ=255)

  gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTg                                                                    >  1:3610266/1‑70 (MQ=255)

  gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTg                                                                    >  1:4004157/1‑70 (MQ=255)

  gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGt                                                                   >  1:1113835/1‑71 (MQ=255)

  gAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGt                                                                   >  1:2506547/1‑71 (MQ=255)

   aGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGtt                                                                  >  1:3375841/1‑71 (MQ=255)

       tGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCg                                                               >  1:3540696/1‑70 (MQ=255)

       tGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCg                                                               >  1:420386/1‑70 (MQ=255)

             gTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGatcat                                                         >  1:5034946/1‑70 (MQ=255)

              ttCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGc                                                                <  1:3169537/62‑1 (MQ=255)

              ttCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtct                                                       <  1:1235344/71‑1 (MQ=255)

              ttCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtct                                                       <  1:4188503/71‑1 (MQ=255)

                cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCgg                                                              >  1:3305240/1‑62 (MQ=255)

                cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCgg                                                              >  1:5314449/1‑62 (MQ=255)

                cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtctc                                                      >  1:3720909/1‑70 (MQ=255)

                cATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCAtctc                                                      >  1:156286/1‑70 (MQ=255)

                 aTCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCgg                                                              >  1:398630/1‑61 (MQ=255)

                 aTCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTcc                                                     <  1:677122/70‑1 (MQ=255)

                    gcgcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                                 <  1:1662411/71‑1 (MQ=255)

                      gcgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                               <  1:4507579/71‑1 (MQ=255)

                       cgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAt                                              <  1:1505373/71‑1 (MQ=255)

                         ttGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCgg                                                              <  1:3189027/53‑1 (MQ=255)

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

                            gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCg                                                               <  1:1617702/49‑1 (MQ=255)

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

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

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

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

                                aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGatcat                                                         >  1:2306733/1‑51 (MQ=255)

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

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

                                  tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTg                                                                    <  1:5458828/38‑1 (MQ=255)

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

                                   ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                                 <  1:4389756/56‑1 (MQ=255)

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

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

                                    tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGtt                                                                  >  1:3198932/1‑38 (MQ=255)

                                    tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGat                                                            >  1:3112716/1‑44 (MQ=255)

                                    tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGc                                                  >  1:5279192/1‑54 (MQ=255)

                                    tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                 >  1:532932/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)

                                      cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                         <  1:2921762/61‑1 (MQ=255)

                                      cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgc                                         <  1:2770303/61‑1 (MQ=255)

                                         aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                                 <  1:2008825/50‑1 (MQ=255)

                                         aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCg                                                 <  1:1369154/50‑1 (MQ=255)

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

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

                                               gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                                   <  1:1751469/42‑1 (MQ=255)

                                               gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                                   <  1:559705/42‑1 (MQ=255)

                                               gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                                   <  1:2043411/42‑1 (MQ=255)

                                               gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                                   <  1:3649309/42‑1 (MQ=255)

                                               gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                                   <  1:3691019/42‑1 (MQ=255)

                                               gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAg                                                   <  1:589468/42‑1 (MQ=255)

                                               gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                               >  1:631560/1‑46 (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:2146850/1‑44 (MQ=255)

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

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

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

                                                      gTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAtt                                             >  1:5068052/1‑41 (MQ=255)

                                                      gTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCCTAtt                                             >  1:4003727/1‑41 (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)

                                                        cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                               >  1:4416215/1‑37 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                               ttgCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCgg                        <  1:1573904/52‑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)

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

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

ATGAGTTTGACCAGTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTG  >  minE/2595153‑2595290

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