Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	547,700	Δ8,125 bp	IS1‑mediated	[nmpC]–[ECB_00513]	[nmpC], ybcR, ybcS, ybcT, ybcU, ECB_00510, nohB, ECB_00512, [ECB_00513]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	546987–547698	555824	8127–8838	22 [20]	[0] 65	[insB‑6]–[ECB_00513]	[insB‑6],insA‑6,nmpC,ybcR,ybcS,ybcT,ybcU,ECB_00510,nohB,ECB_00512,[ECB_00513]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 547699	NA (NA)	64 (0.990)	48/200	0.3	100%	noncoding (1/768 nt)	IS1	repeat region
	?	REL606	555825 =	0 (0.000)					coding (1110/2346 nt)	ECB_00513	conserved hypothetical protein

CCACGCCGAATCCTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCC  >  REL606/546886‑547029                                                                                                     |                                            ccACGCCGAATCCTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACt                                             <  1:989630/101‑1 (MQ=255) ccACGCCGAATCCTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACt                                             <  2:37280/101‑1 (MQ=255)   aCGCCGAATCCTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACttt                                           <  2:952856/101‑1 (MQ=255)     gCCGAATCCTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGt                                         <  1:136944/101‑1 (MQ=255)        gAATCCTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCAt                                      <  1:1463646/101‑1 (MQ=255)        gAATCCTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCAt                                      <  1:172356/101‑1 (MQ=255)            ccTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAg                                  >  2:935782/1‑101 (MQ=255)                 aTACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCa                             >  1:6615/1‑101 (MQ=255)                     tCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGa                         >  1:311434/1‑101 (MQ=255)                      cATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGAt                        >  1:977726/1‑101 (MQ=255)                             agtGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTgaga                 >  2:160388/1‑101 (MQ=255)                              gtGGAGAAACCGAAACCGGTAATGACTCCAACTTAGTGATAGTTTTTTATGTTCAGACAATGCCCGATGACTTTGTCATGCAGCTCGACCGATTTTGAGaa                <  2:1705403/101‑1 (MQ=14)                               tGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAAc               <  1:129810/101‑1 (MQ=255)                                 gagaAACCGAAACCGGTAATGACCCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGa             >  1:1082018/1‑101 (MQ=255)                                     aaccgaaaccgGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGc         <  2:441057/101‑1 (MQ=38)                                      accgaaaccgGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCg        <  1:1455877/101‑1 (MQ=38)                                       ccgaaaccgGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGa       <  2:310448/101‑1 (MQ=35)                                        cgaaaccgGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGAc      <  1:1059233/101‑1 (MQ=35)                                          aaaccgGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACtt    <  2:825431/101‑1 (MQ=32)                                          aaaccgGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCAGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACtt    >  2:435830/1‑101 (MQ=17)                                           aaccgGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTc   >  1:955203/1‑101 (MQ=31)                                            accgGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTcc  <  2:1520752/101‑1 (MQ=11)                                                                                                     |                                            CCACGCCGAATCCTTCATACTCATAAGTAGTGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCC  >  REL606/546886‑547029

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