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

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

CAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTT  >  NC_000913/1299253‑1299729                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        cAAATGAAACGTAGCTGCCAGGCTCCGTTAAATTCGTCATTATTTTTACTTTGTTCCCGCGCAGTGATCAAAATCAAAAGGAAAAGGGAAAAATTTTCTGCTCAAATTACATTTAGTTATAAGTATTTCTTTAATGGTTTCTCTTTCCCCCGACATTATGT‑CGCCACATCAGAATCCTGGTCTTGCTGAAGAATAATGGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTaaataat                                                                                                                                                                                                                                    <  1:58683/250‑6 (MQ=255)       aaaCGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatca                                                                                                                                                                                                                               >  1:7407‑M1/1‑240 (MQ=255)        aaCGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaat                                                                                                                                                                                                                             <  1:393154‑M1/251‑13 (MQ=255)        aaCGGAGCTGCCATGCTCCGTTTACGTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGAGATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaat                                                                                                                                                                                                                             <  1:553297‑M1/251‑13 (MQ=255)          cGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattg                                                                                                                                                                                                                           >  2:70447‑M1/1‑237 (MQ=255)              gCTGCCATGATCCGTTTAACTCGTCGTTATTTTTACTTTGTTCCCGCGCCGTTATCAAAAGCAAACGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgtta                                                                                                                                                                                                                        >  2:498915‑M1/1‑233 (MQ=255)                 gCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaatt                                                                                                                                                                                                                    <  1:70447‑M1/251‑22 (MQ=255)                         ccGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtg                                                                                                                                                                                                              <  1:184483‑M1/249‑28 (MQ=255)                                cTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcact                                                                                                                                                                                                     >  2:34352‑M1/1‑215 (MQ=255)                                   cGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaatgattgtgcatttcactact                                                                                                                                                                                                  <  2:427855‑M1/251‑40 (MQ=255)                                                aCTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATGCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatca                                                                                                                                                                                     <  1:237627‑M1/251‑53 (MQ=255)                                                                      aaaagcaaaagGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTAGTTGCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACGTCCTTTGTTAGATGAAAGTGCTAAgtaagaatcaagtgttaaatgattgggcatttcacgtctggaactgtaatcagaaaagatagacatgcttagcc                                                                                                                                                               <  1:537259‑M1/251‑75 (MQ=255)                                                                       aaagcaaaagGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagcc                                                                                                                                                               >  2:39632‑M1/1‑176 (MQ=255)                                                                              aagGACTAGGTAAACATATTCTTCTCACATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctct                                                                                                                                                        >  1:173916‑M1/1‑169 (MQ=255)                                                                                       gTAAAAATATTCTTCTCAAATTTCAGTTAGTTATAAGGATTTCCTTAACTGATTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGGAGTATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTACataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattg                                                                                                                                               >  1:237353‑M1/1‑160 (MQ=255)                                                                                                  cttctCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagctg                                                                                                                                   >  1:57230‑M1/1‑149 (MQ=255)                                                                                                          aTTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgtta                                                                                                                            >  1:565840‑M1/1‑141 (MQ=255)                                                                                                               agttagttATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcat                                                                                                                      >  1:279745‑M1/1‑136 (MQ=255)                                                                                                                       ataAGGATTTCCTTACCTGCTTCTCCTCACCCTCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgca                                                                                                               >  2:91068‑M1/1‑128 (MQ=255)                                                                                                                               gtcCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGACGAATAATTGACATGATATTATTAATTCCACTGCCTTTGGTAGAGGATAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggtcgcttgcttt                                                                                                         >  2:85240‑M1/2‑120 (MQ=255)                                                                                                                                     aaCTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctga                                                                                                <  1:65590‑M1/251‑138 (MQ=255)                                                                                                                                     aaCTGCTTCTCATCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTCGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtactcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctg                                                                                                 >  2:509571‑M1/1‑114 (MQ=255)                                                                                                                                               cctcACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtca                                                                                        <  1:39632‑M1/249‑146 (MQ=255)                                                                                                                                               ccgcaccATCATGTTATGTTCGCCACATCATAAGCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgggcatttcactactgtaactagaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagaggtttggtaaggcatggatgcaagctatagattcggatacggtc                                                                                         <  2:237353‑M1/245‑145 (MQ=255)                                                                                                                                                         atGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaat                                                                            <  2:173916‑M1/251‑158 (MQ=255)                                                                                                                                                                  tCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcattgatgcaagctatagattctgatacggtcaataaaagataattgcttaac                                                                    <  2:225829‑M1/250‑166 (MQ=255)                                                                                                                                                                      cacaTCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattt                                                               <  1:34352‑M1/251‑171 (MQ=255)                                                                                                                                                                       acaTCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatgtcactactgtaactgtaatcagaaaagatagacatgcttagccaagctctatttgattgaattgaaagatgtttggtaaggcatggatgcaagctatagattctgatacggtcaagaaaagagaattgcttaacaatttt                                                              <  2:57230‑M1/251‑172 (MQ=255)                                                                                                                                                                                    gggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctcgatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtatt                                                  <  1:91068‑M1/250‑184 (MQ=255)                                                                                                                                                                                                             gATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaataagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtact                        >  2:325446‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                                    tGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgggcaaccgccagacgat  <  2:426150‑M1/250‑232 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        CAAATGAAACGGAGCTGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGCCGGCGACCTT  >  NC_000913/1299253‑1299729

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