Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	2,263,508	Δ3,531 bp	IS1‑mediated	rcsC–[atoS]	rcsC, [atoS]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	2262806–2263510	2267038	3529–4233	22 [21]	[0] 59	[insB‑17]–[atoS]	[insB‑17],insA‑17,rcsC,[atoS]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 2263507	NA (NA)	59 (0.920)	46/200	0.4	100%	noncoding (1/768 nt)	IS1	repeat region
	?	REL606	2267039 =	0 (0.000)					coding (297/1827 nt)	atoS	sensory histidine kinase in two‑component regulatory system with AtoC

ACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCC  >  REL606/2262705‑2262837                                                                                                     |                                 aCATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAg                                  >  2:1111412/1‑101 (MQ=255)  cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGc                                 <  2:981630/101‑1 (MQ=255)  cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGc                                 >  1:1628372/1‑101 (MQ=255)   aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCt                                <  1:183070/101‑1 (MQ=255)   aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCt                                >  1:198676/1‑101 (MQ=255)   aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCt                                <  2:1011189/101‑1 (MQ=255)   aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCt                                >  2:1626487/1‑101 (MQ=255)              cAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGAtttt                     >  2:482605/1‑101 (MQ=255)                   aaaTAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGaa                <  2:28936/101‑1 (MQ=255)                    aaTAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAAc               <  2:579284/101‑1 (MQ=255)                     aTAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACg              >  2:800267/1‑101 (MQ=255)                       aGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGAc            >  2:170873/1‑101 (MQ=255)                        atGCAGGCCGTCGGTAATGGCTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACa           <  1:1545311/100‑1 (MQ=255)                          gCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGc         >  1:1219155/1‑101 (MQ=38)                          gCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGc         <  1:667583/101‑1 (MQ=38)                          gCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGc         <  2:1346350/101‑1 (MQ=38)                           cAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCg        >  1:720299/1‑101 (MQ=38)                            aGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGa       >  1:799366/1‑101 (MQ=35)                            aGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGa       <  1:36384/101‑1 (MQ=35)                             ggCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGAc      >  2:311206/1‑101 (MQ=35)                                 gTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTcc  >  2:835196/1‑101 (MQ=11)                                 gTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTcc  >  2:1499133/1‑101 (MQ=11)                                                                                                     |                                 ACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCC  >  REL606/2262705‑2262837

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