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	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	= 4001643	0 (0.000)	41 (0.610)	34/284	1.3	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

ATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCC  >  NC_000913/4001505‑4001781                                                                                                                                           |                                                                                                                                           aTGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                        <  1:735324‑M1/143‑5 (MQ=255) aTGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                        <  1:651676‑M1/143‑5 (MQ=255)   ggATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctt                                                                                                                                      <  1:749589‑M1/143‑7 (MQ=255)     aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                                    >  1:949571‑M1/1‑135 (MQ=255)           ctcgtcgAGCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttct                                                                                                                              <  2:916818‑M1/140‑15 (MQ=255)            ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGCACTGGCAGACATCGAAGCATCGGCCCGTTTCTTTGAAGACGCCGCAGGCCTGCATctccttcttcttctt                                                                                                                             >  1:427089‑M1/1‑128 (MQ=255)              gctgaacCTGACGAGGAGGGGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                           <  2:379541‑M1/140‑18 (MQ=255)                 gAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCCCGTTTCTTTGACGACGACGCCGGCCTGCATcctcttcttctccattcaag                                                                                                                        >  1:135916‑M1/1‑123 (MQ=255)                         cgaggacgcGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCTGCCTGCAACTTGAAGACATCGAAGCATCGGCACGTTTCTGTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                <  2:1018795‑M1/143‑29 (MQ=255)                           atcacgggcgACTGCGTGTACAATCTGAACTTGGCCAGAGCCTGGCACCCCCCCCTGAACTGGAAGACATCGAAGCATCGCCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaga                                                                                                              <  2:191662‑M1/140‑31 (MQ=255)                                gaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacg                                                                                                         >  2:937915‑M1/1‑108 (MQ=255)                                  gCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcc                                                                                                       >  2:913408‑M1/1‑106 (MQ=255)                                                 aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                        <  1:584878‑M1/143‑53 (MQ=255)                                                  tCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtgg                                                                                       >  1:248416‑M1/1‑90 (MQ=255)                                                  tCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGAACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactacactgtgg                                                                                       >  1:10‑M1/1‑90 (MQ=255)                                                        ctgggccaGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattta                                                                                 <  2:221257‑M1/140‑60 (MQ=255)                                                        cgtggCCAGAGCCTGGCAACCCGCCGGGAACGGGAAGACATCGAAGCATCGGCACGTTGCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattta                                                                                 <  2:982218‑M1/141‑60 (MQ=255)                                                        cTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctgcttcttctttgaagatgcggaagatcacgccggtacctccactgtggcattta                                                                                 <  1:941390‑M1/143‑60 (MQ=255)                                                         ttGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgacgctgcggaagctcaccccggtcactccactgcggcatttac                                                                                >  2:268770‑M1/1‑83 (MQ=255)                                                           ggCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttacaa                                                                              >  1:799335‑M1/1‑81 (MQ=255)                                                                 agCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaggatcacgccggtcactcccctgcggcacttcccatccgtg                                                                        >  2:740449‑M1/1‑75 (MQ=255)                                                                 agCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTATTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtg                                                                        >  1:602619‑M1/1‑75 (MQ=255)                                                                  gCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtga                                                                       <  1:141293‑M1/143‑70 (MQ=255)                                                                  gCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtga                                                                       <  1:1133527‑M1/143‑70 (MQ=255)                                                                     tGGCCACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCAGctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgg                                                                    <  2:1002048‑M1/143‑73 (MQ=255)                                                                       gCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCATGCATctacttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcagtgatgggc                                                                  >  1:174950‑M1/1‑69 (MQ=255)                                                                         aaCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgt                                                                <  1:671817‑M1/143‑77 (MQ=255)                                                                            ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgacgatgcggaagatcgcgccggtaactccacactgtctgtttttatccgtgctggtcgtctg                                                             >  1:733149‑M1/1‑64 (MQ=255)                                                                              gcccgGAACTGGAAGACATCGCAGCTTCGGCACGTTTCTGTGAAGACCACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                           <  2:1082697‑M1/143‑82 (MQ=255)                                                                              gcccgGAACTGGAAGACATCGAAGCATCGGCGCGTGTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtt                                                           <  1:911768‑M1/143‑82 (MQ=255)                                                                                          aggacATGGAAGCATCGGTACTTTTCTTTTAAGACTACGACGGCCTGCATctctttcttcttctttgaagatgcggaagatcactccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                               <  2:976473‑M1/141‑94 (MQ=255)                                                                                                 cGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaaattcactgtggcatttcccatccgtgatcgtcgtgtgtttactcggcgtgagcgtgcg                                        >  2:234695‑M1/1‑43 (MQ=255)                                                                                                      cATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcc                                   >  1:113412‑M1/1‑38 (MQ=255)                                                                                                            cccgtTTCTGTGAATACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttt                             <  1:1040392‑M1/141‑112 (MQ=255)                                                                                                               gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcg                          <  1:148832‑M1/143‑115 (MQ=255)                                                                                                                tttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgg                         >  1:1158215‑M1/1‑28 (MQ=255)                                                                                                                 ttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccattcgtgatggtcgtctgattactctgcgtgagcgtgaactccccgcctttcgtc                        >  2:316654‑M1/1‑27 (MQ=255)                                                                                                                  tCTTTGAAGACTACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtct                       <  2:602619‑M1/143‑118 (MQ=255)                                                                                                                      tGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtat                   >  1:566362‑M1/1‑22 (MQ=255)                                                                                                                             gccgacggcCTGCATcgccttcttcttcttggaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgacctgcccgcttttcgtctgtatcgtatgc            <  1:1153284‑M1/141‑129 (MQ=255)                                                                                                                                       tGCTTctccttcttctgctgtgaagatgcggaagttcacgccggtttttccactgtggcatttaccatccgggatggtcgtctgtttactctgcgtgcgcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtag  <  2:758014‑M1/143‑139 (MQ=255)                                                                                                                                           |                                                                                                                                           ATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCC  >  NC_000913/4001505‑4001781

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