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	52 [0]	[0] 53	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)	52 (0.950)	43/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

AGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCC  >  NC_000913/4001502‑4001781                                                                                                                                              |                                                                                                                                           aGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATc                                                                                                                                           <  2:577561‑M1/143‑2 (MQ=255)  gTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATct                                                                                                                                          >  2:421556‑M1/1‑141 (MQ=255)    aTGGTTTGATCTTGTCGAACCGGACGACGGCGATCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                        <  1:439438‑M1/143‑5 (MQ=255)    aTGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                        <  2:362567‑M1/143‑5 (MQ=255)    aTGGATTGAGCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                        <  1:421833‑M1/143‑5 (MQ=255)         ttGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACCATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCGTCGGCACGTCTCTTTGACGACGACGCCGACCTGCAAccccttctc                                                                                                                                   >  1:659797‑M1/1‑134 (MQ=255)         ttGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGCCGACGGCCTGCACctccttctt                                                                                                                                   >  1:696700‑M1/1‑134 (MQ=255)         ctgATCTGGTCGAACCGGAGGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctt                                                                                                                                   <  1:866897‑M1/142‑10 (MQ=255)             tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCCTccccctcttcctc                                                                                                                               >  1:847666‑M1/1‑130 (MQ=255)                  ccgACCCGGACGACGACGAGCGAGTGCGGGTACGATCGGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACCTCGAAGCATCGGCACGTTTCTTTGAAGACGACGCCGGCCTGCATctccttcttcttctttga                                                                                                                          <  2:283181‑M1/142‑19 (MQ=255)                   cGAACCGGACGACGACGAGCGACGGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                         <  2:785246‑M1/143‑20 (MQ=255)                           acgacgacgAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCACCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgga                                                                                                                 <  1:608748‑M1/143‑28 (MQ=255)                           acgacgacgAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgga                                                                                                                 <  1:281475‑M1/143‑28 (MQ=255)                            cgacgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                <  2:545737‑M1/143‑29 (MQ=255)                              gctgcgagcgaCTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGTAGCATCGGCACGTTTCTTTGAAGTCGACGACGGCTTGCATctccttcttcttctttgaagatgcggaaga                                                                                                              <  1:508820‑M1/139‑31 (MQ=255)                              acgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaga                                                                                                              <  2:419460‑M1/143‑31 (MQ=255)                                 tgcagcgGCTGCGCGTGCAAGCTGAAGTTGGCCAGAGTCTGGCAACCCGCCCGGAACTGGAAGACATCGCAGCAGCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatca                                                                                                           <  1:466247‑M1/140‑34 (MQ=255)                                   gaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgacgatgcggtcgctcacg                                                                                                         >  2:43957‑M1/1‑108 (MQ=255)                                          tGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaa                                                                                                  <  2:120544‑M1/143‑43 (MQ=255)                                                tACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccac                                                                                            <  2:44136‑M1/143‑49 (MQ=255)                                                tACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTCTCTTTGAAGACGACGCCGGCGTGCATctccttcttcttctttgaagatgcggaagatcacgcgcgtaactccac                                                                                            >  1:878360‑M1/1‑95 (MQ=255)                                                   aaTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgggaactccactgt                                                                                         >  1:791071‑M1/1‑92 (MQ=255)                                                       tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccagtgtggca                                                                                     >  2:824884‑M1/1‑88 (MQ=255)                                                       tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccacagaggca                                                                                     >  2:816106‑M1/1‑88 (MQ=255)                                                       tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctgctttgaagatgcggaggatcacgccggtaactccactgtggca                                                                                     >  2:592629‑M1/1‑88 (MQ=255)                                                                ccAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatc                                                                            >  2:668998‑M1/1‑79 (MQ=255)                                                                   gagCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgt                                                                         <  2:930821‑M1/143‑68 (MQ=255)                                                                     gCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtga                                                                       <  2:223358‑M1/143‑70 (MQ=255)                                                                      ccTGGCACCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgggat                                                                      >  1:388333‑M1/1‑73 (MQ=255)                                                                       ttgGCAACGCGCCCGTAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatg                                                                     <  2:475142‑M1/142‑72 (MQ=255)                                                                        tGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctgcttgttctttgaagctgcggaagatcacgccggtaactccactacgtcatttaccgcgcgagagag                                                                    >  1:60179‑M1/1‑71 (MQ=255)                                                                               ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcctttaccctccgtgatggtcgtctg                                                             >  2:769541‑M1/1‑64 (MQ=255)                                                                                 gcccgGAACTGGAAGACATCGAAGCATCGGCGCGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                           <  1:729745‑M1/143‑82 (MQ=255)                                                                                   ccgGAACTGGAAGACATCGAAGCATCGGCACGAAAAACTTAAGTCTACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttta                                                         <  2:396896‑M1/143‑84 (MQ=255)                                                                                     ggAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttact                                                       <  2:918950‑M1/143‑86 (MQ=255)                                                                                       aaCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactct                                                     <  1:421556‑M1/143‑88 (MQ=255)                                                                                         cTGGAAGACATCGAAGCAGCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgc                                                   <  2:866025‑M1/143‑90 (MQ=255)                                                                                            gAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttgctctgcgtg                                                >  2:13635‑M1/1‑51 (MQ=255)                                                                                             acgacATCGAATCCTCGGCACGTTTCTTTGCAGACGACCACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                               <  1:473989‑M1/141‑94 (MQ=255)                                                                                               gACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagc                                             >  2:6144‑M1/1‑48 (MQ=255)                                                                                                 cATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgt                                           <  1:611346‑M1/143‑98 (MQ=255)                                                                                                  aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgcgagcgtg                                          >  1:71622‑M1/1‑45 (MQ=255)                                                                                                    cGACGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccgtcttcttcttggaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaa                                        <  1:358785‑M1/143‑101 (MQ=255)                                                                                                           tCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgggtgagagtgaagtcccac                                 >  2:347980‑M1/1‑36 (MQ=255)                                                                                                            cGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgattttcggcggtttactctgcgtgagcgtgaactgcccgc                                <  2:907784‑M1/143‑109 (MQ=255)                                                                                                             ggCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgcgatggtcgtctgtttgctctgcgtgcgcgtgcactgtccgct                               >  2:326563‑M1/1‑34 (MQ=255)                                                                                                                          gAAGAGGACGACGGCCTGCATctccttcttcttcttttaagatgcggaagatcacgccggtaactcctctgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatc                  <  1:392207‑M1/143‑123 (MQ=255)                                                                                                                           atgacgacggcggCGTGCTGctcgttcttcttctttgaggaggcggatgatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaagtgcccgcttttcgtctgtatcg                 <  1:916955‑M1/141‑124 (MQ=255)                                                                                                                             gacgacgacGGCCTGCATctccttcttcttctttgaagatgcggaagataacgccggtagatccactgtggcatttacaatcagtgatggtcgtcagtttgctctgcgggagcgtgaactgcctgcttttggtctgtagcgta               >  2:193577‑M1/1‑18 (MQ=255)                                                                                                                                  cgacgGCCTGCATctccttcttcttcttggaagatgcggaagatcacgccggtaactccactgtggcatgtaccatccgtgttggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgt          <  2:378340‑M1/143‑131 (MQ=255)                                                                                                                                      ggCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccaccgtggcatttaccatccgtgatggtcgtctgttgcctctgcgtgagcgtgaactgcccgcttttcgtctgtgtcgtacgcgtggcc      >  2:63511‑M1/1‑9 (MQ=255)                                                                                                                                          tGCATtttcttcttcttctttgaagatgcggaagatctcgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcttgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtag  <  1:389846‑M1/143‑139 (MQ=255)                                                                                                                                              |                                                                                                                                           AGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCC  >  NC_000913/4001502‑4001781

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