BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,824,410	T→A	pseudogene (543/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,410	0	T	A	100.0%	165.9 / NA	55	pseudogene (543/687 nt)	rph	ECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (46/9); total (46/9)

GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAC  >  W3110S.gb/3824342‑3824475

gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                 >  1:939555/1‑71 (MQ=255)

gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                 >  1:560493/1‑71 (MQ=255)

gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                 >  1:3849542/1‑71 (MQ=255)

gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                 >  1:1435717/1‑71 (MQ=255)

gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGc                                                                 >  1:345405/1‑71 (MQ=255)

                  ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                <  1:707680/67‑1 (MQ=255)

                  ggtaTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGa                                                <  1:702508/67‑1 (MQ=255)

                    taatgtaAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAAc                                              >  1:3130143/4‑70 (MQ=255)

                    taTTGTAAACGGTGAAGGAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:431391/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:4127433/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:3294788/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:3286862/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:1567349/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:307375/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:270667/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:258543/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:2489307/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:2453036/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:705347/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:2091871/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:1974322/2‑71 (MQ=255)

                    taTTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:3998004/2‑71 (MQ=255)

                    taTTGTAAACGG‑TAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:744917/2‑70 (MQ=255)

                        gtaaaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGt                                         >  1:3785973/4‑71 (MQ=255)

                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGAc                                                        >  1:2721455/1‑53 (MQ=255)

                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                  >  1:373845/1‑59 (MQ=255)

                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                  >  1:1574108/1‑59 (MQ=255)

                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                  >  1:2732509/1‑59 (MQ=255)

                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                  >  1:1317123/1‑59 (MQ=255)

                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                  >  1:1230876/1‑59 (MQ=255)

                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                  >  1:402324/1‑59 (MQ=255)

                           aaCGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACAt                                                  >  1:1274545/1‑59 (MQ=255)

                             cGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACGTAGtgatga                                    >  1:4081447/1‑71 (MQ=255)

                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:2274125/51‑1 (MQ=255)

                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:2146383/51‑1 (MQ=255)

                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:1949722/51‑1 (MQ=255)

                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:1818830/51‑1 (MQ=255)

                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:1772586/51‑1 (MQ=255)

                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:1704346/51‑1 (MQ=255)

                                       atttGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             <  1:4065418/51‑1 (MQ=255)

                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:367385/1‑35 (MQ=255)

                                                        cGTTGAAGACTCAGCCGCAGAGACCGACATGAACg                                             >  1:1806610/1‑35 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGa                             >  1:3477937/1‑45 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:62195/1‑57 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:438077/1‑57 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:1174936/1‑57 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:4027878/1‑57 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:3820655/1‑57 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:3409897/1‑57 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:3230101/1‑57 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:2118869/1‑57 (MQ=255)

                                                              aGACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGcatcat                 >  1:1992469/1‑57 (MQ=255)

                                                               gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAc  >  1:3553171/1‑71 (MQ=255)

                                                               gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAc  >  1:1130053/1‑71 (MQ=255)

                                                               gACTCAGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAc  >  1:1183264/1‑71 (MQ=255)

                                                                 ctcaGCCGCAGAGACCGACATGAACGTAGTGATGACCGaa                               <  1:2314859/40‑1 (MQ=255)

GTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACGTAGTGATGACCGAAGACGGGCGCATCATTGAAGTGCAGGGGAC  >  W3110S.gb/3824342‑3824475

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