BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,554,675	G→A	30.9%	S403S (TCG→TCA)	dsdX →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,554,675	0	G	A	30.9%	116.8 / 58.6	123	S403S (TCG→TCA)	dsdX	predicted transporter
Reads supporting (aligned to +/- strand): ref base G (28/56); new base A (18/20); total (46/77)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 123 total aligned reads displayed.

AGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATACTATAC  >  minE/1554610‑1554737

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

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

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

aGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCtt                                                             >  1:1713863/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:2281478/1‑68 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                    gCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc                         <  1:1071691/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:1052547/39‑1 (MQ=37)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                   aTCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                           <  1:1293324/36‑1 (MQ=255)

                                                   aTCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                           <  1:2475221/36‑1 (MQ=255)

                                                   aTCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCaa                                           <  1:488546/36‑1 (MQ=255)

                                                    tCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAAt         >  1:2280981/1‑69 (MQ=255)

                                                    tCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAAt         >  1:3352781/1‑69 (MQ=255)

                                                    tCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAAt         >  1:3274087/1‑69 (MQ=255)

                                                    tCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAAt         >  1:1327541/1‑69 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTaaa          <  1:3050678/66‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTaaa          <  1:1817485/66‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTaaa          <  1:2446928/66‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTaaa          <  1:823617/66‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTaaa          <  1:1118792/66‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTaaa          <  1:2384946/66‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATAc       <  1:3074268/69‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATAc       <  1:2961740/69‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATAc       <  1:2784239/69‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATAc       <  1:1521750/69‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATAc       <  1:1542047/69‑1 (MQ=255)

                                                      gTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGAGCTACGCTCAATGAAACATTTAAATAc       <  1:1099012/69‑1 (MQ=255)

                                                       ttACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATAc       <  1:1648238/68‑1 (MQ=255)

                                                         aCGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTc                         <  1:883154/48‑1 (MQ=255)

                                                         aCGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTa            >  1:2752607/1‑61 (MQ=255)

                                                         aCGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTa            >  1:2446453/1‑61 (MQ=255)

                                                         aCGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTa            >  1:1288112/1‑61 (MQ=255)

                                                         aCGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTa            >  1:1148778/1‑61 (MQ=255)

                                                         aCGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTa            >  1:1072980/1‑61 (MQ=255)

                                                         aCGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTa            >  1:1048618/1‑61 (MQ=255)

                                                            gACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATACTATAc  <  1:784809/68‑1 (MQ=255)

                                                              cTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATAc       <  1:1583915/61‑1 (MQ=255)

AGCCCGGAAATTATTGCGATTGCTATCGGTTCAGGTGCAATTGGCTGCACTATCGTTACGGACTCGCTTTTCTGGCTAGTGAAGCAATATTGCGGCGCTACGCTCAATGAAACATTTAAATACTATAC  >  minE/1554610‑1554737

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