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	34 [0]	[1] 35	insH21	insH21

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

TGCTGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTAATG  >  NC_000913/1300453‑1300943                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      aaatgaatcggagctgccatggtccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcgccatcatgttattttcgccacatcataatccttggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataata                                                                                                                                                                                                                                                  <  2:332915‑M2/6‑1 (MQ=255)       aacggagctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAt                                                                                                                                                                                                                                            <  1:112995‑M2/12‑1 (MQ=255)              ctgccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaa                                                                                                                                                                                                                                      >  1:532965‑M2/233‑250 (MQ=255)                gccatgctccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAatt                                                                                                                                                                                                                                   <  2:218953‑M2/21‑1 (MQ=255)                        ccgtttgcttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccccataataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTtaaattattttg                                                                                                                                                                                                                             <  2:340218‑M2/27‑4 (MQ=255)                        ccgtttacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATGATAATCAATTGTTAAATTATtgtg                                                                                                                                                                                                                             <  1:363903‑M2/27‑1 (MQ=255)                             tacttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTc                                                                                                                                                                                                                       >  2:140106‑M2/218‑250 (MQ=255)                              acttcgtcattatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCa                                                                                                                                                                                                                      <  2:346484‑M2/34‑1 (MQ=255)                                        tatttttactttgttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAAc                                                                                                                                                                                                            <  1:168199‑M2/44‑1 (MQ=255)                                                    gttcccgcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaaa                                                                                                                                                                                                >  1:87701‑M2/195‑250 (MQ=255)                                                          gcgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGAt                                                                                                                                                                                            >  1:85972‑M2/189‑248 (MQ=255)                                                           cgcagttatcaaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGAc                                                                                                                                                                                        <  1:151472‑M2/64‑1 (MQ=255)                                                              agttatcaaaagcaaatggaatagggaaaaatattcttctcaaattaaagttagttataaggatttccttaaatgcttctcctcaccctcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattcttaattccactgcctttggtagaggaaagagctaaCTAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATg                                                                                                                                                                                     >  2:512685‑M2/185‑251 (MQ=255)                                                                 tatcaaaagcaaaaggaataggtaaaaatattcttctcaaattccagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCt                                                                                                                                                                                   <  2:421422‑M2/69‑1 (MQ=255)                                                                     aaaagcaaaagggataggtcaaaatattcttctcaaattacagttcgttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcatttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAACAGATATACATGCTTAGcc                                                                                                                                                                              >  2:454880‑M2/178‑251 (MQ=255)                                                                     aaaagcaaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGcc                                                                                                                                                                              <  2:164515‑M2/74‑1 (MQ=255)                                                                           aaaaggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAtctc                                                                                                                                                                        <  2:419953‑M2/80‑1 (MQ=255)                                                                              aggaataggtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTAt                                                                                                                                                                     >  2:305677‑M2/169‑251 (MQ=255)                                                                                      gtaaaaatattcttctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttatttgcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagagaaaagtgctaaATAATAATCAATTGTTAGATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATttgattg                                                                                                                                                              >  2:335573‑M2/161‑250 (MQ=255)                                                                                         aaaatattcttctcaaattacagttagttataagtattttcttaacttcttctcctcaccaggatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattagtaattccactgcctttggtagaggaaagtgctaaACAATAATCAATTGATAAATTATCGTGCATTTCACTCCTGGAAATGTCATCAGAAAAGATAGATATGCGTCGCTAATCTCTAGTTGattgaatt                                                                                                                                                          <  2:428968‑M2/94‑1 (MQ=255)                                                                                                      caaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtt                                                                                                                                                <  2:509697‑M2/104‑1 (MQ=255)                                                                                                                         aggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCt                                                                                                                          >  2:254526‑M2/126‑251 (MQ=255)                                                                                                                                     actgcgtctcctcaccatcatggtattttcgccacctcataatcctgggcttgctgaagaatcattgaaatgctgttattaattccactgcctttagtagagggaagtactaaATAATGATCAATTGTTAAATTATTGGGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAAATGTTTGTTAAGGCATGGATGCAAGCTATATATTCTGAt                                                                                                              >  2:214661‑M2/114‑251 (MQ=255)                                                                                                                                            ctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTc                                                                                                        >  2:438667‑M2/107‑250 (MQ=255)                                                                                                                                                    catcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATaaa                                                                                                  >  1:30838‑M2/99‑248 (MQ=255)                                                                                                                                                       catgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAgaga                                                                                             >  1:74941‑M2/96‑250 (MQ=255)                                                                                                                                                                     cacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAAttt                                                                              <  1:305677‑M2/170‑1 (MQ=255)                                                                                                                                                                                    ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTgg                                                               <  1:179203‑M2/185‑1 (MQ=255)                                                                                                                                                                                         tctgaagaataattaaaatgatattattaattccactgcctttggtagaggaaagtgctaaAGAATAATCAATTGTTAAATTATTGTGCATTTCTCTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTTAAAGATGTTTGTTAATGCATTGATGCAAGCTATAGATGCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAg                                                           <  1:470337‑M2/189‑1 (MQ=255)                                                                                                                                                                                            gaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATGCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAg                                                       <  1:214661‑M2/193‑1 (MQ=255)                                                                                                                                                                                                        aaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTgg                                           >  1:146256‑M2/47‑251 (MQ=255)                                                                                                                                                                                                                    taattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCggg                                <  1:74330‑M2/216‑1 (MQ=255)                                                                                                                                                                                                                     aattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGc                               >  1:40454‑M2/34‑250 (MQ=255)                                                                                                                                                                                                                           actgcctttggtagaggaaagtgctacATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGcc                        >  1:388178‑M2/28‑251 (MQ=255)                                                                                                                                                                                                                                                   tAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTAATg  >  1:528512/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      TGCTGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAACCGCCAGACGATTCTTTATTGGTAATG  >  NC_000913/1300453‑1300943

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