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	28 [0]	[0] 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)	24 (0.840)	23/144	0.3	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

GAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTC  >  NC_000913/258607‑258763                                                                              |                                                                                gacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                              >  2:56410‑M2/78‑81 (MQ=255)  acgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                              >  2:56411‑M2/77‑80 (MQ=255)       accgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCa                                                                            >  1:56409‑M2/72‑77 (MQ=255)        ccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGAt                                                                       <  1:56268‑M2/11‑1 (MQ=255)          gagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                                                     >  2:56416‑M2/69‑81 (MQ=255)           agtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                                                    <  1:56269‑M2/14‑1 (MQ=255)                acacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGAtt                                                               <  1:56270‑M2/19‑1 (MQ=255)                      gaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttttt                                                         <  1:56178‑M2/25‑1 (MQ=255)                           gcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAtt                                                    >  1:56417‑M2/52‑81 (MQ=255)                               attgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt                                                <  2:56271‑M2/34‑1 (MQ=255)                                      aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                                          <  1:56272‑M2/40‑1 (MQ=255)                                       aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt                                        >  2:56419‑M2/40‑81 (MQ=255)                                       aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt                                        >  2:56414‑M2/40‑81 (MQ=255)                                         aaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGc                                      >  2:56412‑M2/38‑81 (MQ=255)                                            tttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTc                                   <  1:56273‑M2/47‑1 (MQ=255)                                              taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCaa                                 >  2:56438‑M2/33‑81 (MQ=255)                                               accgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAAc                                >  1:56424‑M2/32‑81 (MQ=255)                                                 cgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGt                              >  2:56426‑M2/30‑81 (MQ=255)                                                       aggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc                        >  2:56430‑M2/24‑81 (MQ=255)                                                         gcccgtatattcgttaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCgg                      >  1:56421‑M2/22‑81 (MQ=255)                                                          cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGc                     <  1:56274‑M2/61‑1 (MQ=255)                                                            cgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAc                   >  2:56415‑M2/19‑81 (MQ=255)                                                                 attcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCgg                >  2:56422‑M2/14‑79 (MQ=255)                                                                  ttcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAg             >  1:56428‑M2/13‑81 (MQ=255)                                                                        aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTg       <  1:56383‑M2/75‑1 (MQ=255)                                                                           gtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAt    >  1:56418‑M2/4‑81 (MQ=255)                                                                           gtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAt    >  2:56431‑M2/4‑81 (MQ=255)                                                                             aAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTc  >  1:56427‑M2/2‑81 (MQ=255)                                                                              aGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTc  <  2:56384/80‑1 (MQ=255)                                                                              |                                                                                GAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTC  >  NC_000913/258607‑258763

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑