Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,798,370	(TTC)4→3	coding (124‑126/405 nt)	stpA ←	DNA‑binding transcriptional repressor StpA with RNA chaperone activity

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2798361	2798363	3	27 [0]	[0] 27	stpA	DNA‑binding transcriptional repressor StpA with RNA chaperone activity

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2798364 =	0 (0.000)	18 (0.740)	16/80	0.2	100%	coding (132/405 nt)	stpA	DNA‑binding transcriptional repressor StpA with RNA chaperone activity
	?	NC_000913	= 2798360	0 (0.000)					coding (136/405 nt)	stpA	DNA‑binding transcriptional repressor StpA with RNA chaperone activity

AGGTGCTAATTTTTTCCTGGCGCTCTGCCAGTTCACGCTGCTGCTGTTCTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAATTTTTCGAGCATTT  >  NC_000913/2798315‑2798413                                                  |                                                  gccaggtgctaattttttcctggcgctctgccagttcacgctgctgctgt                                                   >  1:1719359‑M1/50‑50 (MQ=255)   caggtgctaattttttcctggcgctctgccagttcacgctgctgctgttc                                                 >  1:1932654‑M1/48‑50 (MQ=255)   caggtgctaattttttcctggcgctctgccagttcacgctgctgctgttc                                                 >  1:404315‑M1/48‑50 (MQ=255)    aggtgctaattttttcctggcgctctgccagttcacgctgctgctgttct                                                <  1:1587105‑M1/4‑1 (MQ=255)    aggtgctaattttttcctggcgctctgccagttcacgctgctgctgttct                                                >  1:64523‑M1/47‑50 (MQ=255)       tgctaattttttcctggcgctctgccagttcacgctgctgctgttcttct                                             >  1:2106757‑M1/44‑50 (MQ=255)         ctaattttttcctggcgctctgccagttcacgctgctgctgttcttcttc                                           >  1:1104567‑M1/42‑50 (MQ=255)         ctaattttttcctggcgctctgccagttcacgctgctgctgttcttctt                                            >  1:844958‑M1/42‑49 (MQ=255)          taattttttcctggcgctctgccagttcacgctgctgctgTTCTTCTTCa                                          <  1:1540133‑M1/10‑1 (MQ=255)             ttttttcctggcgctctgccagttcacgctgctgctgTTCTTCTTCACGt                                       >  1:45679‑M1/38‑50 (MQ=255)                 ttcctggcgctctgccagttcacgctgctgctgTTCTTCTTCACGtcttt                                   <  1:2009908‑M1/17‑1 (MQ=255)                   cctggcgctctgccagttcacgctgctgctgTTCTTCTTCACGtctttc                                  >  1:371592‑M1/32‑49 (MQ=255)                       gcgctctgccagttcacgctgctgctgTTCTTCTTCACGTCTTTCTTTa                              >  1:2122930‑M1/28‑49 (MQ=255)                       gcgctctgccagttcacgctgctgctgTTCTTCTTCACGTCTTTCTTTa                              >  1:344636‑M1/28‑49 (MQ=255)                         gctctgccagttcacgctgctgctgTTCTTCTTCACGTCTTTCTTTAGTg                           >  1:213358‑M1/26‑50 (MQ=255)                             tgccagttcacgctgctgctgTTCTTCTTCACGTCTTTCTTTAGTGAc                         >  1:1150295‑M1/22‑48 (MQ=255)                               ccagttcacgctgctgctgTTCTTCTTCACGTCTTTCTTTAGTGACa                        >  1:2139290‑M1/20‑47 (MQ=255)                                 agttcacgctgctgctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCt                   >  1:2195749‑M1/18‑50 (MQ=255)                                 agttcacgctgctgctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCt                   >  1:811742‑M1/18‑50 (MQ=255)                                   ttcacgctgctgctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGa                 >  1:2210096‑M1/16‑50 (MQ=255)                                      acgctgctgctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAAtt              >  1:439553‑M1/13‑50 (MQ=255)                                        gctgctgctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAAtttt            <  1:359926‑M1/40‑1 (MQ=255)                                        gctgctgctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAAtttt            >  1:1738613‑M1/11‑50 (MQ=255)                                              gctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAATTTTTCGAg       <  1:277905‑M1/45‑1 (MQ=255)                                               ctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAATTTTTCGAGc      >  1:1976384‑M1/4‑49 (MQ=255)                                               ctgTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAATTTTTCGAGCa     <  1:1719817‑M1/47‑1 (MQ=255)                                                  TTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAATTTTTCGAGCAttt  >  1:2143546‑M1/1‑50 (MQ=255)                                                  |                                                  AGGTGCTAATTTTTTCCTGGCGCTCTGCCAGTTCACGCTGCTGCTGTTCTTCTTCTTCACGTCTTTCTTTAGTGACAACCCTGAATTTTTCGAGCATTT  >  NC_000913/2798315‑2798413

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: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.