BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,595,260	G→A	12.9%	C12C (TGC→TGT)	arsC ←	arsenate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,595,260	0	G	A	12.9%	318.4 / 29.5	124	C12C (TGC→TGT)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base G (55/53); new base A (4/12); total (59/65)
Fisher's exact test for biased strand distribution p-value = 6.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

CCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACG  >  Exported/2595190‑2595328

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

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

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

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

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

                          ttgCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCgg                                                              <  1:1573904/52‑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:4748671/71‑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:2512346/51‑1 (MQ=255)

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

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

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

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

                                     gCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGc                                       >  1:613852/1‑65 (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:136605/52‑1 (MQ=255)

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

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

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

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

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

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

                                          gcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATag                                                    <  1:1504319/46‑2 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                 cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtct                         <  1:2660411/67‑1 (MQ=255)

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

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

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

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

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

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

                                                     gTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACaa                 <  1:1465342/71‑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:5307268/1‑64 (MQ=255)

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

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

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

                                                          aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                       >  1:2854054/1‑60 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                  tGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCat                                   <  1:5025386/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)

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

CCAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACG  >  Exported/2595190‑2595328

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