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

TCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCAC  >  NC_000913/257664‑258144                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              tcaaccgttgacagtaaaTTATTAGCGAAACCCTTTCACCAGGTGTATCGGAATTTTGACAAAGGTCTTCAGGCAATCACCGACTGGATCGACAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGAAGATAGAAATGACGTTACCGAGTGGACAACCGAAGAGCAGATTGATCAAAAACTTTACCGCACTAGGCCCGTATATTCGAGAAGGTAagtgca                                                                                                                                                                                                                                         >  2:293922‑M1/9‑244 (MQ=255)   atcTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaa                                                                                                                                                                                                                                        >  1:277162‑M1/1‑242 (MQ=255)       gctgACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGATATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagata                                                                                                                                                                                                                                   >  2:395844‑M1/1‑238 (MQ=255)          gACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCGTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGGATATGCGGGAAGTTAagtgcaaagataatc                                                                                                                                                                                                                                <  2:88750‑M1/250‑16 (MQ=255)                  attattAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                                                                                                                                                                                                                       <  2:124889‑M1/251‑25 (MQ=255)                          agaTCCCATTTAACCCCGTGTAGCGGAATTTTGACAAAGGTCTTCAGGAAATCACCTTCTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGTAGTTAGCAATGACGTTACCGAGTGGACCCCCGAATAGCAGAGTGATCAAAATATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattg                                                                                                                                                                                                                <  1:303001‑M1/250‑32 (MQ=255)                                    taacctgaagtatcggaatcTTGACAAAGATCTTCAGGAAATCACCGACTGGATCGCCAACCGCTTGTGTTACCAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATAAAAAAATTTACCGCACTAGGCCCGTATATTCGTGACGGTAagtgcaaagataatcgattctttttcgattgtctggctgtat                                                                                                                                                                                                     >  2:259954‑M1/10‑209 (MQ=255)                                                                    ttCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                                                                                                                                                                     >  2:479657‑M1/1‑177 (MQ=255)                                                                        ggAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCCCAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaat                                                                                                                                                                 >  1:80101‑M1/1‑173 (MQ=255)                                                                           aaTCACCGACTGGATCGAAAAACGCGTGTGTTAAAAATTTGCTAAATTGTGACAATTTGGTAAAACAGGTGCATCACAACAGGAGATAGCAATGACGTTACAGAGTGGACAACCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATAGTAGTGAAGGTAagagcaaagataatcgattctttttcgagagtctggctgtatgcgtcaacgtgaaagaggcactggaagtgagagaattct                                                                                                                                                              >  1:118147‑M1/1‑170 (MQ=255)                                                                                ccGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaatgctggggc                                                                                                                                                         <  1:126174‑M1/251‑87 (MQ=255)                                                                                                     tgtgTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgggaattctggggctggtggatggagctggaag                                                                                                                                      <  1:33374‑M1/249‑106 (MQ=255)                                                                                                      gtgtTAAAAATTTGCTAAATTTTTCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTAGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCAGTATATTCGTGAAGGAAagtgaaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgc                                                                                                                                   >  1:146761‑M1/1‑143 (MQ=255)                                                                                                      gtgtTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgc                                                                                                                                   >  1:53269‑M1/1‑143 (MQ=255)                                                                                                           aaaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCCTCACAACAGGAGATAGCAATGACGTTAAAGAGTGGACACCCGAAGAGAAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAaatgaaaagctaatcgattctttttagattgtctggctgtatgcgtcaacgggaaaactgcaccggcagtgcgtgatttctggggatggtgggtgcagctggaagcgcagga                                                                                                                               >  2:399820‑M1/1‑138 (MQ=255)                                                                                                            aaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaat                                                                                                                             <  1:373121‑M1/251‑115 (MQ=255)                                                                                                                                          aaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGAATGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtc                                                                                               <  1:225177‑M1/251‑145 (MQ=255)                                                                                                                                                                           gttacacagtggagaCCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCAGTATATTCGCGAAGGTAagagcaacgataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccagcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttgggctgttcgataaagcaggcgactggaagagtgttc                                                              >  2:240644‑M1/8‑74 (MQ=255)                                                                                                                                                                                  aGTGGACACCCGAAGAGCAGAATGATCAAAAAATTTACCGCACTATGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgttaaagcaggcgactggaagagtgttccggtaaa                                                       <  1:274034‑M1/251‑185 (MQ=255)                                                                                                                                                                                       acacCCGAAGAGCAGATTGATCAAAACATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcacctgaagtgcgtgaattctggtgctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagaca                                                  <  1:431274‑M1/251‑190 (MQ=255)                                                                                                                                                                                            cGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactga                                              >  1:42300‑M1/1‑57 (MQ=255)                                                                                                                                                                                                    aGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgattggaagagtgttccggtaaaagacactgaagtggttg                                      >  1:392649‑M1/1‑49 (MQ=255)                                                                                                                                                                                                         gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgac                                >  1:38605‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                         gCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcg                >  1:339779‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                         gCACTAGGCCCGTATATTCGGGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaacaggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcg                >  2:69181‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                                       tcttagtgaagGTAagtgcaaagatcatcgtttctttttagattgtcttgctgtatgcgtcaacgtgaaaccggcactggaagtgcttgaattctggggctggtggatggagctggaagaacaggaatcacgttttacctacaggtactagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgacagaatggagcacaccatgcgagagtttcacg     >  2:23703‑M1/6‑14 (MQ=255)                                                                                                                                                                                                                                         ttCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttggtctgttcgataaagcaggcgactggaagagtgttccggtaaaagacactgaagtggttgaacgactggagcacaccctgcgtgagtttcacgaga  >  1:254787‑M1/1‑12 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              TCATCTGCTGACGGTAAATTATTAGAGATCCCATTTAACCCGGTGTATCGGAATTTTGACAAAGGTCTTCAGGAAATCACCGACTGGATCGAAAAACGCTTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCAC  >  NC_000913/257664‑258144

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: ‑