Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,361,226	Δ9,811 bp		[gltF]–[nanE]	[gltF], yhcA, yhcD, yhcE, insH10, yhcE, yhcF, yhcG, nanQ, nanK, [nanE]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3361226	3371036	9811	25 [0]	[0] 26	[gltF]–[nanE]	[gltF],yhcA,yhcD,yhcE,insH10,yhcE,yhcF,yhcG,nanQ,nanK,[nanE]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3361225	0 (0.000)	25 (0.410)	23/494	1.0	100%	coding (50/765 nt)	gltF	periplasmic protein GltF
	?	NC_000913	3371037 =	0 (0.000)					intergenic (‑1/+47)	nanE/nanT	putative N‑acetylmannosamine‑6‑phosphate 2‑epimerase/N‑acetylneuraminate:H(+) symporter

GAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCAGATC  >  NC_000913/3370798‑3371287                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           tatcaccacaaaccacctcaccagcacattattattcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCtt                                                                                                                                                                                                                                                 >  1:118125‑M2/240‑251 (MQ=255)  atcaccacaaaccacctcaccagcaccttattattcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCtt                                                                                                                                                                                                                                                 >  1:438910‑M2/239‑250 (MQ=255)                         caccttattattcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTATCGAAAGGCCCGGTACATAGAcc                                                                                                                                                                                                                         >  2:387015‑M2/216‑251 (MQ=255)                                 tcttcacatgatcgtaactataaaattcattaaaaaaatacataccagcaagttggctatttctgaacgcgttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTACCTTAGCGAAATGCCCGGAACCTAGACCGGGCAAc                                                                                                                                                                                                                  >  2:406090‑M2/208‑250 (MQ=255)                                      acatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGAtt                                                                                                                                                                                                            <  2:400355‑M2/49‑1 (MQ=255)                                                 aactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttattaaattcaccaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaaactcacaacagctgctatttgtgcagcactttcACATCTTCCCTTTGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAAACTTTGGtt                                                                                                                                                                                                 >  1:342120‑M2/192‑251 (MQ=255)                                                            cattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATc                                                                                                                                                                                      <  2:434370‑M2/71‑1 (MQ=255)                                                                                      tggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcacg                                                                                                                                                             <  1:220503‑M2/96‑1 (MQ=255)                                                                                                     gagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCCCCGCTGAATGGtt                                                                                                                                             >  2:21815‑M2/140‑251 (MQ=255)                                                                                                     gagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGtt                                                                                                                                             >  1:403748‑M2/140‑251 (MQ=255)                                                                                                                   cattcacaaacagctatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCg                                                                                                                               >  1:96428‑M2/126‑251 (MQ=255)                                                                                                                       cacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCAt                                                                                                                           >  1:234736‑M2/122‑251 (MQ=255)                                                                                                                                 tatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGATGAATGGTTTACCGTCGATAGCGGCATGAGTACGCa                                                                                                                  >  2:237947‑M2/112‑250 (MQ=255)                                                                                                                                        actgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTcc                                                                                                          <  1:383127‑M2/147‑1 (MQ=255)                                                                                                                                                     tttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagcttctacttgtgcagcactttcACATCTTCCCTTAGCGAAAGGTCCGGTACATAGACCGGGCAACAGGATTATCTTTTGGTGTTGATTAAAGCGTTGTTGGCGCTGCCAATCGGCACGGCACCGCTGAATGGTGTACCTTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaaccaa                                                                                              <  2:363575‑M2/159‑1 (MQ=255)                                                                                                                                                                aaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcgcg                                                                                   >  1:28016‑M2/81‑250 (MQ=255)                                                                                                                                                                      atttgttattaaggatatgttcatatgtttttcaaacagaaactcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCAGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGCATGGTTTACAGTCGATAGCGTCATGAGTACGCAACGCTTCCGCGCGCAACCAATGCGGAACGCAAGTAGGc                                                                            >  2:7501‑M2/75‑251 (MQ=255)                                                                                                                                                                                 aggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGGCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGccc                                                                 <  2:118125‑M2/188‑1 (MQ=255)                                                                                                                                                                                    atatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAAt                                                              >  2:499147‑M2/61‑251 (MQ=255)                                                                                                                                                                                           catatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAgg                                                       >  2:225343‑M2/54‑251 (MQ=255)                                                                                                                                                                                                 tttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACGTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATcac                                                  <  2:276827‑M2/203‑1 (MQ=255)                                                                                                                                                                                                      caaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGa                                             >  1:58062‑M2/43‑250 (MQ=255)                                                                                                                                                                                                       aaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAAc                                           >  1:260239‑M2/42‑251 (MQ=255)                                                                                                                                                                                                               cctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACt                                   >  1:144127‑M2/34‑251 (MQ=255)                                                                                                                                                                                                                       cagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACg                              >  2:385419‑M2/26‑248 (MQ=255)                                                                                                                                                                                                                                                acatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCAGATc  <  1:71159/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           GAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCAGATC  >  NC_000913/3370798‑3371287

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.