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	47 [0]	[0] 47	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	45 (1.300)	33/270	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

TAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATA  >  NC_000913/1300566‑1300836                                                                                                                                     |                                                                                                                                             agttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataataat                                                                                                                                       >  1:125032‑M2/133‑140 (MQ=255)  gttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATc                                                                                                                                      >  1:465723‑M2/132‑140 (MQ=255)   ttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCa                                                                                                                                     >  2:217509‑M2/131‑140 (MQ=255)                    actgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                                    <  2:230722‑M2/27‑1 (MQ=255)                    actgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                                    <  1:129624‑M2/27‑1 (MQ=255)                         ttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAttt                                                                                                               >  2:253289‑M2/109‑140 (MQ=255)                            tcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCac                                                                                                            <  2:45917‑M2/35‑1 (MQ=255)                            tcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCac                                                                                                            <  2:465723‑M2/35‑1 (MQ=255)                                    atcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGaa                                                                                                    <  1:195502‑M2/43‑1 (MQ=255)                                      catgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                                  <  1:226457‑M2/45‑1 (MQ=255)                                        tgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGt                                                                                                >  1:561378‑M2/94‑140 (MQ=255)                                            attttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATc                                                                                            <  1:232266‑M2/51‑1 (MQ=255)                                             ttttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCa                                                                                           >  2:204481‑M2/89‑140 (MQ=255)                                               ttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGa                                                                                         <  1:401027‑M2/54‑1 (MQ=255)                                                  gccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaaaa                                                                                      >  1:528434‑M2/84‑140 (MQ=255)                                                       atcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAg                                                                                 >  2:483832‑M2/79‑140 (MQ=255)                                                       atcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAg                                                                                 <  2:372477‑M2/62‑1 (MQ=255)                                                        tcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAGATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑a                                                                              <  1:217509‑M2/63‑1 (MQ=255)                                                                  gggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGc                                                                    <  1:353694‑M2/73‑1 (MQ=255)                                                                  gggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGc                                                                    <  2:430851‑M2/73‑1 (MQ=255)                                                                   ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGcc                                                                   >  2:4850‑M2/67‑140 (MQ=255)                                                                   ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGcc                                                                   >  2:500870‑M2/67‑140 (MQ=255)                                                                        gctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAAtct                                                              >  1:318680‑M2/62‑140 (MQ=255)                                                                        gctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAAtct                                                              >  2:426934‑M2/62‑140 (MQ=255)                                                                        gctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAAtct                                                              >  1:92452‑M2/62‑140 (MQ=255)                                                                                  aattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATttgatt                                                    <  2:528434‑M2/89‑1 (MQ=255)                                                                                            atattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAg                                          >  2:589682‑M2/42‑140 (MQ=255)                                                                                               ttattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtg                                       <  2:595001‑M2/102‑1 (MQ=255)                                                                                                       tccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTaa                               >  2:490518‑M2/31‑140 (MQ=255)                                                                                                           ctgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCa                           >  1:129597‑M2/27‑140 (MQ=255)                                                                                                           ctgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCa                           >  2:309444‑M2/27‑140 (MQ=255)                                                                                                           ctgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCa                           >  1:152307‑M2/27‑140 (MQ=255)                                                                                                            tgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCAt                          >  2:521606‑M2/26‑140 (MQ=255)                                                                                                                tttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGAt                      <  2:526348‑M2/119‑1 (MQ=255)                                                                                                                tttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGAt                      <  2:125032‑M2/119‑1 (MQ=255)                                                                                                                  tggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGc                    <  1:521606‑M2/121‑1 (MQ=255)                                                                                                                   ggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATg                     <  1:204481‑M2/122‑1 (MQ=255)                                                                                                                   ggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCa                   >  2:467216‑M2/19‑140 (MQ=255)                                                                                                                   ggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCa                   <  1:387745‑M2/122‑1 (MQ=255)                                                                                                                        aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGc                    <  2:458380‑M2/121‑1 (MQ=255)                                                                                                                        aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGc                    >  1:458380‑M2/14‑134 (MQ=255)                                                                                                                        aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCta              >  1:351038‑M2/14‑140 (MQ=255)                                                                                                                        aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCta              >  1:62669‑M2/14‑140 (MQ=255)                                                                                                                              gtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGAtt        >  1:144837‑M2/8‑140 (MQ=255)                                                                                                                               tgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTc       >  2:307123‑M2/7‑140 (MQ=255)                                                                                                                                    aATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATa  >  1:214202‑M2/2‑140 (MQ=255)                                                                                                                                    aATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATa  >  1:439830‑M2/2‑140 (MQ=255)                                                                                                                                     |                                                                                                                                             TAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAG‑‑ACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATA  >  NC_000913/1300566‑1300836

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