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] 30	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.550)	26/466	0.7	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

TTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTG  >  NC_000913/258455‑258934                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                           aaattattagagatcccatttaacccggggtatcggaaatttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaaggacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggtccgtatattcgtgaaggtaAGTGCAAAGATAATCGTTTCtt                                                                                                                                                                                                                                       <  2:19608‑M2/22‑1 (MQ=255)                      aacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccattttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATg                                                                                                                                                                                                                  >  2:286602‑M2/209‑251 (MQ=255)                          cggtgtatcggaattttgacaaaggtcttcaggcactcaccgcatggatcgaaaaacgcttgtgttaaaaatttgcaaaattttgccaatttggtaacacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgtagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGGTAATCGATTCTTTTTCGATTGTCGGACTGTATCCGTc                                                                                                                                                                                                              >  2:383810‑M2/205‑251 (MQ=255)                                     aattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTc                                                                                                                                                                                                                            <  2:249677‑M2/33‑1 (MQ=255)                                     aattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTc                                                                                                                                                                                                                            >  1:249677‑M2/194‑226 (MQ=255)                                          tgacaaaggtcttcaggaaatcaccgactggatcgaaatacgcttgtgttaaaaatttgttaaattttgcctatttggtaaaacagttgcatcacaacatgagatatcaatgacgttaccgagtggactcccgaagagtagattgatcaaaaaatttaccgctctaggcccgtatattcgtgatgttaAGTGCAAAGATAAGCGATTCTTTTTCGATTGTCTGTCTGTATGCGTCATCTTGAAACCGGCa                                                                                                                                                                                               <  2:79294‑M2/62‑1 (MQ=255)                                              aaaggtcttcaggaaatcaccgactggatcgacaaacgcttgtgttaaaaatttgctaaattctgccaatttggtaaaacagttgcataacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGGATGAGTCAACGTGAAACCGTCACcgt                                                                                                                                                                                           >  2:98408‑M2/185‑249 (MQ=255)                                                             atcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTg                                                                                                                                                                           >  1:1231‑M2/170‑251 (MQ=255)                                                                 ccgactggagcgacaaacgcttgtgttaaaaatttgctaacttttgcctatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaaacgcagattgataaaaaaatttaccgccataggcccgtatattcgtgaaggtaAGTGCAACGTTAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCggcaccggaagtgcaatcgaggcccac                                                                                                                                                                        >  2:440057‑M2/166‑237 (MQ=255)                                                                                 acgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggccagtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTgg                                                                                                                                                       >  1:227004‑M2/150‑251 (MQ=255)                                                                                   gcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGaa                                                                                                                                                     <  1:286602‑M2/104‑1 (MQ=255)                                                                                      tgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAgcg                                                                                                                                                  >  1:82383‑M2/145‑251 (MQ=255)                                                                                                       aaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTAc                                                                                                                                 >  2:403615‑M2/128‑251 (MQ=255)                                                                                                                      ggtaaaacagttgcatcacaacaggagatagcactgaagtacccgagtcaccccgcgaagagcagattgaccaaaaaacttaccgccctaggcccgtatattcgtgaacgcaAGAGCAAAAATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGCTGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGtt                                                                                                                  >  1:246993‑M2/113‑251 (MQ=255)                                                                                                                           aacagttgcatcacaccaggagatagaaatgacgttaccgagtggacacccgaaaagcagattgatcacaaaatttaccgcactaggaccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGACGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGGAGCGCAGGAGTCCCGTTTTACCTACAGTTACCAGTTTGGTc                                                                                                             >  2:328890‑M2/108‑251 (MQ=255)                                                                                                                            acagatacatcacaacaggagctagcaatgacgttaccgagtggacacccgacgagcagattgatcaaaaactttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCACCGTGAAACCGGCACCAGAAGTGCGTGAATTCTGTGGCTGGTGGATGGAGATGGAAGCGCAGGAATCCCGTTTTACCTAAATTTACCAGTTTGGTCt                                                                                                            >  2:115877‑M2/107‑251 (MQ=255)                                                                                                                             cagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgacgagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGGTTTACCTACAGTTACCAGTTTGGTCTg                                                                                                           >  2:377299‑M2/106‑251 (MQ=255)                                                                                                                                        caacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGc                                                                                                >  1:134874‑M2/95‑251 (MQ=255)                                                                                                                                                                            ccgaagagcagattgaacaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCGGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGa                                                            <  2:430807‑M2/193‑1 (MQ=255)                                                                                                                                                                               aagagcagattgatcaaaaactttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAagg                                                         >  2:96286‑M2/56‑250 (MQ=255)                                                                                                                                                                                       attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGAGCAAAGATAATCGCTTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGATCCGGTAAAAGACACTGAAGTGGTTGAAc                                                 >  2:385976‑M2/48‑251 (MQ=255)                                                                                                                                                                                                              taggcccgtatcttcgtgaaggtcAGTGCAAAGATCATCGATTCTTTTTCGATTGTCTGGCTGCATGCGTCAACGTGAAACAGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTg                            >  2:259490‑M2/25‑249 (MQ=255)                                                                                                                                                                                                               aggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGt                         >  2:94262‑M2/24‑251 (MQ=255)                                                                                                                                                                                                                ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCt                                 <  1:93687‑M2/220‑1 (MQ=255)                                                                                                                                                                                                                ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCt                                 >  2:93687‑M2/23‑242 (MQ=255)                                                                                                                                                                                                                 gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGttt                       <  2:5550‑M2/230‑1 (MQ=255)                                                                                                                                                                                                                             gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAACGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTg           >  1:346767‑M2/10‑251 (MQ=255)                                                                                                                                                                                                                                   gTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCTGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCGGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGaa     <  2:206401/251‑1 (MQ=255)                                                                                                                                                                                                                                     aaGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAAct   >  1:223703/1‑251 (MQ=255)                                                                                                                                                                                                                                      aGTGCAAAGAGAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActg  <  2:1231/251‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                           TTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTG  >  NC_000913/258455‑258934

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