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	32 [0]	[0] 31	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)	25 (1.070)	22/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

AAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTG  >  NC_000913/4626126‑4626224                                                  |                                                  aaGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTa                                                   <  1:288915/50‑1 (MQ=255)  agagCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAg                                                  >  1:1432982‑M2/1‑49 (MQ=255)  agagCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAg                                                  <  1:1784916‑M2/50‑2 (MQ=255)    agCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggc                                                <  1:167141‑M2/50‑4 (MQ=255)    agCTGCCGCGTTTTTCCACCGGATTTAAAGAGATCGACCGCGTACTAggc                                                <  1:2084076‑M2/50‑4 (MQ=255)        gCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcg                                            <  1:404687‑M2/50‑8 (MQ=255)         ccGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggc                                             >  1:34223‑M2/1‑42 (MQ=255)          cgcgTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtg                                          >  1:608868‑M2/1‑41 (MQ=255)            cgTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggt                                        <  1:775956‑M2/50‑12 (MQ=255)                tttCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgc                                      <  1:1189002‑M2/48‑14 (MQ=255)                tttCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgc                                      <  1:357300‑M2/48‑14 (MQ=255)                 ttCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccag                                   <  1:314412‑M2/50‑17 (MQ=255)                 ttCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgcca                                    >  1:1302005‑M2/1‑34 (MQ=255)                  tCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccag                                   >  1:2168910‑M2/1‑33 (MQ=255)                   ccACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccagg                                  <  1:1348415‑M2/49‑18 (MQ=255)                    cACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccagga                                 >  1:2246972‑M2/1‑31 (MQ=255)                         gATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgcc                           >  1:1631401‑M2/1‑26 (MQ=255)                         gATTTAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgc                            <  1:659939‑M2/49‑24 (MQ=255)                             tAAAGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccatt                        >  1:894684‑M2/1‑22 (MQ=255)                              aaaGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattct                      <  1:215773‑M2/50‑30 (MQ=255)                               aaGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattct                      >  1:160106‑M2/1‑20 (MQ=255)                               aaGAGTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattct                      >  1:1432258‑M2/1‑20 (MQ=255)                                  agTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgatt                  >  1:630553‑M2/1‑17 (MQ=255)                                   gTTCGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgattg                 >  1:1895300‑M2/1‑16 (MQ=255)                                      cGACCGCGTACTAggcggcgtggtgccaggaagtgccattctgattggcg              <  1:91574‑M2/50‑38 (MQ=255)                                       gACCGCGTACTAggcggcgtggtgccaggaagtgccattctgattggcgg             <  1:1598614‑M2/50‑39 (MQ=255)                                       gACCGCGTACTAggcggcgtggtgccaggaagtgccattctgattggcgg             >  1:1889717‑M2/1‑12 (MQ=255)                                         ccGCGTACTAggcggcgtggtgccaggaagtgccattctgattggcggta           >  1:1050134‑M2/1‑10 (MQ=255)                                              tACTAggcggcgtggtgccaggaagtgccattctgattggcggtcaccct      <  1:1922647‑M2/50‑46 (MQ=255)                                              tACTAggcggcgtggtgccaggaagtgccattctgattggcggtaaccct      <  1:1198620‑M2/50‑46 (MQ=255)                                                cTAggcggcgtggtgccaggaagtgccattctgattggcggtaaccctgg    >  1:1354650‑M2/1‑3 (MQ=255)                                                  aggcggcgtggtgccaggaagtgccattctgattggcggtaaccctggtg  <  1:691228‑M2/50‑50 (MQ=255)                                                  |                                                  AAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTG  >  NC_000913/4626126‑4626224

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑