BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,595,278	A→G	19.0%	I6I (ATT→ATC)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,595,278	0	A	G	19.0%	302.8 / 63.0	147	I6I (ATT→ATC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base A (56/63); new base G (12/16); total (68/79)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 158 total aligned reads displayed.

TCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGC  >  Exported/2595208‑2595347

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

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

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

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

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

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

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

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

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

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

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

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

         tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCTGGTTGTGATAAATGGTAATg                                                              <  1:68517/71‑1 (MQ=255)

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

          gCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGt                                                             >  1:3184184/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

                    cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATc                                                   <  1:4791176/71‑1 (MQ=255)

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

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

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

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

                         catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                                    >  1:4305419/1‑49 (MQ=255)

                         catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                                    >  1:3673613/1‑49 (MQ=255)

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

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

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

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

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

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

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

                              ccAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                                                 >  1:5326773/1‑63 (MQ=255)

                               cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt                                         >  1:2048545/1‑70 (MQ=255)

                               cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt                                         >  1:2904947/1‑70 (MQ=255)

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

                                 gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                                          <  1:3064936/67‑1 (MQ=255)

                                  cGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATgg                                                                   <  1:1165287/41‑1 (MQ=255)

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

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

                                    tATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa                                     <  1:1366060/69‑1 (MQ=255)

                                    tATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACa                                     <  1:1251121/69‑1 (MQ=255)

                                      ttACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa                                               >  1:4813361/1‑57 (MQ=255)

                                      ttACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAtt                                        >  1:1536322/1‑64 (MQ=255)

                                      ttACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGt                                 >  1:5286726/1‑71 (MQ=255)

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

                                           cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg                                >  1:2678063/1‑67 (MQ=255)

                                           cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                               >  1:1196369/1‑68 (MQ=255)

                                           cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa                               >  1:3237259/1‑68 (MQ=255)

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

                                             aCGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                             >  1:1645505/1‑68 (MQ=255)

                                             aCGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                             >  1:3725178/1‑68 (MQ=255)

                                               gTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACagagag                        <  1:1930385/71‑1 (MQ=255)

                                                tGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCat                                                      <  1:4993596/40‑1 (MQ=255)

                                                tGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAgagaga                       <  1:3901486/71‑1 (MQ=255)

                                                 gCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACagagagc                       >  1:2612185/1‑69 (MQ=255)

                                                  ccGCAGGCGGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACg                     <  1:4105189/71‑1 (MQ=255)

                                                   cacAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCgc                    <  1:3027353/69‑2 (MQ=255)

                                                   cacAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCgc                    >  1:5072877/3‑70 (MQ=255)

                                                    gCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg                                <  1:4014259/58‑1 (MQ=255)

                                                     cAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAg                  <  1:2406121/71‑1 (MQ=255)

                                                      aGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAgc                 >  1:61564/1‑71 (MQ=255)

                                                       ggCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACaaa                                   >  1:5233959/1‑52 (MQ=255)

                                                       ggCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACaaa                                   >  1:4637427/1‑52 (MQ=255)

                                                       ggCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGaca                             <  1:586702/58‑3 (MQ=255)

                                                       ggCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGaca                             <  1:1664765/58‑3 (MQ=255)

                                                       ggCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGaca                             <  1:4191044/58‑3 (MQ=255)

                                                       ggCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt                                                <  1:575553/39‑1 (MQ=255)

                                                       ggCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa                                    <  1:855284/51‑1 (MQ=255)

                                                       ggCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa                                    <  1:5605779/51‑1 (MQ=255)

                                                        gCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAgcgc               <  1:455331/71‑1 (MQ=255)

                                                         ccGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtctc                                           >  1:2681596/1‑42 (MQ=255)

                                                         ccGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                                          >  1:2642026/1‑43 (MQ=255)

                                                         ccGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                                          >  1:1526783/1‑43 (MQ=255)

                                                         ccGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                                          >  1:2551138/1‑43 (MQ=255)

                                                            ggTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAgc           <  1:2453903/71‑1 (MQ=255)

                                                            ggTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga                         >  1:1329552/1‑57 (MQ=255)

                                                             gTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAgc           <  1:3826773/70‑1 (MQ=255)

                                                              ttGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACaa                                    <  1:3386802/44‑1 (MQ=38)

                                                              ttGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACaa                                    <  1:3841108/44‑1 (MQ=38)

                                                              ttGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAgcgc         <  1:340183/71‑1 (MQ=255)

                                                               tgtgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaaa                                   >  1:5132991/1‑44 (MQ=255)

                                                               tgtgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAgcgc         <  1:4702643/70‑1 (MQ=255)

                                                                 tgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                           <  1:954753/50‑1 (MQ=255)

                                                                  gATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                           <  1:44393/49‑1 (MQ=255)

                                                                  gATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                           <  1:1982257/49‑1 (MQ=255)

                                                                   aTAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGt                                 <  1:3648858/42‑1 (MQ=38)

                                                                   aTAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACag                            >  1:3007442/1‑47 (MQ=255)

                                                                   aTAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAg            >  1:3306303/1‑63 (MQ=255)

                                                                    tAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                           <  1:4723629/47‑1 (MQ=255)

                                                                    tAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                           <  1:3409390/47‑1 (MQ=255)

                                                                    tAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccag   <  1:2410854/71‑1 (MQ=255)

                                                                     aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa                              >  1:4359505/1‑43 (MQ=255)

                                                                     aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc  >  1:3469632/1‑71 (MQ=255)

                                                                     aaaTGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagccagc  >  1:3462899/1‑71 (MQ=255)

TCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGC  >  Exported/2595208‑2595347

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