Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,970,427	(ACG)4→3	coding (295‑297/1131 nt)	wecB →	UDP‑N‑acetyl glucosamine‑2‑epimerase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3970418	3970420	3	62 [0]	[0] 62	wecB	UDP‑N‑acetyl glucosamine‑2‑epimerase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3970421 =	0 (0.000)	62 (1.250)	44/218	0.1	100%	coding (289/1131 nt)	wecB	UDP‑N‑acetyl glucosamine‑2‑epimerase
	?	NC_000913	= 3970417	0 (0.000)					coding (285/1131 nt)	wecB	UDP‑N‑acetyl glucosamine‑2‑epimerase

TCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCG  >  NC_000913/3970303‑3970558                                                                                                                       |                                                                                                                                          cgatctcaacataatgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGc                                                                                                                       <  2:307524‑M1/21‑1 (MQ=255)    tctcaacataatgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTg                                                                                                                    <  2:477589‑M1/24‑1 (MQ=255)      tcaacataatgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  >  1:650594‑M1/114‑139 (MQ=255)      tcaacataatgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  <  2:355986‑M1/26‑1 (MQ=255)       caacataatgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggcg                                                                                                                 <  1:407476‑M1/27‑1 (MQ=255)        aacataatgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggcgg                                                                                                                <  2:183830‑M1/28‑1 (MQ=255)           ataatgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGt                                                                                                             <  1:1045053‑M1/31‑1 (MQ=255)             aatgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGttt                                                                                                           <  1:388637‑M1/33‑1 (MQ=255)               tgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAAc                                                                                                                           >  1:267487‑M1/105‑121 (MQ=255)               tgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAAc                                                                                                                           <  2:267487‑M1/17‑1 (MQ=255)               tgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  >  1:879258‑M1/105‑130 (MQ=255)               tgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  <  2:879258‑M1/26‑1 (MQ=255)               tgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTa                                                                                                         <  1:873190‑M1/35‑1 (MQ=255)               tgcagccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTa                                                                                                         <  2:105916‑M1/35‑1 (MQ=255)                  agccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  <  1:898292‑M1/26‑1 (MQ=255)                  agccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  >  2:898292‑M1/102‑127 (MQ=255)                  agccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  <  1:697754‑M1/26‑1 (MQ=255)                  agccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCa                                                                                                      <  2:1018593‑M1/38‑1 (MQ=255)                  agccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCa                                                                                                      <  1:121253‑M1/38‑1 (MQ=255)                  agccaggacagggcctgacagagacaacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  >  2:697754‑M1/102‑127 (MQ=255)                    ccaggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGc                                                                                                    <  2:383360‑M1/40‑1 (MQ=255)                     caggacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCg                                                                                                   <  1:1031481‑M1/41‑1 (MQ=255)                        gacagggcctgacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTAt                                                                                                <  1:217094‑M1/44‑1 (MQ=255)                                  gacagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGt                                                                                      >  1:1104684‑M1/86‑139 (MQ=255)                                   acagagataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTc                                                                                     >  2:497726‑M1/85‑139 (MQ=255)                                       agataacctgtcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCAGTTGGTCACGt                                                                                 >  2:1031481‑M1/81‑139 (MQ=255)                                                 tcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTAt                                                                                                        >  2:732314‑M1/71‑106 (MQ=255)                                                 tcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTAt                                                                                                        <  1:732314‑M1/36‑1 (MQ=255)                                                 tcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGt                                                                       >  2:1045053‑M1/71‑139 (MQ=255)                                                 tcggattctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGt                                                                       <  2:88165‑M1/69‑1 (MQ=255)                                                        ctggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCa                                                                >  1:477589‑M1/64‑139 (MQ=255)                                                         tggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCAc                                                               <  2:771546‑M1/77‑1 (MQ=255)                                                         tggaagggctaaaacctattcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCAc                                                               >  2:407476‑M1/63‑139 (MQ=255)                                                                          attcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGc                                                           >  2:647968‑M1/46‑126 (MQ=255)                                                                          attcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGc                                                           <  1:647968‑M1/81‑1 (MQ=255)                                                                          attcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGc                                              <  2:307020‑M1/94‑1 (MQ=255)                                                                           ttcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGcc                                             <  2:806821‑M1/95‑1 (MQ=255)                                                                           ttcttgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGcc                                             >  1:1018593‑M1/45‑139 (MQ=255)                                                                               tgccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTgg                                         >  1:355986‑M1/41‑139 (MQ=255)                                                                                gccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGAtc                                                       <  2:169247‑M1/85‑1 (MQ=255)                                                                                gccgagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGAtc                                                       >  1:169247‑M1/40‑124 (MQ=255)                                                                                   gagttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCgg                                     >  1:105916‑M1/37‑139 (MQ=255)                                                                                      ttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  <  2:996303‑M1/26‑1 (MQ=255)                                                                                      ttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTggc                                                                                                                  >  1:996303‑M1/34‑59 (MQ=255)                                                                                      ttcaaaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGTGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAag                                  <  1:617548‑M1/106‑1 (MQ=255)                                                                                          aaccagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGc                              >  1:1080854‑M1/30‑139 (MQ=255)                                                                                           accagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAgg                               >  1:1057688‑M1/29‑137 (MQ=255)                                                                                           accagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAgg                               <  2:1057688‑M1/109‑1 (MQ=255)                                                                                           accagacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCt                             >  2:873190‑M1/29‑139 (MQ=255)                                                                                              agacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAAc                          >  2:702076‑M1/26‑139 (MQ=255)                                                                                              agacgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAAc                          >  1:383360‑M1/26‑139 (MQ=255)                                                                                                acgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGc                                                           >  2:155739‑M1/24‑104 (MQ=255)                                                                                                acgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGc                                                           <  1:155739‑M1/81‑1 (MQ=255)                                                                                                acgtcgtgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCg                        <  2:806763‑M1/116‑1 (MQ=255)                                                                                                      tgctggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACAtt                  <  1:1043344‑M1/122‑1 (MQ=255)                                                                                                         tggttcacggcgatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGAc               >  1:307524‑M1/15‑139 (MQ=255)                                                                                                                    gatACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTgg    >  1:461810‑M1/4‑139 (MQ=255)                                                                                                                     atACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTgg                                         >  2:770174‑M1/3‑101 (MQ=255)                                                                                                                     atACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTgg                                         <  1:770174‑M1/99‑1 (MQ=255)                                                                                                                     atACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGc   <  1:548875‑M1/137‑1 (MQ=255)                                                                                                                     atACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGc   >  2:217094‑M1/3‑139 (MQ=255)                                                                                                                      tACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCg  >  1:372629‑M1/2‑139 (MQ=255)                                                                                                                       |                                                                                                                                          TCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTGGCG  >  NC_000913/3970303‑3970558

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