BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,090,347	Δ708 bp	coding (1‑708/708 nt)	endA →	DNA‑specific endonuclease I

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3090347	3091054	708	36 [0]	[0] 36	endA	DNA‑specific endonuclease I

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3090346	0 (0.000)	36 (1.310)	26/96	0.0	100%	intergenic (+94/‑1)	yggI/endA	conserved protein YggI/DNA‑specific endonuclease I
*	?	NC_000913	3091055 =	0 (0.000)	36 (1.310)	26/96	0.0	100%	intergenic (+1/‑79)	endA/rsmE	DNA‑specific endonuclease I/16S rRNA m(3)U1498 methyltransferase

AATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAA  >  NC_000913/3091007‑3091102

caaaacagctttcgctactttgctggctcgttttaacacggagaagtgc                                                 <  1:1743047‑M2/1‑1 (MQ=255)

   aacagctttcgctacgttgctggctcgttttaacacggagaagtgCCTa                                              <  1:957286‑M2/4‑1 (MQ=255)

   aacagctttcgctacgttgctggctcgttttaacacggagaagtgCCTa                                              <  1:2113948‑M2/4‑1 (MQ=255)

   aacagctttcgctacgttgctggctcgttttaacacggagaagtgCCTa                                              <  1:1985117‑M2/4‑1 (MQ=255)

     cagctttcgctacgttgctggctcgttttaacacggagaagtgCCTaca                                            <  1:1039729‑M2/6‑1 (MQ=255)

     cagctttcgctacgttgctggctcgttttaacacggagaagtgCCTaca                                            <  1:699752‑M2/6‑1 (MQ=255)

      agctttcgctacgttgctggctcgttttaacacggagaagtgCCTacac                                           >  1:2176‑M2/43‑49 (MQ=255)

        ctttcgctacgttgctggctcgttttaacacggagaagtgCCTacac                                           >  1:1570082‑M2/41‑47 (MQ=255)

        ctttcgctacgttgctggctcgttttaacacggagaagtgCCTACACTa                                         <  1:992565‑M2/9‑1 (MQ=255)

          ttcgctacgttgctggctcgttttaacacggagaagtgCCTACACTAGc                                       >  1:1735878‑M2/39‑49 (MQ=255)

           tcgctacgttgctggctcgttttaacacggagaagtgCCTACACTAg                                        <  1:493121‑M2/10‑1 (MQ=255)

             gctacgttgctggctcgttttaacacggagaagtgCCTACACTAGCgg                                     >  1:668432‑M2/36‑48 (MQ=255)

             gctacgttgctggctcgttttaacacggagaagtgCCTACACTAGCgg                                     >  1:597282‑M2/36‑48 (MQ=255)

                 cgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGAtt                                 <  1:1296617‑M2/17‑1 (MQ=255)

                 cgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGAt                                  <  1:210612‑M2/16‑1 (MQ=255)

                 cgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTc                                <  1:1076035‑M2/18‑1 (MQ=255)

                 cgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTc                                <  1:733179‑M2/18‑1 (MQ=255)

                 cgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTc                                <  1:950029‑M2/18‑1 (MQ=255)

                 ccttgctggctctttttaacacggagaagtgCCTACACTAGCGGGAtt                                 <  1:1361757‑M2/17‑1 (MQ=255)

                  gttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCt                               >  1:647540‑M2/31‑49 (MQ=255)

                  gttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCt                               <  1:719889‑M2/19‑1 (MQ=255)

                      ctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCttttt                           <  1:1409014‑M2/23‑1 (MQ=255)

                       tggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTg                          >  1:563547‑M2/26‑49 (MQ=255)

                         gctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGtt                        <  1:928520‑M2/26‑1 (MQ=255)

                          ctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGtt                        >  1:442718‑M2/23‑48 (MQ=255)

                          ctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTa                       <  1:141120‑M2/27‑1 (MQ=255)

                           tcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTaa                      >  1:59552‑M2/22‑49 (MQ=255)

                             gttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAc                     <  1:1615100‑M2/29‑1 (MQ=255)

                                ttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAcccc                  <  1:878560‑M2/32‑1 (MQ=255)

                                ttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAcccc                  <  1:528218‑M2/32‑1 (MQ=255)

                                ttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAcccc                  <  1:279183‑M2/32‑1 (MQ=255)

                                 taacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCt                 <  1:478413‑M2/33‑1 (MQ=255)

                                   acacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTAc               <  1:1400220‑M2/35‑1 (MQ=255)

                                         agaagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAcgc        >  1:174617‑M2/8‑49 (MQ=255)

                                             gtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTAc    <  1:735908‑M2/46‑1 (MQ=255)

                                               gCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACaa  >  1:2020228‑M2/2‑49 (MQ=255)

AATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAA  >  NC_000913/3091007‑3091102

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