BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,824,326	T→C	pseudogene (459/687 nt)	rph →	ECK3633:JW3618:b3643; ribonuclease PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,824,326	0	T	C	100.0%	223.3 / NA	62	pseudogene (459/687 nt)	rph	ECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (6/56); total (6/56)

AGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGG  >  W3110S.gb/3824281‑3824373

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:1492926/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:1804242/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:1739901/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:1640704/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:1598840/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:1593854/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:2353426/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:2380515/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:1509149/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:1478457/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:2458587/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:2793008/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:684346/62‑1 (MQ=255)

aGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATgg                                 <  1:892392/62‑1 (MQ=255)

            gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1084173/62‑1 (MQ=255)

            gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:293703/62‑1 (MQ=255)

            gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:57782/62‑1 (MQ=255)

            gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:841734/62‑1 (MQ=255)

            gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1369352/62‑1 (MQ=255)

            gAAGCTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1513242/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2595186/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2755722/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2928261/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2422074/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:1413041/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:500716/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2367135/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2525702/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2598703/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2203948/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:809722/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2162397/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2613026/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:1945893/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2642037/62‑1 (MQ=255)

                cTGGTGGAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:1745705/62‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:2754791/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1071072/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:3008582/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1308971/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1265930/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:803717/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1151551/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:845907/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1749689/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:2457428/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1782057/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:2108633/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:2318730/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:2328334/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:2367467/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:2426083/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:2711827/53‑1 (MQ=255)

                     ggAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:1345394/57‑1 (MQ=255)

                      gAAAACGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGt                 <  1:2892539/56‑1 (MQ=255)

                       aaaaCGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGttt                     <  1:1576550/51‑1 (MQ=255)

                            ggCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCgg              >  1:1889670/1‑53 (MQ=255)

                                aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg  >  1:2965032/1‑61 (MQ=255)

                                aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg  >  1:2396765/1‑61 (MQ=255)

                                aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg  >  1:318934/1‑61 (MQ=255)

                                aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg  >  1:2180372/1‑61 (MQ=255)

                                aGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACgg  >  1:2451466/1‑61 (MQ=255)

AGATGCGCTACAGAAGCTGGTGGAAAACGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGG  >  W3110S.gb/3824281‑3824373

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