BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,554,660	T→G	37.3%	T398T (ACT→ACG)	dsdX →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,554,660	0	T	G	37.3%	57.5 / 91.6	102	T398T (ACT→ACG)	dsdX	predicted transporter
Reads supporting (aligned to +/- strand): ref base T (30/34); new base G (18/20); total (48/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 103 total aligned reads displayed.

GACATCAGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTA  >  minE/1554604‑1554727

gACATCAGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGa                                                          <  1:2151577/68‑1 (MQ=255)

      aGCCCGGAAATTATTGCGATTGGTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCtt                                                   >  1:1632413/1‑69 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGTTGCACTATCGTTACGGACTCGCtt                                                   >  1:1497731/1‑69 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCtt                                                   >  1:1713863/1‑69 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCtt                                                   >  1:550662/1‑69 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCtt                                                   >  1:1999635/1‑69 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCtt                                                   >  1:3435118/1‑69 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCt                                                    >  1:2490925/1‑68 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCt                                                    >  1:203081/1‑68 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCt                                                    >  1:2004980/1‑68 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCt                                                    >  1:515684/1‑68 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCt                                                    >  1:2281478/1‑68 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCt                                                    >  1:1970117/1‑68 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGCACTCGCtt                                                   >  1:1349186/1‑69 (MQ=255)

      aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTAAGGACTCGCt                                                    >  1:2531684/1‑68 (MQ=255)

                 tattGCGATTGCTATCGGTTCGGGTGCAATTGGCTGCACTATCg                                                                 >  1:3270152/1‑44 (MQ=255)

                 tattGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:1859320/1‑45 (MQ=255)

                 tattGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:2954070/1‑45 (MQ=255)

                 tattGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:1707621/1‑45 (MQ=255)

                 tattGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:1917921/1‑45 (MQ=255)

                 tattGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:3049228/1‑45 (MQ=255)

                 tattGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:3396802/1‑45 (MQ=255)

                 tattGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:3399812/1‑45 (MQ=255)

                 tattGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:1227902/1‑45 (MQ=255)

                 tattGCGATTCCTATCGGTTCAGGTGCAATTGGCTGCACTATCGt                                                                >  1:2140081/1‑45 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:2409772/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:206245/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:2213186/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:2256919/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:1801004/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:1016768/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:2853553/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:1746683/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:2631686/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:467597/69‑1 (MQ=255)

                        aTTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                 <  1:364095/69‑1 (MQ=255)

                              aTCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTgc                           <  1:2802810/69‑1 (MQ=255)

                                     cAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGc                                   <  1:1095238/54‑1 (MQ=255)

                                     cAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGc                                   <  1:2163767/54‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCTCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:2287596/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAATCaa                                 <  1:2449068/51‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:3237320/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               >  1:600593/1‑69 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:2290825/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:2405390/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:3354189/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:3042647/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:290237/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               >  1:843967/1‑69 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:2556132/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:2752924/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:218779/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:2057330/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:1071691/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:1128877/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:1530467/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:1450808/69‑1 (MQ=255)

                                          gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCCCTACGCTc               <  1:1718554/69‑1 (MQ=255)

                                           cAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc               <  1:1637104/68‑1 (MQ=255)

                                                ggCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTgc                           >  1:1614906/1‑51 (MQ=255)

                                                ggCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTgc                           >  1:1911052/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:429555/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:3061859/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:1477988/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:3403662/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:13576/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:1052547/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:1079733/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:680091/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:747781/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:754875/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:2274896/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:2552921/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:22473/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:218680/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:2436296/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:2439587/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:1722091/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:2608977/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:1824050/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTg                                       <  1:1893637/39‑1 (MQ=37)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAat                                >  1:2866675/1‑46 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:2033479/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:2896123/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:1649345/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:2126157/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:2047295/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:530270/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:1103758/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:445285/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:2029497/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:3158819/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:1987261/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:3213156/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:3392480/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:3370445/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:1911218/1‑51 (MQ=255)

                                                ggCTGCACGATCGTTACGGACTCACTTTTCTGGCTGGTGAAGCAATATTgc                           >  1:1577678/1‑51 (MQ=255)

                                                       cTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTa  <  1:1137367/69‑1 (MQ=255)

                                                        tATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTa  <  1:2490814/68‑1 (MQ=255)

GACATCAGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTA  >  minE/1554604‑1554727

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