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	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	= 4001643	0 (0.000)	48 (0.770)	42/282	0.5	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

ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCA  >  NC_000913/4001509‑4001789                                                                                                                                             |                                                                                                                                             tatggattgattctgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatctc                                                                                                                                            <  2:960867‑M2/3‑1 (MQ=255)     gattgatcttctcgatccggacgacgacgatcgactgcgcgtacaatctgaacttggccagcgcctggcaatccgcccgtaacttgatgacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                        <  2:111211‑M2/7‑1 (MQ=255)          atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttctt                                                                                                                                   >  1:979153‑M2/132‑143 (MQ=255)           tcttgtcgaaccggacgacgacgagcgactgcgcgcacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                  <  2:425154‑M2/13‑1 (MQ=255)             ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtatctttgaagacgacgacagcctgcatCTCcttcttcatctt                                                                                                                                >  1:738501‑M2/129‑143 (MQ=255)                      cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                       <  2:670979‑M2/24‑1 (MQ=255)                       ggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCg                                                                                                                      >  1:59087‑M2/119‑143 (MQ=255)                            acgacgagcgacggcgcgtacaatgtgaacttggacagagccgggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctggatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                                 <  2:573559‑M2/30‑1 (MQ=255)                                    cgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                         <  1:1012742‑M2/38‑1 (MQ=255)                                    cgactgcgcgtacaatctcaccttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                         <  1:468226‑M2/38‑1 (MQ=255)                                    cgactgcgcgtacaagctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                         <  1:935904‑M2/38‑1 (MQ=255)                                     gactgcgcggacaatcggaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCgg                                                                                                        <  2:887999‑M2/39‑1 (MQ=255)                                             gtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTGCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCa                                                                                                >  1:955407‑M2/97‑143 (MQ=255)                                               acaatctgaacttggctagtgtctttcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                              <  1:16318‑M2/49‑1 (MQ=255)                                                 aatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGCGGACGATCACGCCGGTCACTCCACtgt                                                                                            >  2:794789‑M2/93‑143 (MQ=255)                                                  atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                           >  1:591138‑M2/92‑143 (MQ=255)                                                    ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATGACGCGGGTAACGAGACTGTGGc                                                                                         >  1:710642‑M2/90‑143 (MQ=255)                                                          ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGGAGATCACGCCGGTCACTCCACTGTGGCATTtaa                                                                                   >  1:106110‑M2/84‑142 (MQ=255)                                                          ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCGCTGTGGcattttc                                                                                   >  1:988347‑M2/84‑141 (MQ=255)                                                          gtggccagagccgggcaacccgcccggaactggaagacatcgaagcatcggcacgtctctttgaagacgacgacggcccgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGTATTTAc                                                                                   <  1:432955‑M2/60‑1 (MQ=255)                                                               cagagccttgcaatcctctcgtatcttgaatatatcgaagcatcggcacggttcttggaagacgacgacggcctgcatCTCCTTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATcc                                                                              <  1:742498‑M2/65‑1 (MQ=255)                                                                  agcctggcggtccgcccggacctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGTGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                           <  2:874535‑M2/68‑1 (MQ=255)                                                                   gcccggcaacccgcccggaactggaagacagcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                          <  1:175253‑M2/69‑1 (MQ=255)                                                                     ctcgtggcccgccctgaactggaagacgtcgaagcatcggcacgtttctttgaagacgacgactgcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                        <  2:524171‑M2/71‑1 (MQ=255)                                                                      tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                               >  1:182172‑M2/72‑103 (MQ=255)                                                                      tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                               <  2:182172‑M2/32‑1 (MQ=255)                                                                            cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                 <  1:343205‑M2/78‑1 (MQ=255)                                                                             ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                                >  1:240803‑M2/65‑143 (MQ=255)                                                                             ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                                >  2:33511‑M2/65‑143 (MQ=255)                                                                                 ccggaactggaagacatcgacgcatcggcacgtttgtttgaagacgacgacggccggcatCGCCTTCTTCTTCTTTGAAGAGCCTGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                            <  1:20495‑M2/83‑1 (MQ=255)                                                                                     aactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTGCCATCCGCGATGGTCGCCTGTTTActct                                                        >  2:1071571‑M2/57‑143 (MQ=255)                                                                                       ctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                      <  1:751399‑M2/89‑1 (MQ=255)                                                                                        tggcagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTCACTCTgcg                                                     >  1:87847‑M2/54‑143 (MQ=255)                                                                                        tggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGGTGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGCTGGTCGTCTGTTTACTCTgcg                                                     >  2:855074‑M2/54‑143 (MQ=255)                                                                                          gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                                   <  2:979153‑M2/92‑1 (MQ=255)                                                                                            agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                                 >  2:729724‑M2/50‑143 (MQ=255)                                                                                                      gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                                       <  2:41148‑M2/104‑1 (MQ=255)                                                                                                          cggcacgtttctttgacggcgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGACCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTCCTCTGCGTGAGCGTGAACTGCCcgg                                   >  2:957202‑M2/36‑142 (MQ=255)                                                                                                          cggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCGTCCGTGCTGGTCGTCTGTTTACTCTGCGTGAGCGTGAACtgcccca                                   >  2:792552‑M2/36‑141 (MQ=255)                                                                                                            gcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGGGCGTGAACTGCCGGctg                                 >  1:55950‑M2/34‑142 (MQ=255)                                                                                                                  ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAGGATGCGGAAGGACACGCCGGTAACTCCACTGTGGCATTTACCATCCTTGTTGGTCGGCTGTTTACTCTGAGTGAGCGTGAACGGCCCGCTTTtcgcg                           >  2:176971‑M2/28‑141 (MQ=255)                                                                                                                     tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGGCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTCCATCTGt                        >  1:88608‑M2/25‑143 (MQ=255)                                                                                                                         aagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGCTGGTCGTCTGTTTACTCTCCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCg                    >  2:735850‑M2/21‑143 (MQ=255)                                                                                                                               acgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGACCGTCAAAAGTCCGCTTTTCGTCTGTATCGTATGCg              <  2:955407‑M2/129‑1 (MQ=255)                                                                                                                               acgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCGCTTTACTCTGCGTGAGCGTAACCTGCCCCCTGCTCGTCTGTCTCGTCTGCg              >  2:35236‑M2/15‑143 (MQ=255)                                                                                                                                cgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACGATCCGTGATGGGCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTGTCGTATGCGt             >  1:133907‑M2/14‑143 (MQ=255)                                                                                                                                        tgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGGCTGTTTACTCTGCGTGTGCGTGAACTGCCTGCTTTTCGACTATATCGAATGAGTGCCGGTAg     >  1:333412‑M2/6‑143 (MQ=255)                                                                                                                                           atCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTGTGCGTGCCCGTGgccg  >  1:97388‑M2/3‑142 (MQ=255)                                                                                                                                             |                                                                                                                                             ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCA  >  NC_000913/4001509‑4001789

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