Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_007779	4,310,363	Δ1,199 bp		insH–[alsK]	insH, [alsK]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_007779	4310363–4311557	4311561	5–1199	31 [0]	[0] 31	insH–[alsK]	insH,[alsK]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_007779	= 4310362	0 (0.000)	29 (0.650)	24/192	0.8	100%	intergenic (‑286/‑69)	yjcS/insH	alkyl sulfatase/IS5 element protein
	?	NC_007779	4311562 =	0 (0.000)					coding (916/1029 nt)	alsK	D‑allose kinase

TTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGC  >  NC_007779/4311462‑4311662                                                                                                     |                                                                                                     aatatgagtcgttttatgtgtttattatatatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagc                                                                                                      >  1:231573‑M2/101‑101 (MQ=255)    atgagtcgttttatgtgtttattatatatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagcaca                                                                                                   <  3:1311682‑M2/4‑1 (MQ=255)       agtcgttttatgtgtttattatatatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGaa                                                                                                <  3:86521‑M2/7‑1 (MQ=255)        gtcgttttatgtgtttattatatatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAAc                                                                                               <  3:358481‑M2/8‑1 (MQ=255)            ttttatgtgtttattatatatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACtgtg                                                                                           >  3:382087‑M2/90‑101 (MQ=255)                   tgtttattatatatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCaaaaa                                                                                    <  4:43061‑M2/19‑1 (MQ=255)                   tgtttattatatatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggcgtgagCACAGAACTGTGGCaaaaa                                                                                    <  4:146180‑M2/19‑1 (MQ=255)                       tattatatatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGt                                                                                >  2:31451‑M2/79‑101 (MQ=255)                              tatattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATgcg                                                                         >  2:37455‑M2/72‑101 (MQ=255)                                tattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGcc                                                                       >  3:122199‑M2/70‑101 (MQ=255)                                tattgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGcc                                                                       >  3:643116‑M2/70‑101 (MQ=255)                                 attgtacctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCa                                                                      >  1:300862‑M2/69‑101 (MQ=255)                                        ctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATtgc                                                               <  1:1297201‑M2/40‑1 (MQ=255)                                        ctgtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATtgc                                                               >  3:1081310‑M2/62‑101 (MQ=255)                                          gtcgttggcacagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATtgctg                                                             <  3:1302020‑M2/42‑1 (MQ=255)                                                   acagcgggacattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAg                                                    >  1:1006488‑M2/51‑101 (MQ=255)                                                           acattgcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTa                                            <  1:797801‑M2/59‑1 (MQ=255)                                                                gcgctactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTc                                       >  1:1492327‑M2/38‑101 (MQ=255)                                                                    tactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCagaa                                   >  2:34180‑M2/34‑101 (MQ=255)                                                                    tactgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCagaa                                   >  1:709533‑M2/34‑101 (MQ=255)                                                                      ctgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCagaaga                                 <  3:232975‑M2/70‑1 (MQ=255)                                                                       tgtgccaacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGAt                                <  1:99433‑M2/71‑1 (MQ=255)                                                                             aacgttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGGCAGAAGAGGAggcg                          <  4:34155‑M2/77‑1 (MQ=255)                                                                                gttttactcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCa                       >  1:1090522‑M2/22‑101 (MQ=255)                                                                                      ctcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAg                 >  1:1538635‑M2/16‑101 (MQ=255)                                                                                      ctcatggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAg                 >  1:959397‑M2/16‑101 (MQ=255)                                                                                          tggggctttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCa             <  1:394417‑M2/90‑1 (MQ=255)                                                                                                tttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGacca       <  1:379185‑M2/96‑2 (MQ=255)                                                                                                tttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCt       >  3:495908‑M2/6‑101 (MQ=255)                                                                                                 ttagCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTg      >  1:1003144‑M2/5‑101 (MQ=255)                                                                                                     cacaGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGc  <  1:1388323/101‑1 (MQ=255)                                                                                                     |                                                                                                     TTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCTTAGCACAGAACTGTGGCAAAAAACGTTGATGCGCCAATATTGCTGCACCCTGAGCGCCATTAAAGTCAGAAGATGAGGCGGCAATAAAGCGCACGACCTGATGC  >  NC_007779/4311462‑4311662

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