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	1300697	1–1199	20 [0]	[2] 22	insH21	insH21

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1299498	0 (0.000)	20 (0.540)	18/480	0.8	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

ATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCC  >  NC_000913/1299272‑1299733                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            aTGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttcaattatt                                                                                                                                                                                                                     >  1:211973‑M1/1‑227 (MQ=255)                   cATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCGAAataataatcaattgttaaattattgcgcatttcactactgg                                                                                                                                                                                                    <  2:63102‑M1/250‑42 (MQ=255)                       attTTTACTTTGATCCCGCGCATTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGTATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTTACATGATATTATTAATTCCTCTGCCTTTGGTAGTCGATAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaact                                                                                                                                                                                                <  1:74522‑M1/250‑46 (MQ=255)                             aCTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGTATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAAGTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaatttattgtgcatttcactactggaactgtaatca                                                                                                                                                                                         <  1:372944‑M1/251‑53 (MQ=255)                                ttGTGCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCATGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAatattaatcaattgttaaattattgtgcatttcactactggaactgtaatcaga                                                                                                                                                                                       >  1:355235‑M1/1‑196 (MQ=255)                                    tCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaag                                                                                                                                                                                   <  2:293834‑M1/250‑59 (MQ=255)                                      tcgcgcCGTTATCAATAGCAAAAGGAATAGGTCAAACTATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATACTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTACACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattctgcatttcactacgggaaatgtaatcagaaaaact                                                                                                                                                                                 >  2:257193‑M1/2‑190 (MQ=255)                                                                 tAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgat                                                                                                                                                     >  1:36766‑M1/1‑163 (MQ=255)                                                                                         tACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATAATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaagg                                                                                                                             >  1:130785‑M1/1‑139 (MQ=255)                                                                                                   tataAGGCTTTCCTAAACTGCTTCTCCTCACCATCATGTTATTTTCCCCACATCATCATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGACTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaaactagacatgcttagccaatctctctttgattgaattgaaagatgtttgttaaggaatggatgc                                                                                                                    >  2:249586‑M1/1‑129 (MQ=255)                                                                                                            ttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctat                                                                                                             >  1:374893‑M1/1‑120 (MQ=255)                                                                                                                                        gTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaatt                                                                               <  1:67120‑M1/250‑159 (MQ=255)                                                                                                                                        gTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAAGTCCACTGCCTTTTGTAGAGGAAAGTGCTAAacaataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctattttattgaattgaaagatgtttgtaaaggcatggatgcaagctatagattcttatacggtcaataaaagagaatg                                                                               <  2:221596‑M1/250‑159 (MQ=255)                                                                                                                                           aTTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGCGCTAAataataatcaattgttaaattattgtgcatttcactactggaaatgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgttacggtcaataaaagagaattgctt                                                                           >  1:252052‑M1/1‑89 (MQ=255)                                                                                                                                                cGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgctcagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaaca                                                                       <  1:307722‑M1/250‑167 (MQ=255)                                                                                                                                                            aTCCTGTGCTTGCTGAAGAATAATTGAAATGATATTAATAATTCCACTGCCTGTGGTAGATGAAAGTGCTAAaaaataatcaattgttaaattattttgcagttcactactggaactgtaatcagaaaagatagacatgattagccaatctctatttgattgaatggaaagatttttgttaaggcatggatgcaagctatagatactgatacggtcaatataagataattgcttaacaaatgtacaaaatg                                                          <  2:55224‑M1/251‑180 (MQ=255)                                                                                                                                                                   gCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagatgctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggc                                                   <  1:33936‑M1/251‑187 (MQ=255)                                                                                                                                                                                                         aCTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaacttattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctatatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaccaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaaccgcc             >  1:365495‑M1/1‑27 (MQ=255)                                                                                                                                                                                                           tGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaaacgcc             >  2:112983‑M1/1‑25 (MQ=255)                                                                                                                                                                                                                    tAGAGGAAAGCGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaaccgccagacgattctt  >  1:387288‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            ATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTTATCC  >  NC_000913/1299272‑1299733

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