BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,824,368	G→A	100%	pseudogene (501/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,368	0	G	A	100.0%	80.4 / NA	28	pseudogene (501/687 nt)	rph	ECK3633:JW3618:b3643; ribonuclease PH
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (8/20); total (8/20)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACG  >  W3110S.gb/3824308‑3824432

cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                        >  1:1567859/1‑71 (MQ=255)

cGGCAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAg                                                        >  1:759081/1‑71 (MQ=255)

   cAAGCTGAAAACCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGcaa                                                     <  1:2081560/71‑3 (MQ=255)

           aaaCCAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGAt                                             <  1:2394451/71‑1 (MQ=255)

               cAACCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTgg                                         <  1:739129/71‑1 (MQ=255)

                      aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                      <  1:1374735/51‑2 (MQ=255)

                      aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                      <  1:887256/51‑2 (MQ=255)

                      aTGAAAGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                      <  1:424417/51‑2 (MQ=255)

                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGca                                                      <  1:2493548/48‑2 (MQ=255)

                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:98981/70‑1 (MQ=255)

                         aaaGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:2329738/70‑1 (MQ=255)

                           aGGGATGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATATGCGATCTGGAATACGTTGAAg                             <  1:1683518/71‑1 (MQ=255)

                            gggATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                <  1:630153/35‑1 (MQ=37)

                            gggATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                <  1:622924/35‑1 (MQ=37)

                            gggATGGTAGCCGCAGTTTCTGTCGGTATTGTaaa                                                                <  1:165997/35‑1 (MQ=37)

                                tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:1885686/63‑1 (MQ=255)

                                tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:1758700/63‑1 (MQ=255)

                                tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTg                                <  1:2628388/63‑1 (MQ=255)

                                tGGTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActca                        <  1:2489363/71‑2 (MQ=255)

                                 ggTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGActca                        <  1:887380/70‑2 (MQ=255)

                                  gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              <  1:1245086/63‑1 (MQ=255)

                                  gTAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              <  1:1756343/63‑1 (MQ=255)

                                  gTAGCC‑CAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGaa                              <  1:1088448/62‑1 (MQ=255)

                                   tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                   >  1:1718350/1‑57 (MQ=255)

                                   tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                   >  1:1611925/1‑57 (MQ=255)

                                   tAGCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACg                                   >  1:1329422/1‑57 (MQ=255)

                                     gCCGCAGTTTCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCg                    <  1:2147498/70‑1 (MQ=255)

                                            tttCTGTCTGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCag                 <  1:2460481/66‑1 (MQ=255)

                                              tCTGTCGGTATTGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCagag               >  1:1864132/1‑66 (MQ=255)

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

                                                        ttGTAAACGGTGAAGCAATTTGCGATCTGGAATACGTTGAAGACTCAGCCGCAGAGACCGACATGAACg  >  1:773132/1‑69 (MQ=255)

CGGCAAGCTGAAAACCAATCCGATGAAAGGGATGGTAGCCGCAGTTTCTGTCGGAATTGTGAACGGCGAAGCGGTTTGCGATCTGGAATACGTTGAAGACTCTGCCGCAGAGACCGACATGAACG  >  W3110S.gb/3824308‑3824432

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