Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp	100%		insB9–[crl]	insB9, insA9, [crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908	258683	776	24 [0]	[0] 24	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	23 (1.030)	18/248	NT	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTT  >  NC_000913/258544‑258803                                                                                                                                             |                                                                                                                        ttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatataccgcactaggcccgtatattcgtgacggtaAGt                                                                                                                       >  2:201384‑M2/141‑143 (MQ=255)      aatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaaa                                                                                                                  <  1:423947‑M2/8‑1 (MQ=255)                       caatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                                                                 <  1:251215‑M2/25‑1 (MQ=255)                          tttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaacaaatttacggcactaggcccgtatatgcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTcgg                                                                                              >  2:251041‑M2/116‑142 (MQ=255)                            tggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacgcctgaccagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAtt                                                                                            <  1:216072‑M2/30‑1 (MQ=255)                            tggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagaggacacccgaagagcagattgatcaaaaaatttaccgcacaaggcccgtaaatgcgtgaaggtaAGTGCAAAGATAAACGATTCTTTTCCGatg                                                                                            >  2:237582‑M2/114‑142 (MQ=255)                                aaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatcttcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt                                                                                        >  1:356410‑M2/110‑143 (MQ=255)                                aaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt                                                                                        >  1:36594‑M2/110‑143 (MQ=255)                                          catcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatcttcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTCTATGc                                                                              >  1:300025‑M2/100‑143 (MQ=255)                                               caacaggagatagcaatgactttaccgagtgtccacccgtagagcagattgatcaaaaaatttaccgctcttggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCGTTTTCGATTGTCTGGCTGTATGCGTCaa                                                                         <  1:112519‑M2/49‑1 (MQ=255)                                                  cgggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGt                                                                      <  1:250861‑M2/52‑1 (MQ=255)                                                  caggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGt                                                                      <  1:205003‑M2/52‑1 (MQ=255)                                                  caggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGt                                                                      <  1:158859‑M2/52‑1 (MQ=255)                                                  caggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgaatatgcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGt                                                                      >  2:164201‑M2/92‑143 (MQ=255)                                                  caggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatctaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGt                                                                      <  1:239271‑M2/52‑1 (MQ=255)                                                          agcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtctattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCgg                                                              >  1:126940‑M2/84‑143 (MQ=255)                                                                   gttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccggatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCGGGCTGTATGCGTCAACGTGAAACCGCGACCGGAAg                                                     >  2:195984‑M2/75‑143 (MQ=255)                                                                        ggagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGt                                                <  1:56295‑M2/74‑1 (MQ=255)                                                                                        gagcagatttatcaaaaaatttaccgcactgggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCtggt                                <  1:361726‑M2/90‑1 (MQ=255)                                                                                        gagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACGGGAAGTGCGTGAATTCTGGGGCtggt                                >  1:163642‑M2/54‑143 (MQ=255)                                                                                              attgatcaaaaaatttaccgcacaaggcccgtatattcgtgaaggtaAGTGCAAAGAT‑ATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAGGTGCGTGAATTCTGGGGCTGGgggaggga                         >  2:466846‑M2/48‑143 (MQ=255)                                                                                                     aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAATGTCTGGTTGTATGCGTCAACGTGAAACCGGCACCGGAAGTACGTGAAGTCTGGGCGTGGTGGATGGAGGTGGa                   >  2:334335‑M2/41‑143 (MQ=255)                                                                                                     aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAGTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGCAGTGCGTGAGTTGTGGGGCTGGTGgctggagctgga                   >  2:31619‑M2/41‑143 (MQ=255)                                                                                                                      aggcccgaatattcgggaaggtaAGTGCAAAGATAAGCGAGTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGtt  <  1:32312‑M2/120‑1 (MQ=255)                                                                                                                                             |                                                                                                                        CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTT  >  NC_000913/258544‑258803

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