BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,737,305	G→A	47.4%	L231L (CTG→TTG)	lepB ←	leader peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,737,305	0	G	A	47.4%	5.5 / 66.9	78	L231L (CTG→TTG)	lepB	leader peptidase
Reads supporting (aligned to +/- strand): ref base G (17/24); new base A (21/16); total (38/40)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

ATACATCCCCACCTGATCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCGG  >  minE/1737243‑1737359

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:2495385/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:385233/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:1782817/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:3120121/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:3024292/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:1637313/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:1616610/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:735801/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:1573570/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:2815860/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:445390/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:615513/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:250913/1‑69 (MQ=255)

aTACATCCCCACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGtctc                                                  >  1:1104774/1‑69 (MQ=255)

aTACATCCCCACCTGATCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAGTGtctc                                                  >  1:347521/1‑69 (MQ=255)

aTACATCCCCACCTGATCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAGTGtctc                                                  >  1:1983900/1‑69 (MQ=255)

aTACATCCCCACCTGATCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAGTGtctc                                                  >  1:3184844/1‑69 (MQ=255)

aTACATCCCCACCTGATCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAGTGtctc                                                  >  1:290297/1‑69 (MQ=255)

      ccccACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTAc                                             >  1:379492/1‑68 (MQ=255)

      ccccACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACg                                            >  1:1668875/1‑69 (MQ=255)

      ccccACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACg                                            >  1:924429/1‑69 (MQ=255)

      ccccACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACg                                            >  1:346798/1‑69 (MQ=255)

      ccccACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACg                                            >  1:1248670/1‑69 (MQ=255)

      ccccACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACg                                            >  1:2378971/1‑69 (MQ=255)

      ccccACCTGGTCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACg                                            >  1:329260/1‑69 (MQ=255)

       cccACCTGATCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGc                                           >  1:3196787/1‑69 (MQ=255)

       cccACCTGATCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGc                                           >  1:3201205/1‑69 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGTATTCCATTTTCt                   <  1:2547566/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:513543/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:2747011/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:2540080/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:2175400/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:2089700/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:2005511/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:965582/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:1839792/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:928739/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:1237861/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:1378203/69‑1 (MQ=255)

                               aCTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:773844/69‑1 (MQ=255)

                                cTGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:2901284/68‑1 (MQ=255)

                                 tGTCAGAATGCGGTGCGTCACATCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:457109/67‑1 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTctct                                                 >  1:1830205/1‑36 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                <  1:1814366/53‑1 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                <  1:998610/53‑1 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                <  1:3246826/53‑1 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                <  1:3338765/53‑1 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                >  1:1674789/1‑53 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                <  1:665622/53‑1 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                >  1:658257/1‑53 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                <  1:3114887/53‑1 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                <  1:1015941/53‑1 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACg                                >  1:256803/1‑53 (MQ=255)

                                  gTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTg                >  1:1869485/1‑69 (MQ=255)

                                     aGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCAttt                      >  1:2905233/1‑60 (MQ=255)

                                     aGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGtt             >  1:2845490/1‑69 (MQ=255)

                                     aGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGtt             >  1:2433344/1‑69 (MQ=255)

                                     aGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGtt             >  1:207883/1‑69 (MQ=255)

                                        aTGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCg          >  1:1209910/1‑69 (MQ=255)

                                        aTGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCg          >  1:1800668/1‑69 (MQ=255)

                                              tGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTc    <  1:1750367/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:1035541/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:1204423/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:873417/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:853050/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:750575/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:1609132/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:534930/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:2264610/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:513431/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:2668275/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:3115260/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:3128340/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:3149242/69‑1 (MQ=255)

                                                cGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:3156750/69‑1 (MQ=255)

                                                      aTCACCCAATGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCt                   <  1:1543014/46‑1 (MQ=255)

                                                        cACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:3173782/61‑1 (MQ=255)

                                                         aCCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCgg  <  1:1526124/60‑1 (MQ=255)

ATACATCCCCACCTGATCCTGCGCAATCGGCACTGTCAGAATGCGGTGCGTCACATCACCCAGTGTCTCTTTACGCTCGGAAAGACGAATTCCATTTTCTTTGGTTTCGTTTTTCGG  >  minE/1737243‑1737359

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