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

CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTA  >  NC_000913/4001502‑4001769                                                                                                                                               |                                                                                                                              cAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGCAGCATCGGCACGTTTCTTTGCAGACGACGACGGCCTGCat                                                                                                                               >  2:1056352/1‑143 (MQ=255)  aGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATc                                                                                                                              <  2:797348‑M1/143‑2 (MQ=255)    tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctc                                                                                                                            <  2:703842‑M1/143‑4 (MQ=255)    tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTGTCTTTGAAGACGACGACGGCCTGCATctc                                                                                                                            >  1:7522‑M1/1‑140 (MQ=255)     aTGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                           <  1:753509‑M1/143‑5 (MQ=255)        gATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCTGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                        <  2:669482‑M1/143‑8 (MQ=255)        gATTGATCTTGTCGAACCGGAAGACGACGAGCGACTGCGAGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                        >  1:414708‑M1/1‑136 (MQ=255)         aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGACCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                       >  1:684204‑M1/1‑135 (MQ=255)         aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTCCTTTGAAGAGGACGACGGCCTGCGTgtccttct                                                                                                                       >  2:607854‑M1/1‑135 (MQ=255)         aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGACCGGTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                       >  1:924126‑M1/1‑135 (MQ=255)             aTCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctt                                                                                                                   <  1:219160‑M1/143‑13 (MQ=255)               cTTGTCGAACCGGCCGACGACGCGCGACTGCGCGTACAATCTGACCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACACCGAAGCATCGTCAAGTTTATTCGTAGACGACGACGGCCTGCATctccttcttcttct                                                                                                                 <  2:972715‑M1/143‑15 (MQ=255)                ttGTCGAGCCGGACGACGCCCAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                <  2:1162386‑M1/143‑16 (MQ=255)                ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                >  1:1114386‑M1/1‑128 (MQ=255)                tggtcGAACCGGGTGACGACGAGCGACTGCGCGTACAATCTGAGCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACAGCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                <  2:1260104‑M1/141‑16 (MQ=255)                atgTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACATGGCCAGAGCATGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                <  1:16606‑M1/142‑16 (MQ=255)                 tGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                               <  1:1158522‑M1/143‑17 (MQ=255)                  ggcgaACCGGACGACGACGAGCGCCTGCGCGTAAAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACGGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                              <  1:1058914‑M1/141‑18 (MQ=255)                   tCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTCCCTGTGAAGACGACGCCGGCCTGCCTctccttcttcttcttcga                                                                                                             >  2:739189‑M1/1‑125 (MQ=255)                        ccGGACGACGACGAGCGACTGCGCGTACAATCTGACCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatg                                                                                                        <  2:10623‑M1/143‑24 (MQ=255)                              gacgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaag                                                                                                  >  1:139305‑M1/1‑114 (MQ=255)                                 gacgaGGGACGTCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatc                                                                                               <  2:1072793‑M1/143‑33 (MQ=255)                                   cgcgcgaCGGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcac                                                                                             <  1:705609‑M1/140‑35 (MQ=255)                                       cgcctgcTCGTACACTCTGAAGGTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccg                                                                                         <  2:511678‑M1/140‑39 (MQ=255)                                         aCTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggt                                                                                       <  1:704153‑M1/143‑41 (MQ=255)                                                  aCAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                              >  1:1327760‑M1/1‑94 (MQ=255)                                                      tCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtgg                                                                          >  2:1059143‑M1/1‑90 (MQ=255)                                                        tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggca                                                                        <  1:182723‑M1/143‑56 (MQ=255)                                                          aaCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcctcttcagattcgatgtggaagatcacgccgctaactcctctatggcatt                                                                      >  1:371365‑M1/1‑86 (MQ=255)                                                            cctggCCCGTGCCTGGCCACCCGCCTGGAACTGGAAGACACCGAAGCATCGGCAAGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattta                                                                    <  1:244694‑M1/141‑60 (MQ=255)                                                             ttGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccacggtggcgtttac                                                                   >  1:228635‑M1/1‑83 (MQ=255)                                                                 ccAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtcactccactgtggcatttaccatc                                                               >  2:1229010‑M1/1‑79 (MQ=255)                                                                 ccAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGACTGCATctccttcttcttgtttgaagatgcggaagatcacgccggtagctccactgtggcatttaccatg                                                               >  1:523122‑M1/1‑79 (MQ=255)                                                                  cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcc                                                              >  1:1200453‑M1/1‑78 (MQ=255)                                                                   agagCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGATTCTTTGAAGACGACGACGGCCTGCATctccttcttctgctttgaagatgcggaagatcacgccggtaactccactgtggcatctaccacgta                                                             >  1:359141‑M1/1‑77 (MQ=255)                                                                     agCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtg                                                           <  1:461307‑M1/143‑69 (MQ=255)                                                                       ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtacctccactgtggaatttaccatccgtgat                                                         >  2:189539‑M1/1‑73 (MQ=255)                                                                        ctcgcaaCCCGCCCGGAACTGCAAGACATCGAAGCAGCGGCACGTTCCTTCGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatg                                                        <  1:268109‑M1/140‑72 (MQ=255)                                                                           caaaccCGCACGGGACTGGAAGACATCGAAGCATCGGCACGTGTCTTTGAAGACGACGACGGCCGGCATctccttcttctcctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtc                                                     <  1:12521‑M1/141‑75 (MQ=255)                                                                               cccgcccgGGACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactctactgtggcattgaccatccgtgatggtcgtct                                                 <  2:711964‑M1/143‑79 (MQ=255)                                                                                   cccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGACTGCTActccttcttcttctttgaagatgcggaagatcacgccggtcactccactgtggcatttaccatccgtgatggtcgtctcact                                             >  1:572663‑M1/1‑61 (MQ=255)                                                                                      cgaACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttact                                          <  1:1202524‑M1/142‑86 (MQ=255)                                                                                           tGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                     <  1:517277‑M1/143‑91 (MQ=255)                                                                                           tGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                     >  2:536562‑M1/1‑53 (MQ=255)                                                                                             gAAGACATCGATGCATGGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtg                                   <  2:816080‑M1/143‑93 (MQ=255)                                                                                              aaGACCTCGAAGCATCGGCTCGTCTCTTTGAGGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                  <  2:182512‑M1/143‑94 (MQ=255)                                                                                               aGACATCGAAGCTGCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactcctctgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgag                                 <  2:1258390‑M1/143‑95 (MQ=255)                                                                                                  gcgcgaaGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgt                              <  1:1112018‑M1/140‑98 (MQ=255)                                                                                                   agcgaAGCCTTGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtg                             <  1:1324895‑M1/141‑99 (MQ=255)                                                                                                   aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgggcgtctgtttactctgcgtgagcgtg                             >  1:458918‑M1/1‑45 (MQ=255)                                                                                                     cGAATCCTCGGCATGTTTCTTTGAAGACGACGACGGCCTGCAGctccttcttcttctttgaagatgcggcagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                           <  2:785674‑M1/143‑101 (MQ=255)                                                                                                     cGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                           <  1:1074868‑M1/143‑101 (MQ=255)                                                                                                     cGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                           <  1:1176843‑M1/143‑101 (MQ=255)                                                                                                         gCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgacctgc                       >  1:1115111‑M1/1‑39 (MQ=255)                                                                                                          cATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggcagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcc                      >  2:1291949‑M1/1‑38 (MQ=255)                                                                                                            tCGGCACGTTTCTTTGATGGCGCCGACGGCCTGCATctccttcttcttcgttgaagatgcggaagatcacgccggtaactccactgtggcatttcccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccg                    <  2:1284100‑M1/143‑108 (MQ=255)                                                                                                                cACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgatgatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgacttccccgcttt                <  2:236405‑M1/143‑112 (MQ=255)                                                                                                                      tCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtct          <  2:452638‑M1/143‑118 (MQ=255)                                                                                                                        tttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctgttcgtctgt        >  1:91750‑M1/1‑24 (MQ=255)                                                                                                                         tgtaagacGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcgcgccggtaactcctctgttgcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgta       <  2:449025‑M1/140‑121 (MQ=255)                                                                                                                              gacgacgacGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactggggcatgtcccatccgtggtggtcgtcggtttactctgcgtgagcgtgcactgcccgctggtcgtctgtgtcgta  >  1:515519‑M1/1‑18 (MQ=255)                                                                                                                              gacgacgacGGCCTGCACctccttcttcttctttgacgatgcggaggatcacgccggtaactccactgtggcctttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgta  <  2:367387‑M1/143‑126 (MQ=255)                                                                                                                                               |                                                                                                                              CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTA  >  NC_000913/4001502‑4001769

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