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

AAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438079‑4438169

aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                            <  1:314720/49‑1 (MQ=255)

aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           >  1:4095249/1‑50 (MQ=255)

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

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

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

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

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

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

      cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      >  1:598733/1‑49 (MQ=255)

       gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:3851646/1‑50 (MQ=255)

          aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:2368285/1‑47 (MQ=255)

          aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 <  1:1051014/50‑1 (MQ=255)

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

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

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

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

               ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             <  1:2421317/49‑1 (MQ=255)

                cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:3844267/1‑50 (MQ=255)

                cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:2067475/1‑50 (MQ=255)

                cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:1591316/1‑50 (MQ=255)

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

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

                  ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACa                         <  1:2404923/50‑1 (MQ=255)

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

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

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

                        cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:856727/1‑50 (MQ=255)

                        cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   <  1:1914624/50‑1 (MQ=255)

                         aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:194845/1‑49 (MQ=255)

                         aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:3970457/1‑49 (MQ=255)

                          aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 <  1:2658474/50‑1 (MQ=255)

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

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

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

                             cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              >  1:4030478/1‑50 (MQ=255)

                              aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:2934641/1‑50 (MQ=255)

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

                               gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:3836980/50‑1 (MQ=255)

                               gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:1536319/50‑1 (MQ=255)

                                aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            >  1:2002110/1‑49 (MQ=255)

                                 cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          <  1:3687818/50‑1 (MQ=255)

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

                                    ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       <  1:3603597/50‑1 (MQ=255)

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

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

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

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

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

AAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438079‑4438169

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