BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,458,605	G→A	52.2%	pseudogene (3350/4605 nt) pseudogene (3350/4104 nt)	yfaS ← yfaS ←	ECK2220:JW2222+JW2221:b4500; hypothetical protein ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,458,605	0	G	A	52.2%	‑5.1 / 16.9	23	pseudogene (3350/4605 nt) pseudogene (3350/4104 nt)	yfaS yfaS	ECK2220:JW2222+JW2221:b4500; hypothetical protein ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment
Reads supporting (aligned to +/- strand): ref base G (7/4); new base A (5/7); total (12/11)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTT  >  minE/1458546‑1458669

atgcccatgccTTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                                                         <  1:691474/69‑1 (MQ=255)

atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                                                         <  1:1637416/69‑1 (MQ=255)

atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                                                         <  1:1945024/69‑1 (MQ=255)

atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                                                         <  1:1987051/69‑1 (MQ=255)

atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                                                         <  1:2869734/69‑1 (MQ=255)

atgcccatgccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGcc                                                         <  1:2704278/69‑1 (MQ=255)

atgcccatgccGTAAGGAAGGCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGc                                                          <  1:687114/68‑1 (MQ=255)

atgcccataccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTACGGGcc                                                         <  1:1308139/69‑1 (MQ=255)

        gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCgg                                                            >  1:2196791/1‑58 (MQ=255)

        gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGc                                                          >  1:1214899/1‑60 (MQ=255)

        gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGc                                                          >  1:3434234/1‑60 (MQ=255)

        gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGc                                                          >  1:522200/1‑60 (MQ=255)

        gccGTGAGGAAGGCGTCGCCATTGCCATCTTCACCCCACCAGGTAAAGCGCACTCCGGGc                                                          >  1:1326567/1‑60 (MQ=255)

                         aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAg                                   >  1:2226407/1‑66 (MQ=255)

                         aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACg                                >  1:922468/1‑69 (MQ=255)

                         aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACg                                >  1:2273580/1‑69 (MQ=255)

                         aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACg                                >  1:1133598/1‑69 (MQ=255)

                         aCCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACg                                >  1:3085351/1‑69 (MQ=255)

                                     ttCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAat                    >  1:511329/1‑69 (MQ=255)

                                     ttCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAat                    >  1:937614/1‑69 (MQ=255)

                                         ccccACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAatca                  <  1:456730/67‑1 (MQ=255)

                                         ccccACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAatca                  <  1:3056700/67‑1 (MQ=255)

                                              ccaGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACg           <  1:408929/69‑1 (MQ=255)

                                                       gcgcgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGtt  >  1:3175150/1‑69 (MQ=255)

                                                       gcgcgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGtt  <  1:3142969/69‑1 (MQ=255)

                                                       gcgcgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGtt  >  1:2903472/1‑69 (MQ=255)

                                                       gcgcgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGtt  <  1:908381/69‑1 (MQ=255)

                                                       gcgcgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGtt  >  1:2399353/1‑69 (MQ=255)

                                                       gcgcgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGtt  >  1:1786403/1‑69 (MQ=255)

                                                       gcgcgcTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGt   >  1:901431/1‑68 (MQ=255)

ATGCCCATGCCGTAAGGAAGGCGTCACCATTGCCATCTTCACCCCACCAGGTAAAGCGCGCTCCGGGCCCCGCCAGTTGCATCAGCCGCAGACGGTTATCCTGAATCATCTGACGAATGTCGTT  >  minE/1458546‑1458669

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