Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,626,183	(GCG)3→2	coding (272‑274/1383 nt)	radA →	DNA recombination protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4626176	4626178	3	15 [0]	[0] 16	radA	DNA recombination protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4626179 =	0 (0.000)	12 (0.380)	12/82	0.7	100%	coding (268/1383 nt)	radA	DNA recombination protein
	?	NC_000913	= 4626175	0 (0.000)					coding (264/1383 nt)	radA	DNA recombination protein

GAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGT  >  NC_000913/4626128‑4626217                                                |                                           gagCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAg                                           <  1:1279227‑M2/49‑2 (MQ=255)   gCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggc                                         >  1:2353126‑M2/1‑46 (MQ=255)     tGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcgg                                       <  1:2715088‑M2/49‑6 (MQ=255)           gTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggt                                 <  1:1838247‑M2/49‑12 (MQ=255)                tCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgcca                             >  1:2724526‑M2/1‑33 (MQ=255)                   aCCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccagg                           <  1:1780533‑M2/47‑18 (MQ=255)                     cGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagt                       <  1:1449795‑M2/49‑22 (MQ=255)                     cGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagt                       >  1:2329626‑M2/1‑28 (MQ=255)                       gATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgc                     >  1:2231991‑M2/1‑26 (MQ=255)                            aaaGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattc                <  1:203770‑M2/49‑29 (MQ=255)                                agTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgat            >  1:338042‑M2/1‑17 (MQ=255)                                        cgcgTACTAggcggcgtggtgccaggaagtgccattctgattggcggta    >  1:1035980‑M2/1‑9 (MQ=255)                                         gcgTACTAggcggcgtggtgccaggaagtgccattctgattggcggtaa   <  1:1559052‑M2/49‑42 (MQ=255)                                         gcgTACTAggcggcgtggtgccaggaagtgccattctgattggcggtaa   >  1:2458392‑M2/1‑8 (MQ=255)                                          cgTACTAggcggcgtggtgccaggaagtgccattctgattggcggtaac  <  1:910303‑M2/49‑43 (MQ=255)                                                |                                           GAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGT  >  NC_000913/4626128‑4626217

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.