Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,001,644	Δ5 bp	coding (220‑224/951 nt)	corA →	magnesium/nickel/cobalt transporter

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4001644	4001648	5	62 [0]	[0] 62	corA	magnesium/nickel/cobalt transporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4001643	0 (0.000)	62 (0.860)	53/284	0.2	100%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001649 =	0 (0.000)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

TGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGT  >  NC_000913/4001506‑4001784                                                                                                                                          |                                                                                                                                              tGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACACTCTGAACTTGGCCAGCGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGCCGGCCCGCAActcct                                                                                                                                          >  2:683111‑M1/1‑138 (MQ=255)  ggATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctt                                                                                                                                         >  2:1212291‑M1/1‑137 (MQ=255)    aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                                       >  2:715673‑M1/1‑135 (MQ=255)      tGATCGTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGACGACGACGACGGCCTGCATctccttcttc                                                                                                                                     <  2:386231‑M1/143‑11 (MQ=255)         tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCGGAACTTGGCCAGAGCCGGGCAACCCGCCCGGAACTGGAACACATCGAAGCATCGGCGCGTTTCTTTGGAGACGACGACGGCCTGCATcgccttctgcctc                                                                                                                                  >  2:409272‑M1/1‑130 (MQ=255)          cTTCTCGAACCCGACGACGCCGAGCGACTGCCCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGACGCATCGGCACGTTGCTTTGAAGACGACGACGGCCTGCATctccttcttcttct                                                                                                                                 <  1:215262‑M1/143‑15 (MQ=255)           ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGCAGCAGCGGCACGGTTGGTTGAAGACGGCGACGGCCTGCAGctccctctccggcct                                                                                                                                >  1:1018528‑M1/1‑128 (MQ=255)               cGAAACGGCCGAGGACGAGCGACTGCGCGTACAATCTGAACTGGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                            <  1:661605‑M1/143‑20 (MQ=255)                 aaCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaaga                                                                                                                          <  2:1052535‑M1/143‑22 (MQ=255)                   ccGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCTGGAACTGGAAGACATCGAAGCATCGGCACGTCTCTTTGAAGACGACGACGGCCTGCATctgcgtgttcttctttgacgatg                                                                                                                        >  1:1133154‑M1/1‑120 (MQ=255)                   ccGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatg                                                                                                                        <  2:1186589‑M1/143‑24 (MQ=255)                   ccGCACCACGACGACCGACCGCGCGCACAATCTCAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTCCATctccttcttcttctttgaagatg                                                                                                                        <  2:1127176‑M1/143‑24 (MQ=255)                    cgcacgacgacgagcGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                       <  1:773212‑M1/140‑25 (MQ=255)                    cGGACGACGGCGAGCGACAGCGCGTACAATCTGAACTTGGCCAGAGCCCGGCAACCCGCCCGGAACTGGAAGACAGCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                       <  2:495776‑M1/143‑25 (MQ=255)                    cGGACGACGACGAGCGACTGCGCGGACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                       <  2:666242‑M1/143‑25 (MQ=255)                     ggACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctccttcttctttggcgatgcg                                                                                                                      >  2:202617‑M1/1‑118 (MQ=255)                      gacgacgacGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgggg                                                                                                                     >  1:964486‑M1/1‑117 (MQ=255)                      gacgacgacGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCGGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgg                                                                                                                     <  1:877516‑M1/143‑27 (MQ=255)                                 gCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcc                                                                                                          >  2:29037‑M1/1‑106 (MQ=255)                                    cctGCGCGTACAATCTGAACTTGGCCATAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggt                                                                                                       <  2:731358‑M1/142‑41 (MQ=255)                                      tGCGCGGGCAATCTGAACTGGGCCAGAGGCTGGCGACCCGCGCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctgtgaagatgcggaagatcacgccggtaa                                                                                                     <  1:68455‑M1/143‑43 (MQ=255)                                       gcgcCTACAATCGGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCAGCGGCACGTTTCTTGGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaac                                                                                                    <  1:214243‑M1/143‑44 (MQ=255)                                          cgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactcc                                                                                                 <  2:969513‑M1/143‑47 (MQ=255)                                            tACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagcacgccggtaactgcac                                                                                               >  2:947607‑M1/1‑95 (MQ=255)                                                ttcttgatttgtccagagCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                           <  1:1070201‑M1/135‑53 (MQ=255)                                                cgcggaccttggccaGAGCCTGGCAACGCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                           <  1:423544‑M1/136‑53 (MQ=255)                                                  cTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggc                                                                                         <  1:943743‑M1/143‑55 (MQ=255)                                                   tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgcagatgcggaagataacgccggtaactccactgtggca                                                                                        >  1:993528‑M1/1‑88 (MQ=255)                                                      aCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgccgaagcacacggccgtgcttccaccgtggcattt                                                                                     >  1:1143805‑M1/1‑85 (MQ=255)                                                       cTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttc                                                                                    >  2:471647‑M1/1‑84 (MQ=255)                                                       cTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACAACAACGGACTGCATctccttcttcttatttgaagatgcggaagatcacgccggtaactccactgtggcattta                                                                                    >  2:1100110‑M1/1‑84 (MQ=255)                                                            ccAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcttttaccatc                                                                               >  2:613739‑M1/1‑79 (MQ=255)                                                             cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaaccccactgtggcatttaccatcc                                                                              >  1:1145525‑M1/1‑78 (MQ=255)                                                               gagCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgt                                                                            <  1:947938‑M1/143‑68 (MQ=255)                                                                 gCCTGGCAACCCGCCCGGAACGGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGGCGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtga                                                                          <  1:29250‑M1/143‑70 (MQ=255)                                                                       cAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcacggatggtaa                                                                    <  2:1228903‑M1/143‑76 (MQ=255)                                                                          cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                                 >  2:500053‑M1/1‑65 (MQ=255)                                                                          cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                                 <  2:805134‑M1/143‑79 (MQ=255)                                                                            cgcccgGAACTGGAAGACATCGAAGCACCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgt                                                               <  2:784713‑M1/143‑81 (MQ=255)                                                                             gcccgGAACTGCAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcgtcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                              <  2:358622‑M1/143‑82 (MQ=255)                                                                              cccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccattcgtgatggtcgtctgttt                                                             >  2:272605‑M1/1‑61 (MQ=255)                                                                               ctagaacTGGAAGACATCTAAGCATCCGCACGTTTCTGTGAAGACGACGACGGCCTGCATctccttctacttctttgaagatgctgaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttta                                                            <  2:1004528‑M1/140‑84 (MQ=255)                                                                               ccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttta                                                            >  1:787202‑M1/1‑60 (MQ=255)                                                                                  gAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttggaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactc                                                         <  2:792873‑M1/143‑87 (MQ=255)                                                                                      tGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctgtgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgggatggtcgtctgtttactctgcg                                                     >  1:607598‑M1/1‑53 (MQ=255)                                                                                       ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                    >  2:554249‑M1/1‑52 (MQ=255)                                                                                       ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                    <  1:715673‑M1/143‑92 (MQ=255)                                                                                                ggaggcatcGTCGCGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                           <  2:901367‑M1/139‑101 (MQ=255)                                                                                                cGAAGCATCGGCACGTTTCTTTGGAGACGGGGACCGCCTGCATccgcggcttcggctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccttgatggtcgtctgtttactctgcgtgagcgtgaa                                           <  2:170599‑M1/143‑101 (MQ=255)                                                                                                           cccgtTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttt                                <  2:1091216‑M1/141‑112 (MQ=255)                                                                                                              gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcg                             >  2:344019‑M1/1‑29 (MQ=255)                                                                                                                ttCTTTGAAGACGACCACGTCTTGCATctccttcttcttctttgaagatgcggcagatcacgccggtaactccactgccgcatttaccatccgtgatgttcgtctctttactctgcgtgagcgtgaactgcccgctcttcgtc                           <  2:189481‑M1/143‑117 (MQ=255)                                                                                                                    ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctgttcgtgtgta                       >  2:500580‑M1/1‑23 (MQ=255)                                                                                                                          acgacgccgGCCTGCATctccttcttcttcgttgcagatgcggaagattacgccggtcactccactgtggcatttaccatccgtgatggtcttctgttttctctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtat                 <  1:1240535‑M1/143‑127 (MQ=255)                                                                                                                          acgacgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgcactgcccgcttttcgtctgtatcgtat                 >  1:805428‑M1/1‑17 (MQ=255)                                                                                                                          acgacgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtat                 >  2:194045‑M1/1‑17 (MQ=255)                                                                                                                            gactacgGTCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc               <  2:1156118‑M1/143‑129 (MQ=255)                                                                                                                            gacgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc               >  1:644507‑M1/1‑15 (MQ=255)                                                                                                                                 cgGCCTGCTCcctctccttcttctttgaagctgcggaagttcacgcccgtaactccactgtggcacttaccacccgtgatggtcgtctgtttaccctgcgtgaccgtgaactgcccgcttttcgtctgtatcgtatgcgtgcc          <  2:804767‑M1/143‑134 (MQ=255)                                                                                                                                 cgGCCTGCATctcctgcttctgctttgaagatgcggaagatcccgagggtaactccactgcggcagttaccttccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcc          <  2:341200‑M1/143‑134 (MQ=255)                                                                                                                                  ggCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgccc         >  2:594386‑M1/1‑9 (MQ=255)                                                                                                                                         aTctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttctcgtctgtatcgtatgcgtgcccgtagcca  >  1:516361‑M1/1‑2 (MQ=255)                                                                                                                                          |                                                                                                                                              TGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGT  >  NC_000913/4001506‑4001784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑