Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,872,030	Δ1,437 bp		[dgcJ]	[dgcJ]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1872030	1873466	1437	30 [0]	[0] 30	[dgcJ]	[dgcJ]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1872029	0 (0.000)	30 (1.010)	27/96	0.0	100%	intergenic (+169/‑12)	cdgI/dgcJ	putative c‑di‑GMP binding protein CdgI/putative diguanylate cyclase DgcJ
	?	NC_000913	1873467 =	0 (0.000)					coding (1427/1491 nt)	dgcJ	putative diguanylate cyclase DgcJ

GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACA  >  NC_000913/1873419‑1873514                                                 |                                                gccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagAt                                                <  1:2622826‑M2/2‑1 (MQ=255)  ccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagATg                                               >  1:27985‑M2/48‑50 (MQ=255)    gcaatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                             >  1:1587829‑M2/46‑50 (MQ=255)      aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                           >  1:1390587‑M2/44‑50 (MQ=255)       atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtat                                          <  1:1567560‑M2/8‑1 (MQ=255)       atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                           >  1:470261‑M2/43‑49 (MQ=255)       atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                           >  1:132712‑M2/43‑49 (MQ=255)         cgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                         >  1:930425‑M2/41‑49 (MQ=255)           aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                      <  1:1938772‑M2/12‑1 (MQ=255)            actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                     <  1:2549356‑M2/13‑1 (MQ=255)                tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                 >  1:349626‑M2/34‑50 (MQ=255)                 atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                 <  1:1817914‑M2/17‑1 (MQ=255)                   ctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                              >  1:884021‑M2/31‑50 (MQ=255)                    tggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                              >  1:2014547‑M2/30‑49 (MQ=255)                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGa                            <  1:1988871‑M2/22‑1 (MQ=255)                        tttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                         >  1:2076982‑M2/26‑50 (MQ=255)                        tttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                         >  1:1123868‑M2/26‑50 (MQ=255)                         ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                         <  1:1080644‑M2/25‑1 (MQ=255)                         ttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGt                        >  1:2349208‑M2/25‑50 (MQ=255)                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                     <  1:1407335‑M2/29‑1 (MQ=255)                               tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                  >  1:754477‑M2/19‑50 (MQ=255)                                ttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                 >  1:2617605‑M2/18‑50 (MQ=255)                                  cactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                >  1:1250293‑M2/16‑49 (MQ=255)                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc               >  1:2091203‑M2/15‑49 (MQ=255)                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc               >  1:1557987‑M2/15‑49 (MQ=255)                                     taaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa             >  1:1773661‑M2/13‑49 (MQ=255)                                      aaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt            >  1:1007298‑M2/12‑49 (MQ=255)                                            aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg      >  1:1069394‑M2/6‑49 (MQ=255)                                              gagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGaa    >  1:891890‑M2/4‑49 (MQ=255)                                               agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa  <  1:1318055‑M2/48‑1 (MQ=255)                                                 |                                                GCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACA  >  NC_000913/1873419‑1873514

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