Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	729,365	Δ9,598 bp		rhsC–ybfD	rhsC, ybfB, ybfO, ybfC, ybfQ, ybfL, ybfD

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	729365	738962	9598	38 [0]	[0] 36	rhsC–ybfD	rhsC,ybfB,ybfO,ybfC,ybfQ,ybfL,ybfD

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 729364	0 (0.000)	33 (1.150)	27/94	0.0	100%	intergenic (+24/‑219)	ybfA/rhsC	DUF2517 domain‑containing protein YbfA/rhs element protein RhsC
	?	NC_000913	738963 =	0 (0.000)					intergenic (+110/‑38)	ybfD/ybgA	H repeat‑associated putative transposase YbfD/DUF1722 domain‑containing protein YbgA

ATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACA  >  NC_000913/738916‑739011                                                |                                                 aggcaactggtgatttttattgattgctgagttttgaacacattgaaCGt                                                >  1:1859350‑M2/48‑50 (MQ=255)   gcaactggtgatttttattgattgctgagttttgaatacattgaaCGTCt                                              <  1:152676‑M2/5‑1 (MQ=255)     aactggtgatttttattgattgctgagttttgaatacattgaaCGTCTac                                            <  1:119853‑M2/7‑1 (MQ=255)         ggtgatttttattgattgctgagttttgaatacattgaaCGTCTACACt                                         >  1:2138943‑M2/40‑49 (MQ=255)          gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACtt                                        <  1:1695975‑M2/11‑1 (MQ=255)          gtgatttttattgattgctgagttttgaatacattgaaCGTCTACACtt                                        <  1:1662019‑M2/11‑1 (MQ=255)           tgatttttattgattgctgagttttgaatacattgaaCGTCTACACTTAc                                      <  1:347054‑M2/13‑1 (MQ=255)             atttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTc                                    >  1:1426788‑M2/36‑50 (MQ=255)               ttttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCt                                  >  1:224596‑M2/34‑50 (MQ=255)                tttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCt                                  <  1:1554798‑M2/17‑1 (MQ=255)                 ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                 >  1:1799192‑M2/32‑49 (MQ=255)                 ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                 >  1:962075‑M2/32‑49 (MQ=255)                 ttattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTg                                 >  1:1007535‑M2/32‑49 (MQ=255)                  tattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGt                                <  1:1005308‑M2/19‑1 (MQ=255)                  tattgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGt                                <  1:733978‑M2/19‑1 (MQ=255)                     tgattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAc                            >  1:493143‑M2/28‑50 (MQ=255)                      gattgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAAcc                           >  1:223303‑M2/27‑50 (MQ=255)                       attgctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCg                          >  1:1337400‑M2/26‑50 (MQ=255)                          gctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCt                        >  1:2170403‑M2/23‑49 (MQ=255)                           ctgagttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCa                      <  1:1413336‑M2/29‑1 (MQ=255)                              agttttgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAg                     <  1:1814632‑M2/30‑1 (MQ=255)                                   tgaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAg               >  1:466035‑M2/14‑49 (MQ=255)                                    gaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGt              >  1:1012958‑M2/13‑49 (MQ=255)                                    gaatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGt              >  1:47071‑M2/13‑49 (MQ=255)                                     aatacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTaa            <  1:1287706‑M2/39‑1 (MQ=255)                                       tacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAAt           >  1:1127380‑M2/10‑49 (MQ=255)                                       tacattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATg          <  1:1608978‑M2/41‑1 (MQ=255)                                        acattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATg          >  1:1302320‑M2/9‑49 (MQ=255)                                        acattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATg          >  1:922856‑M2/9‑49 (MQ=255)                                         cattgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGa         <  1:1660898‑M2/42‑1 (MQ=255)                                          attgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt       >  1:934184‑M2/7‑50 (MQ=255)                                          attgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAAt       <  1:100192‑M2/44‑1 (MQ=255)                                            tgaaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATc      >  1:770175‑M2/5‑49 (MQ=255)                                              aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTa    <  1:764644‑M2/47‑1 (MQ=255)                                              aaCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTAc   <  1:2137447‑M2/48‑1 (MQ=255)                                               aCGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACa  <  1:1719615‑M2/49‑1 (MQ=255)                                                |                                                 ATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCGCTCAGGAGCAGTAATGAATCTACA  >  NC_000913/738916‑739011

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: ‑