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–1300696	1300697	2–1199	39 [0]	[0] 39	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	39 (1.380)	31/256	0.1	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

GCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTAT  >  NC_000913/1300572‑1300825                                                                                                                               |                                                                                                                                aaggatttccttaactgcttgtcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaggaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataataa                                                                                                                           >  1:94309‑M2/127‑133 (MQ=255) aaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataataa                                                                                                                           >  2:493727‑M2/127‑133 (MQ=255)                tgcttctcctcaccatcatgttattttcggcagatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgtctttggtagaggaaagtgctaaATAATAATCAATTGTTAAatta                                                                                                            <  2:94309‑M2/22‑1 (MQ=255)                  cttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAattatt                                                                                                          >  2:97565‑M2/110‑133 (MQ=255)                     ctcctcaccatcatgttattttggccacataataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                       <  1:249603‑M2/27‑1 (MQ=255)                        ctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAt                                                                                                    <  2:73248‑M2/30‑1 (MQ=255)                          caccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAttt                                                                                                  <  2:444357‑M2/32‑1 (MQ=255)                          caccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAttt                                                                                                  >  2:189384‑M2/102‑133 (MQ=255)                                  tgttattttcgccacatcatagtcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTc                                                                                                 <  2:197577‑M2/33‑1 (MQ=255)                                  tgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTc                                                                                                 >  1:197577‑M2/94‑126 (MQ=255)                                     tattttcgccacatcttaatccagggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGaa                                                                                       <  1:252030‑M2/43‑1 (MQ=255)                                           cgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTCAATTATTAGGCATTTCACTACTTGAACTGTaa                                                                                 >  2:252030‑M2/85‑133 (MQ=255)                                                 atcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaa                                                                           >  1:337763‑M2/79‑133 (MQ=255)                                                          ctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGAc                                                                  <  2:50050‑M2/64‑1 (MQ=255)                                                             ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCAGTACTGGACCTGTAATAAGAAATGATAGACATg                                                               >  2:168582‑M2/67‑133 (MQ=255)                                                             ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATg                                                               >  2:399957‑M2/67‑133 (MQ=255)                                                             ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATg                                                               >  2:305149‑M2/67‑133 (MQ=255)                                                             ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATg                                                               >  1:459132‑M2/67‑133 (MQ=255)                                                               cttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCt                                                             <  1:305149‑M2/69‑1 (MQ=255)                                                                  gctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAg                                                          >  2:49386‑M2/62‑133 (MQ=255)                                                                   ctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGc                                                         <  1:273691‑M2/73‑1 (MQ=255)                                                                            aattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTa                                                >  1:393856‑M2/52‑133 (MQ=255)                                                                                 aaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttga                                           <  1:399957‑M2/87‑1 (MQ=255)                                                                                       tattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGattgaat                                     >  1:329627‑M2/41‑133 (MQ=255)                                                                                                    actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttg                        >  2:396938‑M2/28‑133 (MQ=255)                                                                                                    actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTACATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttg                        >  1:406223‑M2/28‑133 (MQ=255)                                                                                                     ctgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgt                       >  2:483644‑M2/27‑133 (MQ=255)                                                                                                      tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGaaa                                <  2:360729‑M2/98‑1 (MQ=255)                                                                                                      tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGaaa                                >  1:360729‑M2/26‑123 (MQ=255)                                                                                                      tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgtt                      <  2:459132‑M2/108‑1 (MQ=255)                                                                                                      tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgtt                      <  1:97565‑M2/108‑1 (MQ=255)                                                                                                      tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgtt                      >  1:67086‑M2/26‑133 (MQ=255)                                                                                                           ttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTACATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGc                 <  2:406223‑M2/113‑1 (MQ=255)                                                                                                             ggtagaggaaagtgctaaATAAGAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCAt               <  1:189384‑M2/115‑1 (MQ=255)                                                                                                                  aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGAt           >  2:402416‑M2/14‑132 (MQ=255)                                                                                                                  aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGAt           <  1:402416‑M2/119‑1 (MQ=255)                                                                                                                  aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATg          <  1:396938‑M2/120‑1 (MQ=255)                                                                                                                  aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATg          >  1:71060‑M2/14‑133 (MQ=255)                                                                                                                          gctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCtat  <  2:67086‑M2/128‑1 (MQ=255)                                                                                                                               |                                                                                                                                GCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTAT  >  NC_000913/1300572‑1300825

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