Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,720,634	Δ1,442 bp	IS150‑mediated	insJ–insK	insJ, insK

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3720634	3722075	1442	38 [0]	[0] 34	insJ–insK	insJ,insK

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3720633	0 (0.000)	34 (1.130)	25/98	0.1	100%	noncoding (1/1443 nt)	IS150	repeat region
	?	NC_000913	3722076 =	0 (0.000)					intergenic (+27/+252)	insK/glyS	IS150 conserved protein InsB/glycine‑‑tRNA ligase subunit beta

CCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGG  >  NC_000913/3722030‑3722122                                               |                                               tatttctgcgcattcccagcctcctgtgttgatttccaacgagtgtAACt                                             >  1:2102394‑M2/47‑50 (MQ=255) tatttctgcgcattcccagcctcctgtgttgatttccaacgagtata                                                <  1:1827738‑M2/1‑1 (MQ=255) tatttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACt                                             >  1:1046644‑M2/47‑50 (MQ=255)   tttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACttt                                           <  1:1202222‑M2/6‑1 (MQ=255)   tttctgcgcattcccagcctcctgtgttgatttccaacgagtatAACttt                                           <  1:2230046‑M2/6‑1 (MQ=255)         cgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGAt                                        <  1:1744596‑M2/9‑1 (MQ=255)         cgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGAt                                        <  1:1791461‑M2/9‑1 (MQ=255)         cgcattcccagcctcctgtgttgatttccaacgagtatAACTTTGAt                                        <  1:646274‑M2/9‑1 (MQ=255)          gcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTtat                                    <  1:5980‑M2/13‑1 (MQ=255)          gcattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTta                                     >  1:1882578‑M2/38‑49 (MQ=255)           cattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTtata                                   <  1:1687166‑M2/14‑1 (MQ=255)           cattcccagcctcctgtgttgatttccaacgagtatAACTTTGATTtata                                   <  1:2138883‑M2/14‑1 (MQ=255)            attcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                  >  1:193414‑M2/36‑50 (MQ=255)            attcccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAg                                  >  1:696318‑M2/36‑50 (MQ=255)                ccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCa                               <  1:1364689‑M2/18‑1 (MQ=255)                ccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCa                               <  1:317093‑M2/18‑1 (MQ=255)                ccagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAg                              <  1:1924094‑M2/19‑1 (MQ=255)                 cagcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAgg                             <  1:1822306‑M2/20‑1 (MQ=255)                  agcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGt                            <  1:2055397‑M2/21‑1 (MQ=255)                   gcctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTg                           <  1:1932323‑M2/22‑1 (MQ=255)                    cctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTgg                          >  1:106898‑M2/28‑50 (MQ=255)                    cctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTgg                          >  1:1973719‑M2/28‑50 (MQ=255)                     ctcctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTggg                         <  1:1358913‑M2/24‑1 (MQ=255)                       cctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGc                       >  1:842773‑M2/25‑50 (MQ=255)                        ctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTgggg                        <  1:2171790‑M2/25‑1 (MQ=255)                        ctgtgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCt                      <  1:1994999‑M2/27‑1 (MQ=255)                           tgttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCttt                    <  1:600032‑M2/29‑1 (MQ=255)                             ttgatttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTtct                 <  1:1251046‑M2/32‑1 (MQ=255)                                  ttccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtc            >  1:974486‑M2/14‑50 (MQ=255)                                    ccaacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgtct           <  1:567037‑M2/38‑1 (MQ=255)                                     caacgagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGc         >  1:920871‑M2/11‑50 (MQ=255)                                         gagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt     >  1:1212307‑M2/7‑50 (MQ=255)                                         gagtatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt     <  1:1186298‑M2/44‑1 (MQ=255)                                            tatAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg  <  1:1674306‑M2/47‑1 (MQ=255)                                               |                                               CCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGG  >  NC_000913/3722030‑3722122

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