BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,813,364	T→C	V250A (GTC→GCC)	emrB →	multidrug efflux pump membrane subunit EmrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,813,364	0	T	C	100.0%	194.2 / NA	50	V250A (GTC→GCC)	emrB	multidrug efflux pump membrane subunit EmrB
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (27/23); total (27/23)

TCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGTT  >  NC_000913/2813319‑2813409

tcatcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCt                                             <  1:2808678/48‑1 (MQ=255)

  atcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                           >  1:1047069/1‑48 (MQ=255)

   tcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                           >  1:650841/1‑47 (MQ=255)

   tcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                           >  1:1175433/1‑47 (MQ=255)

   tcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                           >  1:2782285/1‑47 (MQ=255)

   tcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTgg                                           >  1:2653374/1‑47 (MQ=255)

   tcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGa                                         <  1:726139/49‑1 (MQ=255)

   tcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGa                                         <  1:2490226/49‑1 (MQ=255)

     cTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAg                                        <  1:1249470/48‑1 (MQ=255)

     cTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAg                                        <  1:1424899/48‑1 (MQ=255)

     cTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAg                                        <  1:247133/48‑1 (MQ=255)

     cTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGc                                       >  1:25847/1‑49 (MQ=255)

        aCCGTGGTGGCCGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGa                                    >  1:1718957/1‑49 (MQ=39)

         ccGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAc                                   >  1:904987/1‑49 (MQ=255)

         ccGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAc                                   >  1:1779210/1‑49 (MQ=255)

          cGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAcc                                  >  1:1521137/1‑49 (MQ=255)

          cGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGAcc                                  >  1:2366333/1‑49 (MQ=255)

                ggCGGTGGTGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACcgacga                             >  1:2174942/1‑48 (MQ=255)

                   ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAc                         <  1:573783/49‑1 (MQ=255)

                   ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAc                         >  1:382035/1‑49 (MQ=255)

                   ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAc                         >  1:290974/1‑49 (MQ=255)

                   ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAc                         >  1:1371676/1‑49 (MQ=255)

                   ggtggtGGCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAc                         >  1:1651562/1‑49 (MQ=255)

                     tggtggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAAccc                       >  1:1784771/1‑49 (MQ=255)

                       gtggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCg                      <  1:1392697/48‑1 (MQ=255)

                         ggCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGa                     >  1:2452586/1‑47 (MQ=255)

                          gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATa                   <  1:747986/48‑1 (MQ=255)

                          gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATa                   <  1:358179/48‑1 (MQ=255)

                          gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAg                  <  1:1942010/49‑1 (MQ=255)

                          gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAg                  <  1:2101355/49‑1 (MQ=255)

                          gCTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAg                  <  1:2132940/49‑1 (MQ=255)

                           cTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGt                 <  1:1072922/49‑1 (MQ=255)

                           cTATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGt                 <  1:577970/49‑1 (MQ=255)

                            tATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTc                >  1:1677865/1‑49 (MQ=255)

                            tATCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTc                >  1:1393582/1‑49 (MQ=255)

                             aTCTGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTc                >  1:208588/1‑48 (MQ=255)

                                tGCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATc            >  1:1044466/1‑49 (MQ=255)

                                 gCTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCt           >  1:783113/1‑49 (MQ=255)

                                  cTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGAt             <  1:2159245/46‑1 (MQ=255)

                                  cTTCCTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCt           <  1:1022205/48‑1 (MQ=255)

                                     ccTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGt         <  1:2849923/47‑1 (MQ=255)

                                     ccTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCg       >  1:1657407/1‑49 (MQ=255)

                                      cTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgt      >  1:844929/1‑49 (MQ=255)

                                      cTGATTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCg       <  1:2567260/48‑1 (MQ=255)

                                        gATTGCCTGGGATCTGACCGACGATAACCCGATAGTCGATCTGTCgttg    >  1:2132878/1‑49 (MQ=38)

                                         aTTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgttgt   <  1:1149822/49‑1 (MQ=255)

                                         aTTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgttgt   <  1:1898212/49‑1 (MQ=255)

                                         aTTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgttgt   >  1:794970/1‑49 (MQ=255)

                                         aTTGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgttg    <  1:395053/48‑1 (MQ=255)

                                          ttGCCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCgttgtt  <  1:2432856/49‑1 (MQ=255)

TCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGTT  >  NC_000913/2813319‑2813409

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