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

CTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGT  >  NC_000913/4626131‑4626226                                                 |                                                gagctgccgcgtttttccaccggatttaaagagttcgaccgcgtactag                                                 <  1:1279227‑M1/1‑1 (MQ=255)   gctgccgcgtttttccaccggatttaaagagttcgaccgcgtactaggc                                               >  1:2353126‑M1/47‑49 (MQ=255)     tgccgcgtttttccaccggatttaaagagttcgaccgcgtactaggcgg                                             <  1:2715088‑M1/5‑1 (MQ=255)           gtttttccaccggatttaaagagttcgaccgcgtactaGGCGGCgtggt                                       <  1:1838247‑M1/11‑1 (MQ=255)                tccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCa                                   >  1:2724526‑M1/34‑48 (MQ=255)                   accggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAgg                                 <  1:1780533‑M1/17‑1 (MQ=255)                     cggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGt                             <  1:1449795‑M1/21‑1 (MQ=255)                     cggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGt                             >  1:2329626‑M1/29‑49 (MQ=255)                       gatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGc                           >  1:2231991‑M1/27‑49 (MQ=255)                            aaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTc                      <  1:203770‑M1/28‑1 (MQ=255)                                agttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGAt                  >  1:338042‑M1/18‑49 (MQ=255)                                        cgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTa          >  1:1035980‑M1/10‑49 (MQ=255)                                         gcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTaa         <  1:1559052‑M1/41‑1 (MQ=255)                                         gcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTaa         >  1:2458392‑M1/9‑49 (MQ=255)                                          cgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAc        <  1:910303‑M1/42‑1 (MQ=255)                                                 GGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGt  >  1:1766317‑M1/1‑48 (MQ=255)                                                 |                                                CTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGT  >  NC_000913/4626131‑4626226

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.