Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,102,944	Δ1,424 bp	IS5‑mediated	[wbbL]–[wbbK]	[wbbL], [wbbK]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2101983–2102948	2104367	1420–2385	5 [4]	[0] 18	[insH7]–[wbbK]	[insH7],wbbL,[wbbK]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2104368 =	0 (0.000)	18 (0.450)	17/482	1.1	100%	coding (142/1119 nt)	wbbK	putative glycosyltransferase WbbK
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

TTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGACT  >  NC_000913/2104143‑2104617                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          gccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAAt                                                                                                                                                                                                                                  >  1:91270‑M1/226‑251 (MQ=255)   cggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGc                                                                                                                                                                                                                                >  2:296743‑M1/224‑251 (MQ=255)    ggccattaccagccttggggtaaaacgtctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgagctcatatcagggacttgttcgcaccttccGTAACTCTTTAGCAGAATGGCCTAATGCg                                                                                                                                                                                                                               >  2:260402‑M1/223‑251 (MQ=255)                ccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGAc                                                                                                                                                                                                                    <  2:328370‑M1/40‑1 (MQ=255)                          aaacggctcgacgacttccacaatgttttgccatggcagaatctgctccatgcgggaaaagacaatctcttttctggtctgccggcgcttactgctgacttcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTttat                                                                                                                                                                                                         <  2:275731‑M1/51‑1 (MQ=255)                           aacggctcgatggcttccacaatgttttgtcctggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaatttactgtcggcgaaggtaagttgatgactcctgatgaaccctgttctatggctccagatgacaaaaatgatctcatatcaggggcttgttcgcaccttccTTAACTCTTTTGCAGAATGGACTAATGCGATaaa                                                                                                                                                                                                                          <  2:28778‑M1/34‑1 (MQ=255)                           aacggctcgatgacttccacaatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATaaa                                                                                                                                                                                                                          >  1:28778‑M1/200‑233 (MQ=255)                                     tgacttccaccatgttttgccatggcagaatctgctccacgcgcgacaagaaaatctcttttctggtctgacggcgcttacggccgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttatatggctcaagatgaaaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTtgct                                                                                                                                                                                              <  1:260462‑M1/62‑1 (MQ=255)                                                  gttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAAtttt                                                                                                                                                                                  >  1:339653‑M1/177‑250 (MQ=255)                                                                         ccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCGGCCAAAAATTTTTTCAAAATGGTAAATGGACCGcc                                                                                                                                                           >  1:328370‑M1/154‑250 (MQ=255)                                                                                ggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTa                                                                                                                                                    <  1:120227‑M1/104‑1 (MQ=255)                                                                                        aaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaggggaagttgatgactcatgatgaaccctgttctatggctcaagatgacaaacatgatcacatatcagggactttttcgcaccttccTTAACTTTTTAGCAGAATTAACTAATGCGATAAAACTGTCATTTTCTTTATTATTAGTTGTTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAattgcac                                                                                                                                           >  1:221553‑M1/139‑248 (MQ=255)                                                                                                             ggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTtgct                                                                                                                                                                                              >  1:70604‑M1/118‑179 (MQ=255)                                                                                                             ggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTtgct                                                                                                                                                                                              <  2:70604‑M1/62‑1 (MQ=255)                                                                                                                                              taagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCt                                                                                     <  2:364809‑M1/167‑1 (MQ=255)                                                                                                                                                                        gttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGc                                                            <  1:35141‑M1/192‑1 (MQ=255)                                                                                                                                                                                                                         gcaccttccTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCg           <  2:412508‑M1/241‑1 (MQ=255)                                                                                                                                                                                                                                 cTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGACt  <  1:74660‑M1/250‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          TTTCACGGAATAAAATTCCTTTATAAAAAGGTGCAGGGTTATGACAATACACATATCTTTTTTTAGTGACGACATTGGCCGTAATATCATGCAGACAAATCCAATGCGTAGCATTCAGCTCTTTTGAAAGTTTTTTACAAACTACATATTCAAAGTGCAAACGTTTTAGCCACGACCCTTTAACCTCAGGAAACTCAATGAATTTAACCCATGGATAACTTTCTTTTAACTCTTTAGCAGAATGGACTAATGCGATAAAACTGACATTTTCTTTATTATTAGTTGCTGCCAAAAATTTTTTCAAAATGGTAAATGGACCGCCAGTGGTAAAATTGACCGCAGAAACAACGACTATGCTTTTTCCCATAATTTGGTCTCATGATTGTATTTCTTTATGATTTTTGCTGGTACTCCCGCAATGACAGTATTTTCGGGAATAGAACCTCTAACAACAGAATTGGCGCCGACTACGACT  >  NC_000913/2104143‑2104617

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