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	59 [0]	[0] 59	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)	59 (0.850)	53/284	0.2	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

TGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCA  >  NC_000913/4001512‑4001790                                                                                                                                           |                                                                                                                                             tggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacactctgaacttggccagcgcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgccggcccgcaaCTCct                                                                                                                                          >  1:683111‑M2/139‑143 (MQ=255)  ggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCctt                                                                                                                                         >  1:1212291‑M2/138‑143 (MQ=255)    attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                       >  1:715673‑M2/136‑143 (MQ=255)      tgatcgtgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgacgacgacgacggcctgcatCTCcttcttc                                                                                                                                     <  1:386231‑M2/10‑1 (MQ=255)         tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatcggaacttggccagagccgggcaacccgcccggaactggaacacatcgaagcatcggcgcgtttctttggagacgacgacggcctgcatcgccttctgcctc                                                                                                                                  >  1:409272‑M2/131‑138 (MQ=255)          cttctcgaacccgacgacgccgagcgactgcccgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgacgcatcggcacgttgctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                                 <  2:215262‑M2/14‑1 (MQ=255)           ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgcagcagcggcacggttggttgaagacggcgacggcctgcagctccctctccggcct                                                                                                                                >  2:1018528‑M2/129‑136 (MQ=255)               cgaaacggccgaggacgagcgactgcgcgtacaatctgaactgggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                            <  2:661605‑M2/19‑1 (MQ=255)                 aaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGa                                                                                                                          <  1:1052535‑M2/21‑1 (MQ=255)                   ccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcctggaactggaagacatcgaagcatcggcacgtctctttgaagacgacgacggcctgcatCTGCGTGTTCTTCTTTGACGATg                                                                                                                        >  2:1133154‑M2/121‑143 (MQ=255)                   ccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATg                                                                                                                        <  1:1186589‑M2/23‑1 (MQ=255)                   ccgcaccacgacgaccgaccgcgcgcacaatctcaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctccatCTCCTTCTTCTTCTTTGAAGATg                                                                                                                        <  1:1127176‑M2/23‑1 (MQ=255)                    cggacgacggcgagcgacagcgcgtacaatctgaacttggccagagcccggcaacccgcccggaactggaagacagcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                       <  1:495776‑M2/24‑1 (MQ=255)                    cggacgacgacgagcgactgcgcggacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                       <  1:666242‑M2/24‑1 (MQ=255)                    cgcacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                       <  2:773212‑M2/24‑1 (MQ=255)                     ggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTCCTTCTTCTTTGGCGATGCg                                                                                                                      >  1:202617‑M2/119‑143 (MQ=255)                      gacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATgggg                                                                                                                     >  2:964486‑M2/118‑143 (MQ=255)                      gacgacgacgagcgactgcgcgtacaatctgaacttggccagagccgggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCgg                                                                                                                     <  2:877516‑M2/26‑1 (MQ=255)                                 gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGcc                                                                                                          >  1:29037‑M2/107‑143 (MQ=255)                                    cctgcgcgtacaatctgaacttggccatagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGt                                                                                                       <  1:731358‑M2/40‑1 (MQ=255)                                      tgcgcgggcaatctgaactgggccagaggctggcgacccgcgcggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTGTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                     <  2:68455‑M2/42‑1 (MQ=255)                                       gcgcctacaatcggaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcagcggcacgtttcttggaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAc                                                                                                    <  2:214243‑M2/43‑1 (MQ=255)                                          cgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTcc                                                                                                 <  1:969513‑M2/46‑1 (MQ=255)                                            tacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAGCACGCCGGTAACTGCAc                                                                                               >  1:947607‑M2/96‑143 (MQ=255)                                                cgcggaccttggccagagcctggcaacgcgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                           <  2:423544‑M2/52‑1 (MQ=255)                                                  ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGc                                                                                         <  2:943743‑M2/54‑1 (MQ=255)                                                   tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGCGGAAGATAACGCCGGTAACTCCACTGTGGCa                                                                                        >  2:993528‑M2/89‑143 (MQ=255)                                                      acttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCCGAAGCACACGGCCGTGCTTCCACCGTGGCAttt                                                                                     >  2:1143805‑M2/86‑143 (MQ=255)                                                       cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCatttc                                                                                    >  1:471647‑M2/85‑142 (MQ=255)                                                       cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacaacaacggactgcatCTCCTTCTTCTTATTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                    >  1:1100110‑M2/85‑143 (MQ=255)                                                            ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCTTTTACCATc                                                                               >  1:613739‑M2/80‑143 (MQ=255)                                                             cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACCCCACTGTGGCATTTACCATcc                                                                              >  2:1145525‑M2/79‑143 (MQ=255)                                                               gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGt                                                                            <  2:947938‑M2/67‑1 (MQ=255)                                                                 gcctggcaacccgcccggaacgggaagacatcgaagcatcggcacgtttctttgaagacgacggcggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                          <  2:29250‑M2/69‑1 (MQ=255)                                                                          cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                 <  1:805134‑M2/78‑1 (MQ=255)                                                                          cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                 >  1:500053‑M2/66‑143 (MQ=255)                                                                            cgcccggaactggaagacatcgaagcaccggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                               <  1:784713‑M2/80‑1 (MQ=255)                                                                             gcccggaactgcaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCGTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                              <  1:358622‑M2/81‑1 (MQ=255)                                                                              cccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATTCGTGATGGTCGTCTGttt                                                             >  1:272605‑M2/62‑143 (MQ=255)                                                                               ctagaactggaagacatctaagcatccgcacgtttctgtgaagacgacgacggcctgcatCTCCTTCTACTTCTTTGAAGATGCTGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                            <  1:1004528‑M2/83‑1 (MQ=255)                                                                               ccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                            >  2:787202‑M2/61‑143 (MQ=255)                                                                                  gaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActc                                                         <  1:792873‑M2/86‑1 (MQ=255)                                                                                      tggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTGTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGGGATGGTCGTCTGTTTACTCTgcg                                                     >  2:607598‑M2/54‑143 (MQ=255)                                                                                       ggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                    >  1:554249‑M2/53‑143 (MQ=255)                                                                                       ggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                    <  2:715673‑M2/91‑1 (MQ=255)                                                                                                ggaggcatcgtcgcgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                           <  1:901367‑M2/100‑1 (MQ=255)                                                                                                cgaagcatcggcacgtttctttggagacggggaccgcctgcatCCGCGGCTTCGGCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCTTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                           <  1:170599‑M2/100‑1 (MQ=255)                                                                                                           cccgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCttt                                <  1:1091216‑M2/111‑1 (MQ=255)                                                                                                              gtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCg                             >  1:344019‑M2/30‑143 (MQ=255)                                                                                                                ttctttgaagacgaccacgtcttgcatCTCCTTCTTCTTCTTTGAAGATGCGGCAGATCACGCCGGTAACTCCACTGCCGCATTTACCATCCGTGATGTTCGTCTCTTTACTCTGCGTGAGCGTGAACTGCCCGCTCTTCGTc                           <  1:189481‑M2/116‑1 (MQ=255)                                                                                                                    ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTGTTCGTGTGTa                       >  1:500580‑M2/24‑143 (MQ=255)                                                                                                                          acgacgccggcctgcatCTCCTTCTTCTTCGTTGCAGATGCGGAAGATTACGCCGGTCACTCCACTGTGGCATTTACCATCCGTGATGGTCTTCTGTTTTCTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAt                 <  2:1240535‑M2/126‑1 (MQ=255)                                                                                                                          acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGCACTGCCCGCTTTTCGTCTGTATCGTAt                 >  2:805428‑M2/18‑143 (MQ=255)                                                                                                                          acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAt                 >  1:194045‑M2/18‑143 (MQ=255)                                                                                                                            gactacggtctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGc               <  1:1156118‑M2/128‑1 (MQ=255)                                                                                                                            gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGc               >  2:644507‑M2/16‑143 (MQ=255)                                                                                                                                 cggcctgcatCTCCTGCTTCTGCTTTGAAGATGCGGAAGATCCCGAGGGTAACTCCACTGCGGCAGTTACCTTCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGcc          <  1:341200‑M2/133‑1 (MQ=255)                                                                                                                                  ggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGccc         >  1:594386‑M2/10‑143 (MQ=255)                                                                                                                                         atCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTCTCGTCTGTATCGTATGCGTGCCCGTAGCCa  >  2:516361‑M2/3‑143 (MQ=255)                                                                                                                                           |                                                                                                                                             TGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCA  >  NC_000913/4001512‑4001790

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