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	40 [1]	[0] 39	corA	magnesium/nickel/cobalt transporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4001643	1 (0.020)	39 (0.930)	30/284	0.4	98.7%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001649 =	0 (0.000)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001509‑4001786                                                                                                                                             |                                                                                                                                          tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatctc                                                                                                                                         <  2:46913‑M2/3‑1 (MQ=255) tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatctc                                                                                                                                         >  1:537968‑M2/141‑143 (MQ=255)      attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                    >  2:597435‑M2/136‑143 (MQ=255)       gttatcctgtcgatccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatctgcacgtttctttgaagacgacgacggcctgcatCTCcttctt                                                                                                                                   <  1:254910‑M2/9‑1 (MQ=255)             ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                             <  1:281451‑M2/15‑1 (MQ=255)             ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctcgcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                             <  1:280055‑M2/15‑1 (MQ=255)              tgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                            <  1:618503‑M2/16‑1 (MQ=255)               gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                           <  1:170553‑M2/17‑1 (MQ=255)               gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                           <  2:687021‑M2/17‑1 (MQ=255)                      cggacgacgatgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                    >  2:525956‑M2/120‑143 (MQ=255)                      cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGc                                                                                                                    >  1:91587‑M2/120‑143 (MQ=255)                      cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCCCCTTCTTCTTCTTTGAAGATGc                                                                                                                    >  1:624024‑M2/120‑143 (MQ=255)                       ggccgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccccccggaaccggaagacatcgaagcatcggcacgttcctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCg                                                                                                                   <  2:359418‑M2/25‑1 (MQ=255)                          cgacgacgggcgactgcgcgtacaatctgaagttggccggagcctggcagcccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                <  1:543339‑M2/28‑1 (MQ=255)                             cgacgctctgcttcgcgtagaatcttaacttggccagagcctggcacgcctcccggatctggaagacatctaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAt                                                                                                             <  1:504068‑M2/31‑1 (MQ=255)                              gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggcctgcatCTCCTTCTTCTTCTTTGAAGCTGCGGAAGATc                                                                                                            >  2:608514‑M2/112‑143 (MQ=255)                              gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgccggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGCTc                                                                                                            >  2:50861‑M2/112‑143 (MQ=255)                              gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                            >  1:98952‑M2/112‑143 (MQ=255)                                cgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAGCAc                                                                                                          >  2:249067‑M2/110‑143 (MQ=255)                                                caatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCAGGTAACTCCACtg                                                                                          >  2:210691‑M2/94‑143 (MQ=255)                                                  atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGCAACTCCACcgtg                                                                                        >  1:70410‑M2/92‑143 (MQ=255)                                                    ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGCAACTCCACTGTGGc                                                                                      >  1:650227‑M2/90‑143 (MQ=255)                                                      gaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCCTCTTCTTCTTTGCAGTCGCGGAAGATCACGCCGGTAACTCCACTGTGGCAt                                                                                    <  2:503222‑M2/56‑1 (MQ=255)                                                         cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                 <  2:101908‑M2/59‑1 (MQ=255)                                                                   gcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                       <  2:599273‑M2/69‑1 (MQ=255)                                                                      tggccacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgaccgcctgcatCTCCTTCTTCTTCTGTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATgg                                                                    <  2:479781‑M2/72‑1 (MQ=255)                                                                            cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTATCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                              <  2:473190‑M2/78‑1 (MQ=255)                                                                            cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                              <  2:537968‑M2/78‑1 (MQ=255)                                                                                 ccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAACATCCGTGATGGTCGTCTGTTTa                                                         >  1:71909‑M2/61‑143 (MQ=255)                                                                                         ggaggacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCGTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                 <  1:689045‑M2/91‑1 (MQ=255)                                                                                         ggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                 <  1:472971‑M2/91‑1 (MQ=255)                                                                                           aagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                               <  2:596327‑M2/93‑1 (MQ=255)                                                                                                atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtg                                          >  2:117238‑M2/46‑143 (MQ=255)                                                                                                      gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                                    >  2:226041‑M2/40‑143 (MQ=255)                                                                                                           ggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGCAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCCTGAGCGTGAGCTGCCCGCt                               >  1:510507‑M2/35‑143 (MQ=255)                                                                                                                   tcttgggagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGGTCTCGCCGGTAACACAACTGTGCCCTTTTCCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                       <  1:10034‑M2/117‑1 (MQ=255)                                                                                                                       agaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGTCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                   <  2:77251‑M2/121‑1 (MQ=255)                                                                                                                           gtcgaccacgccctgcatCTCCTTCTTCTGCTTTTGAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTa               <  1:37404‑M2/125‑1 (MQ=255)                                                                                                                                        tgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGtac  >  2:707918‑M2/6‑142 (MQ=255)                                                                                                                                             |                                                                                                                                          ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001509‑4001786

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑