Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,761,993	Δ8,250 bp		rhsA–yibV	rhsA, yibA, yibJ, yibG, yibS, yibV

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3761993	3770242	8250	37 [0]	[0] 38	rhsA–yibV	rhsA,yibA,yibJ,yibG,yibS,yibV

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3761992	0 (0.000)	36 (1.360)	33/94	0.0	100%	intergenic (‑37/‑191)	yibF/rhsA	glutathione transferase‑like protein YibF/rhs element protein RhsA
	?	NC_000913	3770243 =	0 (0.000)					coding (1137/1137 nt)	yibH	inner membrane protein YibH

GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3761944‑3761992 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAG  >  NC_000913/3770243‑3770287                                                                                                GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                               <  1:1267738/49‑1  ACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCT                                              >  1:363884/1‑49   CGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTC                                             <  1:1697730/49‑1    GAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCA                                            <  1:1025681/49‑1     AGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAA                                           >  1:1196848/1‑49     AGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAA                                           <  1:1438526/49‑1      GTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAAT                                          <  1:36347/49‑1       TTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAAT                                          >  1:1425177/1‑48         TCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGA                                        <  1:339138/48‑1           ATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATC                                      >  1:1484055/1‑48            TGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCA                                    <  1:286712/49‑1            TGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCC                                     >  1:1952003/1‑48                 AGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAA                               <  1:932747/49‑1                 AGTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAA                               <  1:1426893/49‑1                  GTCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAA                              <  1:72708/49‑1                   TCCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAG                             <  1:491815/49‑1                    CCTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATA                                >  1:1963002/1‑45                     CTTTGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGA                            <  1:1005276/48‑1                        TGTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAA                         <  1:268759/48‑1                         GTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATG                       >  1:120261/1‑49                         GTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATG                       >  1:370726/1‑49                         GTGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAAT                        <  1:420422/48‑1                          TGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATGC                      <  1:1172526/49‑1                          TGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATG                       >  1:1965684/1‑48                          TGCGAGGAAAAATATCAGTATGGCTCAATGATCCAAATAAAGATAATG                       <  1:1017499/48‑1                                    AATATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGG            <  1:1608905/49‑1                                     ATATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGT           <  1:574668/49‑1                                       ATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCA         >  1:154736/1‑49                                       ATCAGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCA         >  1:1793857/1‑49                                        TCAGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCAT        >  1:129895/1‑49                                          AGTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTC      >  1:358447/1‑49                                           GTATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCT     >  1:462325/1‑49                                            TATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTT    <  1:155906/49‑1                                             ATGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCT     >  1:348495/1‑47                                              TGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAG  >  1:1003804/1‑49                                              TGGCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTA   <  1:445356/48‑1                                                GCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAG  >  1:2004733/1‑47                                                                                                GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3761944‑3761992 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAG  >  NC_000913/3770243‑3770287

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