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

GCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTG  >  NC_000913/4626130‑4626227                                                  |                                                 agagctgccgcgtttttccaccggatttaaagagttcgaccgcgtactag                                                  <  1:1784916‑M1/1‑1 (MQ=255) agagctgccgcgtttttccaccggatttaaagagttcgaccgcgtactag                                                  >  1:1432982‑M1/50‑50 (MQ=255)   agctgccgcgtttttccaccggatttaaagagttcgaccgcgtactaggc                                                <  1:167141‑M1/3‑1 (MQ=255)   agctgccgcgtttttccaccggatttaaagagatcgaccgcgtactaggc                                                <  1:2084076‑M1/3‑1 (MQ=255)       gccgcgtttttccaccggatttaaagagttcgaccgcgtactaGgcggcg                                            <  1:404687‑M1/7‑1 (MQ=255)        ccgcgtttttccaccggatttaaagagttcgaccgcgtactaggcggc                                             >  1:34223‑M1/43‑48 (MQ=255)         cgcgtttttccaccggatttaaagagttcgaccgcgtactaGGCGGCgtg                                          >  1:608868‑M1/42‑50 (MQ=255)           cgtttttccaccggatttaaagagttcgaccgcgtactaGGCGGCgtggt                                        <  1:775956‑M1/11‑1 (MQ=255)               tttccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGc                                      <  1:357300‑M1/13‑1 (MQ=255)               tttccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGc                                      <  1:1189002‑M1/13‑1 (MQ=255)                ttccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCa                                    >  1:1302005‑M1/35‑49 (MQ=255)                ttccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAg                                   <  1:314412‑M1/16‑1 (MQ=255)                 tccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAg                                   >  1:2168910‑M1/34‑49 (MQ=255)                  ccaccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAgg                                  <  1:1348415‑M1/17‑1 (MQ=255)                   caccggatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGa                                 >  1:2246972‑M1/32‑49 (MQ=255)                        gatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGcc                           >  1:1631401‑M1/27‑50 (MQ=255)                        gatttaaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGc                            <  1:659939‑M1/23‑1 (MQ=255)                            taaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCAtt                        >  1:894684‑M1/23‑49 (MQ=255)                             aaagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCt                      <  1:215773‑M1/29‑1 (MQ=255)                              aagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCt                      >  1:160106‑M1/21‑49 (MQ=255)                              aagagttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCt                      >  1:1432258‑M1/21‑49 (MQ=255)                                 agttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGAtt                  >  1:630553‑M1/18‑50 (MQ=255)                                  gttcgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTg                 >  1:1895300‑M1/17‑50 (MQ=255)                                     cgaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCg              <  1:91574‑M1/37‑1 (MQ=255)                                      gaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCgg             >  1:1889717‑M1/13‑50 (MQ=255)                                      gaccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCgg             <  1:1598614‑M1/38‑1 (MQ=255)                                        ccgcgtactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTa           >  1:1050134‑M1/11‑50 (MQ=255)                                             tactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTCACCCt      <  1:1922647‑M1/45‑1 (MQ=255)                                             tactaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCt      <  1:1198620‑M1/45‑1 (MQ=255)                                               ctaGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTgg    >  1:1354650‑M1/4‑50 (MQ=255)                                                 aGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTg  <  1:691228‑M1/49‑1 (MQ=255)                                                  |                                                 GCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTG  >  NC_000913/4626130‑4626227

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