BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	39,505	T→C	29.2%	A35A (GCA→GCG)	caiT ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	39,505	0	T	C	29.2%	137.3 / 66.2	140	A35A (GCA→GCG)	caiT	predicted transporter
Reads supporting (aligned to +/- strand): ref base T (59/40); new base C (10/31); total (69/71)
Fisher's exact test for biased strand distribution p-value = 1.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

CAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGGAAAGAAAAC  >  minE/39444‑39570

cAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCa                                                                >  1:1738517/1‑65 (MQ=255)

cAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCa                                                                >  1:1602710/1‑65 (MQ=255)

cAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCa                                                                >  1:2344014/1‑65 (MQ=255)

cAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCa                                                                >  1:692996/1‑65 (MQ=255)

cAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCa                                                                >  1:2192942/1‑65 (MQ=255)

cAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCa                                                                >  1:932241/1‑65 (MQ=255)

cAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCa                                                            <  1:466198/69‑1 (MQ=255)

cAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCa                                                            <  1:3392477/69‑1 (MQ=255)

            cccATACATTGGTGACGTAACTGAATACAGCATTAATAACGCCATTCGCTGCATCCAGATCTCTGACt                                                 >  1:2321189/1‑68 (MQ=255)

            cccATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACt                                                 >  1:1309821/1‑68 (MQ=255)

            cccATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACt                                                 >  1:1422885/1‑68 (MQ=255)

            cccATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACt                                                 >  1:865316/1‑68 (MQ=255)

            cccATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACt                                                 >  1:1241962/1‑68 (MQ=255)

            cccATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACt                                                 >  1:2090482/1‑68 (MQ=255)

            cccATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACt                                                 >  1:1064167/1‑68 (MQ=255)

                  cATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGcc                                          >  1:73413/1‑69 (MQ=255)

                  cATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGcc                                          >  1:2559710/1‑69 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:1593311/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:2886091/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:3400943/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:2163188/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:325359/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:2099485/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:1653909/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:556766/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:1557159/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:2721007/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:2497986/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:1127037/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:1095698/69‑1 (MQ=255)

                         gACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:965995/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:509473/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:2813597/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:237127/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:2612765/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:2859051/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:2777171/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:2527478/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:254089/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:935779/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:3259647/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:1330938/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:358271/69‑1 (MQ=255)

                                         gCATTAATAACGACATTCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGattatt                   <  1:3425367/69‑1 (MQ=255)

                                          cATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  <  1:2009333/69‑1 (MQ=255)

                                          cATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  <  1:1286919/69‑1 (MQ=255)

                                          cATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  <  1:1505933/69‑1 (MQ=255)

                                          cATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  <  1:566558/69‑1 (MQ=255)

                                          cATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  <  1:1846540/69‑1 (MQ=255)

                                          cATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  <  1:2949911/69‑1 (MQ=255)

                                           aTTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  <  1:1181608/68‑1 (MQ=255)

                                            ttAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACaaat                                   <  1:1826623/50‑2 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  >  1:1174923/1‑65 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:32637/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:3088229/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:3108305/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:2715841/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:2464520/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:100792/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:1232960/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:175776/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:2277375/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcgg              >  1:2133053/1‑69 (MQ=255)

                                              aataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAAAAAAGTATGCCGACGATTATTAAcgg              >  1:1611729/1‑69 (MQ=255)

                                               ataaCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAg                                   <  1:826601/47‑1 (MQ=255)

                                                  aCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcggcg            >  1:3012213/1‑67 (MQ=255)

                                                  aCGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAAcggc             >  1:2206602/1‑66 (MQ=255)

                                                    gACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGAttatta                  <  1:1012594/59‑1 (MQ=255)

                                                     aCATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaag       <  1:1968543/69‑1 (MQ=255)

                                                     aCATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaag       <  1:656023/69‑1 (MQ=255)

                                                      cATTCGCCGCATCCAGATCTCTGACTGTAAGCcaacaa                                     >  1:3330196/1‑38 (MQ=255)

                                                      cATTCGCCGCATCCAGATCTCTGACTGTAAGCcaacaa                                     >  1:1863721/1‑38 (MQ=255)

                                                      cATTCGCCGCATCCAGATCTCTGACTGTAAGCcaacaa                                     >  1:2974500/1‑38 (MQ=255)

                                                      cATTCGCCGCATCCAGATCTCTGACTGTAAGCcaacaa                                     >  1:420819/1‑38 (MQ=255)

                                                      cATTCGCCGCATCCAGATCTCTGACTGTAAGCcaacaa                                     >  1:667789/1‑38 (MQ=255)

                                                      cATTCGCCGCATCCAGATCTCTGACTGTAAGCcaaca                                      >  1:2631564/1‑37 (MQ=255)

                                                      cATTCGCCGCATCCAGATATCTGACTGTAAGCcaacaa                                     >  1:2283258/1‑38 (MQ=38)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGTCGgaaagaaa    <  1:178766/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgcaagaaa    <  1:1507/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    <  1:1988910/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    <  1:2069450/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    <  1:3219101/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    <  1:1371756/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    <  1:1735399/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    <  1:1680195/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    <  1:117333/68‑1 (MQ=255)

                                                         tCGCCGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    <  1:2482836/68‑1 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACaaa                                    >  1:3386917/1‑35 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTa                                 >  1:2343871/1‑38 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTa                                 >  1:2849409/1‑38 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTa                                 >  1:1523156/1‑38 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGCCGATTATTAACGGCGGAAAGAAAAc  >  1:1602979/1‑69 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    >  1:762375/1‑67 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    >  1:2685628/1‑67 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    >  1:189573/1‑67 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    >  1:3196698/1‑67 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    >  1:2382281/1‑67 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    >  1:2804133/1‑67 (MQ=255)

                                                          cGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGgaaagaaa    >  1:2577545/1‑67 (MQ=255)

                                                          cGCTGCATCCAGATCTATGACTGTAAGCCAAAAAAGTATGCCGACGATTATTAACGGCGGAAAGAAAAc  >  1:3339461/1‑69 (MQ=255)

CAAATGCCCATCCCCATACATTGGTGACGTAACTGAATACAGCATTAATAACGACATTCGCTGCATCCAGATCTCTGACTGTAAGCCAACAAAGTATGCCGACGATTATTAACGGCGGAAAGAAAAC  >  minE/39444‑39570

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