BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,595,232	T→G	12.2%	I22L (ATC→CTC) ‡	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,595,232	0	T	G	12.2%	233.9 / 27.5	106	I22L (ATC→CTC) ‡	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base T (33/60); new base G (3/10); total (36/70)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

GTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTA  >  Exported/2595166‑2595302

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

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

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

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

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

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

          cgTTGGCGGAGTTTCCAGATAATGGATAATAGTCTGTTCTGTGCCGCTGTTGCGGatc                                                                       <  1:5663473/58‑1 (MQ=255)

          cgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGatc                                                                       <  1:1422093/58‑1 (MQ=255)

          cgTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTAt                                                          <  1:1505373/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)

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

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

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

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

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

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

                       tCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGcc                                             >  1:532932/1‑71 (MQ=255)

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

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

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

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

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

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

                            aTAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCg                                                             <  1:4125509/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                  gATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTa                                                           >  1:1638362/1‑46 (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:2990478/67‑1 (MQ=255)

                                      attgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                               >  1:765989/4‑70 (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:5451187/55‑1 (MQ=255)

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

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

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

                                          tCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                          <  1:930599/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:1609557/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                 tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCg                                     <  1:2059901/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)

                                                     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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

GTTCATCGCGCGTTGGCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTA  >  Exported/2595166‑2595302

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