Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp		[crl]	[crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908–258674	258683	10–776	19 [0]	[0] 19	[crl]	[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	17 (1.020)	13/116	0.2	100%	intergenic (+8/‑769)	crl/crl	pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers/pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	pseudogene, sigma factor‑binding protein, RNA polymerase holoenzyme formation stimulator;regulator; Surface structures; transcriptional regulator of cryptic csgA gene for curli surface fibers

CGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/257833‑257907 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGG  >  NC_000913/258684‑258749                                                                                                                                               CGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGT                                                                           >  2:237555/1‑68     ACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAG                                                                       >  1:471368/1‑68         AGTGGACACCCGAAGAGCAGATTGATCTAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAA                                                                   >  1:662399/1‑68           TGGTCACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGT                                                                 <  2:700057/68‑1           TGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGT                                                                 <  2:70673/68‑1           TGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGT                                                                 <  2:471368/68‑1            GGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTG                                                                >  1:519114/1‑68            GGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTG                                                                >  2:548294/1‑68            GGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTG                                                                <  2:519114/68‑1                           AGATTGATCAAAAAATTTTCCGCACTAGGCCCGTATATTCGTGATGGTAAGTGCAAAGATAATCGATT                                                 <  2:678219/68‑1                              TTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTT                                              >  2:573023/1‑68                                                  ACTAGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAT                          <  1:237555/68‑1                                                     AGGCCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCG                       <  1:624655/68‑1                                                        CCCGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCA                    >  2:7045/1‑68                                                          CGTATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAAC                  <  2:178760/68‑1                                                            TATATTCGAGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGT                >  1:274280/1‑68                                                             ATATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTG               >  1:99293/1‑68                                                              TATTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGA              >  2:242635/1‑68                                                                cTCGTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAA            >  2:509158/2‑68                                                                    TGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGG        <  1:146279/68‑1                                                                          TAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGG  >  2:424007/1‑68                                                                                                                                               CGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/257833‑257907 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGG  >  NC_000913/258684‑258749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.