BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,595,278	0	A	G	29.4%	142.1 / 77.9	102	I6I (ATT→ATC)	arsC	arsenate reductase
Reads supporting (aligned to +/- strand): ref base A (42/30); new base G (13/17); total (55/47)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

CGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCG  >  minE/2595209‑2595348

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)

        tgCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATg                                                               >  1:2599270/1‑71 (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:2881427/1‑67 (MQ=255)

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

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

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

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

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

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

                    ggATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATaa                                                                        <  1:3132729/50‑1 (MQ=255)

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

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

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

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

                      atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt                                                 <  1:1984321/71‑1 (MQ=255)

                      atcatcTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGt                                                 <  1:3667590/71‑1 (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:1636529/1‑71 (MQ=255)

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

                          tctcCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTAtct                                             >  1:2415416/1‑71 (MQ=255)

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

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

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

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

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

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

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

                              cAGCGTGTTGCGTGAGGTGCCACAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTAtctaa                                           <  1:3380061/69‑3 (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)

                                       aCGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTa                                                                  >  1:2242104/1‑37 (MQ=255)

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

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

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

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

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

                                          cgACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAag                             <  1:1810013/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)

                                               tGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAaga                            <  1:716150/67‑1 (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:2170751/1‑71 (MQ=255)

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

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

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

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

                                                  cacAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAg                                                  >  1:3812191/3‑42 (MQ=39)

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

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

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

                                                   acAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCgcc                    <  1:1770789/70‑3 (MQ=255)

                                                    cAGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACaga                            >  1:3912182/1‑62 (MQ=255)

                                                     aGGCCGGGTTGTGATAGATGGTAATGTTGCTCATATCa                                                   <  1:3493372/38‑1 (MQ=255)

                                                     aGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAgc                  <  1:3577484/71‑1 (MQ=255)

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

                                                        ccGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAagaga                          >  1:1028909/1‑60 (MQ=255)

                                                           ggTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAgc            >  1:3680107/1‑71 (MQ=255)

                                                           ggTTGTGATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAg             >  1:2377166/1‑70 (MQ=255)

                                                            gTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                              <  1:1473824/52‑1 (MQ=255)

                                                            gTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGaaa                              <  1:1481245/52‑1 (MQ=255)

                                                            gTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAgcg           <  1:3544600/71‑1 (MQ=255)

                                                            gTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAgcg           <  1:523028/71‑1 (MQ=255)

                                                              tgtgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAgcg           >  1:1493478/1‑69 (MQ=255)

                                                               gtgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcag        >  1:2165543/1‑71 (MQ=255)

                                                                tgATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTGACAGAGAGCCGCCACGCCAGCGcggc       <  1:3812469/71‑1 (MQ=255)

                                                                tgATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAg             >  1:88660/1‑65 (MQ=255)

                                                                 gATAGATGGTAATGTTGCTCATATCAGTATCTAATTACAAAGTg                                 <  1:297475/44‑1 (MQ=38)

                                                                  aTAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGcagcc      <  1:2661785/70‑1 (MQ=255)

                                                                     aaTGGTAATGTTGCTCATAACAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCg  <  1:2762614/71‑1 (MQ=255)

CGGTTCTGTGCCGCTGTTGCGGATCATCTCCAGCGTATTACGCGACGTGCCGCAGGCCGGGTTGTGATAAATGGTAATGTTGCTCATATCAGTATCTCATTACAAAGTGAAAGAGAGACGTAGCGCCAGCGCAGCCAGCG  >  minE/2595209‑2595348

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