BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,459,166	G→A	A53V (GCT→GTT)	secE ←	preprotein translocase membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,459,166	0	G	A	100.0%	178.5 / NA	57	A53V (GCT→GTT)	secE	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (26/31); total (26/31)

GCGCTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAGCAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAATGTCGCGATAAAGATAGTTGCCGACA  >  W3110S.gb/3459098‑3459234

gcgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTgcaaca                                                                    <  1:221410/71‑1 (MQ=255)

gcgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTgcaaca                                                                    <  1:3473882/71‑1 (MQ=255)

 cgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGcaacaa                                                                   <  1:1036695/71‑1 (MQ=255)

 cgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGcaacaa                                                                   <  1:152252/71‑1 (MQ=255)

 cgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGcaacaa                                                                   <  1:4013989/71‑1 (MQ=255)

 cgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGcaacaa                                                                   <  1:3887108/71‑1 (MQ=255)

 cgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGcaacaa                                                                   <  1:3840430/71‑1 (MQ=255)

 cgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGcaacaa                                                                   <  1:3710881/71‑1 (MQ=255)

 cgcTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGcaacaa                                                                   <  1:3807939/71‑1 (MQ=255)

       aCGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGa                                                              <  1:2315759/70‑1 (MQ=255)

       aCGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGa                                                              <  1:1963386/70‑1 (MQ=255)

       aCGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGa                                                              <  1:3550054/70‑1 (MQ=255)

       aCGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGa                                                              <  1:1656310/70‑1 (MQ=255)

       aCGGGAAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGa                                                              <  1:2958535/70‑1 (MQ=255)

          ggCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATc                                                                 >  1:2286082/1‑64 (MQ=255)

          ggCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATta                                                          >  1:165086/1‑71 (MQ=255)

          ggCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATta                                                          >  1:3342076/1‑71 (MQ=255)

          ggCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATta                                                          >  1:820088/1‑71 (MQ=255)

          ggCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGAGACACCCCCCGCTGCAACAATCAGAATta                                                          >  1:181729/1‑71 (MQ=255)

               aaGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGa                                                              <  1:4020530/62‑1 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATc                                                                 >  1:1290176/1‑57 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:3352596/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:1692413/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:1655173/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:1882710/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:3659755/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:2015647/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:3460587/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:3022416/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:2763502/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:2558912/1‑71 (MQ=255)

                 gCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGc                                                   >  1:67200/1‑71 (MQ=255)

                    aCGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGcaacaa                                                                   <  1:2627427/52‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCTCGACACCACCCGCTGCAACAATCAg                                                               <  1:3035351/55‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCGACAATCAg                                                               <  1:580044/55‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAg                                                               <  1:985921/55‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAg                                                               <  1:2932052/55‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAg                                                               <  1:3807362/55‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAg                                                               <  1:2412992/55‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAg                                                               <  1:1349016/55‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAg                                                               <  1:1551799/55‑1 (MQ=255)

                     cGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAg                                                               <  1:544012/55‑1 (MQ=255)

                       gTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGGTGCAACAATCAGaa                                                             >  1:3432752/1‑55 (MQ=255)

                       gTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGaa                                                             >  1:3957098/1‑55 (MQ=255)

                       gTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGaa                                                             >  1:2726669/1‑55 (MQ=255)

                       gTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGaa                                                             >  1:3936348/1‑55 (MQ=255)

                        tAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGcgca                                            <  1:1885803/71‑1 (MQ=255)

                                 cTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGgc                                   <  1:3692777/71‑1 (MQ=255)

                                 cTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGgc                                   <  1:4041538/71‑1 (MQ=255)

                                 cTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGgc                                   <  1:3819676/71‑1 (MQ=255)

                                 cTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGgc                                   <  1:2749390/71‑1 (MQ=255)

                                       tCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCAt                              <  1:2437342/70‑1 (MQ=255)

                                       tCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCAt                              <  1:3531890/70‑1 (MQ=255)

                                       tCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCAt                              <  1:3971159/70‑1 (MQ=255)

                                       tCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCAt                              <  1:978972/70‑1 (MQ=255)

                                       tCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCAt                              <  1:20304/70‑1 (MQ=255)

                                       tCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCAt                              <  1:264674/70‑1 (MQ=255)

                                       tCGTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCAt                              <  1:2268683/70‑1 (MQ=255)

                                         gTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAAt                           >  1:1820922/1‑71 (MQ=255)

                                         gTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAAt                           >  1:863559/1‑71 (MQ=255)

                                         gTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAAt                           >  1:710766/1‑71 (MQ=255)

                                         gTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAAt                           >  1:3173895/1‑71 (MQ=255)

                                         gTTAACAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAAt                           >  1:1437847/1‑71 (MQ=255)

                                              cAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAATGt                         <  1:602402/68‑1 (MQ=255)

                                              cAGCGCGACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAATGt                         <  1:498209/68‑1 (MQ=255)

                                                   cgACACCACCCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAATGt                         <  1:4139630/63‑1 (MQ=255)

                                                        ccaccCGCTGCAACAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCAt                              <  1:3385061/53‑1 (MQ=255)

                                                                  cagcaATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAATGTCGCGATAAAGATAGTTGCCGACa  <  1:3680370/71‑1 (MQ=255)

GCGCTTCACGGGCAAAAGCAACGGTAGCTTTACCTTTTGTCGTTAACAGCGCGACACCACCCGCTGCAGCAATCAGAATTACTACGGCCAGCGCACGCAGCGGCAGCATAATGTCGCGATAAAGATAGTTGCCGACA  >  W3110S.gb/3459098‑3459234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`