BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,090,282	C→T	intergenic (+30/‑65)	yggI → / → endA	conserved protein YggI/DNA‑specific endonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,090,282	0	C	T	100.0%	461.0 / NA	113	intergenic (+30/‑65)	yggI/endA	conserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (44/69); total (44/69)

Only 100 of 113 total aligned reads displayed.

AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090236‑3090328

aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             >  1:3952112/1‑50 (MQ=255)

aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             <  1:2645076/50‑1 (MQ=255)

aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             <  1:643105/50‑1 (MQ=255)

aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             <  1:3310309/50‑1 (MQ=255)

  cTGGTTGCGAAATAACTATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  1:3786131/1‑50 (MQ=39)

  cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           <  1:2465888/50‑1 (MQ=255)

   tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  1:5125499/1‑49 (MQ=255)

    ggTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                         >  1:2927287/1‑50 (MQ=255)

     gTTGCGAAATAACCATCTGGACTATCAGGAACTTTCCTGATTTGGCTg                                          <  1:2876407/48‑1 (MQ=39)

     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                          <  1:3720145/48‑1 (MQ=255)

     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  1:3633866/50‑1 (MQ=255)

     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        >  1:2903539/1‑50 (MQ=255)

     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        >  1:254896/1‑50 (MQ=255)

     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  1:1290600/50‑1 (MQ=255)

     gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        >  1:1432798/1‑50 (MQ=255)

      ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  1:2627582/50‑1 (MQ=255)

      ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  1:2422853/50‑1 (MQ=255)

        gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      >  1:602288/1‑49 (MQ=255)

        gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      >  1:755942/1‑49 (MQ=255)

         cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                     <  1:153098/49‑1 (MQ=255)

         cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    <  1:38186/50‑1 (MQ=255)

          gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:5610575/1‑49 (MQ=255)

          gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:963899/1‑49 (MQ=255)

          gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                   <  1:3801803/50‑1 (MQ=255)

           aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  >  1:1502127/1‑50 (MQ=255)

           aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  <  1:2194715/50‑1 (MQ=255)

           aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  <  1:1561487/50‑1 (MQ=255)

           aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  <  1:4361873/50‑1 (MQ=255)

            aaTATCCATCTGAACTATGAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 <  1:2349993/50‑1 (MQ=255)

            aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  1:1394407/1‑50 (MQ=255)

            aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 <  1:4370893/50‑1 (MQ=255)

             aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                                <  1:3795919/50‑1 (MQ=255)

             aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  1:2765667/1‑49 (MQ=255)

              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                                >  1:2438043/1‑49 (MQ=255)

              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                                >  1:6057004/1‑49 (MQ=255)

              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               <  1:1222814/50‑1 (MQ=255)

              tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               <  1:1887275/50‑1 (MQ=255)

               aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGTTTGCATACCa                               >  1:5140323/1‑49 (MQ=255)

               aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               >  1:1390782/1‑49 (MQ=255)

                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                             <  1:3208891/50‑1 (MQ=255)

                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                             <  1:3799234/50‑1 (MQ=255)

                aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              <  1:1167029/49‑1 (MQ=255)

                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  1:2182956/1‑50 (MQ=255)

                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            <  1:2102853/50‑1 (MQ=255)

                 ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            <  1:4651639/50‑1 (MQ=255)

                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATGCCAAAAc                           <  1:2649198/50‑1 (MQ=255)

                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  1:3508437/1‑49 (MQ=255)

                  cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:1578649/50‑1 (MQ=255)

                   aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAcg                          <  1:3094010/50‑2 (MQ=255)

                   aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACa                          <  1:4349789/50‑1 (MQ=255)

                      tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       >  1:1753063/1‑50 (MQ=255)

                       gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      >  1:6332881/1‑50 (MQ=255)

                       gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      >  1:5314346/1‑50 (MQ=255)

                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:3952511/50‑1 (MQ=255)

                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:2780581/50‑1 (MQ=255)

                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:801792/50‑1 (MQ=255)

                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     >  1:2712901/1‑50 (MQ=255)

                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     >  1:2675760/1‑50 (MQ=255)

                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      <  1:307227/49‑1 (MQ=255)

                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       <  1:2956165/48‑1 (MQ=255)

                         aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    >  1:2573249/1‑50 (MQ=255)

                          cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                   <  1:1693986/50‑1 (MQ=255)

                          cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                   >  1:2598755/1‑50 (MQ=255)

                            aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:4033547/1‑49 (MQ=255)

                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  <  1:4396656/48‑1 (MQ=255)

                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  <  1:2580697/48‑1 (MQ=255)

                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt                 <  1:3443355/49‑1 (MQ=255)

                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  1:5480242/50‑1 (MQ=255)

                             tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                >  1:1315709/1‑50 (MQ=255)

                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCtgc               <  1:6393172/49‑3 (MQ=39)

                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               <  1:5403919/49‑1 (MQ=255)

                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               <  1:6120036/49‑1 (MQ=255)

                                ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt             <  1:1084200/50‑1 (MQ=255)

                                 gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            >  1:5163036/1‑50 (MQ=255)

                                 gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:969874/50‑1 (MQ=255)

                                 gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:472804/50‑1 (MQ=255)

                                 gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:197687/50‑1 (MQ=255)

                                 gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:4702588/50‑1 (MQ=255)

                                  aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           >  1:4013380/1‑50 (MQ=255)

                                   aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          <  1:1543658/50‑1 (MQ=255)

                                   aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          >  1:3552108/1‑50 (MQ=255)

                                    cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         >  1:3219914/1‑50 (MQ=255)

                                    cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         <  1:718536/50‑1 (MQ=255)

                                    cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         <  1:2963990/50‑1 (MQ=255)

                                    cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         <  1:3876862/50‑1 (MQ=255)

                                    cTTTCCTGATTTGGCTGATTGCATACAAAAACAGCTTTCGCTACGTTGCt         <  1:3002935/50‑1 (MQ=255)

                                      ttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg       >  1:3413692/1‑50 (MQ=255)

                                      ttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg       >  1:2227935/1‑50 (MQ=255)

                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        <  1:5164407/48‑1 (MQ=255)

                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        <  1:4075759/48‑1 (MQ=255)

                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        <  1:4188712/48‑1 (MQ=255)

                                       tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  1:5573995/1‑50 (MQ=255)

                                         cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt     <  1:1597741/49‑1 (MQ=255)

                                         cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt     <  1:4492654/49‑1 (MQ=255)

                                         cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:369145/50‑1 (MQ=255)

                                         cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:6245430/50‑1 (MQ=255)

                                         cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:166832/50‑1 (MQ=255)

                                          tGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg   >  1:2472696/1‑50 (MQ=255)

                                           gATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt  >  1:5291796/1‑50 (MQ=255)

                                           gATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt  >  1:851395/1‑50 (MQ=255)

AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090236‑3090328

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