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	33 [0]	[0] 35	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)	29 (1.120)	23/84	0.0	100%	coding (268/1383 nt)	radA	DNA recombination protein
	?	NC_000913	= 4626175	0 (0.000)					coding (264/1383 nt)	radA	DNA recombination protein

GAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACC  >  NC_000913/4626128‑4626221                                                |                                               gagCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAgg                                              <  1:477890‑M2/50‑3 (MQ=255) gagCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAgg                                              >  1:1678656‑M2/1‑48 (MQ=255)      gCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcg                                         <  1:2675796‑M2/50‑8 (MQ=255)      gCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcg                                         <  1:2395097‑M2/50‑8 (MQ=255)           gTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtg                                    <  1:1327673‑M2/50‑13 (MQ=255)            tttttCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgc                                   <  1:833108‑M2/50‑14 (MQ=255)             ttttCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgc                                   <  1:1827837‑M2/49‑14 (MQ=255)                 ccACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccagg                               <  1:723257‑M2/49‑18 (MQ=255)                  cACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccagga                              >  1:94782‑M2/1‑31 (MQ=255)                  cACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccagga                              >  1:2278392‑M2/1‑31 (MQ=255)                  cACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccagg                               <  1:1710698‑M2/48‑18 (MQ=255)                   aCCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaag                            <  1:787443‑M2/50‑21 (MQ=255)                    ccGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagt                           <  1:1401736‑M2/50‑22 (MQ=255)                     cGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtg                          <  1:923153‑M2/50‑23 (MQ=255)                        aTTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgcca                       <  1:855466‑M2/50‑26 (MQ=255)                           tAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattc                    <  1:1330309‑M2/50‑29 (MQ=255)                           tAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattc                    <  1:2005711‑M2/50‑29 (MQ=255)                            aaaGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattct                   <  1:325705‑M2/50‑30 (MQ=255)                            aaaGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattc                    >  1:1599267‑M2/1‑21 (MQ=255)                            aaaGAGTTCGACCGCGTACTAggcggcatggtgccaggaagtgccattc                    >  1:2157088‑M2/1‑21 (MQ=255)                              agagTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctga                 <  1:605252‑M2/50‑32 (MQ=255)                                agTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgatt               >  1:2619866‑M2/1‑17 (MQ=255)                                agTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgatt               >  1:1461890‑M2/1‑17 (MQ=255)                                 gTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgattg              >  1:2568726‑M2/1‑16 (MQ=255)                                 gTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgattg              <  1:1665593‑M2/50‑35 (MQ=255)                                  ttCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgattgg             >  1:1376500‑M2/1‑15 (MQ=255)                                   tCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgattggc            <  1:726173‑M2/50‑37 (MQ=255)                                      aCCGCGTACTAggcggcgtggtgccaggaagtgccattctgattggcggt         >  1:726713‑M2/1‑11 (MQ=255)                                         gcgTACTAggcggcgtggtgccaggaagtgccattctgattggcggtaac      >  1:1958707‑M2/1‑8 (MQ=255)                                            tACTAggcggcgtggtgccaggaagtgccattctgattggcggtaaccct   >  1:2236092‑M2/1‑5 (MQ=255)                                            tACTAggcggcgtggtgccaggaagtgccattctgattggcggtaaccct   >  1:2043819‑M2/1‑5 (MQ=255)                                            tACTAggcggcgtggtgccaggaagtgccattctgattggcggtaaccct   >  1:1101688‑M2/1‑5 (MQ=255)                                              cTAggcggcgtggtgccaggaagtgccattctgattggcggtaaccctg  >  1:2591303‑M2/1‑3 (MQ=255)                                                |                                               GAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACC  >  NC_000913/4626128‑4626221

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ 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.