Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,320,144	Δ16,489 bp		[atoS]–[yfaA]	[atoS], atoC, atoD, atoA, atoE, atoB, yfaP, yfaQ, yfaS, yfaT, [yfaA]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2320144	2336632	16489	35 [0]	[0] 33	[atoS]–[yfaA]	[atoS],atoC,atoD,atoA,atoE,atoB,yfaP,yfaQ,yfaS,yfaT,[yfaA]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2320143	0 (0.000)	32 (1.190)	28/94	0.0	100%	coding (101/1827 nt)	atoS	sensory histidine kinase AtoS
	?	NC_000913	2336633 =	0 (0.000)					coding (12/1689 nt)	yfaA	DUF2138 domain‑containing protein YfaA

GCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAC  >  NC_000913/2336586‑2336680                                                |                                                atggcaatcctgatggtcattgtcccaacgcttactattggttatattt                                                >  1:587796‑M2/48‑49 (MQ=255)  tggcaatcctgatggtcattgtcccaacgcttactattggttatatTTc                                               <  1:1543257‑M2/3‑1 (MQ=255)   ggcaatcctgatggtcattgtcccaacgcttactattggttatatTTca                                              >  1:852653‑M2/46‑49 (MQ=255)     caatcctgatggtcattgtcccaacgcttactattggttatatTTcac                                             >  1:945583‑M2/44‑48 (MQ=255)       atcctgatggtcattgtcccaacgcttactattggttatatTTcaccac                                          >  1:1233163‑M2/42‑49 (MQ=255)          ctgatggtcattgtcccaacgcttactattggttatatTTCACCACTCa                                       >  1:442265‑M2/39‑49 (MQ=255)             atggtcattgtcccaacgcttactattggttatatTTCACCACTCAt                                      <  1:828536‑M2/12‑1 (MQ=255)             atggtcattgtcccaacgcttactattggttatatTTCACCACTCATTg                                    <  1:1384723‑M2/14‑1 (MQ=255)             atggtcattgtcccaacgcttactattggttatatTTCACCACTCATTg                                    >  1:1866190‑M2/36‑49 (MQ=255)              tggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                   <  1:1941366‑M2/15‑1 (MQ=255)               ggtcattgtcccaacgcttactattggttatatTTCACCACTCATTGc                                   >  1:464638‑M2/34‑48 (MQ=255)                gtcattgtcccaacgcttactattggttatatTTCACCACTCATTGcc                                  <  1:955956‑M2/16‑1 (MQ=255)                 tcattgtcccaacgcttactattggttatatTTCACCACTCATTGCcac                                <  1:722055‑M2/18‑1 (MQ=255)                    ttgtcccaacgcttactattggttatatTTCACCACTCATTGCCACAtt                             >  1:1516513‑M2/29‑49 (MQ=255)                          caacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgt                       <  1:1766779‑M2/27‑1 (MQ=255)                          caacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgt                       >  1:667669‑M2/23‑49 (MQ=255)                          caacgcttactattggttatatTTCACCACTCATTGCCACATTCCTtgt                       <  1:1924821‑M2/27‑1 (MQ=255)                              gcttactattggttatatTTCACCACTCATTGCCACATTCCTTGTGta                    >  1:1059969‑M2/19‑48 (MQ=255)                                 tactattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                >  1:1865813‑M2/16‑49 (MQ=255)                                   ctattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagcc               <  1:23905‑M2/35‑1 (MQ=255)                                   ctattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                >  1:1887750‑M2/14‑47 (MQ=255)                                   ctattggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagc                >  1:782398‑M2/14‑47 (MQ=255)                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagcc           <  1:1390967‑M2/39‑1 (MQ=255)                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagcc           <  1:1233732‑M2/39‑1 (MQ=255)                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc            <  1:1726898‑M2/38‑1 (MQ=255)                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc            <  1:1551046‑M2/38‑1 (MQ=255)                                      ttggttatatTTCACCACTCATTGCCACATTCCTTGTGTATagccagc            >  1:476444‑M2/11‑48 (MQ=255)                                        ggttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAt         >  1:691733‑M2/9‑49 (MQ=255)                                          ttatatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttt       <  1:1739521‑M2/43‑1 (MQ=255)                                             tatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttttt    <  1:1296314‑M2/46‑1 (MQ=255)                                             tatTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAttttt     <  1:213274‑M2/45‑1 (MQ=255)                                              atTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTa   >  1:1872254‑M2/3‑49 (MQ=255)                                               tTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAc  >  1:699895‑M2/2‑49 (MQ=255)                                                |                                                GCGCCAAAACCTACAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAC  >  NC_000913/2336586‑2336680

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ 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.