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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	10 [0]	[0] 10	insB5–insA5	insB5,insA5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	10 (0.200)	8/468	2.2	100%	intergenic (‑305/+16)	flhD/insB5	DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA5/uspC	IS1 family protein InsA/universal stress protein C

CTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATC  >  NC_000913/1978308‑1978734                                                                                                                                                                                                   |                                                                                                                                                                                                                                         cTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGACCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGCCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                                                                                                                   >  1:465915‑M1/1‑195 (MQ=255)                                          ccaAATAGAATTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAAGGGGACTGCGATTTTCAATAATGCGTGATGCATATCACACAAAACACTCCATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagctagcaaaaatgcaacaggcataatcgaaatgggg                                                                                                                                         >  2:17032‑M1/2‑154 (MQ=255)                                              aTAGAAATGGGTCTTTACACTTATGTCAGATTTTTCCTAAATCGACGCAAATGTACTCGTCACTACACGCACATACAAAGGAGGGGGGATGCGATTATCAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttgggtgaatttcaatctgttaaaaagtttttcgctacgctagcaagaaaaaatgaaacagaaataatcgaaatgggatgt                                                                                                                                      >  2:509004‑M1/1‑150 (MQ=255)                                              aTAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                                                                                                                     >  1:16044‑M1/1‑150 (MQ=255)                                                                               ttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataat                                                                                                     <  1:562205‑M1/250‑134 (MQ=255)                                                                                 ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATACTGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat                                                                                                  >  1:23254‑M1/1‑115 (MQ=255)                                                                                                                acGCACATACAACGGAGGGGTGCTGCGAGTTTCAATAATTCTTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtattgggtgtatttcaatctgttaaaaagtttttcgctacgcaagcaagcaaaaatgaaacaggaataatcgaaatgggatttttcgcacagtcaaaataactcaccgtaaataatcatccgctataaataatccctttcatgcaatacc                                                                   <  2:465915‑M1/251‑168 (MQ=255)                                                                                                                                        gCGATATTCAATAATGCGGGAGGCAGATCACAAAAAACACTCAATTACTCAACATAAATGgataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcccagtcaaaataactcaccgtaaataatcatctgcgataaataatcacttgcatgcaataccagataatctattatcaaacaga                                             <  2:23254‑M1/249‑190 (MQ=255)                                                                                                                                                                             cacTCAATTACTTAACATAAATGtttaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaataggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcactacaaactaattaaatacaaagaaagg      >  2:407023‑M1/1‑23 (MQ=255)                                                                                                                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcacc  >  1:105126‑M1/1‑18 (MQ=255)                                                                                                                                                                                                   |                                                                                                                                                                                                                                         CTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATC  >  NC_000913/1978308‑1978734

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