BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	827,597	C→T	P26P (CCG→CCA)	ybhS ←	predicted transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	827,597	0	C	T	100.0%	186.9 / NA	58	P26P (CCG→CCA)	ybhS	predicted transporter subunit
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (34/24); total (34/24)

GTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGG  >  W3110S.gb/827540‑827660

gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGa                                                    >  1:347230/1‑71 (MQ=255)

gTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGa                                                    >  1:825842/1‑71 (MQ=255)

 tAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGa                                                              <  1:739322/60‑1 (MQ=255)

 tAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGa                                                              <  1:2227265/60‑1 (MQ=255)

 tAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGa                                                              <  1:3015348/60‑1 (MQ=255)

 tAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGa                                                              <  1:2148746/60‑1 (MQ=255)

 tAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGa                                                              <  1:2735025/60‑1 (MQ=255)

 tAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGa                                                              <  1:2156615/60‑1 (MQ=255)

   aCCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCt                                                 >  1:2530928/1‑71 (MQ=255)

   aCCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCt                                                 >  1:2695609/1‑71 (MQ=255)

   aCCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCt                                                 >  1:1957976/1‑71 (MQ=255)

   aCCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCt                                                 >  1:1232792/1‑71 (MQ=255)

                    gTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACg                                <  1:4032529/71‑1 (MQ=255)

                    gTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACg                                <  1:698299/71‑1 (MQ=255)

                       gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGa                                                        <  1:3589836/44‑1 (MQ=255)

                       gcagcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGca                              >  1:1827415/1‑70 (MQ=255)

                         agcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAAc                                                      >  1:2457498/1‑44 (MQ=255)

                         agcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAAc                                                      >  1:384793/1‑44 (MQ=255)

                         agcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAAc                                                      >  1:1886203/1‑44 (MQ=255)

                         agcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAAc                                                      >  1:3702428/1‑44 (MQ=255)

                         agcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAAc                                                      >  1:1694416/1‑44 (MQ=255)

                         agcGGGATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAAc                                                      >  1:3019900/1‑44 (MQ=255)

                            ccgATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAgcg                         <  1:3413156/68‑1 (MQ=255)

                            ccgATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAgcg                         <  1:3987182/68‑1 (MQ=255)

                            ccgATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAgcg                         <  1:1830782/68‑1 (MQ=255)

                            ccgATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAgcg                         <  1:1402336/68‑1 (MQ=255)

                            ccgATCACTACCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAgcg                         <  1:1171170/68‑1 (MQ=255)

                                   cTACCGCAATCAGCCAGCTACTTGGATCGCGAACg                                                     <  1:200885/35‑1 (MQ=255)

                                     aCCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                >  1:1968858/1‑70 (MQ=255)

                                     aCCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                >  1:1756130/1‑70 (MQ=255)

                                     aCCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                >  1:1547193/1‑70 (MQ=255)

                                     aCCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                >  1:2940001/1‑70 (MQ=255)

                                     aCCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                >  1:3520801/1‑70 (MQ=255)

                                     aCCGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                >  1:3118933/1‑70 (MQ=255)

                                      cctcAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                <  1:2241699/66‑1 (MQ=255)

                                      ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                <  1:377667/69‑1 (MQ=255)

                                      ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                <  1:718534/69‑1 (MQ=255)

                                      ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                <  1:2754638/69‑1 (MQ=255)

                                      ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                <  1:1217346/69‑1 (MQ=255)

                                      ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTAcgc                <  1:174303/69‑1 (MQ=255)

                                      ccGCAATCAGCCAGCTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCCCCCGTAcgc                <  1:2496256/69‑1 (MQ=255)

                                             cagccagcTACTTGGATCGCGAACGATCTGCCGCGtc                                         >  1:3297456/1‑37 (MQ=255)

                                             cagccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcag        >  1:3235277/1‑70 (MQ=255)

                                             cagccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcag        >  1:2792771/1‑70 (MQ=255)

                                             cagccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCca         >  1:1384106/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:3881924/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:920344/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:1489307/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:160267/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:2342743/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:1075469/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:3532475/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:2949721/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:3195591/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcagga      >  1:3303737/1‑69 (MQ=255)

                                                ccagcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGAAGCcagga      >  1:2160880/1‑69 (MQ=255)

                                                  agcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcaggacagg  <  1:4074071/71‑1 (MQ=255)

                                                  agcTACTTGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCcaggacagg  <  1:729553/71‑1 (MQ=255)

GTAACCAAAAATAAACAGCAGTAGCAGCGGGATCACTACCGCAATCAGCCAGCTACTCGGATCGCGAACGATCTGCCGCGTCTCTTTAACGCACAGCGCCCGTACGCGACGCCAGGACAGG  >  W3110S.gb/827540‑827660

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