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	20 [0]	[0] 20	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)	18 (0.770)	17/82	0.1	100%	coding (268/1383 nt)	radA	DNA recombination protein
	?	NC_000913	= 4626175	0 (0.000)					coding (264/1383 nt)	radA	DNA recombination protein

GCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGG  >  NC_000913/4626133‑4626225                                               |                                               gctgccgcgtttttccaccggatttaaagagttcgaccgcgtactaggc                                              <  1:1113199‑M1/3‑1 (MQ=255)  ctgccgcgtttttccaccggatttaaagagttcgaccgcgtactaggcg                                             <  1:2037416‑M1/4‑1 (MQ=255)         gtttttccaccggatttaaagagttcgaccgcgtactaGGCGGCgtggt                                      <  1:1130489‑M1/11‑1 (MQ=255)            tttccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGc                                    <  1:1432491‑M1/13‑1 (MQ=255)            tttccaccggatttaaagagttcgaccgcgtactaGGCGGCATGGTGc                                    <  1:1853335‑M1/13‑1 (MQ=255)             ttccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGcc                                   >  1:109349‑M1/35‑48 (MQ=255)              tccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCa                                  >  1:1124138‑M1/34‑48 (MQ=255)                caccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAgg                                <  1:1619773‑M1/17‑1 (MQ=255)                    ggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTg                           >  1:2120297‑M1/28‑49 (MQ=255)                      atttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGc                          <  1:206339‑M1/23‑1 (MQ=255)                         taaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCAtt                      <  1:1152774‑M1/27‑1 (MQ=255)                           aagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCt                    <  1:225207‑M1/29‑1 (MQ=255)                            agagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCt                    <  1:2024287‑M1/29‑1 (MQ=255)                             gagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                  <  1:1971945‑M1/31‑1 (MQ=255)                                   gaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGc             >  1:1948203‑M1/13‑48 (MQ=255)                                        cgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAc       >  1:1067204‑M1/8‑49 (MQ=255)                                         gtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAcc      >  1:2118258‑M1/7‑49 (MQ=255)                                         gtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAcc      <  1:128439‑M1/43‑1 (MQ=255)                                           actaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCt    <  1:226006‑M1/45‑1 (MQ=255)                                             taGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTgg  <  1:694086‑M1/47‑1 (MQ=255)                                               |                                               GCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGG  >  NC_000913/4626133‑4626225

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