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	37 [0]	[0] 37	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.910)	24/268	0.5	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

TTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTT  >  NC_000913/1299379‑1299637                                                                                                                        |                                                                                                                                            ttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTTGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaat                                                                                                                         >  1:495118‑M1/1‑120 (MQ=255)        aCTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtg                                                                                                                  <  1:12787‑M1/140‑28 (MQ=255)        aCTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtg                                                                                                                  <  1:330621‑M1/140‑28 (MQ=255)            cTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatt                                                                                                              >  2:635194‑M1/1‑109 (MQ=255)              tctcCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttc                                                                                                            >  2:162944‑M1/1‑107 (MQ=255)                      ccATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgg                                                                                                    >  2:15106‑M1/1‑99 (MQ=255)                      ccATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactgg                                                                                                    >  1:332043‑M1/1‑99 (MQ=255)                          catGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtg                                                                                                                  <  1:191255‑M1/122‑28 (MQ=255)                          catGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtg                                                                                                                  >  2:191255‑M1/1‑95 (MQ=255)                                   ttCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcaga                                                                                       >  1:603858‑M1/1‑86 (MQ=255)                                            tcatAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagataga                                                                              <  2:495118‑M1/140‑64 (MQ=255)                                            tcatAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagataga                                                                              <  2:303270‑M1/140‑64 (MQ=255)                                            tcatAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagataga                                                                              <  1:431636‑M1/140‑64 (MQ=255)                                                aaTCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatg                                                                          <  1:162944‑M1/140‑68 (MQ=255)                                                  tCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgct                                                                        <  2:112099‑M1/140‑70 (MQ=255)                                                       ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagcc                                                                   >  1:10438‑M1/1‑66 (MQ=255)                                                       ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagcc                                                                   >  2:78250‑M1/1‑66 (MQ=255)                                                       ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagcc                                                                   >  2:98524‑M1/1‑66 (MQ=255)                                                              tGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctct                                                            >  1:87590‑M1/1‑59 (MQ=255)                                                                                  attattAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagat                                        >  1:146918‑M1/1‑39 (MQ=255)                                                                                             cACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaagg                             >  2:624437‑M1/1‑28 (MQ=255)                                                                                                tGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgtt                                 <  2:358419‑M1/133‑109 (MQ=255)                                                                                                tGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgtt                                 >  1:358419‑M1/1‑25 (MQ=255)                                                                                                  ccTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggca                           >  1:379720‑M1/1‑23 (MQ=255)                                                                                                  ccTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggca                           <  2:379720‑M1/137‑115 (MQ=255)                                                                                                    tttGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggat                      <  2:146918‑M1/140‑120 (MQ=255)                                                                                                    tttGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggat                      <  1:635194‑M1/140‑120 (MQ=255)                                                                                                    tttGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcat                          >  1:502773‑M1/1‑21 (MQ=255)                                                                                                    tttGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcat                          <  2:502773‑M1/136‑116 (MQ=255)                                                                                                     ttGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatg                     >  1:127680‑M1/1‑20 (MQ=255)                                                                                                            agGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagcta              >  1:118968‑M1/1‑13 (MQ=255)                                                                                                            agGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagcta              >  1:539667‑M1/1‑13 (MQ=255)                                                                                                             ggAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctat             <  2:332043‑M1/140‑129 (MQ=255)                                                                                                             ggAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctat             <  1:78250‑M1/140‑129 (MQ=255)                                                                                                                       aaataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgat   <  2:539667‑M1/140‑139 (MQ=255)                                                                                                                        aataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgata  >  2:38915‑M1/1‑1 (MQ=255)                                                                                                                        aataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgata  <  2:127680‑M1/140‑140 (MQ=255)                                                                                                                        |                                                                                                                                            TTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTT  >  NC_000913/1299379‑1299637

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