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%	122.8 / NA	34	W140R (TGG→CGG)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (16/18); total (16/18)

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

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

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

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

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

   aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCACGGGata                                         <  1:1612929/50‑1 (MQ=255)

    aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        <  1:1443012/50‑1 (MQ=255)

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

      ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      <  1:1137702/50‑1 (MQ=255)

       cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:2427054/50‑1 (MQ=255)

         tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:2650756/1‑49 (MQ=255)

         tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                   <  1:2367218/50‑1 (MQ=255)

          gAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:556143/1‑48 (MQ=255)

           aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                  <  1:1000910/49‑1 (MQ=255)

            aTTTCCGTTAGAACAAAACGGACGATTGTATGCCCGGGATATTCTACAGa                                <  1:546836/50‑1 (MQ=39)

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

               tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              <  1:2133449/49‑1 (MQ=255)

                ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:1671359/1‑49 (MQ=255)

                 cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaa                            <  1:125784/49‑1 (MQ=255)

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

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

                     aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       >  1:2108198/1‑50 (MQ=255)

                       aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     >  1:2185591/1‑50 (MQ=255)

                        aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     >  1:1496670/1‑49 (MQ=255)

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

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

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

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               <  1:762291/50‑1 (MQ=255)

                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               <  1:1397061/50‑1 (MQ=255)

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

                                    aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          >  1:435761/1‑48 (MQ=255)

                                    aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          >  1:2683092/1‑48 (MQ=255)

                                     ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:2643322/1‑50 (MQ=255)

                                          tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   >  1:864675/1‑49 (MQ=255)

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

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

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