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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908–258676	258683	8–776	27 [0]	[2] 29	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)	27 (0.570)	21/478	1.0	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 family transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

TGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTG  >  NC_000913/258440‑258929                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      tcaaccgttgacagtaaattattagcgaaaccctttcaccaggtgtatcggaattttgacaaaggtcttcaggcaatcaccgactggatcgacaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacagaagatagaaatgacgttaccgagtggacaaccgaagagcagattgatcaaaaactttaccgcactaggcccgtatattcgagaaggtaAGTGCa                                                                                                                                                                                                                                                  >  2:293922‑M2/245‑250 (MQ=255)   atctgctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaa                                                                                                                                                                                                                                                 >  1:277162‑M2/243‑249 (MQ=255)       gctgacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggatatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATa                                                                                                                                                                                                                                            >  2:395844‑M2/239‑250 (MQ=255)          gacggtaaattattagagatcccatttaacccggtgtatcggaattttgacaaaggtcgtcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccggatatgcgggaagttaAGTGCAAAGATAATc                                                                                                                                                                                                                                         <  2:88750‑M2/15‑1 (MQ=255)                  attattagagatcccatttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCtttt                                                                                                                                                                                                                                <  2:124889‑M2/24‑1 (MQ=255)                          agatcccatttaaccccgtgtagcggaattttgacaaaggtcttcaggaaatcaccttctggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacagtagttagcaatgacgttaccgagtggacccccgaatagcagagtgatcaaaatatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTg                                                                                                                                                                                                                         <  1:303001‑M2/31‑1 (MQ=255)                                    taacctgaagtatcggaatcttgacaaagatcttcaggaaatcaccgactggatcgccaaccgcttgtgttaccaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgataaaaaaatttaccgcactaggcccgtatattcgtgacggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt                                                                                                                                                                                                              >  2:259954‑M2/210‑251 (MQ=255)                                                                    ttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGt                                                                                                                                                                              >  2:479657‑M2/178‑251 (MQ=255)                                                                        ggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcccaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAt                                                                                                                                                                          >  1:80101‑M2/174‑251 (MQ=255)                                                                           aatcaccgactggatcgaaaaacgcgtgtgttaaaaatttgctaaattgtgacaatttggtaaaacaggtgcatcacaacaggagatagcaatgacgttacagagtggacaaccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatagtagtgaaggtaAGAGCAAAGATAATCGATTCTTTTTCGAGAGTCTGGCTGTATGCGTCAACGTGAAAGAGGCACTGGAAGTGAGAGAATTCt                                                                                                                                                                       >  1:118147‑M2/171‑251 (MQ=255)                                                                                ccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATGCTGGGGc                                                                                                                                                                  <  1:126174‑M2/86‑1 (MQ=255)                                                                                                     tgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGGGAATTCTGGGGCTGGTGGATGGAGCTGGAAg                                                                                                                                               <  1:33374‑M2/105‑1 (MQ=255)                                                                                                      gtgttaaaaatttgctaaatttttccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtagacacccgaagagcagattgatcaaaaaatttaccgcactaggccagtatattcgtgaaggaaAGTGAAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgcgc                                                                                                                                            >  1:146761‑M2/144‑251 (MQ=255)                                                                                                      gtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgcgc                                                                                                                                            >  1:53269‑M2/144‑251 (MQ=255)                                                                                                           aaaaatttgctaaattttgccaatttggtaaaacagttgcctcacaacaggagatagcaatgacgttaaagagtggacacccgaagagaagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAATGAAAAGCTAATCGATTCTTTTTAGATTGTCTGGCTGTATGCGTCAACGGGAAAACTGCACCGGCAGTGCGTGATTTCTGGGGATGGTGGGTGCAGCTGGAAGCGCAGGa                                                                                                                                        >  2:399820‑M2/139‑250 (MQ=255)                                                                                                            aaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAAt                                                                                                                                      <  1:373121‑M2/114‑1 (MQ=255)                                                                                                                                          aacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagaatgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTc                                                                                                        <  1:225177‑M2/144‑1 (MQ=255)                                                                                                                                                                           gttacacagtggagacccgaagagcagattgatcaaaaaatttaccgcactaggccagtatattcgcgaaggtaAGAGCAACGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCAGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGGCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTc                                                                       >  2:240644‑M2/75‑251 (MQ=255)                                                                                                                                                                                  agtggacacccgaagagcagaatgatcaaaaaatttaccgcactatgcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGTTAAAGCAGGCGACTGGAAGAGTGTTCCGGTaaa                                                                <  1:274034‑M2/184‑1 (MQ=255)                                                                                                                                                                                       acacccgaagagcagattgatcaaaacatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCTGAAGTGCGTGAATTCTGGTGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGaca                                                           <  1:431274‑M2/189‑1 (MQ=255)                                                                                                                                                                                            cgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGa                                                       >  1:42300‑M2/58‑250 (MQ=255)                                                                                                                                                                                                    agattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGATTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTg                                               >  1:392649‑M2/50‑250 (MQ=255)                                                                                                                                                                                                         gatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGAc                                         >  1:38605‑M2/45‑251 (MQ=255)                                                                                                                                                                                                                         gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCg                         >  1:339779‑M2/29‑251 (MQ=255)                                                                                                                                                                                                                         gcactaggcccgtatattcgggaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACAGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCg                         >  2:69181‑M2/29‑251 (MQ=255)                                                                                                                                                                                                                                       tcttagtgaaggtaAGTGCAAAGATCATCGTTTCTTTTTAGATTGTCTTGCTGTATGCGTCAACGTGAAACCGGCACTGGAAGTGCTTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGAACAGGAATCACGTTTTACCTACAGGTACTAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGACAGAATGGAGCACACCATGCGAGAGTTTCACg              >  2:23703‑M2/15‑248 (MQ=255)                                                                                                                                                                                                                                         ttcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACgaga           >  1:254787‑M2/13‑249 (MQ=255)                                                                                                                                                                                                                                               aaGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCgg   >  1:242712/1‑250 (MQ=255)                                                                                                                                                                                                                                                 ggTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTg  >  1:278765/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTG  >  NC_000913/258440‑258929

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