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	63 [0]	[0] 64	[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)	63 (1.230)	48/98	0.0	100%	pseudogene (2515/2525 nt)	yejO	adhesin‑like autotransporter YejO
	?	NC_000913	2290181 =	0 (0.000)					intergenic (‑1/‑319)	yejO/narP	adhesin‑like autotransporter YejO/DNA‑binding transcriptional dual regulator NarP

AGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATT  >  NC_000913/2286353‑2286449                                                |                                                  aGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCt                                                  >  1:1827939‑M1/1‑48 (MQ=255)  gACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtga                                                <  1:1780650‑M1/50‑4 (MQ=255)  gACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCt                                                  >  1:2475280‑M1/1‑47 (MQ=255)    ccGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagc                                              >  1:780116‑M1/1‑45 (MQ=255)    ccGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagc                                              <  1:1424684‑M1/50‑6 (MQ=255)    ccGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagc                                              >  1:2718044‑M1/1‑45 (MQ=255)      gTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctt                                            <  1:2648084‑M1/50‑8 (MQ=255)      gTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctt                                            >  1:310701‑M1/1‑43 (MQ=255)         acacATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctttca                                         <  1:1686752‑M1/50‑11 (MQ=255)            caTCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcc                                      >  1:3380235‑M1/1‑37 (MQ=255)             aTCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcc                                      >  1:2470305‑M1/1‑36 (MQ=255)             aTCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcc                                      >  1:2674542‑M1/1‑36 (MQ=255)             aTCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcc                                      >  1:1070175‑M1/1‑36 (MQ=255)                acaGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcctat                                   <  1:3242653‑M1/49‑17 (MQ=255)                 caGCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcctatg                                  >  1:312258‑M1/1‑32 (MQ=255)                   gCCTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaa                               <  1:1554841‑M1/50‑21 (MQ=255)                    ccTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaat                              <  1:2846877‑M1/50‑22 (MQ=255)                    ccTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaat                              >  1:3833302‑M1/1‑29 (MQ=255)                     cTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaatt                             <  1:2003337‑M1/50‑23 (MQ=255)                     cTGTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaatt                             <  1:62114‑M1/50‑23 (MQ=255)                       gTTTATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaa                           >  1:2012499‑M1/1‑26 (MQ=255)                        tttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaat                          >  1:3234678‑M1/1‑25 (MQ=255)                        tttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaat                          >  1:2654175‑M1/1‑25 (MQ=255)                        tttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaat                          >  1:1075271‑M1/1‑25 (MQ=255)                        tttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaat                          <  1:2169453‑M1/50‑26 (MQ=255)                        tttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaat                          >  1:159655‑M1/1‑25 (MQ=255)                        tttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaa                           <  1:3890416‑M1/49‑25 (MQ=255)                        tttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaatta                            <  1:61927‑M1/48‑24 (MQ=255)                        tttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaatta                            <  1:1826374‑M1/48‑24 (MQ=255)                         ttATTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaatt                         >  1:180771‑M1/1‑24 (MQ=255)                           aTTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgc                       >  1:3836446‑M1/1‑22 (MQ=255)                           aTTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgc                       >  1:2792091‑M1/1‑22 (MQ=255)                           aTTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattg                        <  1:2557273‑M1/49‑28 (MQ=255)                           aTTTTCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattg                        <  1:2041162‑M1/49‑28 (MQ=255)                            ttttCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgct                      >  1:1334527‑M1/1‑21 (MQ=255)                             tttCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctg                     >  1:2800‑M1/1‑20 (MQ=255)                              ttCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctgt                    >  1:334218‑M1/1‑19 (MQ=255)                              ttCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctg                     <  1:2132322‑M1/49‑31 (MQ=255)                              ttCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctg                     <  1:1360876‑M1/49‑31 (MQ=255)                               tCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctgt                    >  1:2894403‑M1/1‑18 (MQ=255)                               tCTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctgt                    <  1:1350100‑M1/49‑32 (MQ=255)                                cTGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctgtt                   >  1:501519‑M1/1‑17 (MQ=255)                                 tGTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctgtta                  >  1:1560794‑M1/1‑16 (MQ=255)                                  gTTATCAGAACGTCCtgagctttcatcctatgaaattaattgctgttaa                 >  1:3818412‑M1/1‑15 (MQ=255)                                   ttATCAGAACGTCCtgagctttcatcctatgaaattaattgctgttaaaa               <  1:740634‑M1/50‑37 (MQ=255)                                   ttATCAGAACGTCCtgagctttcatcctatgaaattaattgctgttaaaa               <  1:2459941‑M1/50‑37 (MQ=255)                                   ttATCAGAACGTCCtgagctttcatcctatgaaattaattgctgttaaa                <  1:3191107‑M1/49‑36 (MQ=255)                                   ttATCAGAACGTCCtgagctttcatcctatgaaattaattgctgttaa                 <  1:3932774‑M1/48‑35 (MQ=255)                                      tCAGAACGTCCtgagctttcatcctatgaaattaattgctgttaaaagca            >  1:3134879‑M1/1‑11 (MQ=255)                                      tCAGAACGTCCtgagctttcatcctatgaaattaattgctgttaaaagca            <  1:452870‑M1/50‑40 (MQ=255)                                       cAGAACGTCCtgagctttcatcctatgaaattaattgctgttaaaagcat           >  1:964177‑M1/1‑10 (MQ=255)                                       cAGAACGTCCtgagctttcatcctatgaaattaattgctgttaaaagcat           >  1:298822‑M1/1‑10 (MQ=255)                                        aGAACGTCCtgagctttcatcctatgaaattaattgctgttaaaagcat           <  1:3875389‑M1/49‑41 (MQ=255)                                         gAACGTCCtgagctttcatcctatgaaattaattgctgttaaaagcattg         <  1:1629609‑M1/50‑43 (MQ=255)                                          aaCGTCCtgagctttcatcctatgaaattaattgctgttaaaagcattgg        <  1:2663563‑M1/50‑44 (MQ=255)                                           aCGTCCtgagctttcatcctatgaaattaattgctgttaaaagcattggg       >  1:779706‑M1/1‑6 (MQ=255)                                             gTCCtgagctttcatcctatgaaattaattgctgttaaaagcattgggta     <  1:3273019‑M1/50‑47 (MQ=255)                                             gTCCtgagctttcatcctatgaaattaattgctgttaaaagcattgggta     >  1:2494971‑M1/1‑4 (MQ=255)                                             gTCCtgagctttcatcctatgaaattaattgctgttaaaagcattgggta     <  1:1080442‑M1/50‑47 (MQ=255)                                             gTCCtgagctttcatcctatgaaattaattgctgttaaaagcattgggta     >  1:842460‑M1/1‑4 (MQ=255)                                              tCCtgagctttcatcctatgaaattaattgctgttaaaagcattgggtac    <  1:2308800‑M1/50‑48 (MQ=255)                                                ctgagctttcatcctatgaaattaattgctgttaaaagcattgggtacag  <  1:928906‑M1/50‑50 (MQ=255)                                                ctgagctttcatcctatgaaattaattgctgttaaaagcattgggtaca   <  1:1533894‑M1/49‑49 (MQ=255)                                                |                                                  AGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATT  >  NC_000913/2286353‑2286449

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