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] 61	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)	61 (0.810)	54/284	0.3	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

GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTA  >  NC_000913/4001508‑4001774                                                                                                                                               |                                                                                                                             agtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatc                                                                                                                              <  2:797348‑M2/1‑1 (MQ=255)   tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatctc                                                                                                                            <  2:703842‑M2/3‑1 (MQ=255)   tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtgtctttgaagacgacgacggcctgcatctc                                                                                                                            >  1:7522‑M2/141‑143 (MQ=255)    atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                           <  1:753509‑M2/4‑1 (MQ=255)       gattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggcctgagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                        <  2:669482‑M2/7‑1 (MQ=255)       gattgatcttgtcgaaccggaagacgacgagcgactgcgagtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                        >  1:414708‑M2/137‑143 (MQ=255)        attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggacctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                       >  1:684204‑M2/136‑143 (MQ=255)        attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacggtcctttgaagaggacgacggcctgcgtGTCcttct                                                                                                                       >  2:607854‑M2/136‑143 (MQ=255)        attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggaccggttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                       >  1:924126‑M2/136‑143 (MQ=255)            atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttctt                                                                                                                   <  1:219160‑M2/12‑1 (MQ=255)              cttgtcgaaccggccgacgacgcgcgactgcgcgtacaatctgaccttggccagagcctggcaacccgcccggaactggaagacaccgaagcatcgtcaagtttattcgtagacgacgacggcctgcatCTCcttcttcttct                                                                                                                 <  2:972715‑M2/14‑1 (MQ=255)               ttgtcgagccggacgacgcccagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                <  2:1162386‑M2/15‑1 (MQ=255)               ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                >  1:1114386‑M2/129‑143 (MQ=255)               tggtcgaaccgggtgacgacgagcgactgcgcgtacaatctgagcttggccagagcctggcaacccgcccggaactggaagacagcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                <  2:1260104‑M2/15‑1 (MQ=255)               atgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacatggccagagcatggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                <  1:16606‑M2/15‑1 (MQ=255)                tgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                               <  1:1158522‑M2/16‑1 (MQ=255)                 ggcgaaccggacgacgacgagcgcctgcgcgtaaaatctgaacttggccagagcctggcaacccgcccggaacgggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                              <  1:1058914‑M2/17‑1 (MQ=255)                  tcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtccctgtgaagacgacgccggcctgcctCTCCTTCTTCTTCTTCGa                                                                                                             >  2:739189‑M2/126‑143 (MQ=255)                       ccggacgacgacgagcgactgcgcgtacaatctgaccttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATg                                                                                                        <  2:10623‑M2/23‑1 (MQ=255)                             gacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAg                                                                                                  >  1:139305‑M2/115‑143 (MQ=255)                                gacgagggacgtcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                               <  2:1072793‑M2/32‑1 (MQ=255)                                  cgcgcgacggcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCAc                                                                                             <  1:705609‑M2/34‑1 (MQ=255)                                      cgcctgctcgtacactctgaaggtggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                         <  2:511678‑M2/38‑1 (MQ=255)                                        actgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGt                                                                                       <  1:704153‑M2/40‑1 (MQ=255)                                                 acaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                              >  1:1327760‑M2/95‑143 (MQ=255)                                                     tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTgg                                                                          >  2:1059143‑M2/91‑143 (MQ=255)                                                       tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCa                                                                        <  1:182723‑M2/55‑1 (MQ=255)                                                         aacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCCTCTTCAGATTCGATGTGGAAGATCACGCCGCTAACTCCTCTATGGCAtt                                                                      >  1:371365‑M2/87‑143 (MQ=255)                                                           cctggcccgtgcctggccacccgcctggaactggaagacaccgaagcatcggcaagtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                    <  1:244694‑M2/59‑1 (MQ=255)                                                            ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACGGTGGCGTTTAc                                                                   >  1:228635‑M2/84‑143 (MQ=255)                                                                ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTCACTCCACTGTGGCATTTACCATc                                                               >  2:1229010‑M2/80‑143 (MQ=255)                                                                ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggactgcatCTCCTTCTTCTTGTTTGAAGATGCGGAAGATCACGCCGGTAGCTCCACTGTGGCATTTACCatg                                                               >  1:523122‑M2/80‑142 (MQ=255)                                                                 cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATcc                                                              >  1:1200453‑M2/79‑143 (MQ=255)                                                                  agagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgattctttgaagacgacgacggcctgcatCTCCTTCTTCTGCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATCtaccacgta                                                             >  1:359141‑M2/78‑139 (MQ=255)                                                                    agcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                           <  1:461307‑M2/68‑1 (MQ=255)                                                                      cctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTACCTCCACTGTGGAATTTACCATCCGTGAt                                                         >  2:189539‑M2/74‑143 (MQ=255)                                                                       ctcgcaacccgcccggaactgcaagacatcgaagcagcggcacgttccttcgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                        <  1:268109‑M2/71‑1 (MQ=255)                                                                          caaacccgcacgggactggaagacatcgaagcatcggcacgtgtctttgaagacgacgacggccggcatCTCCTTCTTCTCCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtc                                                     <  1:12521‑M2/74‑1 (MQ=255)                                                                              cccgcccgggactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCTACTGTGGCATTGACCATCCGTGATGGTCGTCt                                                 <  2:711964‑M2/78‑1 (MQ=255)                                                                                  cccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggactgctaCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTCACTCCACTGTGGCATTTACCATCCGTGATGGTcgtctcact                                             >  1:572663‑M2/62‑139 (MQ=255)                                                                                     cgaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAct                                          <  1:1202524‑M2/85‑1 (MQ=255)                                                                                          tggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcg                                     <  1:517277‑M2/90‑1 (MQ=255)                                                                                          tggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcg                                     >  2:536562‑M2/54‑143 (MQ=255)                                                                                            gaagacatcgatgcatgggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                   <  2:816080‑M2/92‑1 (MQ=255)                                                                                             aagacctcgaagcatcggctcgtctctttgaggacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                  <  2:182512‑M2/93‑1 (MQ=255)                                                                                              agacatcgaagctgcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCTCTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                 <  2:1258390‑M2/94‑1 (MQ=255)                                                                                                 gcgcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgt                              <  1:1112018‑M2/97‑1 (MQ=255)                                                                                                  atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGGCGTCTGTTTACTCTgcgtgagcgtg                             >  1:458918‑M2/46‑143 (MQ=255)                                                                                                  agcgaagccttggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtg                             <  1:1324895‑M2/98‑1 (MQ=255)                                                                                                    cgaatcctcggcatgtttctttgaagacgacgacggcctgcagCTCCTTCTTCTTCTTTGAAGATGCGGCAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                           <  2:785674‑M2/100‑1 (MQ=255)                                                                                                    cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                           <  1:1074868‑M2/100‑1 (MQ=255)                                                                                                    cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                           <  1:1176843‑M2/100‑1 (MQ=255)                                                                                                        gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGACCTGc                       >  1:1115111‑M2/40‑143 (MQ=255)                                                                                                         catcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGCAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                      >  2:1291949‑M2/39‑143 (MQ=255)                                                                                                           tcggcacgtttctttgatggcgccgacggcctgcatCTCCTTCTTCTTCGTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTCCCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCg                    <  2:1284100‑M2/107‑1 (MQ=255)                                                                                                               cacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGATGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGACTTCCCCGCttt                <  2:236405‑M2/111‑1 (MQ=255)                                                                                                                     tctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt          <  2:452638‑M2/117‑1 (MQ=255)                                                                                                                       tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTGTTCGTCTGt        >  1:91750‑M2/25‑143 (MQ=255)                                                                                                                        tgtaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCGCGCCGGTAACTCCTCTGTTGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTa       <  2:449025‑M2/120‑1 (MQ=255)                                                                                                                             gacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGGGGCATGTCCCATCCGTGGTGGTCGTCGGTTTACTCTGCGTGAGCGTGCACTGCCCGCTGGTCGTCTGTGTCGTa  >  1:515519‑M2/19‑143 (MQ=255)                                                                                                                             gacgacgacggcctgcacCTCCTTCTTCTTCTTTGACGATGCGGAGGATCACGCCGGTAACTCCACTGTGGCCTTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTa  <  2:367387‑M2/125‑1 (MQ=255)                                                                                                                                               |                                                                                                                             GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTA  >  NC_000913/4001508‑4001774

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