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	36 [0]	[0] 35	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)	35 (0.870)	32/284	0.3	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

TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCC  >  NC_000913/4001511‑4001789                                                                                                                                            |                                                                                                                                            gcagtttgctctcgtcgaagctgacaacgccgagcgactgctcgtacaatctgaacttggccagagcctggcaccccgcccggaactggaatacatcgaagcatcggtacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                          <  2:212600‑M2/4‑1 (MQ=255)     attactcttgtctacccggacgacgacgagcgactgcgcgtccaacccgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                      <  1:345550‑M2/8‑1 (MQ=255)          tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                 >  1:424931‑M2/131‑143 (MQ=255)               ccgaaccggacgacggcgagcgactgcgcgtacaatctgaacctggccagagcctggcagcccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                            <  1:464990‑M2/18‑1 (MQ=255)                cgacccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                           <  1:81461‑M2/19‑1 (MQ=255)                cgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                           <  1:581910‑M2/19‑1 (MQ=255)                  aaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGa                                                                                                                         <  1:238631‑M2/21‑1 (MQ=255)                         cgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                  <  1:352261‑M2/28‑1 (MQ=255)                                gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTCCTTTGAAGATGCGGAAGATCACg                                                                                                           >  2:466608‑M2/109‑143 (MQ=255)                                gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTCCTTCTCTGAACATGCGGAACGTCACg                                                                                                           >  2:584716‑M2/109‑143 (MQ=255)                                    gactgcgcggacaagctgaagttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCgg                                                                                                       <  1:316637‑M2/39‑1 (MQ=255)                                          gcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGACGCTGCGGACGGTCACGCCGGTAGctg                                                                                                 >  1:296069‑M2/99‑142 (MQ=255)                                                 atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                          <  1:237793‑M2/52‑1 (MQ=255)                                                   ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGTGGAAGATCACGCCGGTAACTCCACTGTGGc                                                                                        >  1:395295‑M2/90‑143 (MQ=255)                                                    tgaacttggccagagcctggcaacccgcccggaactggaagacatcgcagcatcggcacgtttcttggaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGAGGCGGAAGATCACGCCGGTAAAgccactgagcgc                                                                                       >  1:670761‑M2/89‑138 (MQ=255)                                                             cccgagcctggcaacccgcccggaactggacgccatcgacgcatcggcacgtttctttgaagacgacgaccgcctgcatCTCCTTCTTCTTCTTTGCAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATc                                                                              <  1:176229‑M2/64‑1 (MQ=255)                                                              cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATcc                                                                             >  2:517973‑M2/79‑143 (MQ=255)                                                                gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGt                                                                           <  1:176889‑M2/67‑1 (MQ=255)                                                                         aacccgcccggaactggaagacatcgaagcatcggcacggttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTCTGGCATTTACCATCCGTGATGgtcgt                                                                  <  2:694786‑M2/76‑1 (MQ=255)                                                                           cccgcccggagcgggaagacatcgaagcatcggcacgtttcttttaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                <  2:402298‑M2/78‑1 (MQ=255)                                                                                   gaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGCGATGGTCGTCTGTTTActc                                                        >  2:531430‑M2/58‑143 (MQ=255)                                                                                    aactggaagacagcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActct                                                       <  1:111158‑M2/87‑1 (MQ=255)                                                                                        ggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                   <  1:13540‑M2/91‑1 (MQ=255)                                                                                        ggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCGTCTTTGAAGATGCGGAAGATCACGCCGGTACCTCCGCTGTGGCATTTGCCATCCGTGATGGTCGGCTGTTTACTCTgcgt                                                   >  1:549282‑M2/53‑143 (MQ=255)                                                                                            gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                               >  1:252733‑M2/49‑143 (MQ=255)                                                                                               atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTC‑‑‑CTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATACGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAAc                                         >  1:177761‑M2/46‑143 (MQ=255)                                                                                                           gcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGTGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtt                                <  2:395295‑M2/110‑1 (MQ=255)                                                                                                           gcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtt                                <  2:675406‑M2/110‑1 (MQ=255)                                                                                                                   ctttgaagacgacgccggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTACCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTg                        <  1:315235‑M2/118‑1 (MQ=255)                                                                                                                    tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTCTCGGCTGt                       >  1:177128‑M2/25‑143 (MQ=255)                                                                                                                           acgacgtcggcctgcatCTCCTTTTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTATTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAt                <  1:696586‑M2/126‑1 (MQ=255)                                                                                                                             tgcgtcggcctgcatCTCCTTCTTCTGCCTTGACGTTGCGGCAGCTCTCGCCGGTACCTCCTCTGTGGCACTTACCCTCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGc              <  2:722643‑M2/128‑1 (MQ=255)                                                                                                                                 gcggcatgcatCTCCTTCTTCTTCTCTGAAGATGCGGAAGATCACGCCGGTAACTCCCCTGTGGCATTTACCATCCGTGATGGTCGTTTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGc          <  1:145406‑M2/132‑1 (MQ=255)                                                                                                                                        gcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGCGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGCGCCCGTAGc   >  2:615002‑M2/5‑143 (MQ=255)                                                                                                                                         catCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAAATATTCGCTGCCCGTCTGTATCGTATGCGTGCCCGTTGcc  >  1:302224‑M2/4‑143 (MQ=255)                                                                                                                                            |                                                                                                                                            TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCC  >  NC_000913/4001511‑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: ‑