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	50 [0]	[0] 51	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)	50 (0.950)	42/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

GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001508‑4001787                                                                                                                                               |                                                                                                                                          agtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatc                                                                                                                                           <  1:577561‑M2/1‑1 (MQ=255)  gtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatct                                                                                                                                          >  1:421556‑M2/142‑143 (MQ=255)    atggtttgatcttgtcgaaccggacgacggcgatcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                        <  2:439438‑M2/4‑1 (MQ=255)    atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                        <  1:362567‑M2/4‑1 (MQ=255)    atggattgagcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                        <  2:421833‑M2/4‑1 (MQ=255)         ttgatcttgtcgaaccggacgacgacgagcgactgcgcgtaccatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcgtcggcacgtctctttgacgacgacgccgacctgcaaCCccttctc                                                                                                                                   >  2:659797‑M2/135‑142 (MQ=255)         ttgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggcctgcacCTCcttctt                                                                                                                                   >  2:696700‑M2/135‑143 (MQ=255)         ctgatctggtcgaaccggaggacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttctt                                                                                                                                   <  2:866897‑M2/9‑1 (MQ=255)                  ccgacccggacgacgacgagcgagtgcgggtacgatcggaacttggccagagcctggcaacccgcccggaactggaagacctcgaagcatcggcacgtttctttgaagacgacgccggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                          <  1:283181‑M2/18‑1 (MQ=255)                   cgaaccggacgacgacgagcgacggcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                         <  1:785246‑M2/19‑1 (MQ=255)                           acgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaccccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGa                                                                                                                 <  2:608748‑M2/27‑1 (MQ=255)                           acgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGa                                                                                                                 <  2:281475‑M2/27‑1 (MQ=255)                            cgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                <  1:545737‑M2/28‑1 (MQ=255)                              gctgcgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgtagcatcggcacgtttctttgaagtcgacgacggcttgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                              <  2:508820‑M2/30‑1 (MQ=255)                              acgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                              <  1:419460‑M2/30‑1 (MQ=255)                                 tgcagcggctgcgcgtgcaagctgaagttggccagagtctggcaacccgcccggaactggaagacatcgcagcagcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCa                                                                                                           <  2:466247‑M2/33‑1 (MQ=255)                                   gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGACGATGCGGTCGCTCACg                                                                                                         >  1:43957‑M2/109‑143 (MQ=255)                                          tgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                  <  1:120544‑M2/42‑1 (MQ=255)                                                tacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCAc                                                                                            <  1:44136‑M2/48‑1 (MQ=255)                                                tacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtctctttgaagacgacgccggcgtgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCGCGTAACTCCAc                                                                                            >  2:878360‑M2/96‑143 (MQ=255)                                                   aatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGGAACTCCACtgt                                                                                         >  2:791071‑M2/93‑143 (MQ=255)                                                       tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCAGTGTGGCa                                                                                     >  1:824884‑M2/89‑143 (MQ=255)                                                       tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACAGAGGCa                                                                                     >  1:816106‑M2/89‑143 (MQ=255)                                                       tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTGCTTTGAAGATGCGGAGGATCACGCCGGTAACTCCACTGTGGCa                                                                                     >  1:592629‑M2/89‑143 (MQ=255)                                                                ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATc                                                                            >  1:668998‑M2/80‑143 (MQ=255)                                                                   gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGt                                                                         <  1:930821‑M2/67‑1 (MQ=255)                                                                     gcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                       <  1:223358‑M2/69‑1 (MQ=255)                                                                      cctggcaccccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGGGAt                                                                      >  2:388333‑M2/74‑143 (MQ=255)                                                                       ttggcaacgcgcccgtaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                     <  1:475142‑M2/71‑1 (MQ=255)                                                                               ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCCTTTACCCTCCGTGATGGTCGTCTg                                                             >  1:769541‑M2/65‑143 (MQ=255)                                                                                 gcccggaactggaagacatcgaagcatcggcgcgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                           <  2:729745‑M2/81‑1 (MQ=255)                                                                                   ccggaactggaagacatcgaagcatcggcacgaaaaacttaagtctacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                         <  1:396896‑M2/83‑1 (MQ=255)                                                                                     ggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAct                                                       <  1:918950‑M2/85‑1 (MQ=255)                                                                                       aactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActct                                                     <  2:421556‑M2/87‑1 (MQ=255)                                                                                         ctggaagacatcgaagcagcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                   <  1:866025‑M2/89‑1 (MQ=255)                                                                                            gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTGCTCTgcgtg                                                >  1:13635‑M2/52‑143 (MQ=255)                                                                                             acgacatcgaatcctcggcacgtttctttgcagacgaccacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                               <  2:473989‑M2/93‑1 (MQ=255)                                                                                               gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                             >  1:6144‑M2/49‑143 (MQ=255)                                                                                                 catcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgt                                           <  2:611346‑M2/97‑1 (MQ=255)                                                                                                  atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgcgagcgtg                                          >  2:71622‑M2/46‑143 (MQ=255)                                                                                                    cgacgcatcggcacgtttctttgaagacgacgacggcctgcatCTCCGTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                        <  2:358785‑M2/100‑1 (MQ=255)                                                                                                           tcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGGGTgagagtgaagtcccac                                 >  1:347980‑M2/37‑136 (MQ=255)                                                                                                            cggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATTTTCGGCGGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                <  1:907784‑M2/108‑1 (MQ=255)                                                                                                             ggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGCGATGGTCGTCTGTTTGCTCTGCGTGCGCGTGCACTGTCCGCt                               >  1:326563‑M2/35‑143 (MQ=255)                                                                                                                          gaagaggacgacggcctgcatCTCCTTCTTCTTCTTTTAAGATGCGGAAGATCACGCCGGTAACTCCTCTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATc                  <  2:392207‑M2/122‑1 (MQ=255)                                                                                                                           atgacgacggcggcgtgctgCTCGTTCTTCTTCTTTGAGGAGGCGGATGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAAGTGCCCGCTTTTCGTCTGTATCg                 <  2:916955‑M2/123‑1 (MQ=255)                                                                                                                             gacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATAACGCCGGTAGATCCACTGTGGCATTTACAATCAGTGATGGTCGTCAGTTTGCTCTGCGGGAGCGTGAACTGCCTGCTTTTGGTCTGTAGCGTa               >  1:193577‑M2/19‑143 (MQ=255)                                                                                                                                  cgacggcctgcatCTCCTTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATGTACCATCCGTGTTGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGt          <  1:378340‑M2/130‑1 (MQ=255)                                                                                                                                      ggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACCGTGGCATTTACCATCCGTGATGGTCGTCTGTTGCCTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTGTCGTACGCGTGgcc      >  1:63511‑M2/10‑143 (MQ=255)                                                                                                                                         cttcaccTCCCTCTTCATCTTTGGACAGGCGGAATATCACGCCTGTACCTCCACGGTGGCATTTACCATCCGTGATTGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTa   <  2:824884/137‑1 (MQ=255)                                                                                                                                          tgcatTTTCTTCTTCTTCTTTGAAGATGCGGAAGATCTCGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCTTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAg  <  2:389846‑M2/138‑1 (MQ=255)                                                                                                                                               |                                                                                                                                          GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001508‑4001787

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