BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,595,245	0	A	G	14.8%	303.0 / 37.6	128	N17N (AAT→AAC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base A (50/59); new base G (4/15); total (54/74)
Fisher's exact test for biased strand distribution p-value = 4.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

GCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAA  >  Exported/2595175‑2595313

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

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

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

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

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

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

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

             ttCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGAc                                                             <  1:32204/67‑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:3918693/69‑2 (MQ=255)

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

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

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

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

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

                         gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                                                      >  1:631560/1‑46 (MQ=255)

                         gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                                                      >  1:1638362/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:3382434/55‑1 (MQ=255)

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

                                gTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAtt                                                                    >  1:5068052/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)

                                  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:5081422/50‑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:5660910/53‑1 (MQ=255)

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

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

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

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

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

                                           gCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGt                           >  1:3184184/1‑71 (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:2260799/51‑1 (MQ=255)

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

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

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

                                                  ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCat                    <  1:3272102/71‑1 (MQ=255)

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

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

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

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

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

                                                        atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATgg                                 <  1:136605/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)

                                                         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)

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

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

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

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

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

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

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

                                                              tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa        >  1:1777915/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)

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

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

                                                                   cGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa   <  1:2069308/71‑1 (MQ=255)

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

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

GCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAA  >  Exported/2595175‑2595313

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