BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,595,245	A→G	21.8%	N17N (AAT→AAC)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,245	0	A	G	21.8%	232.0 / 62.7	119	N17N (AAT→AAC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base A (54/39); new base G (12/14); total (66/53)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 119 total aligned reads displayed.

TGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACA  >  minE/2595179‑2595312

tGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAc                                                                 >  1:3429468/1‑71 (MQ=255)

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

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

  gCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGATGTTGCGGATCATCTCCAGCGTATTAcgc                                                               <  1:572767/71‑1 (MQ=255)

   cGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTAcgcg                                                              <  1:3269313/71‑1 (MQ=255)

      aGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                                  >  1:833701/1‑64 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

             aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                                  >  1:3419162/1‑57 (MQ=255)

             aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTa                                                                  >  1:3889924/1‑57 (MQ=255)

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

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

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

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

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

                   tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAt                                                                    <  1:3019104/49‑1 (MQ=255)

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

                    ggATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAt                                                                    >  1:1299870/1‑48 (MQ=255)

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

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

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

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

                       taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                                             >  1:3084155/1‑52 (MQ=255)

                       taataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGa                                                             >  1:2391373/1‑52 (MQ=255)

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

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

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

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

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

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

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

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

                          taGTCGGTTCTGTGCAGCTGTTGCGGATCATCTCCAGCGTATTAcg                                                                >  1:1114239/1‑46 (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)

                               ggTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGa                                                             >  1:1047466/1‑44 (MQ=255)

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

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

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

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

                                              ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                 <  1:29745/57‑1 (MQ=255)

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

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

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

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

                                               tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCata                  >  1:2579112/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)

                                                    atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt             <  1:3667590/71‑1 (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)

                                                      catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa            >  1:2241728/1‑70 (MQ=255)

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

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

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

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

                                                        tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtct         >  1:2415416/1‑71 (MQ=255)

                                                        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:2656478/68‑1 (MQ=255)

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

                                                            cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa       <  1:3209063/69‑1 (MQ=255)

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

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

                                                              gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtctaa       >  1:307983/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)

TGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACA  >  minE/2595179‑2595312

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