BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,595,251	C→T	15.8%	S15S (TCG→TCA)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,595,251	0	C	T	15.8%	299.2 / 39.3	126	S15S (TCG→TCA)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base C (53/53); new base T (4/16); total (57/69)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 137 total aligned reads displayed.

GCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAA  >  Exported/2595181‑2595320

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

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

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

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

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

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

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

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

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

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

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

                aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCa                                                            >  1:4934723/1‑66 (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)

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

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

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

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

                            cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa                                           >  1:951161/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)

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

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

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

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

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

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

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

                                  tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCg                                                       <  1:5660910/53‑1 (MQ=255)

                                  tgtgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCg                                                       <  1:2059901/53‑1 (MQ=255)

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

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

                                    tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGAt                                              <  1:166696/60‑1 (MQ=255)

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

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

                                    tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                   <  1:5096028/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:2512346/51‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                  agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                                 <  1:4023910/41‑1 (MQ=255)

                                                  agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                                 <  1:2933133/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:3673613/1‑49 (MQ=255)

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

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

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

                                                        tcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAt                   >  1:4740944/1‑67 (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:1323531/1‑70 (MQ=255)

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

                                                          cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt              >  1:557149/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)

                                                            gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa               >  1:5641697/1‑67 (MQ=255)

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

                                                             cGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGTTAATGTTGCTCATATCAGTATCTCATTAc           >  1:2334989/1‑70 (MQ=255)

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

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

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

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

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

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

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

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

                                                                   aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGt                                       <  1:316472/36‑1 (MQ=255)

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

                                                                    cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa   <  1:2640712/71‑1 (MQ=255)

                                                                    cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaa   <  1:4988950/71‑1 (MQ=255)

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

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

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

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

GCGGAGTTTCCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAA  >  Exported/2595181‑2595320

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