Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,755,661	Δ2,703 bp		pykF–[ldtE]	pykF, lpp, [ldtE]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1755661	1758363	2703	49 [0]	[1] 50	pykF–[ldtE]	pykF, lpp, [ldtE]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1755660	0 (0.000)	49 (1.140)	45/296	0.1	99.0%	intergenic (‑519/‑38)	ydhZ/pykF	uncharacterized protein/pyruvate kinase I
	?	NC_000913	1758364 =	1 (0.020)					coding (362/1005 nt)	ldtE	murein L,D‑transpeptidase

CAGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCG  >  NC_000913/1758218‑1758510                                                                                                                                                   |                                                                                                                                                   gactcaaaaacaggaactcacgttgggctgagacacaagcaaaaattcgtatgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgagcacagaaagcaagtGGAt                                                                                                                                                 >  1:422850‑M2/147‑150 (MQ=255)  actcaaaaacaggcactcacgttgggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATa                                                                                                                                                >  2:50395‑M2/146‑150 (MQ=255)         aacaggcactcacgttgggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGca                                                                                                                                        <  1:19812‑M2/13‑1 (MQ=255)                 ctcacgttgggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATAttt                                                                                                                                >  2:379606‑M2/131‑151 (MQ=255)                   cacgttgggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATattg                                                                                                                                <  2:421424‑M2/21‑2 (MQ=255)                   cacgttgggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATAtttt                                                                                                                               <  2:195941‑M2/22‑1 (MQ=255)                   cacgttgggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTtc                                                                                                                              <  2:624600‑M2/23‑1 (MQ=255)                      gttgggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTcc                                                                                                                           <  1:346056‑M2/26‑1 (MQ=255)                        cgggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagcGGATAAACCTGCACAATATTTTCTCCcg                                                                                                                         >  2:76296‑M2/124‑151 (MQ=255)                         gggctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCcg                                                                                                                         <  1:511426‑M2/28‑1 (MQ=255)                           gctgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCcggcg                                                                                                                      >  1:166377‑M2/121‑151 (MQ=255)                             tgagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCcggcgg                                                                                                                     >  2:192444‑M2/119‑150 (MQ=255)                              gagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCcggcgg                                                                                                                     <  1:231340‑M2/32‑1 (MQ=255)                              gagacacaagcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGa                                                                                                                    <  2:170615‑M2/33‑1 (MQ=255)                                      agcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGAtaataataa                                                                                                           <  1:36034‑M2/42‑1 (MQ=255)                                       gcacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATaaa                                                                                                          <  2:452798‑M2/43‑1 (MQ=255)                                        cacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCGCAATATTTTCTCCCGGCGGATAATAATAAAg                                                                                                         <  1:538595‑M2/44‑1 (MQ=255)                                        cacacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAg                                                                                                         <  2:601592‑M2/44‑1 (MQ=255)                                        acctcTGACGAATGCCTGCCGTATGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCTCCGGCGTTTCCAGCCCCTGCAACCCAATACCCATTCGATACACCTGCACACTACTTTCTCCCGGCGGatcatattt                                                                                                            <  2:313825/144‑4 (MQ=255)                                          cacattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGgc                                                                                                       <  1:611173‑M2/46‑1 (MQ=255)                                            cattcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGgcgc                                                                                                     <  1:138886‑M2/48‑1 (MQ=255)                                             attcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCa                                                                                                    <  2:340008‑M2/49‑1 (MQ=255)                                              ttcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGgcgc                                                                                                     >  1:543304‑M2/102‑149 (MQ=255)                                               tcctctgcacgctttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAg                                                                                                   <  2:355349‑M2/50‑1 (MQ=255)                                                        cgctttttcgatgtcacctatccttagagctaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTCTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAgg                                                                                         <  2:334627‑M2/60‑1 (MQ=255)                                                            ttttcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTa                                                                                      >  1:6901‑M2/88‑150 (MQ=255)                                                               tcgatgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACg                                                                                   >  2:343211‑M2/85‑150 (MQ=255)                                                                   tgtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAAt                                                                              >  2:390209‑M2/81‑151 (MQ=255)                                                                    gtcacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAAt                                                                              <  1:577118‑M2/71‑1 (MQ=255)                                                                      cacctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATccc                                                                           >  2:606451‑M2/78‑151 (MQ=255)                                                                       acctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATcccc                                                                          >  2:587704‑M2/77‑151 (MQ=255)                                                                        cctatccttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCt                                                                         <  2:190526‑M2/76‑1 (MQ=255)                                                                              cttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGAcgcg                                                                   >  1:424115‑M2/70‑151 (MQ=255)                                                                              cttagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGAcgcg                                                                   >  2:276765‑M2/70‑151 (MQ=255)                                                                                tagagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGAcgcg                                                                   <  1:226833‑M2/82‑1 (MQ=255)                                                                                 agagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGt                                                                 <  1:305845‑M2/84‑1 (MQ=255)                                                                                 agagcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTg                                                                <  1:272630‑M2/85‑1 (MQ=255)                                                                                    gcgaggcaccaccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCAt                                                             <  1:390220‑M2/88‑1 (MQ=255)                                                                                             caccactttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGtaat                                                     >  1:447334‑M2/55‑150 (MQ=255)                                                                                                 actttcgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAg                                                <  1:394240‑M2/101‑1 (MQ=255)                                                                                                      cgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTtgtg                                           <  2:402944‑M2/106‑1 (MQ=255)                                                                                                      cgtaataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTtgtg                                           >  2:120237‑M2/46‑151 (MQ=255)                                                                                                          ataccggattcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAgg                                       <  2:92858‑M2/110‑1 (MQ=255)                                                                                                                 attcgctttccggcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGtt                                <  1:665113‑M2/117‑1 (MQ=255)                                                                                                                             gcagtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTAcc                     >  2:579062‑M2/23‑150 (MQ=255)                                                                                                                               agtgcgcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAgg                  >  1:649396‑M2/21‑151 (MQ=255)                                                                                                                                    gcccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTg             >  2:334038‑M2/16‑151 (MQ=255)                                                                                                                                     cccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTgg            >  1:96054‑M2/15‑151 (MQ=255)                                                                                                                                     cccagaaagcaagtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTgg            <  1:513633‑M2/137‑1 (MQ=255)                                                                                                                                                agtGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCg  >  2:280927‑M2/4‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CAGGAACGACTGGCGGTAATTTAATGCCACGCTCCAGCGAACGCTGACGAATGCCTGCCGTAGGCGTCCAGGTTGGGTTAGGGATTTTCTGCCCAACACGCGTTTCCATCACCGGCGTTTCCAGCCCCTGCAATCCAATACCTATTGGATAAACCTGCACAATATTTTCTCCCGGCGGATAATAATAAAGGCGCAGCTCTGCAAGGTTAACGATAATCCCCTGACGCGGTGCATCAGGTAATAACAGTTGTGAAGGAATAGTTATCGTCGTACCAGGTTTTGGCACCGGGGCG  >  NC_000913/1758218‑1758510

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