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

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

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

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