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–258674	258683	10–776	73 [0]	[0] 73	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)	68 (1.450)	23/118	0.2	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

ACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGG  >  NC_000913/258619‑258749                                                                  |                                                                  tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:1062594‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:365099‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:2074447‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:2074228‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:1828587‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:1403374‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:82492‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:983760‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:910825‑M2/3‑1 (MQ=255) tggacacccgaagagcagattgagcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGt                                                                 <  1:184134‑M2/3‑1 (MQ=255)  ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                >  1:1343692‑M2/65‑68 (MQ=255)  ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                <  1:2537718‑M2/4‑1 (MQ=255)  ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                >  1:1164725‑M2/65‑68 (MQ=255)  ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                <  1:3910368‑M2/4‑1 (MQ=255)  ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTg                                                                >  1:2015761‑M2/65‑68 (MQ=255)       cccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAg                                                           >  1:3643820‑M2/60‑68 (MQ=255)        ccgaagagcagattgatcaaaaaatttaccgcaccaggcccgtatattcgtgaaggtaAGTGCAAAGa                                                          <  1:2922256‑M2/10‑1 (MQ=255)         cgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGAt                                                         <  1:3344892‑M2/11‑1 (MQ=255)          gaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATa                                                        >  1:305328‑M2/57‑68 (MQ=255)            agagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAACGATAAt                                                      >  1:3891485‑M2/55‑68 (MQ=255)            agagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAAt                                                      <  1:684128‑M2/14‑1 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccttatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:2725303‑M2/53‑68 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:1766008‑M2/53‑68 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:584948‑M2/53‑68 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:2190239‑M2/53‑68 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:2205171‑M2/53‑68 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:2838144‑M2/53‑68 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:2729151‑M2/53‑68 (MQ=255)              agcagattgatcaaaaaatttaccgcactaggcccatatattcgtgaaggtaAGTGCAAAGATAATCg                                                    >  1:2747589‑M2/53‑68 (MQ=255)                  gattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTc                                                >  1:2877216‑M2/49‑68 (MQ=255)                              aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                    <  1:3529176‑M2/32‑1 (MQ=255)                              aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                    <  1:2996994‑M2/32‑1 (MQ=255)                              aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                    <  1:593925‑M2/32‑1 (MQ=255)                              aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                    <  1:3552033‑M2/32‑1 (MQ=255)                              aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                    <  1:2384186‑M2/32‑1 (MQ=255)                              aatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                    <  1:1738849‑M2/32‑1 (MQ=255)                              aagttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGt                                    <  1:3779500‑M2/32‑1 (MQ=255)                                  taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                                <  1:620724‑M2/36‑1 (MQ=255)                                  taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                                <  1:2700766‑M2/36‑1 (MQ=255)                                  taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                                <  1:2789104‑M2/36‑1 (MQ=255)                                  taccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTgg                                <  1:1335659‑M2/36‑1 (MQ=255)                                    ccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCt                              <  1:984218‑M2/38‑1 (MQ=255)                                      gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                            <  1:52993‑M2/40‑1 (MQ=255)                                      gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                            <  1:1423156‑M2/40‑1 (MQ=255)                                      gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                            <  1:1180931‑M2/40‑1 (MQ=255)                                      gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                            <  1:336465‑M2/40‑1 (MQ=255)                                      gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                            <  1:3905520‑M2/40‑1 (MQ=255)                                      gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                            <  1:3896601‑M2/40‑1 (MQ=255)                                      gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                            <  1:218486‑M2/40‑1 (MQ=255)                                      gcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGt                            <  1:3709518‑M2/40‑1 (MQ=255)                                          taggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGc                        <  1:121912‑M2/44‑1 (MQ=255)                                              cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGtcc                    >  1:1174528‑M2/21‑67 (MQ=255)                                              cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                    >  1:3045885‑M2/21‑68 (MQ=255)                                              cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                    >  1:1780248‑M2/21‑68 (MQ=255)                                              cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                    >  1:3776047‑M2/21‑68 (MQ=255)                                              cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                    >  1:355844‑M2/21‑68 (MQ=255)                                              cccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                    >  1:328794‑M2/21‑68 (MQ=255)                                              cccgtatattcgtgaagggaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCa                    >  1:745683‑M2/21‑68 (MQ=255)                                                cgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAAc                  <  1:2275655‑M2/50‑1 (MQ=255)                                                cgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAAc                  <  1:1348980‑M2/50‑1 (MQ=255)                                                cgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAAc                  <  1:1808009‑M2/50‑1 (MQ=255)                                                 gtgtattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACg                 <  1:3467752‑M2/51‑1 (MQ=255)                                                 gtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACg                 <  1:2423185‑M2/51‑1 (MQ=255)                                                  tatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGt                >  1:110551‑M2/17‑68 (MQ=255)                                                     attcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGAGTGTCTGGCTGTATGCGTCAACGTGaa             >  1:3272558‑M2/14‑68 (MQ=255)                                                       tcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAc           >  1:3948148‑M2/12‑68 (MQ=255)                                                        cgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc          >  1:1759411‑M2/11‑68 (MQ=255)                                                        cgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAAcc          >  1:2277256‑M2/11‑68 (MQ=255)                                                            aaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCa      <  1:2690064‑M2/62‑1 (MQ=255)                                                              ggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc    <  1:1554577‑M2/64‑1 (MQ=255)                                                              ggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCAcc    <  1:300717‑M2/64‑1 (MQ=255)                                                               gtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCg   >  1:2725036‑M2/4‑68 (MQ=255)                                                                taAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCgg  <  1:3291286‑M2/66‑1 (MQ=255)                                                                  |                                                                  ACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGG  >  NC_000913/258619‑258749

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