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%	426.2 / NA	102	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 (51/51); total (51/51)

Only 100 of 102 total aligned reads displayed.

AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090236‑3090328

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

 gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            <  1:486438/50‑1 (MQ=255)

 gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            >  1:7034042/1‑50 (MQ=255)

 gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            <  1:3656169/50‑1 (MQ=255)

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

  cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  1:1626225/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                   aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  1:4237000/1‑48 (MQ=255)

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

                    tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         >  1:3254836/1‑50 (MQ=255)

                    tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         >  1:4693995/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

                        aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         <  1:4063350/46‑1 (MQ=255)

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

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

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

                         aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  1:3172170/50‑1 (MQ=255)

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

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

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

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

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

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

                           tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:3666064/1‑50 (MQ=255)

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

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

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

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

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

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

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

                              cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  1:1328613/1‑50 (MQ=255)

                              cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  1:6381674/1‑50 (MQ=255)

                               aGGAACTTTCCTGATTTGGCTGATTGCATGCCAAAACAGCTTTCGCTAc               <  1:1740923/49‑1 (MQ=39)

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

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

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

                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              <  1:6473728/50‑1 (MQ=255)

                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              <  1:1446331/50‑1 (MQ=255)

                               aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              >  1:4884308/1‑50 (MQ=255)

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

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

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

                                  aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  1:4881424/50‑1 (MQ=255)

                                  aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  1:4958172/50‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                        ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt     >  1:2471094/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090236‑3090328

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