BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,871,189	A→G	14.1%	P372P (CCA→CCG)	yqcE →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,871,189	0	A	G	14.1%	108.5 / 7.0	63	P372P (CCA→CCG)	yqcE	predicted transporter
Reads supporting (aligned to +/- strand): ref base A (28/26); new base G (2/7); total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCAT  >  minE/1871125‑1871257

ttATGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAAtt                                                                   <  1:1159342/68‑1 (MQ=255)

ttATGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCCGATGTCTTCGTTTACCCAAtt                                                                   <  1:301206/68‑1 (MQ=255)

  aTGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATc                                                                <  1:614691/69‑1 (MQ=255)

  aTGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATc                                                                <  1:527327/69‑1 (MQ=255)

  aTGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATc                                                                <  1:3385099/69‑1 (MQ=255)

  aTGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATc                                                                <  1:893475/69‑1 (MQ=255)

  aTGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATc                                                                <  1:3038513/69‑1 (MQ=255)

  aTGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATc                                                                <  1:897433/69‑1 (MQ=255)

  aTGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATc                                                                <  1:2088021/69‑1 (MQ=255)

     ggCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGc                                                             <  1:1806270/69‑1 (MQ=255)

              gtaggtaTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCaa                                                    <  1:331206/69‑1 (MQ=255)

               taggtaTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCaa                                                    <  1:2817573/68‑1 (MQ=255)

                   taTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               <  1:2418956/69‑1 (MQ=255)

                   taTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               <  1:3277553/69‑1 (MQ=255)

                   taTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               <  1:1192275/69‑1 (MQ=255)

                   taTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               <  1:1311258/69‑1 (MQ=255)

                   taTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               <  1:1400264/69‑1 (MQ=255)

                   taTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               <  1:3436240/69‑1 (MQ=255)

                    aTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACAcc                                              >  1:460475/1‑69 (MQ=255)

                    aTTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACAcc                                              >  1:3067338/1‑69 (MQ=255)

                       tGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCt                                            >  1:978843/1‑68 (MQ=255)

                       tGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCt                                            >  1:1944572/1‑68 (MQ=255)

                       tGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTg                                           >  1:1963812/1‑69 (MQ=255)

                       tGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTg                                           >  1:871240/1‑69 (MQ=255)

                       tGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTg                                           >  1:2663024/1‑69 (MQ=255)

                        gTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATCGGCTACTGGCAAGACACCCTGc                                          >  1:3429765/1‑69 (MQ=255)

                        gTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCGATTATCGGCTACTGGCAAGACACCCTGc                                          >  1:1345773/1‑69 (MQ=255)

                                  tGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               <  1:966576/54‑1 (MQ=255)

                                  tGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGa                                >  1:664920/1‑69 (MQ=255)

                                  tGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGa                                >  1:1650559/1‑69 (MQ=255)

                                  tGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGa                                >  1:1779045/1‑69 (MQ=255)

                                  tGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATAGGCCACTGGCAAGacac                                               <  1:16393/54‑1 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               >  1:3020600/1‑53 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               >  1:662073/1‑53 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               >  1:2696193/1‑53 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               >  1:452084/1‑53 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               >  1:1959370/1‑53 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGaca                                                >  1:1839578/1‑52 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGaca                                                >  1:1928305/1‑52 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGaca                                                >  1:289859/1‑52 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGaca                                                >  1:1313380/1‑52 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGaca                                                >  1:2110848/1‑52 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGaa                               >  1:135686/1‑69 (MQ=255)

                                   ggATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGa                                >  1:2748624/1‑68 (MQ=255)

                                     aTTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAg                              <  1:3096866/68‑1 (MQ=255)

                                      ttCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGcc                            <  1:1233904/69‑1 (MQ=255)

                                      ttCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGcc                            <  1:255237/69‑1 (MQ=255)

                                      ttCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGcc                            <  1:965567/69‑1 (MQ=255)

                                      ttCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGcc                            <  1:2213451/69‑1 (MQ=255)

                                      ttCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGcc                            <  1:473146/69‑1 (MQ=255)

                                      ttCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGcc                            <  1:1163980/69‑1 (MQ=255)

                                                 tGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGacac                                               <  1:1823398/39‑1 (MQ=255)

                                                   tCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATAtgt               >  1:1111560/1‑69 (MQ=255)

                                                   tCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATAtgt               >  1:979524/1‑69 (MQ=255)

                                                     ttCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTgg             >  1:64992/1‑69 (MQ=255)

                                                     ttCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTgg             >  1:486074/1‑69 (MQ=255)

                                                     ttCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTgg             >  1:930324/1‑69 (MQ=255)

                                                     ttCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTgg             >  1:2181047/1‑69 (MQ=255)

                                                         tttACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCt           <  1:3275400/67‑1 (MQ=255)

                                                         tttACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCt           <  1:242854/67‑1 (MQ=255)

                                                              ccAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATggg     <  1:469570/68‑1 (MQ=255)

                                                              ccAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATggg     <  1:3144903/68‑1 (MQ=255)

                                                              ccAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATggg     <  1:204763/68‑1 (MQ=255)

                                                                aaTTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCAt  <  1:2836684/69‑1 (MQ=255)

                                                                aaTTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCAt  <  1:90306/69‑1 (MQ=255)

                                                                aaTTATCGGCCACTGGCAAGACACCATGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCAt  <  1:2302637/69‑1 (MQ=255)

TTATGGGCACCACGGTAGGTATTTGTTCGGTGATTGGATTCCTGCCGGATGTCTTCGTTTACCCAATTATCGGCCACTGGCAAGACACCCTGCCCGCTGCAGAAGCCTACCGCAATATGTGGCTGATGGGCAT  >  minE/1871125‑1871257

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