BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,286,401	Δ3,780 bp		[yejO]–[yejO]	[yejO], insH‑8, [yejO]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2286401	2290180	3780	48 [0]	[0] 48	[yejO]–[yejO]	[yejO],insH‑8,[yejO]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2286400	0 (0.000)	48 (1.280)	33/98	0.0	100%	pseudogene (2515/2525 nt)	yejO	adhesin‑like autotransporter YejO
*	?	NC_000913	2290181 =	0 (0.000)	48 (1.280)	33/98	0.0	100%	intergenic (‑1/‑319)	yejO/narP	adhesin‑like autotransporter YejO/DNA‑binding transcriptional dual regulator NarP

GAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACA  >  NC_000913/2290132‑2290228

cagaccgtgacacatcacagcctgtttattttctgttatcagaacgtcct                                                 <  1:2317826‑M2/1‑1 (MQ=255)

cagaccgtgacacatcacagcctgtttattttctgttatcagaacgtcct                                                 <  1:860487‑M2/1‑1 (MQ=255)

cagaccgtgacacatcacagcctgtttattttctgttatcagaacgtcct                                                 <  1:617590‑M2/1‑1 (MQ=255)

cagaccgtgacacatcacagcctgtttattttctgttatcagaacgtcct                                                 <  1:1414254‑M2/1‑1 (MQ=255)

 agaccgtgacacatcacagcctgtttattttctgttatcagaacgtcct                                                 >  1:513765‑M2/49‑49 (MQ=255)

 agaccgtgacacatcacagcctgtttattttctgttatcagaacgtcct                                                 >  1:46835‑M2/49‑49 (MQ=255)

       tgacacatcacagcctgtttattttctgttatcagaacgtccTGAGCttt                                          <  1:1401634‑M2/8‑1 (MQ=255)

       tgacacatcacagcctgtttattttctgttatcagaacgtccTGAGCttt                                          <  1:1370399‑M2/8‑1 (MQ=255)

       tgacacatcacagcctgtttattttctgttatcagaacgtccTGAGCttt                                          <  1:2062815‑M2/8‑1 (MQ=255)

        gacacatcacagcctgtttattttctgttatcagaacgtccTGAGCttt                                          >  1:347343‑M2/42‑49 (MQ=255)

         acacatcacagcctgtttattttctgttatcagaacgtccTGAGCTTTc                                         <  1:295709‑M2/9‑1 (MQ=255)

         acacatcacagcctgtttattttctgttatcagaacgtccTGAGCTTTCa                                        >  1:1026428‑M2/41‑50 (MQ=255)

          cacatcacagcctgtttattttctgttatcagaacgtccTGAGCTTTCAt                                       <  1:1701520‑M2/11‑1 (MQ=255)

                acagcctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTAt                                  <  1:613290‑M2/16‑1 (MQ=255)

                acagcctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTAt                                  >  1:1063918‑M2/34‑49 (MQ=255)

                acagcctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTAt                                  >  1:1018298‑M2/34‑49 (MQ=255)

                acagcctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTAt                                  <  1:2586713‑M2/16‑1 (MQ=255)

                 cagcctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATg                                 >  1:2717792‑M2/33‑49 (MQ=255)

                 cagcctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATGa                                <  1:2164492‑M2/18‑1 (MQ=255)

                   gcctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATGaaa                              <  1:509477‑M2/20‑1 (MQ=255)

                   gcctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATGaaa                              <  1:1736025‑M2/20‑1 (MQ=255)

                    cctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaat                             <  1:1475507‑M2/21‑1 (MQ=255)

                    cctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaat                             <  1:1629325‑M2/21‑1 (MQ=255)

                    cctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaat                             <  1:344871‑M2/21‑1 (MQ=255)

                     ctgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaat                             >  1:1596107‑M2/29‑49 (MQ=255)

                      tgtttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaatt                            >  1:880319‑M2/28‑49 (MQ=255)

                        tttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaattaat                         <  1:1421656‑M2/25‑1 (MQ=255)

                        tttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaattaat                         <  1:40617‑M2/25‑1 (MQ=255)

                        tttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaattaa                          <  1:2741096‑M2/24‑1 (MQ=255)

                        tttattttctgttatcagaacgtccTGAGCTTTCATCCTATGAaatta                           <  1:2283446‑M2/23‑1 (MQ=255)

                          tattttctgttatcagaacgtccTGAGCTTTCATCCTATGAAATTAATTg                       >  1:2669894‑M2/24‑50 (MQ=255)

                           attttctgttatcagaacgtccTGAGCTTTCATCCTATGAAATTAATTg                       <  1:2749649‑M2/27‑1 (MQ=255)

                             tttctgttatcagaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTg                    >  1:2197056‑M2/21‑50 (MQ=255)

                              ttctgttatcagaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTg                    <  1:808532‑M2/30‑1 (MQ=255)

                                 ttttatcagaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTa                 >  1:361193‑M2/17‑49 (MQ=255)

                                   ttatcagaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTaaaa              <  1:576711‑M2/36‑1 (MQ=255)

                                   ttatcagaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTaaaa              <  1:717143‑M2/36‑1 (MQ=255)

                                   ttatcagaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTaa                >  1:919455‑M2/15‑48 (MQ=255)

                                       cagaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAt          >  1:1808090‑M2/11‑50 (MQ=255)

                                        agaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAtt         <  1:833814‑M2/41‑1 (MQ=255)

                                        agaacgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAtt         >  1:2329563‑M2/10‑50 (MQ=255)

                                           acgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTggg      <  1:389005‑M2/44‑1 (MQ=255)

                                            cgtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGt     <  1:128221‑M2/45‑1 (MQ=255)

                                             gtccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTa    <  1:951461‑M2/46‑1 (MQ=255)

                                              tccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTa    >  1:601613‑M2/4‑49 (MQ=255)

                                               ccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTa    >  1:2115060‑M2/3‑48 (MQ=255)

                                               ccTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTa    >  1:998750‑M2/3‑48 (MQ=255)

                                                cTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACa  <  1:2030809‑M2/48‑1 (MQ=255)

GAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACA  >  NC_000913/2290132‑2290228

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