Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,299,499	Δ1,199 bp		insH21	insH21

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1299499	1300697	1199	45 [0]	[0] 45	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	42 (1.690)	34/88	0.0	100%	intergenic (+253/‑33)	ychE/insH21	putative inner membrane protein/IS5 transposase and trans‑activator
	?	NC_000913	1300698 =	0 (0.000)					intergenic (+151/‑484)	insH21/oppA	IS5 transposase and trans‑activator/oligopeptide ABC transporter periplasmic binding protein

ACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTA  >  NC_000913/1300650‑1300745                                                 |                                                tgaaatgatattattaattccactgcctttggtagaggaaagtgctaaa                                                 >  1:1860401‑M2/49‑49 (MQ=255)  gaaatgatattattaattccactgcctttggtagaggaaagtgctaaa                                                 >  1:1580184‑M2/48‑48 (MQ=255)   aaatgatattattaattccactgcctttggtagaggaaagtgctaaat                                                <  1:1528409‑M2/2‑1 (MQ=255)    aatgatattattaattccactgcctttggtagaggaaagtgctaaataa                                              >  1:1794989‑M2/46‑49 (MQ=255)    aatgatattattaattccactgcctttggtagaggaaagtgctaaataa                                              >  1:1553601‑M2/46‑49 (MQ=255)     atgatattattaattccactgcctttggtagaggaaagtgctaaataa                                              >  1:919579‑M2/45‑48 (MQ=255)     atgatattattaattccactgcctttggtagaggaaagtgctaaataa                                              >  1:355035‑M2/45‑48 (MQ=255)      tgatattattaattccactgcctttggcagaggaaagtgctaaataata                                            <  1:2033855‑M2/6‑1 (MQ=255)        atattattaattccactgcctttggtagaggaaagtgctaaataataa                                           <  1:101302‑M2/7‑1 (MQ=255)           ttattaattccactgcctttggtagaggaaagtgctaaATAATAATCaa                                       >  1:544256‑M2/39‑49 (MQ=255)           ttattaattccactgcctttggtagaggaaagtgctaaATAATAATCa                                        <  1:1675608‑M2/10‑1 (MQ=255)            tattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAt                                      >  1:1316674‑M2/38‑49 (MQ=255)            tattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAt                                      >  1:1160577‑M2/38‑49 (MQ=255)            tattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAt                                      >  1:590132‑M2/38‑49 (MQ=255)            tattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAt                                      >  1:1473570‑M2/38‑49 (MQ=255)             attaattccactgcctttggtagaggaaagtgctaaATAATAATCAAtt                                     >  1:1629209‑M2/37‑49 (MQ=255)              ttaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTg                                    >  1:1597020‑M2/36‑49 (MQ=255)               taattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGt                                   >  1:1352370‑M2/35‑49 (MQ=255)               taattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGt                                   >  1:1356247‑M2/35‑49 (MQ=255)                 attccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTa                                 <  1:1935781‑M2/17‑1 (MQ=255)                 attccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTa                                 <  1:2041811‑M2/17‑1 (MQ=255)                  ttccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTa                                 >  1:55094‑M2/32‑48 (MQ=255)                    ccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaaa                               <  1:25715‑M2/19‑1 (MQ=255)                    ccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAat                              >  1:1535759‑M2/30‑49 (MQ=255)                    ccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAat                              >  1:1673780‑M2/30‑49 (MQ=255)                     cactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAat                              >  1:1429256‑M2/29‑48 (MQ=255)                     cactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAat                              <  1:1456436‑M2/20‑1 (MQ=255)                      actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAatta                            >  1:125705‑M2/28‑49 (MQ=255)                      actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAatta                            <  1:1197137‑M2/22‑1 (MQ=255)                        tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAattatt                          >  1:87252‑M2/26‑49 (MQ=255)                          cctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgt                        >  1:1580227‑M2/24‑49 (MQ=255)                          cctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgt                        <  1:1438219‑M2/26‑1 (MQ=255)                            tttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                       >  1:1566594‑M2/22‑48 (MQ=255)                             ttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCa                     >  1:472062‑M2/21‑49 (MQ=255)                              tggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCa                     <  1:556261‑M2/29‑1 (MQ=255)                                gtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAttt                  >  1:1582773‑M2/18‑49 (MQ=255)                                  agaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCa                >  1:442336‑M2/16‑49 (MQ=255)                                      gaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCactac            <  1:1139050‑M2/38‑1 (MQ=255)                                        aagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTg          >  1:849319‑M2/10‑49 (MQ=255)                                        aagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTg          <  1:1134631‑M2/40‑1 (MQ=255)                                          gtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTgg         >  1:91115‑M2/8‑48 (MQ=255)                                           tgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGaa       <  1:1205375‑M2/43‑1 (MQ=255)                                             ctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt     >  1:492432‑M2/5‑49 (MQ=255)                                             ctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt     <  1:1506346‑M2/45‑1 (MQ=255)                                                aATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTa  <  1:1178256‑M2/48‑1 (MQ=255)                                                 |                                                ACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTA  >  NC_000913/1300650‑1300745

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