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	60 [0]	[1] 62	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)	60 (0.970)	49/284	0.2	99.2%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001650 =	1 (0.020)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGT  >  NC_000913/4001511‑4001792                                                                                                                                            |                                                                                                                                                agggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                              <  1:290235‑M2/4‑1 (MQ=255)    gattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                           <  2:869210‑M2/7‑1 (MQ=255)     attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggcagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                          >  1:471042‑M2/136‑143 (MQ=255)          tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggccacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagccgacgacggcctgcatCTCcttcttcttc                                                                                                                                     >  2:350184‑M2/131‑143 (MQ=255)           cttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                                    <  1:782085‑M2/14‑1 (MQ=255)             tgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                  <  2:635499‑M2/16‑1 (MQ=255)             ggtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                  <  1:167707‑M2/16‑1 (MQ=255)             ggtcgaaccggacgacgacgagcgactgcgcggcaaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                  <  1:564907‑M2/16‑1 (MQ=255)               tcggaccggacgacgacgagcgactgcgcgtacaatctgagctttttttgagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                                <  2:1085797‑M2/18‑1 (MQ=255)               tcgaaccggccgccggctagcgacggcgcgtacgatctgaacttggccagagcctggcatcccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                                <  2:786071‑M2/18‑1 (MQ=255)                cgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                               <  1:550571‑M2/19‑1 (MQ=255)                 gaaccggacgacgacgagcgactgcgcgtacactctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcacCACCTTCTTCTTCTTCGAAg                                                                                                                              >  2:619096‑M2/124‑143 (MQ=255)                 gaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTCCTTTGAAg                                                                                                                              >  1:1032765‑M2/124‑143 (MQ=255)                 gaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgcttctttgaagtcgacgacggcccgcatCTCCTTCTTCTTCTTTgaac                                                                                                                              >  2:732412‑M2/124‑142 (MQ=255)                           acgacgagcgactgcgcgtacaatctgaacttgcccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                                    <  1:1107564‑M2/30‑1 (MQ=255)                            cgacaaacgactgctcataaaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAt                                                                                                                   <  1:681955‑M2/31‑1 (MQ=255)                             gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                                  >  1:96126‑M2/112‑143 (MQ=255)                              acgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCa                                                                                                                 <  1:535261‑M2/33‑1 (MQ=255)                              acgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCa                                                                                                                 <  1:650433‑M2/33‑1 (MQ=255)                                          gcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcccgtctctttgaagacgacgacggcccgcatCCCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTc                                                                                                     >  1:468271‑M2/99‑143 (MQ=255)                                            gtacaggctgcactacgccagagccgggcaacccgcccggaactggaagacatcgaagcagctgcacggttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCa                                                                                                   <  2:301566‑M2/47‑1 (MQ=255)                                            atacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcatcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTCGAAGATGCGGAAGATCACGCCGGTAACTCCa                                                                                                   >  1:158355‑M2/97‑143 (MQ=255)                                             tacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgcttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATACGGAAGATCACGCCGGTAACTCCAc                                                                                                  >  2:236582‑M2/96‑143 (MQ=255)                                              acaatctgagcttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                                 <  2:8606‑M2/49‑1 (MQ=255)                                              acaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                                 <  1:822388‑M2/49‑1 (MQ=255)                                                aatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgt                                                                                               <  1:932797‑M2/51‑1 (MQ=255)                                                 atctgaacttggccatagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                              <  2:647392‑M2/52‑1 (MQ=255)                                                 atctgaacttggccagagcctggcaacccgcccggaactggaagacgtcgaagcatcggcacgtttcgttgacgacgaagacggactgcatCGACTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTCACTCCACcgtg                                                                                              >  2:16381‑M2/92‑143 (MQ=255)                                                  tctcaactgggccagagcctgtcaacgcgcccgggactggaagacatcgaagcatcggcacgtttctttggagacgacgacggcctgcatCTCCTTCTTCGTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTgg                                                                                             <  2:503324‑M2/53‑1 (MQ=255)                                                    tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCa                                                                                           <  1:857103‑M2/55‑1 (MQ=255)                                                    tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTACTTCTTCTTTGAAGCTGCGAAAGCTCACGCCGCTAACTccactgtcgaa                                                                                           >  2:140536‑M2/89‑139 (MQ=255)                                                         ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAGCACGCCGGTAACTCCACTGTGGCATTGAc                                                                                      >  2:345463‑M2/84‑143 (MQ=255)                                                           ggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCGTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCa                                                                                    <  1:211273‑M2/62‑1 (MQ=255)                                                            gccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTtaccct                                                                                   >  1:89447‑M2/81‑141 (MQ=255)                                                              cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTaccacga                                                                                 >  2:108990‑M2/79‑140 (MQ=255)                                                                  gccgggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                             <  2:1011073‑M2/69‑1 (MQ=255)                                                                     tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATgg                                                                          >  1:526905‑M2/72‑143 (MQ=255)                                                                      ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGACGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTggtggt                                                                         >  1:430063‑M2/71‑143 (MQ=255)                                                                            ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGAGATGGTCGTCTg                                                                   >  1:771678‑M2/65‑143 (MQ=255)                                                                             cgccctgaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGCGATGGTCGTCTGt                                                                  <  1:1063115‑M2/80‑1 (MQ=255)                                                                              gcccggaactggacgacatcgaagcagcggcacgtgtctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                                 <  1:654107‑M2/81‑1 (MQ=255)                                                                                  ggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAct                                                             <  2:691513‑M2/85‑1 (MQ=255)                                                                                         gaagacagcgacgcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                                      <  1:745430‑M2/92‑1 (MQ=255)                                                                                            gacatcgaggcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag‑c                                                  <  2:593012‑M2/95‑1 (MQ=255)                                                                                              catcgctgcatcgtctcgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGCTCACGCCGGTAACTCCTCTGTGGCAGTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag‑cgt                                                <  2:291035‑M2/97‑1 (MQ=255)                                                                                              catcgaagcatcggcacgtttctttgaaggcgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag‑cgt                                                <  2:969721‑M2/97‑1 (MQ=255)                                                                                              catagaagcatcggcacgtgtcgttgaagacgacgacggccggcatCGCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag‑cgt                                                <  1:336089‑M2/97‑1 (MQ=255)                                                                                                tcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag‑cgtga                                              <  1:849043‑M2/99‑1 (MQ=255)                                                                                                         cggcacgtttacttgaaggcgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTATCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGc                                     <  1:227632‑M2/108‑1 (MQ=255)                                                                                                              cgtttcttttcggtcgacgacggcctgcgtCTCCTTCTTCTTCTTTTGAGGTTCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTc                                <  2:838323‑M2/113‑1 (MQ=255)                                                                                                                tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTCGt                              <  1:785472‑M2/115‑1 (MQ=255)                                                                                                                 ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTCGTc                             >  1:919370‑M2/28‑143 (MQ=255)                                                                                                                 ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTCGTc                             <  1:60277‑M2/116‑1 (MQ=255)                                                                                                                    tttgaagacgacgacggcctgcatCTCCTTCTGCTTCTGTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAGCATCCGTGATGGGCGTCTTTTTACTCTGCG‑GAGAGGTGAACTGCACGCTTTTCGTCtgg                          >  2:1063441‑M2/25‑142 (MQ=255)                                                                                                                     ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGCACTGCCCGCTGTTCCTCTgtc                         >  2:1120695‑M2/24‑142 (MQ=255)                                                                                                                        aagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCGTTTCGTCGGTATCg                      >  2:1021762‑M2/21‑143 (MQ=255)                                                                                                                        aagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTACACTGTGGCATTTACCATCCGTAATGGTCGCCGGTTTAGTCTGCGTAAG‑CGTGAACTACGATCTTTTCGTCTGTGtcc                      >  2:831936‑M2/21‑142 (MQ=255)                                                                                                                                gacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTg              <  1:1042811‑M2/131‑1 (MQ=255)                                                                                                                                 acggccttcacCTCCTTCGTCTTCCTCGAAGATGCGCAAGATCACGCCGGTAACTCCACTGTGGCCTTTACCATCCGGGCTGGTCGTCTGTTTACTCTGCGTGAG‑CGTCAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGc             <  2:300842‑M2/132‑1 (MQ=255)                                                                                                                                      ctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTa        >  1:619416‑M2/7‑143 (MQ=255)                                                                                                                                        gcaaCTGCTTGTTCTGCGGTGAAGTTGCGGCAGCTCACGCCCGTACCTCCATTGTTGCATTTACCATCCGTGCTGGTCGTCTGTTTACTCTGCGTGTG‑CGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGc      <  2:658898/140‑1 (MQ=255)                                                                                                                                            cTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTCCCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTCCTGCGTATCGTATGCGTGCCCGTAGGCAgc  >  2:637886/1‑142 (MQ=255)                                                                                                                                            |                                                                                                                                                TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG‑CGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGT  >  NC_000913/4001511‑4001792

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