BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,260	0	G	A	24.5%	198.7 / 65.7	110	C12C (TGC→TGT)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base G (46/37); new base A (10/17); total (56/54)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 116 total aligned reads displayed.

CAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACG  >  minE/2595191‑2595328

cAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                                                     <  1:545457/71‑1 (MQ=255)

 aGATAATTAATTATTGTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCcaca                                                                    <  1:3064865/71‑1 (MQ=255)

 aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                                                      >  1:325666/1‑69 (MQ=255)

 aGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                                                      >  1:2546880/1‑69 (MQ=255)

   aTAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAgg                                                                  >  1:31737/1‑71 (MQ=255)

    tAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGc                                                                 >  1:2706900/1‑71 (MQ=255)

     aaTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCa                                                                    >  1:2572014/1‑67 (MQ=255)

      aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGGGACGTGCCGCAGGcc                                                                <  1:1000847/70‑1 (MQ=255)

      aTGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCg                                                                      <  1:3754353/64‑1 (MQ=255)

       tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCgg                                                              >  1:696829/1‑71 (MQ=255)

       tGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCgg                                                              >  1:693838/1‑71 (MQ=255)

          ataataGTCGGTTCTGTGCCGCTGTTGCTGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                           <  1:1677810/71‑1 (MQ=255)

          ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                           <  1:817088/71‑1 (MQ=255)

          ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                           <  1:1060286/71‑1 (MQ=255)

          ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGtt                                                           <  1:1867950/71‑1 (MQ=255)

          ataataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGt                                                            >  1:3857410/1‑70 (MQ=255)

            aataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgt                                                         >  1:1556364/1‑71 (MQ=255)

             attgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgt                                                         <  1:2587251/67‑1 (MQ=255)

             ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                                                        <  1:3801900/71‑1 (MQ=255)

             ataGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                                                        <  1:676431/71‑1 (MQ=255)

              ttgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGa                                                       <  1:672465/69‑1 (MQ=255)

              taGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGc                                                                     >  1:2458016/1‑57 (MQ=255)

               tgTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgtg                                                        >  1:3721245/2‑69 (MQ=255)

                gTCGGTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTtgtg                                                        >  1:217582/1‑68 (MQ=255)

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

                  cGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                                    >  1:1503227/1‑70 (MQ=255)

                   ggTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaaa                                                   <  1:3376541/70‑1 (MQ=255)

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

                   ggTTCGTTGCCGTTGTTATGAAGCATCTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGAt                                                  >  1:3355037/1‑71 (MQ=255)

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

                          tgCCGCTGTTGCGTATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTtgtg                                                        <  1:1139418/58‑1 (MQ=255)

                          tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                           >  1:2599270/1‑71 (MQ=255)

                            ccGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                             >  1:3530304/1‑67 (MQ=255)

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

                                  ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATg                                                 <  1:29745/57‑1 (MQ=255)

                                  ttGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTT‑CTCat                                   <  1:2393005/70‑1 (MQ=255)

                                   tGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCt                                      >  1:1855616/1‑67 (MQ=255)

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

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

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

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

                                     cGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                                    <  1:1213714/51‑1 (MQ=255)

                                      ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTaa                                             <  1:2820987/57‑1 (MQ=255)

                                      ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCat                                   <  1:3618691/67‑1 (MQ=255)

                                       gATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                              <  1:2906178/71‑1 (MQ=255)

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

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

                                        agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATc                                <  1:471410/66‑1 (MQ=255)

                                        agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAg                              >  1:3036614/3‑70 (MQ=255)

                                        agcatcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGt                             <  1:2259287/69‑1 (MQ=255)

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

                                          catcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa                            >  1:2241728/1‑70 (MQ=255)

                                           atcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGt                                          >  1:3055752/1‑55 (MQ=255)

                                           atcTCCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGt                                          >  1:2378487/1‑55 (MQ=255)

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

                                            tctcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGtt                                                           <  1:1749654/37‑1 (MQ=16)

                                            tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTAATATCAGTAtct                         >  1:3822883/1‑71 (MQ=255)

                                            tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGT‑ATGTTGCTCATATCAGTAtct                         >  1:2076000/1‑70 (MQ=255)

                                              tcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGAt                                                  <  1:1714810/44‑1 (MQ=255)

                                              tcCAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGAt                                                  <  1:201980/44‑1 (MQ=255)

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

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

                                               ccAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCAt                      >  1:740956/1‑71 (MQ=255)

                                                cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATCGTAATGTTGCTCATATCAGTAtct                         <  1:3723317/67‑1 (MQ=255)

                                                cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtctc                        <  1:1560621/68‑1 (MQ=255)

                                                cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtctc                        <  1:2656478/68‑1 (MQ=255)

                                                cAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                       <  1:3209063/69‑1 (MQ=255)

                                                  gCGTGTTGCGTGAGGTGCTACAGGCCGGGTTGTGATAGATGGTAATGTTGATCATATCAGTATCTAATTa                    >  1:3424739/1‑70 (MQ=255)

                                                  gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtctaa                       >  1:1868750/1‑65 (MQ=255)

                                                  gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtctaa                       >  1:307983/1‑65 (MQ=255)

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

                                                  gCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTAc                   <  1:2753213/71‑1 (MQ=255)

                                                  gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCatat                                 <  1:2495028/57‑1 (MQ=255)

                                                  gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTa                            >  1:2126527/1‑62 (MQ=255)

                                                  gCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTa                    >  1:1590420/1‑70 (MQ=255)

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

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

                                                         aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg             <  1:754711/70‑1 (MQ=255)

                                                          cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTg             <  1:1221495/69‑1 (MQ=255)

                                                          cgcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGa            >  1:168189/1‑70 (MQ=255)

                                                           gcgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAg                              >  1:275117/1‑51 (MQ=255)

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

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

                                                              aggTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACaga        <  1:2437832/68‑1 (MQ=255)

                                                               ggTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTg             <  1:529118/63‑1 (MQ=255)

                                                               ggTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATAAGTATCTAATTACAAAGTg             <  1:813768/63‑1 (MQ=255)

                                                               cGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga      <  1:2427159/71‑1 (MQ=255)

                                                               cGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga        >  1:697413/1‑69 (MQ=255)

                                                               cGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga        >  1:1999142/1‑69 (MQ=255)

                                                                 tGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATc                                <  1:3841555/43‑1 (MQ=38)

                                                                 tGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATc                                <  1:2302452/43‑1 (MQ=38)

                                                                 tGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACa                  <  1:3883971/57‑1 (MQ=255)

                                                                  gCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGAc   >  1:655747/1‑71 (MQ=255)

                                                                  gCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGAc   >  1:2170751/1‑71 (MQ=255)

                                                                   ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACg  >  1:2712208/1‑71 (MQ=255)

                                                                   ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACg  >  1:2716847/1‑71 (MQ=255)

                                                                   ccGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACg  >  1:1950270/1‑71 (MQ=255)

                                                                    cGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCa                       <  1:949235/49‑1 (MQ=255)

                                                                     gCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt                             <  1:1489750/42‑1 (MQ=255)

                                                                     gCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt                             <  1:1409828/42‑1 (MQ=255)

CAGATAATGGATAATAGTCGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACG  >  minE/2595191‑2595328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`