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%	160.8 / NA	40	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 (18/22); total (18/22)

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCT  >  NC_000913/3090235‑3090325

cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                           <  1:1356064/50‑1 (MQ=255)

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

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

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

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

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

         gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                   <  1:2730194/49‑1 (MQ=255)

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

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

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

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

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

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

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

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

                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                           >  1:3668409/1‑49 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt          >  1:1771478/1‑50 (MQ=255)

                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt          >  1:1894579/1‑50 (MQ=255)

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

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

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

                                      tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg     >  1:1446960/1‑50 (MQ=255)

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

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

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

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCT  >  NC_000913/3090235‑3090325

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: `‑`