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	258683	776	33 [0]	[0] 33	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)	28 (1.230)	22/80	0.0	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

TCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGT  >  NC_000913/258642‑258729                                           |                                              tcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCaa                                         <  1:484733‑M2/7‑1 (MQ=255) tcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCa                                          >  1:53179‑M2/43‑48 (MQ=255)   aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCa                                          >  1:66212‑M2/41‑46 (MQ=255)   aaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCa                                          >  1:306480‑M2/41‑46 (MQ=255)    aaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCa                                          >  1:2122125‑M2/40‑45 (MQ=255)       aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                   <  1:1318462‑M2/13‑1 (MQ=255)       aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATaa                                   <  1:1481853‑M2/13‑1 (MQ=255)        atttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                  <  1:73733‑M2/14‑1 (MQ=255)          ttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATc                                 >  1:1462307‑M2/34‑48 (MQ=255)              cgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTc                            >  1:77599‑M2/30‑49 (MQ=255)               gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCt                           >  1:1447510‑M2/29‑49 (MQ=255)                 actaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCttt                         >  1:1564315‑M2/27‑49 (MQ=255)                  ctaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCtttt                        <  1:496774‑M2/24‑1 (MQ=255)                     ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTc                      <  1:1135159‑M2/26‑1 (MQ=255)                     ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTc                      <  1:2115672‑M2/26‑1 (MQ=255)                     ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTc                      <  1:1643114‑M2/26‑1 (MQ=255)                     ggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCg                     <  1:226800‑M2/27‑1 (MQ=255)                      gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGa                    >  1:557816‑M2/22‑49 (MQ=255)                      gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGa                    <  1:1240712‑M2/28‑1 (MQ=255)                       cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAt                   <  1:409651‑M2/29‑1 (MQ=255)                        ccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAtt                  >  1:1233613‑M2/20‑49 (MQ=255)                        ccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAtt                  >  1:295085‑M2/20‑49 (MQ=255)                        ccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAtt                  >  1:1761350‑M2/20‑49 (MQ=255)                          gtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                >  1:605899‑M2/18‑49 (MQ=255)                             tattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCt              <  1:1060504‑M2/34‑1 (MQ=255)                             tattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTg             >  1:546837‑M2/15‑49 (MQ=255)                                  gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTg         <  1:1872422‑M2/39‑1 (MQ=255)                                  gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt        <  1:1837174‑M2/40‑1 (MQ=255)                                   tgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTa       <  1:1944911‑M2/41‑1 (MQ=255)                                   tgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTa       <  1:1138196‑M2/41‑1 (MQ=255)                                     aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTAt      <  1:1965829‑M2/42‑1 (MQ=255)                                      aggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGc    >  1:1978068‑M2/6‑49 (MQ=255)                                         taAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGt  >  1:32431‑M2/3‑48 (MQ=255)                                           |                                              TCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGT  >  NC_000913/258642‑258729

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