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	48 [0]	[0] 48	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)	48 (0.800)	42/282	0.4	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

TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGT  >  NC_000913/4001505‑4001785                                                                                                                                            |                                                                                                                                              tATGGATTGATTCTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctc                                                                                                                                            <  2:960867‑M1/143‑4 (MQ=255)     gATTGATCTTCTCGATCCGGACGACGACGATCGACTGCGCGTACAATCTGAACTTGGCCAGCGCCTGGCAATCCGCCCGTAACTTGATGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttc                                                                                                                                        <  2:111211‑M1/143‑8 (MQ=255)          aTCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctt                                                                                                                                   >  1:979153‑M1/1‑131 (MQ=255)           tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGCACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttc                                                                                                                                  <  2:425154‑M1/143‑14 (MQ=255)             ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTATCTTTGAAGACGACGACAGCCTGCATctccttcttcatctt                                                                                                                                >  1:738501‑M1/1‑128 (MQ=255)                      cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                       <  2:670979‑M1/143‑25 (MQ=255)                       ggACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcg                                                                                                                      >  1:59087‑M1/1‑118 (MQ=255)                            acgacgaGCGACGGCGCGTACAATGTGAACTTGGACAGAGCCGGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGGATctccttcttcttctttgaagatgcggaaga                                                                                                                 <  2:573559‑M1/143‑31 (MQ=255)                                    cGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccg                                                                                                         <  1:1012742‑M1/143‑39 (MQ=255)                                    cGACTGCGCGTACAATCTCACCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccg                                                                                                         <  1:468226‑M1/143‑39 (MQ=255)                                    cGACTGCGCGTACAAGCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccg                                                                                                         <  1:935904‑M1/143‑39 (MQ=255)                                     gACTGCGCGGACAATCGGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgg                                                                                                        <  2:887999‑M1/143‑40 (MQ=255)                                             gTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctgctttgaagatgcggaagatcacgccggtaactcca                                                                                                >  1:955407‑M1/1‑96 (MQ=255)                                               aCAATCTGAACTTGGCTAGTGTCTTTCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccact                                                                                              <  1:16318‑M1/143‑50 (MQ=255)                                                 aaTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgcagatgcggacgatcacgccggtcactccactgt                                                                                            >  2:794789‑M1/1‑92 (MQ=255)                                                  aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                           >  1:591138‑M1/1‑91 (MQ=255)                                                    cTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatgacgcgggtaacgagactgtggc                                                                                         >  1:710642‑M1/1‑89 (MQ=255)                                                          ttGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgggagatcacgccggtcactccactgtggcatttaa                                                                                   >  1:106110‑M1/1‑83 (MQ=255)                                                          ttGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccgctgtggcattttc                                                                                   >  1:988347‑M1/1‑83 (MQ=255)                                                          gtgGCCAGAGCCGGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTCTCTTTGAAGACGACGACGGCCCGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggtatttac                                                                                   <  1:432955‑M1/142‑61 (MQ=255)                                                               cAGAGCCTTGCAATCCTCTCGTATCTTGAATATATCGAAGCATCGGCACGGTTCTTGGAAGACGACGACGGCCTGCATctccttctgcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcc                                                                              <  1:742498‑M1/143‑66 (MQ=255)                                                                  agCCTGGCGGTCCGCCCGGACCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgtggaagatcacgccggtaactccactgtggcatttaccatccgtg                                                                           <  2:874535‑M1/143‑69 (MQ=255)                                                                   gCCCGGCAACCCGCCCGGAACTGGAAGACAGCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtga                                                                          <  1:175253‑M1/143‑70 (MQ=255)                                                                     ctcgtggcccgccctgaactggAAGACGTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACTGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatg                                                                        <  2:524171‑M1/136‑72 (MQ=255)                                                                      tGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatc                                                                                                               >  1:182172‑M1/1‑71 (MQ=255)                                                                      tGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatc                                                                                                               <  2:182172‑M1/103‑33 (MQ=255)                                                                            cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                                 <  1:343205‑M1/143‑79 (MQ=255)                                                                             ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctg                                                                >  1:240803‑M1/1‑64 (MQ=255)                                                                             ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctg                                                                >  2:33511‑M1/1‑64 (MQ=255)                                                                                 ccgGAACTGGAAGACATCGACGCATCGGCACGTTTGTTTGAAGACGACGACGGCCGGCATcgccttcttcttctttgaagagcctgaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttta                                                            <  1:20495‑M1/143‑84 (MQ=255)                                                                                     aaCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttgccatccgcgatggtcgcctgtttactct                                                        >  2:1071571‑M1/1‑56 (MQ=255)                                                                                       cTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgc                                                      <  1:751399‑M1/143‑90 (MQ=255)                                                                                        tGGCAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgcagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttcactctgcg                                                     >  1:87847‑M1/1‑53 (MQ=255)                                                                                        tGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaaggtgcggaagatcacgccggtaactccactgtggcatttaccatccgtgctggtcgtctgtttactctgcg                                                     >  2:855074‑M1/1‑53 (MQ=255)                                                                                          gAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtg                                                   <  2:979153‑M1/143‑93 (MQ=255)                                                                                            aGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgag                                                 >  2:729724‑M1/1‑49 (MQ=255)                                                                                                      gCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgc                                       <  2:41148‑M1/143‑105 (MQ=255)                                                                                                          cGGCACGTTTCTTTGACGGCGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagaccacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttcctctgcgtgagcgtgaactgcccgg                                   >  2:957202‑M1/1‑35 (MQ=255)                                                                                                          cGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgcagatgcggaagatcacgccggtaactccactgtggcatttaccgtccgtgctggtcgtctgtttactctgcgtgagcgtgaactgcccca                                   >  2:792552‑M1/1‑35 (MQ=255)                                                                                                            gCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgggcgtgaactgccggctg                                 >  1:55950‑M1/1‑33 (MQ=255)                                                                                                                  ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaggatgcggaaggacacgccggtaactccactgtggcatttaccatccttgttggtcggctgtttactctgagtgagcgtgaacggcccgcttttcgcg                           >  2:176971‑M1/1‑27 (MQ=255)                                                                                                                     tttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcggctgtttactctgcgtgagcgtgaactgcccgctttccatctgt                        >  1:88608‑M1/1‑24 (MQ=255)                                                                                                                         aaGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgctggtcgtctgtttactctccgtgagcgtgaactgcccgcttttcgtctgtatcg                    >  2:735850‑M1/1‑20 (MQ=255)                                                                                                                               acgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgaccgtcaaaagtccgcttttcgtctgtatcgtatgcg              <  2:955407‑M1/143‑130 (MQ=255)                                                                                                                               acgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtcgctttactctgcgtgagcgtaacctgccccctgctcgtctgtctcgtctgcg              >  2:35236‑M1/1‑14 (MQ=255)                                                                                                                                cgacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttacgatccgtgatgggcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtgtcgtatgcgt             >  1:133907‑M1/1‑13 (MQ=255)                                                                                                                                        tGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcggctgtttactctgcgtgtgcgtgaactgcctgcttttcgactatatcgaatgagtgccggtag     >  1:333412‑M1/1‑5 (MQ=255)                                                                                                                                           aTctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtgtgcgtgcccgtggccg  >  1:97388‑M1/1‑2 (MQ=255)                                                                                                                                            |                                                                                                                                              TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGT  >  NC_000913/4001505‑4001785

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