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	41 [0]	[0] 41	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)	41 (0.630)	33/284	1.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

TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001511‑4001787                                                                                                                                            |                                                                                                                                          atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                        <  1:735324‑M2/4‑1 (MQ=255) atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                        <  1:651676‑M2/4‑1 (MQ=255)   ggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCctt                                                                                                                                      <  1:749589‑M2/6‑1 (MQ=255)     attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                    >  1:949571‑M2/136‑143 (MQ=255)           ctcgtcgagccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                              <  2:916818‑M2/14‑1 (MQ=255)            ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggcactggcagacatcgaagcatcggcccgtttctttgaagacgccgcaggcctgcatCTCcttcttcttctt                                                                                                                             >  1:427089‑M2/129‑143 (MQ=255)              gctgaacctgacgaggaggggcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                           <  2:379541‑M2/17‑1 (MQ=255)                 gaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcccgtttctttgacgacgacgccggcctgcatCCTCTTCTTCTCCATTCAAg                                                                                                                        >  1:135916‑M2/124‑143 (MQ=255)                         cgaggacgcgcgactgcgcgtacaatctgaacttggccagagcctggcaacctgcctgcaacttgaagacatcgaagcatcggcacgtttctgtgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                <  2:1018795‑M2/28‑1 (MQ=255)                           atcacgggcgactgcgtgtacaatctgaacttggccagagcctggcacccccccctgaactggaagacatcgaagcatcgccacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                              <  2:191662‑M2/30‑1 (MQ=255)                                gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                         >  2:937915‑M2/109‑143 (MQ=255)                                  gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGcc                                                                                                       >  2:913408‑M2/107‑143 (MQ=255)                                                 atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                        <  1:584878‑M2/52‑1 (MQ=255)                                                  tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTgg                                                                                       >  1:248416‑M2/91‑143 (MQ=255)                                                  tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggaacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTACACTGTgg                                                                                       >  1:10‑M2/91‑143 (MQ=255)                                                        cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTGCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTACCTCCACTGTGGCATTTa                                                                                 <  1:941390‑M2/59‑1 (MQ=255)                                                        ctgggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                 <  2:221257‑M2/59‑1 (MQ=255)                                                        cgtggccagagcctggcaacccgccgggaacgggaagacatcgaagcatcggcacgttgctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                 <  2:982218‑M2/59‑1 (MQ=255)                                                         ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGACGCTGCGGAAGCTCACCCCGGTCACTCCACTGCGGCATTTAc                                                                                >  2:268770‑M2/84‑143 (MQ=255)                                                           ggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTtacaa                                                                              >  1:799335‑M2/82‑141 (MQ=255)                                                                 agcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAGGATCACGCCGGTCACTCCCCTGCGGCACTTCCCATCCGTg                                                                        >  2:740449‑M2/76‑143 (MQ=255)                                                                 agcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttatttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                        >  1:602619‑M2/76‑143 (MQ=255)                                                                  gcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                       <  1:141293‑M2/69‑1 (MQ=255)                                                                  gcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                       <  1:1133527‑M2/69‑1 (MQ=255)                                                                     tggccacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcagCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATgg                                                                    <  2:1002048‑M2/72‑1 (MQ=255)                                                                       gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcatgcatCTACTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCAGTGatgggc                                                                  >  1:174950‑M2/70‑141 (MQ=255)                                                                         aacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcgt                                                                <  1:671817‑M2/76‑1 (MQ=255)                                                                            ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGACGATGCGGAAGATCGCGCCGGTAACTCCACACTGTCTGTTTTTATCCGTGCTGGTCGTCTg                                                             >  1:733149‑M2/65‑143 (MQ=255)                                                                              gcccggaactggaagacatcgcagcttcggcacgtttctgtgaagaccacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                           <  2:1082697‑M2/81‑1 (MQ=255)                                                                              gcccggaactggaagacatcgaagcatcggcgcgtgtctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                           <  1:911768‑M2/81‑1 (MQ=255)                                                                                          aggacatggaagcatcggtacttttcttttaagactacgacggcctgcatCTCTTTCTTCTTCTTTGAAGATGCGGAAGATCACTCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                               <  2:976473‑M2/93‑1 (MQ=255)                                                                                                 cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAATTCACTGTGGCATTTCCCATCCGTGATCGTCGTGTGTTTACTcggcgtgagcgtgcg                                        >  2:234695‑M2/44‑141 (MQ=255)                                                                                                      catcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                   >  1:113412‑M2/39‑143 (MQ=255)                                                                                                            cccgtttctgtgaatacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCttt                             <  1:1040392‑M2/111‑1 (MQ=255)                                                                                                               gtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCg                          <  1:148832‑M2/114‑1 (MQ=255)                                                                                                                tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTcgg                         >  1:1158215‑M2/29‑142 (MQ=255)                                                                                                                 ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATTCGTGATGGTCGTCTGATTACTCTGCGTGAGCGTGAACTCCCCGCCTTTCGTc                        >  2:316654‑M2/28‑143 (MQ=255)                                                                                                                  tctttgaagactacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                       <  2:602619‑M2/117‑1 (MQ=255)                                                                                                                      tgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                   >  1:566362‑M2/23‑143 (MQ=255)                                                                                                                             gccgacggcctgcatCGCCTTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGACCTGCCCGCTTTTCGTCTGTATCGTATGc            <  1:1153284‑M2/128‑1 (MQ=255)                                                                                                                                       tgcttCTCCTTCTTCTGCTGTGAAGATGCGGAAGTTCACGCCGGTTTTTCCACTGTGGCATTTACCATCCGGGATGGTCGTCTGTTTACTCTGCGTGCGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAg  <  2:758014‑M2/138‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001511‑4001787

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