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

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4001643	0 (0.000)	62 (0.860)	53/284	0.2	100%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001649 =	0 (0.000)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

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