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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_pae_pgi	257908	258683	776	11 [0]	[0] 11	[crl]	[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	= 257907	0 (0.000)	11 (0.490)	9/270	1.3	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_3_pae_pgi	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

ACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACA  >  NC_000913_3_pae_pgi/258564‑258812                                                                                                                         |                                                                                                                                 gccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAt                                                                                                      <  2:65150‑M2/29‑1 (MQ=255)        tggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                                                                                               <  1:192080‑M2/36‑1 (MQ=255)        tggtaaaacagttgcaacacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                                                                                               <  2:240488‑M2/36‑1 (MQ=255)                             acaggagatagcaatgacgttaccgagtggacacccgaagggcagattgattaaaaaatttaccgcattaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGGCAACGTGAAAc                                                                          >  2:248626‑M2/93‑149 (MQ=255)                                              cgttaccgggtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGt                                                         <  1:125730‑M2/74‑1 (MQ=255)                                                ttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcatgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGa                                                       >  1:357262‑M2/74‑149 (MQ=255)                                                ttaccgagtggacacccgaagagcagagtgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGa                                                       >  2:153390‑M2/74‑149 (MQ=255)                                                                                           ccgcactaggcacgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCGGGAAGCGCAGGAATCCCGt            >  1:253499‑M2/31‑149 (MQ=255)                                                                                                     cccgtatttttgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACa  <  1:75430‑M2/129‑1 (MQ=255)                                                                                                                gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATccc              >  1:46717‑M2/10‑126 (MQ=255)                                                                                                                gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATccc              <  2:46717‑M2/117‑1 (MQ=255)                                                                                                                         |                                                                                                                                 ACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACA  >  NC_000913_3_pae_pgi/258564‑258812

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