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

Only 100 of 121 total aligned reads displayed.

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

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

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

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

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

 aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           <  1:827852/50‑1 (MQ=255)

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

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

   aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         <  1:2799132/50‑1 (MQ=255)

   aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         <  1:6232872/50‑1 (MQ=255)

   aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         <  1:191525/50‑1 (MQ=255)

    aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        >  1:1747031/1‑50 (MQ=255)

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

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

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

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

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

     aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       >  1:5135483/1‑50 (MQ=255)

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

      ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      >  1:4344899/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

               tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              >  1:4929713/1‑49 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                     aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        <  1:1069479/49‑1 (MQ=255)

                     aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        <  1:1273214/49‑1 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       <  1:14928/49‑1 (MQ=255)

                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       <  1:3583258/49‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                               aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:3533373/50‑1 (MQ=255)

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

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

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

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

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

                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:2876996/50‑1 (MQ=255)

                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:5800725/50‑1 (MQ=255)

                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:4707487/50‑1 (MQ=255)

                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          >  1:3529358/1‑49 (MQ=255)

                                    aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa        <  1:2939530/50‑1 (MQ=255)

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

                                       gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg      <  1:2562837/49‑1 (MQ=255)

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

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

                                        tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    <  1:4109807/50‑1 (MQ=255)

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

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

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

                                          tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  <  1:3992097/50‑1 (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` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`