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	258683	776	20 [0]	[2] 22	[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 (0.690)	14/116	0.4	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

ACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCG  >  NC_000913/258619‑258748                                                                  |                                                                 tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                <  1:913255‑M2/3‑1 (MQ=255)  ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                               >  1:559706‑M2/65‑68 (MQ=255)         cgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGAGCAAAGAt                                                        >  2:470241‑M2/58‑68 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCg                                                   >  1:690209‑M2/53‑68 (MQ=255)                    ttgatctaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCtt                                             >  2:452958‑M2/47‑68 (MQ=255)                    ttgatcaaaaaatttaccgcactaggcccttatattcgtgaaggtaAGCGCAAAGATAATCGATTCtt                                             <  2:441701‑M2/22‑1 (MQ=255)                         caaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCg                                        <  2:609150‑M2/27‑1 (MQ=255)                              aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                   <  2:43493‑M2/32‑1 (MQ=255)                              aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                   <  2:53661‑M2/32‑1 (MQ=255)                                  taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCtt                                             <  1:608549‑M2/22‑1 (MQ=255)                                  taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTCTTTCGATTGTCTgg                               <  1:100807‑M2/36‑1 (MQ=255)                                  taccgcactaggaccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCtt                                             >  2:608554‑M2/33‑54 (MQ=255)                                   accgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGc                              <  2:849479‑M2/37‑1 (MQ=255)                                    ccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCt                             <  2:637571‑M2/38‑1 (MQ=255)                                       cactaggaccgtatattcgtgaaggtaAGAGCAAAGATAATCGATTCTTTTTCGATAGTCTGGCTGTa                          >  2:934576‑M2/28‑68 (MQ=255)                                              cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                   >  1:877754‑M2/21‑68 (MQ=255)                                              cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                   >  2:71572‑M2/21‑68 (MQ=255)                                                        cgagAAGGTAAGAGAAAAGATAATCGATACTTTATCGATTGTCTGGCTGTATCAGACAACGTGAAAcc         >  2:391152/4‑68 (MQ=255)                                                         gtGATGGTTAGTGCAAAGATAATCGATTCTTATTCGATTTTCTGTCTGAATGCGTCATCGTGAAACCg        <  2:765651/67‑1 (MQ=255)                                                            aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTTTGGCTGTATGCGTCAACGTGAAACCGGCa     <  1:293910‑M2/62‑1 (MQ=255)                                                            aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCa     <  2:548657‑M2/62‑1 (MQ=255)                                                               gtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCg  >  1:588588‑M2/4‑68 (MQ=255)                                                                  |                                                                 ACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCG  >  NC_000913/258619‑258748

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