BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,438,124	T→C	W140R (TGG→CGG)	ytfI →	protein YtfI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,438,124	0	T	C	100.0%	162.8 / NA	42	W140R (TGG→CGG)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (22/20); total (22/20)

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             >  1:3884512/1‑49 (MQ=255)

cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             <  1:128926/49‑1 (MQ=255)

  aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:349755/50‑1 (MQ=255)

  aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           >  1:2415456/1‑49 (MQ=255)

   aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         >  1:1287158/1‑50 (MQ=255)

     aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       <  1:1083200/50‑1 (MQ=255)

        gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:1415197/50‑1 (MQ=255)

            aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                                <  1:857985/50‑1 (MQ=255)

             tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               <  1:3602725/50‑1 (MQ=255)

              ttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               >  1:2438766/1‑49 (MQ=255)

                acGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:2139054/2‑49 (MQ=255)

                 cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           <  1:394754/50‑1 (MQ=255)

                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:1393292/1‑49 (MQ=255)

                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:3419789/1‑49 (MQ=255)

                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                          <  1:3790739/50‑1 (MQ=255)

                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                          <  1:3352964/50‑1 (MQ=255)

                   ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:1087402/1‑48 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:3286608/50‑1 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      >  1:1673445/1‑50 (MQ=255)

                       aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:2233584/50‑1 (MQ=255)

                       aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:370464/50‑1 (MQ=255)

                        aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                    <  1:263010/50‑1 (MQ=255)

                         cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAATAAACAGTCTGa                   >  1:3741316/1‑50 (MQ=255)

                          aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                  <  1:48456/50‑1 (MQ=255)

                           aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 >  1:2250632/1‑50 (MQ=255)

                            aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 >  1:2069063/1‑49 (MQ=255)

                            aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt                <  1:388641/50‑1 (MQ=255)

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               >  1:306210/1‑50 (MQ=255)

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               >  1:1487592/1‑50 (MQ=255)

                              cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              <  1:503027/50‑1 (MQ=255)

                              cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              <  1:158597/50‑1 (MQ=255)

                                gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:3276423/1‑49 (MQ=255)

                                 aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa           <  1:1133229/50‑1 (MQ=255)

                                  cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa           >  1:3871226/1‑49 (MQ=255)

                                  cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa           >  1:2556899/1‑49 (MQ=255)

                                      tGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:1478674/1‑49 (MQ=255)

                                       gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa     >  1:597581/1‑50 (MQ=255)

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:3069074/1‑50 (MQ=255)

                                         aTGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:297050/1‑49 (MQ=255)

                                          tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:1308420/49‑1 (MQ=255)

                                          tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  <  1:1733955/50‑1 (MQ=255)

                                           gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  >  1:1614186/1‑49 (MQ=255)

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

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