Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503	1979278	776	12 [0]	[0] 11	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	10 (0.440)	7/116	1.0	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

AAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAG  >  NC_000913/1979229‑1979344                                                   |                                                                  ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                  <  1:10478‑M2/18‑1 (MQ=255)        gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                           <  1:851770‑M2/25‑1 (MQ=255)             gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                       >  1:222190‑M2/39‑51 (MQ=255)             gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                       <  2:222190‑M2/13‑1 (MQ=255)                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:151957‑M2/34‑68 (MQ=255)                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:332785‑M2/34‑68 (MQ=255)                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:544242‑M2/34‑68 (MQ=255)                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  2:412511‑M2/34‑68 (MQ=255)                      acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                             >  1:844620‑M2/30‑68 (MQ=255)                            cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                       >  1:471915‑M2/24‑68 (MQ=255)                                                 tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  2:673724‑M2/3‑68 (MQ=255)                                                   |                                                                  AAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAG  >  NC_000913/1979229‑1979344

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