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%	139.9 / NA	35	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 (20/15); total (20/15)

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCG  >  NC_000913/3090235‑3090327

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

cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                             >  1:2543581/1‑50 (MQ=255)

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

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

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

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

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

          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                   >  1:395442/1‑50 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg      <  1:2518275/49‑1 (MQ=255)

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

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

                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc   >  1:1101060/1‑50 (MQ=255)

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

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCG  >  NC_000913/3090235‑3090327

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