Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	1,304,524	Δ367 bp		[yciU]–[cls]	[yciU], [cls]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	1304524	1304890	367	33 [0]	[0] 33	[yciU]–[cls]	[yciU],[cls]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 1304523	0 (0.000)	32 (0.590)	26/148	1.3	100%	coding (222/330 nt)	yciU	dsDNA‑mimic protein
	?	REL606	1304891 =	0 (0.000)					coding (1349/1461 nt)	cls	cardiolipin synthetase

CGACACGTTGCCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAATCGTCCTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTC  >  REL606/1304817‑1304965                                                                           |                                                                           ctttttcgcggcataacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTg                                                                           >  1:1440223‑M2/75‑76 (MQ=255)     ttcgcggcataacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAAcg                                                                       <  2:1132638‑M2/6‑1 (MQ=255)        gcggcataacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAAcggcg                                                                    <  1:1286316‑M2/9‑1 (MQ=255)         cggcataacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAAcggcgg                                                                   <  1:1088470‑M2/10‑1 (MQ=255)            cataacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCaa                                                                <  1:1327148‑M2/13‑1 (MQ=255)             ataacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAg                                                               >  2:1374079‑M2/63‑76 (MQ=255)             ataacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAg                                                               <  1:1353605‑M2/14‑1 (MQ=255)              taacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAgg                                                              <  2:607747‑M2/15‑1 (MQ=255)              taacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAgg                                                              <  1:631888‑M2/15‑1 (MQ=255)               aacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAaggg                                                             >  2:1315014‑M2/61‑75 (MQ=255)               aacagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGt                                                             >  2:1670853‑M2/61‑76 (MQ=255)                 cagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCt                                                           <  1:308651‑M2/18‑1 (MQ=255)                 cagaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCt                                                           <  1:1267610‑M2/18‑1 (MQ=255)                  agaatgcgggcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTg                                                          <  2:320731‑M2/19‑1 (MQ=255)                           gcgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAAc                                                 >  2:1205744‑M2/49‑76 (MQ=255)                            cgaaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAAcc                                                >  2:1335579‑M2/48‑76 (MQ=255)                              aaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCtt                                              >  2:1256481‑M2/46‑76 (MQ=255)                              aaaacatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGGCTGCACCAAAACCtt                                              >  2:338654‑M2/46‑76 (MQ=255)                                 acatcgttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTAt                                           <  2:75431‑M2/34‑1 (MQ=255)                                      gttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTAtcgtcg                                      <  1:372402‑M2/39‑1 (MQ=255)                                      gttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTAtcgtcg                                      <  1:372423‑M2/39‑1 (MQ=255)                                       ttgatttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGa                                     >  1:175736‑M2/37‑76 (MQ=255)                                          gtttcaccatcttcgctgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTg                                  <  2:684239‑M2/43‑1 (MQ=255)                                           tttcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCCAAACCTTTATCGTCGATTGc                                 >  1:329602‑M2/33‑76 (MQ=255)                                             tcaccatcttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCa                               <  1:935284‑M2/46‑1 (MQ=255)                                                    cttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAAt                        <  2:1441552‑M2/53‑1 (MQ=255)                                                    cttcactgtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAAt                        <  1:1478434‑M2/53‑1 (MQ=255)                                                           gtccgccagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGaaa                 >  1:1420511‑M2/17‑76 (MQ=255)                                                                 cagaccaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGc           >  2:982413‑M2/11‑76 (MQ=255)                                                                    accaatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCcac        <  1:126378‑M2/69‑1 (MQ=255)                                                                       aatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGa     <  2:1306731‑M2/72‑1 (MQ=255)                                                                       aatcTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGa     >  1:1484381‑M2/5‑76 (MQ=255)                                                                          cTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTc  >  1:702141‑M2/2‑76 (MQ=255)                                                                           |                                                                           CGACACGTTGCCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAATCGTCCTGAACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTC  >  REL606/1304817‑1304965

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