BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,309,484	T→C	K199E (AAA→GAA)	yjcS ←	predicted alkyl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,309,484	0	T	C	100.0%	208.2 / NA	60	K199E (AAA→GAA)	yjcS	predicted alkyl sulfatase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (29/31); total (29/31)

CAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAT  >  W3110S.gb/4309416‑4309550

cAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCaa                                                                  <  1:1250065/71‑1 (MQ=255)

cAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCaa                                                                  <  1:1753957/71‑1 (MQ=255)

 aGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCaa                                                                  >  1:3301716/1‑70 (MQ=255)

  gCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACa                                                                <  1:2476451/71‑1 (MQ=255)

  gCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACa                                                                <  1:4508737/71‑1 (MQ=255)

    aCGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATc                                                              >  1:2842276/1‑71 (MQ=255)

    aCGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATc                                                              >  1:2320796/1‑71 (MQ=255)

         ttCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCtt                                                         >  1:4027828/1‑71 (MQ=255)

                aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:1035576/71‑1 (MQ=255)

                aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:2386492/71‑1 (MQ=255)

                aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:1344792/71‑1 (MQ=255)

                aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:1259112/71‑1 (MQ=255)

                aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:1121054/71‑1 (MQ=255)

                aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:80835/71‑1 (MQ=255)

                 tGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGatat                                                  <  1:2066094/70‑1 (MQ=255)

                   gCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAAt                                               >  1:1263194/1‑71 (MQ=255)

                   gCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAAt                                               >  1:3759854/1‑71 (MQ=255)

                       cGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGc                                             <  1:1089766/69‑1 (MQ=255)

                         tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:499513/62‑1 (MQ=255)

                         tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCttt                                         <  1:1443515/71‑1 (MQ=255)

                         tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCttt                                         <  1:3593345/71‑1 (MQ=255)

                         tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCttt                                         <  1:3414421/71‑1 (MQ=255)

                         tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCttt                                         <  1:1986306/71‑1 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  >  1:2944276/1‑56 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTc                                        >  1:1280232/1‑66 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGc                                    >  1:557820/1‑70 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGc                                    >  1:281677/1‑70 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGc                                    >  1:4308401/1‑70 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGc                                    >  1:2945705/1‑70 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGc                                    >  1:1356585/1‑70 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGc                                    >  1:3896227/1‑70 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGc                                    >  1:2487251/1‑70 (MQ=255)

                               aaGCCTGCAGGGGCAATCACCTGAACTTTCCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGc                                    >  1:3857258/1‑70 (MQ=255)

                                  ccTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATg                                              <  1:3004025/57‑1 (MQ=255)

                                   cTGCAGGGGCAATCACCTGAACTTTGCCGGATTCaa                                                                  <  1:1032841/36‑1 (MQ=255)

                                   cTGCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGccacca                               <  1:2775803/71‑1 (MQ=255)

                                     gCAGGGGCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:1571574/50‑1 (MQ=255)

                                         gggCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGAt                                                  <  1:692473/46‑1 (MQ=255)

                                         gggCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTg                          <  1:1503637/70‑1 (MQ=255)

                                          ggCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGt                        >  1:3188970/1‑71 (MQ=255)

                                           gCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTc                       >  1:4268244/1‑71 (MQ=255)

                                           gCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTc                       >  1:4151170/1‑71 (MQ=255)

                                           gCAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACAATAGTGGTc                       >  1:2980199/1‑71 (MQ=255)

                                            cAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATaa                                                >  1:1075705/1‑45 (MQ=255)

                                            cAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATaa                                                >  1:3946751/1‑45 (MQ=255)

                                            cAATCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATaa                                                >  1:712172/1‑45 (MQ=255)

                                            cAATCACCTGAACTTTGCCGGATTCAACATCGGCTTAAGAGATaa                                                >  1:510626/1‑45 (MQ=255)

                                             aaTCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCgg                     <  1:756233/71‑1 (MQ=255)

                                             aaTCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCgg                     <  1:2731414/71‑1 (MQ=255)

                                              aTCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGt                           <  1:2635994/64‑1 (MQ=255)

                                               tCACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGgt                    <  1:1207562/70‑1 (MQ=255)

                                                cACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATaa                                                >  1:1782854/1‑41 (MQ=255)

                                                cACCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATa                             >  1:678065/1‑60 (MQ=255)

                                                 aCCTGAACTTTGCCGTATACAACATCGGCTACAGAGATAATGCCTTTCACGCCACCATAg                            >  1:2662858/1‑60 (MQ=255)

                                                 aCCTGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGtgtg                 >  1:366459/1‑71 (MQ=255)

                                                    tGAACTTTGCCGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCt              >  1:2672756/1‑71 (MQ=255)

                                                             ccGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTa     <  1:4681006/71‑1 (MQ=255)

                                                             ccGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTa     <  1:3353868/71‑1 (MQ=255)

                                                             ccGGATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTa     <  1:888814/71‑1 (MQ=255)

                                                                gATTCAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAt  <  1:3146230/71‑1 (MQ=255)

CAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAT  >  W3110S.gb/4309416‑4309550

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