BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,254	C→G	22.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	22.3%	220.2 / 59.3	112	T14T (ACG→ACC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base C (48/39); new base G (11/14); total (59/53)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 115 total aligned reads displayed.

AGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGA  >  minE/2595185‑2595324

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

aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACg                                                                       <  1:879926/71‑1 (MQ=255)

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

  tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                                                     >  1:2816590/1‑71 (MQ=255)

  tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                                                     >  1:3519560/1‑71 (MQ=255)

  tttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTg                                                                     >  1:2440572/1‑71 (MQ=255)

   ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                                                                    >  1:1903062/1‑71 (MQ=255)

    tcgaGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGc                                                                    <  1:210215/67‑1 (MQ=255)

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

    tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGc                                                                    <  1:1160763/70‑1 (MQ=255)

     ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                                                   <  1:3564941/70‑1 (MQ=255)

     ccAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                                                   <  1:1578844/70‑1 (MQ=255)

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

      cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                                                 <  1:545457/71‑1 (MQ=255)

      cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                                                 <  1:3867493/71‑1 (MQ=255)

       aGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCcaca                                                                <  1:3064865/71‑1 (MQ=255)

       aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                                                  >  1:2546880/1‑69 (MQ=255)

       aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                                                  >  1:325666/1‑69 (MQ=255)

         aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAgg                                                              >  1:31737/1‑71 (MQ=255)

          tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                                             >  1:2706900/1‑71 (MQ=255)

          tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                                             <  1:2711013/71‑1 (MQ=255)

            aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGGGACGTGCCGCAGGcc                                                            <  1:1000847/70‑1 (MQ=255)

            aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                                                  <  1:3754353/64‑1 (MQ=255)

             tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCgg                                                          >  1:693838/1‑71 (MQ=255)

                ataataGTCGGTTCTGTGCCGCTGTTGCTGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                       <  1:1677810/71‑1 (MQ=255)

                ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                       <  1:817088/71‑1 (MQ=255)

                ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                       <  1:1060286/71‑1 (MQ=255)

                ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                       <  1:1867950/71‑1 (MQ=255)

                ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGt                                                        >  1:3857410/1‑70 (MQ=255)

                  aataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                                     >  1:1556364/1‑71 (MQ=255)

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

                   ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                                                    <  1:3801900/71‑1 (MQ=255)

                   ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                                                    <  1:57239/71‑1 (MQ=255)

                   ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                                                    <  1:676431/71‑1 (MQ=255)

                    ttgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGa                                                   <  1:672465/69‑1 (MQ=255)

                    taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                                                 >  1:2458016/1‑57 (MQ=255)

                     tgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgtg                                                    >  1:3721245/2‑69 (MQ=255)

                      gTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgtg                                                    >  1:217582/1‑68 (MQ=255)

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

                        cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                                >  1:1503227/1‑70 (MQ=255)

                         ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa                                               <  1:3376541/70‑1 (MQ=255)

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

                         ggTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGAt                                              >  1:3355037/1‑71 (MQ=255)

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

                              tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCggg                                                         <  1:1398041/55‑1 (MQ=255)

                                tgCCGCTGTTGCGTATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                                                    <  1:1139418/58‑1 (MQ=255)

                                tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                       >  1:2599270/1‑71 (MQ=255)

                                  ccGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                         >  1:3530304/1‑67 (MQ=255)

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

                                        ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTT‑CTCat                               <  1:2393005/70‑1 (MQ=255)

                                         tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt                                  >  1:2881427/1‑67 (MQ=255)

                                         tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt                                  >  1:1855616/1‑67 (MQ=255)

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

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

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

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

                                           cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                                <  1:1213714/51‑1 (MQ=255)

                                            ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                                <  1:3132729/50‑1 (MQ=255)

                                            ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                         <  1:2820987/57‑1 (MQ=255)

                                            ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCat                               <  1:3618691/67‑1 (MQ=255)

                                             gATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                          <  1:2906178/71‑1 (MQ=255)

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

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

                                              agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATc                            <  1:471410/66‑1 (MQ=255)

                                              agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAg                          >  1:3036614/3‑70 (MQ=255)

                                              agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGt                         <  1:2259287/69‑1 (MQ=255)

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

                                                 atcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGt                                      >  1:3055752/1‑55 (MQ=255)

                                                 atcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGt                                      >  1:2378487/1‑55 (MQ=255)

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

                                                  tctcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGtt                                                       <  1:1749654/37‑1 (MQ=16)

                                                  tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTAATATCAGTAtct                     >  1:3822883/1‑71 (MQ=255)

                                                  tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGT‑ATGTTGCTCATATCAGTAtct                     >  1:2076000/1‑70 (MQ=255)

                                                    tcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGAt                                              <  1:201980/44‑1 (MQ=255)

                                                    tcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGAt                                              <  1:1714810/44‑1 (MQ=255)

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

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

                                                     ccAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt                  >  1:740956/1‑71 (MQ=255)

                                                      cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtctaa                   <  1:3380061/69‑3 (MQ=255)

                                                      cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATCGTAATGTTGCTCATATCAGTAtct                     <  1:3723317/67‑1 (MQ=255)

                                                      cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtctc                    <  1:1560621/68‑1 (MQ=255)

                                                      cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtctc                    <  1:2656478/68‑1 (MQ=255)

                                                        gCGTGTTGCGTGAGGTGCTACAGGCCGGGTTGTGATAGATGGTAATGTTGATCATATCAGTATCTAATTa                >  1:3424739/1‑70 (MQ=255)

                                                        gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtctaa                   >  1:307983/1‑65 (MQ=255)

                                                        gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtctaa                   >  1:1868750/1‑65 (MQ=255)

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

                                                        gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTAc               <  1:2753213/71‑1 (MQ=255)

                                                        gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCatat                             <  1:2495028/57‑1 (MQ=255)

                                                        gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa                        >  1:2126527/1‑62 (MQ=255)

                                                        gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTa                >  1:1590420/1‑70 (MQ=255)

                                                               aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTa                                          >  1:2242104/1‑37 (MQ=255)

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

                                                               aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg         <  1:754711/70‑1 (MQ=255)

                                                                cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa        >  1:168189/1‑70 (MQ=255)

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

                                                                  cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag     <  1:1810013/71‑1 (MQ=255)

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

                                                                     cGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga  <  1:2427159/71‑1 (MQ=255)

                                                                     cGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga    >  1:697413/1‑69 (MQ=255)

                                                                     cGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga    >  1:1999142/1‑69 (MQ=255)

AGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGA  >  minE/2595185‑2595324

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