Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,001,645	Δ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	4001645	4001649	5	62 [0]	[0] 63	corA	magnesium/nickel/cobalt transporter

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

GATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTA  >  NC_000913/4001509‑4001786                                                                                                                                        |                                                                                                                                               tatTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                                         <  2:92682‑M1/142‑8 (MQ=255) ggttgaTCGTGTTGAGCCGGAGGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAGCTGGAATACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                                         <  2:248047‑M1/141‑8 (MQ=255)   ttGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctt                                                                                                                                       <  1:738146‑M1/143‑10 (MQ=255)   ttGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctt                                                                                                                                       <  1:166836‑M1/143‑10 (MQ=255)   ttGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctt                                                                                                                                       >  1:469154‑M1/1‑134 (MQ=255)       tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCGTctccttcttctgc                                                                                                                                   >  1:12536‑M1/1‑130 (MQ=255)         ttcttgaaccgCGCGCCGACGAGCGACTGCGCGTACAACCTGAACTTGGCCACAGCCTGGCAACCCGCCCGGAACTGGAAGCCATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                                 <  1:334634‑M1/138‑16 (MQ=255)         ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGCGCCTGGCACCTCGCCCGGAACTGGAAGACATCGAAGCCTAGGCACGCTTCCTTGAAGACGACGACCGCCTGCATctccttcttcttctt                                                                                                                                 >  2:1129615‑M1/1‑128 (MQ=255)         gtgTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                                 >  1:437244‑M1/2‑128 (MQ=255)          tgccgaaCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                                <  1:40541‑M1/140‑17 (MQ=255)           gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                               <  1:97018‑M1/143‑18 (MQ=255)              gAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaag                                                                                                                            <  1:1097078‑M1/143‑21 (MQ=255)               aaCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGGAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGGCGGCCTGCATctccttcttcttctttgaaga                                                                                                                           >  1:91576‑M1/1‑122 (MQ=255)                  cGGTCGACGATGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                        <  1:733947‑M1/143‑25 (MQ=255)                    gacgacgacGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgg                                                                                                                      <  1:292823‑M1/143‑27 (MQ=255)                     acgacgacgAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgga                                                                                                                     <  1:791499‑M1/143‑28 (MQ=255)                       gacgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaag                                                                                                                   >  1:316040‑M1/1‑114 (MQ=255)                         cgacggGCGACTGCGCGTGCAAGATCACCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGTAAGACATCGTAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagat                                                                                                                 <  2:880705‑M1/143‑32 (MQ=255)                         cgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagat                                                                                                                 <  1:873257‑M1/143‑32 (MQ=255)                         cgacgaGCGACTGCGCGTACAAGCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagat                                                                                                                 <  1:1044806‑M1/143‑32 (MQ=255)                            cgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggcagctcac                                                                                                              >  2:262216‑M1/1‑109 (MQ=255)                              aGCGCCCGCGCGCACAGTGTGAACGTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTTAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgc                                                                                                            <  2:1082652‑M1/143‑37 (MQ=255)                               gCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcc                                                                                                           <  1:760288‑M1/143‑38 (MQ=255)                               gCGACTGCGCGCACAACCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCAGCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcc                                                                                                           <  1:809619‑M1/143‑38 (MQ=255)                                 gcctccgcgtacaaTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgg                                                                                                         <  2:1107449‑M1/138‑40 (MQ=255)                                 gACTGCGCGTACAATCTGAACGTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgg                                                                                                         <  2:1002277‑M1/143‑40 (MQ=255)                                 gACTGCGCCTACACTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgg                                                                                                         <  2:587582‑M1/143‑40 (MQ=255)                                   cTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggta                                                                                                       >  1:989551‑M1/1‑102 (MQ=255)                                    tGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaa                                                                                                      >  1:629023‑M1/1‑101 (MQ=255)                                      ctcgtaCAATCTGAGCTTGGGCAGCGCCTGGCAACCCGCCCGGAACTTGAAGACATCGAAGCATCGGCACGTGTCTTTGCAGACGACGACGGCCTGCATctcctttttcttctttgaagatgcggaagatcacgccggtaact                                                                                                    <  2:406173‑M1/141‑45 (MQ=255)                                         gcacaAGCTGACCGTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGCAGCCTCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactcca                                                                                                 <  2:241304‑M1/141‑48 (MQ=255)                                         gTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGCCGGCCCGCATctccttcttcggctttggcgttgcggaagatcacgccggtaactccc                                                                                                 >  1:210168‑M1/1‑96 (MQ=255)                                            cAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactg                                                                                              <  1:737410‑M1/143‑51 (MQ=255)                                              aTCTGAACGGGGCCAGAGCCTGGCAACCCGCCCGGAGCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                            <  1:64180‑M1/143‑53 (MQ=255)                                                 tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggca                                                                                         >  2:999442‑M1/1‑88 (MQ=255)                                                  gAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcat                                                                                        >  2:123980‑M1/1‑87 (MQ=255)                                                      ttGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAATACGATGTCGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttac                                                                                    >  1:162631‑M1/1‑83 (MQ=255)                                                             gagCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgt                                                                             >  1:470598‑M1/1‑76 (MQ=255)                                                                    gCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtc                                                                      <  2:537043‑M1/143‑75 (MQ=255)                                                                           gcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                               <  1:496787‑M1/143‑82 (MQ=255)                                                                              cgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttac                                                            <  1:999442‑M1/143‑85 (MQ=255)                                                                              cgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttac                                                            <  1:468128‑M1/143‑85 (MQ=255)                                                                               ggAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaacaaccgtgatggtcgtctgtttact                                                           >  2:975288‑M1/1‑58 (MQ=255)                                                                                      gAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcgtttaccatccgtgatggtggtctgtttactctgggtg                                                    >  2:721972‑M1/1‑51 (MQ=255)                                                                                       gagcacgcgaagcatCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                                   <  2:1057180‑M1/136‑94 (MQ=255)                                                                                       cagACACCGTATCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATccccttcttcttcgttgacgatgcggaagatcacgccggtaactccactgtggcacttaccatccgtgatggtcgtctgtttactctgcgtga                                                   <  2:1056252‑M1/142‑94 (MQ=255)                                                                                        aGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgag                                                  <  1:123980‑M1/143‑95 (MQ=255)                                                                                        aGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccacggtggcatttaccatccgtgatggtcgtctgtttactctgcgtgag                                                  >  2:578221‑M1/1‑49 (MQ=255)                                                                                             tCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtga                                             >  1:711501‑M1/1‑44 (MQ=255)                                                                                             tCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcggctgtttgctctgcgtgagcgtga                                             >  1:573774‑M1/1‑44 (MQ=255)                                                                                              cGAACCATCGGCACGTTCCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                            <  1:754892‑M1/143‑101 (MQ=255)                                                                                                aaGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttcctctgcgtgagcgtgacct                                          >  1:555598‑M1/1‑41 (MQ=255)                                                                                                   cATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtttgtttactctgcgtgagcgtgagctgcc                                       >  1:533692‑M1/1‑38 (MQ=255)                                                                                                      cGGCACGTTTCTTTGAAGACGACGCCGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcctgtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgc                                    <  2:346231‑M1/143‑109 (MQ=255)                                                                                                        gCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgagggtcgtctgtttactctgcgtgtgcgtgagctgccggctg                                  >  1:1082065‑M1/1‑33 (MQ=255)                                                                                                            gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtcactccactgtggcatttaccatccgtgacggtcgtctgtttactcggcgtcagcgagaaatgcccgcctttcg                              >  2:201923‑M1/1‑29 (MQ=255)                                                                                                             tttCTGTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgt                             <  1:611572‑M1/143‑116 (MQ=255)                                                                                                             ttgggttgaagaCGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgt                             <  1:146815‑M1/138‑116 (MQ=255)                                                                                                                    gAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaaccgcccgcgtttcgtctgtatc                      >  2:1056854‑M1/1‑21 (MQ=255)                                                                                                                       gacgacgacGGCCTGCATctccttcttcttctttgaagatgcggaagatgccgttattacgcccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgta                   <  1:262216‑M1/143‑126 (MQ=255)                                                                                                                                ggCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgctggtcgtctgtttactctgcgtgagcgtgaattgcccgcttttcgtatgtatcgtatgcgtgacc          >  2:960992‑M1/1‑9 (MQ=255)                                                                                                                                        tctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttgtcgtctgtatcgtatgcatgcccggacccag  >  2:334713‑M1/1‑1 (MQ=255)                                                                                                                                        |                                                                                                                                               GATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTA  >  NC_000913/4001509‑4001786

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