Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,618,992	Δ6,688 bp		rhsB–[yhhJ]	rhsB, yhhH, yhhI, [yhhJ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3618992	3625679	6688	34 [0]	[0] 33	rhsB–[yhhJ]	rhsB,yhhH,yhhI,[yhhJ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3618991	0 (0.000)	30 (0.820)	24/96	0.2	100%	intergenic (+2/‑201)	nikR/rhsB	DNA‑binding transcriptional repressor NikR/rhs element protein RhsB
	?	NC_000913	3625680 =	NA (NA)					coding (1124/1125 nt)	yhhJ	ABC transporter family protein YhhJ

TGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625632‑3625728                                                 |                                                 cggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTa                                                 >  1:861136‑M2/49‑50 (MQ=255) cggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTa                                                 <  1:805909‑M2/2‑1 (MQ=255)    cgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:963133‑M2/46‑50 (MQ=255)     gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGc                                              >  1:1105742‑M2/45‑49 (MQ=255)     gggcggcacgggcatttgcagtgcttgccgaaggaagattgagcTACGc                                              >  1:1693972‑M2/45‑49 (MQ=255)         ggcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTg                                         >  1:2794323‑M2/41‑50 (MQ=255)         gacacgggcatttgcagtgcttgccgaaggtagattgagtTACGCCATTg                                         >  1:2024917‑M2/41‑50 (MQ=255)          gcacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGt                                        >  1:2496496‑M2/40‑50 (MQ=255)           cacgggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTc                                       >  1:3959‑M2/39‑50 (MQ=255)              gggcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTccc                                     >  1:688811‑M2/36‑49 (MQ=255)                gcatttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAAt                                  <  1:1275359‑M2/17‑1 (MQ=255)                 catttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATc                                 <  1:1953056‑M2/18‑1 (MQ=255)                  atttgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCg                                <  1:1142511‑M2/19‑1 (MQ=255)                     tgcagtgcttgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCt                             >  1:265934‑M2/29‑50 (MQ=255)                              tgccgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATc                    >  1:2254465‑M2/20‑50 (MQ=255)                                 cgaaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgca                 <  1:479657‑M2/34‑1 (MQ=255)                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagc               <  1:1865683‑M2/36‑1 (MQ=255)                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:1109947‑M2/15‑49 (MQ=255)                                   aaggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcag                >  1:2636948‑M2/15‑49 (MQ=255)                                    aggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCgcagca              <  1:1336067‑M2/37‑1 (MQ=255)                                     ggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGcagcag             >  1:1753871‑M2/13‑50 (MQ=255)                                     ggaagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGcagcag             >  1:1510161‑M2/13‑50 (MQ=255)                                       aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           <  1:2302409‑M2/40‑1 (MQ=255)                                       aagattgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAgcg           <  1:949759‑M2/40‑1 (MQ=255)                                          attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        <  1:581286‑M2/43‑1 (MQ=255)                                          attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        >  1:1757826‑M2/8‑50 (MQ=255)                                          attgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCaa        >  1:1012160‑M2/8‑50 (MQ=255)                                            tgagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATg      >  1:2443489‑M2/6‑50 (MQ=255)                                             gagtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATgg     >  1:1566729‑M2/5‑50 (MQ=255)                                               gtTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTg   <  1:1124781‑M2/48‑1 (MQ=255)                                                tTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTg   >  1:2178998‑M2/2‑49 (MQ=255)                                                tTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTg   >  1:183981‑M2/2‑49 (MQ=255)                                                tTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGa  <  1:2644140‑M2/49‑1 (MQ=255)                                                 |                                                 TGCACCGATGCGTAGTCCCTATAAGCTTACGCATCGGGCAATTGTATTTACGCCATTGTCCCAATCGTCTTCCTGAATCGCAGCAGCGCAATGGTGA  >  NC_000913/3625632‑3625728

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