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

AAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAA  >  NC_000913/257679‑258149                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   aaaTTATTAGAGATCCCATTTAACCCGGGGTATCGGAAATTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAAGGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGTCCGTATATTCGTGAAGGTAagtgcaaagataatcgtttctt                                                                                                                                                                                                                              <  2:19608‑M1/251‑23 (MQ=255)                      aaCCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCATTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatg                                                                                                                                                                                                         >  2:286602‑M1/1‑208 (MQ=255)                          cGGTGTATCGGAATTTTGACAAAGGTCTTCAGGCACTCACCGCATGGATCGAAAAACGCTTGTGTTAAAAATTTGCAAAATTTTGCCAATTTGGTAACACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGTAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaaggtaatcgattctttttcgattgtcggactgtatccgtc                                                                                                                                                                                                     >  2:383810‑M1/1‑204 (MQ=255)                                     aaTTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                                                                                                                                                                                                                   >  1:249677‑M1/1‑193 (MQ=255)                                     aaTTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                                                                                                                                                                                                                   <  2:249677‑M1/226‑34 (MQ=255)                                          tGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAATACGCTTGTGTTAAAAATTTGTTAAATTTTGCCTATTTGGTAAAACAGTTGCATCACAACATGAGATATCAATGACGTTACCGAGTGGACTCC‑CGAAGAGTAGATTGATCAAAAAATTTACCGCTCTAGGCCCGTATATTCGTGATGTTAagtgcaaagataagcgattctttttcgattgtctgtctgtatgcgtcatcttgaaaccggca                                                                                                                                                                                      <  2:79294‑M1/250‑63 (MQ=255)                                              aaaGGTCTTCAGGAAATCACCGACTGGATCGACAAACGCTTGTGTTAAAAATTTGCTAAATTCTGCCAATTTGGTAAAACAGTTGCATAACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctggatgagtcaacgtgaaaccgtcaccgt                                                                                                                                                                                  >  2:98408‑M1/1‑184 (MQ=255)                                                             aTCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctg                                                                                                                                                                  >  1:1231‑M1/1‑169 (MQ=255)                                                                 ccGACTGGAGCGACAAACGCTTGTGTTAAAAATTTGCTAACTTTTGCCTATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAACGCAGATTGATAAAAAAATTTACCGCCATAGGCCCGTATATTCGTGAAGGTAagtgcaacgttaatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcaatcgaggcccac                                                                                                                                                               >  2:440057‑M1/1‑165 (MQ=255)                                                                                 aCGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCAGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctgg                                                                                                                                              >  1:227004‑M1/1‑149 (MQ=255)                                                                                   gCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaa                                                                                                                                            <  1:286602‑M1/251‑105 (MQ=255)                                                                                      tgtgTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcg                                                                                                                                         >  1:82383‑M1/1‑144 (MQ=255)                                                                                                       aaaTTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttac                                                                                                                        >  2:403615‑M1/1‑127 (MQ=255)                                                                                                                      ggTAAAACAGTTGCATCACAACAGGAGATAGCACTGAAGTACCCGAGT‑CACCCCGCGAAGAGCAGATTGACCAAAAAACTTACCGCCCTAGGCCCGTATATTCGTGAACGCAagagcaaaaataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggctggagctggaagcgcaggaatcccgttttacctacagttaccagtt                                                                                                         >  1:246993‑M1/1‑112 (MQ=255)                                                                                                                           aaCAGTTGCATCACACCAGGAGATAGAAATGACGTTACCGAGTGGACACC‑CGAAAAGCAGATTGATCACAAAATTTACCGCACTAGGACCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggacgtgcgtgaattctggggctggtggatggagctgggagcgcaggagtcccgttttacctacagttaccagtttggtc                                                                                                    >  2:328890‑M1/1‑107 (MQ=255)                                                                                                                            acagatacatcacaaCAGGAGCTAGCAATGACGTTACCGAGTGGACACC‑CGACGAGCAGATTGATCAAAAACTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaccgtgaaaccggcaccagaagtgcgtgaattctgtggctggtggatggagatggaagcgcaggaatcccgttttacctaaatttaccagtttggtct                                                                                                   >  2:115877‑M1/8‑106 (MQ=255)                                                                                                                             cAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGACGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccggtttacctacagttaccagtttggtctg                                                                                                  >  2:377299‑M1/1‑105 (MQ=255)                                                                                                                                        caacaGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagc                                                                                       >  1:134874‑M1/1‑94 (MQ=255)                                                                                                                                                                            c‑cGAAGAGCAGATTGAACAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtcggttcgataaagcaggcgactggaagagtgttccggtaaaagacactga                                                   <  2:430807‑M1/251‑194 (MQ=255)                                                                                                                                                                                aaGAGCAGATTGATCAAAAACTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagg                                                >  2:96286‑M1/1‑55 (MQ=255)                                                                                                                                                                                        aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagagcaaagataatcgcttctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgatccggtaaaagacactgaagtggttgaac                                        >  2:385976‑M1/1‑47 (MQ=255)                                                                                                                                                                                                               tAGGCCCGTATCTTCGTGAAGGTCagtgcaaagatcatcgattctttttcgattgtctggctgcatgcgtcaacgtgaaacaggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtg                   >  2:259490‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                aGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagt                >  2:94262‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                 ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccct                        >  2:93687‑M1/1‑22 (MQ=255)                                                                                                                                                                                                                 ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccct                        <  1:93687‑M1/242‑221 (MQ=255)                                                                                                                                                                                                                  gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagttt              <  2:5550‑M1/251‑231 (MQ=255)                                                                                                                                                                                                                              gTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataacgcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacgagaagctg  >  1:346767‑M1/1‑9 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   AAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACC‑CGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAA  >  NC_000913/257679‑258149

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