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

CGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTG  >  NC_000913/3970300‑3970555                                                                                                                      |                                                                                                                                           cGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagc                                                                                                                       <  2:307524‑M2/139‑22 (MQ=255)    tctcAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctg                                                                                                                    <  2:477589‑M2/139‑25 (MQ=255)      tcAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  >  1:650594‑M2/1‑113 (MQ=255)      tcAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  <  2:355986‑M2/139‑27 (MQ=255)       cAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcg                                                                                                                 <  1:407476‑M2/139‑28 (MQ=255)        aaCATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcgg                                                                                                                <  2:183830‑M2/139‑29 (MQ=255)           aTAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgt                                                                                                             <  1:1045053‑M2/139‑32 (MQ=255)             aaTGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttt                                                                                                           <  1:388637‑M2/139‑34 (MQ=255)               tGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgtttta                                                                                                         <  2:105916‑M2/139‑36 (MQ=255)               tGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgtttta                                                                                                         <  1:873190‑M2/139‑36 (MQ=255)               tGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  <  2:879258‑M2/130‑27 (MQ=255)               tGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  >  1:879258‑M2/1‑104 (MQ=255)               tGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaac                                                                                                                           <  2:267487‑M2/121‑18 (MQ=255)               tGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaac                                                                                                                           >  1:267487‑M2/1‑104 (MQ=255)                  aGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatca                                                                                                      <  2:1018593‑M2/139‑39 (MQ=255)                  aGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatca                                                                                                      <  1:121253‑M2/139‑39 (MQ=255)                  aGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  >  2:898292‑M2/1‑101 (MQ=255)                  aGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  <  1:898292‑M2/127‑27 (MQ=255)                  aGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  <  1:697754‑M2/127‑27 (MQ=255)                  aGCCAGGACAGGGCCTGACAGAGACAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  >  2:697754‑M2/1‑101 (MQ=255)                    ccAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagc                                                                                                    <  2:383360‑M2/139‑41 (MQ=255)                     cAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcg                                                                                                   <  1:1031481‑M2/139‑42 (MQ=255)                        gACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtat                                                                                                <  1:217094‑M2/139‑45 (MQ=255)                                  gACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggt                                                                                      >  1:1104684‑M2/1‑85 (MQ=255)                                   aCAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtc                                                                                     >  2:497726‑M2/1‑84 (MQ=255)                                       agaTAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattccagttggtcacgt                                                                                 >  2:1031481‑M2/1‑80 (MQ=255)                                                 tCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggt                                                                       <  2:88165‑M2/139‑70 (MQ=255)                                                 tCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggt                                                                       >  2:1045053‑M2/1‑70 (MQ=255)                                                 tCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttat                                                                                                        <  1:732314‑M2/106‑37 (MQ=255)                                                 tCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttat                                                                                                        >  2:732314‑M2/1‑70 (MQ=255)                                                        cTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgca                                                                >  1:477589‑M2/1‑63 (MQ=255)                                                         tGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcac                                                               <  2:771546‑M2/139‑78 (MQ=255)                                                         tGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcac                                                               >  2:407476‑M2/1‑62 (MQ=255)                                                                          aTTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgc                                              <  2:307020‑M2/139‑95 (MQ=255)                                                                          aTTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggc                                                           <  1:647968‑M2/126‑82 (MQ=255)                                                                          aTTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggc                                                           >  2:647968‑M2/1‑45 (MQ=255)                                                                           ttCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgcc                                             <  2:806821‑M2/139‑96 (MQ=255)                                                                           ttCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgcc                                             >  1:1018593‑M2/1‑44 (MQ=255)                                                                               tGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtgg                                         >  1:355986‑M2/1‑40 (MQ=255)                                                                                gCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatc                                                       >  1:169247‑M2/1‑39 (MQ=255)                                                                                gCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatc                                                       <  2:169247‑M2/124‑86 (MQ=255)                                                                                   gAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccgg                                     >  1:105916‑M2/1‑36 (MQ=255)                                                                                      ttCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggtggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaag                                  <  1:617548‑M2/139‑107 (MQ=255)                                                                                      ttCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  <  2:996303‑M2/59‑27 (MQ=255)                                                                                      ttCAAACCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggc                                                                                                                  >  1:996303‑M2/1‑33 (MQ=255)                                                                                          aaCCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggc                              >  1:1080854‑M2/1‑29 (MQ=255)                                                                                           aCCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggct                             >  2:873190‑M2/1‑28 (MQ=255)                                                                                           aCCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagagg                               <  2:1057688‑M2/137‑110 (MQ=255)                                                                                           aCCAGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagagg                               >  1:1057688‑M2/1‑28 (MQ=255)                                                                                              aGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaac                          >  2:702076‑M2/1‑25 (MQ=255)                                                                                              aGACGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaac                          >  1:383360‑M2/1‑25 (MQ=255)                                                                                                aCGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaaccg                        <  2:806763‑M2/139‑117 (MQ=255)                                                                                                aCGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggc                                                           >  2:155739‑M2/1‑23 (MQ=255)                                                                                                aCGTCGTGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggc                                                           <  1:155739‑M2/104‑82 (MQ=255)                                                                                                      tGCTGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaaccgtacatt                  <  1:1043344‑M2/139‑123 (MQ=255)                                                                                                         tGGTTCACGGCGATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaaccgtacattgac               >  1:307524‑M2/1‑14 (MQ=255)                                                                                                                    gATacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaaccgtacattgaccgggcatctgg    >  1:461810‑M2/1‑3 (MQ=255)                                                                                                                     aTacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaaccgtacattgaccgggcatctggc   <  1:548875‑M2/139‑138 (MQ=255)                                                                                                                     aTacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaaccgtacattgaccgggcatctggc   >  2:217094‑M2/1‑2 (MQ=255)                                                                                                                     aTacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtgg                                         >  2:770174‑M2/1‑2 (MQ=255)                                                                                                                     aTacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtgg                                         <  1:770174‑M2/101‑100 (MQ=255)                                                                                                                      tacgacgacgctggcaaccagcctggcggcgttttatcagcgtattcctgttggtcacgttgaggctggtctgcgcacgggcgatctctattcgccgtggccggaagaggctaaccgtacattgaccgggcatctggcg  >  1:372629‑M2/1‑1 (MQ=255)                                                                                                                      |                                                                                                                                           CGATCTCAACATAATGCAGCCAGGACAGGGCCTGACAGAGATAACCTGTCGGATTCTGGAAGGGCTAAAACCTATTCTTGCCGAGTTCAAACCAGACGTCGTGCTGGTTCACGGCGATACGACGACGACGCTGGCAACCAGCCTGGCGGCGTTTTATCAGCGTATTCCTGTTGGTCACGTTGAGGCTGGTCTGCGCACGGGCGATCTCTATTCGCCGTGGCCGGAAGAGGCTAACCGTACATTGACCGGGCATCTG  >  NC_000913/3970300‑3970555

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: ‑