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

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4001644	1 (0.020)	38 (0.940)	30/284	0.4	98.7%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001650 =	0 (0.000)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCC  >  NC_000913/4001505‑4001782                                                                                                                                            |                                                                                                                                           tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctc                                                                                                                                         <  1:46913‑M1/143‑4 (MQ=255) tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctc                                                                                                                                         >  2:537968‑M1/1‑140 (MQ=255)  aTGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTTCTTTGAAGACGACGACGGCCTGCatatat                                                                                                                                        >  1:265367/1‑141 (MQ=255)      aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                                    >  1:597435‑M1/1‑135 (MQ=255)       gttatccTGTCGATCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCTGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctt                                                                                                                                   <  2:254910‑M1/140‑10 (MQ=255)             ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                             <  2:281451‑M1/143‑16 (MQ=255)             ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTCGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctt                                                                                                                             <  2:280055‑M1/143‑16 (MQ=255)              tGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttt                                                                                                                            <  2:618503‑M1/143‑17 (MQ=255)               gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                           <  2:170553‑M1/143‑18 (MQ=255)               gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                           <  1:687021‑M1/143‑18 (MQ=255)                      cGGACGACGATGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                    >  1:525956‑M1/1‑119 (MQ=255)                      cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGACGGCCTGCATctccttcttcttctttgcagatgc                                                                                                                    >  2:91587‑M1/1‑119 (MQ=255)                      cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATccccttcttcttctttgaagatgc                                                                                                                    >  2:624024‑M1/1‑119 (MQ=255)                       ggccgacgacgagcgACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCCCCCGGAACCGGAAGACATCGAAGCATCGGCACGTTCCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcg                                                                                                                   <  1:359418‑M1/140‑26 (MQ=255)                          cgacgacggGCGACTGCGCGTACAATCTGAAGTTGGCCGGAGCCTGGCAGCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                <  2:543339‑M1/143‑29 (MQ=255)                              gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGCCGACGGCCTGCATctccttcttcttctttgaagctgcggaagatc                                                                                                            >  1:608514‑M1/1‑111 (MQ=255)                              gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGCCGGCCTGCATctccttcttcttctttgaagatgcggaagctc                                                                                                            >  1:50861‑M1/1‑111 (MQ=255)                              gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatc                                                                                                            >  2:98952‑M1/1‑111 (MQ=255)                                cgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagagcac                                                                                                          >  1:249067‑M1/1‑109 (MQ=255)                                                cAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcaggtaactccactg                                                                                          >  1:210691‑M1/1‑93 (MQ=255)                                                  aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggcaactccaccgtg                                                                                        >  2:70410‑M1/1‑91 (MQ=255)                                                    cTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggcaactccactgtggc                                                                                      >  2:650227‑M1/1‑89 (MQ=255)                                                      gAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccctcttcttctttgcagtcgcggaagatcacgccggtaactccactgtggcat                                                                                    <  1:503222‑M1/143‑57 (MQ=255)                                                         cTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattta                                                                                 <  1:101908‑M1/143‑60 (MQ=255)                                                                   gCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtga                                                                       <  1:599273‑M1/143‑70 (MQ=255)                                                                      tGGCCACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACCGCCTGCATctccttcttcttctgtgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgg                                                                    <  1:479781‑M1/143‑73 (MQ=255)                                                                            cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtatctccactgtggcatttaccatccgtgatggtcgtct                                                              <  1:473190‑M1/143‑79 (MQ=255)                                                                            cccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                              <  1:537968‑M1/143‑79 (MQ=255)                                                                                 ccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaacatccgtgatggtcgtctgttta                                                         >  2:71909‑M1/1‑60 (MQ=255)                                                                                         ggaggaCATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcgtcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                 <  2:689045‑M1/143‑92 (MQ=255)                                                                                         ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                 <  2:472971‑M1/143‑92 (MQ=255)                                                                                           aaGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                               <  1:596327‑M1/143‑94 (MQ=255)                                                                                                aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtg                                          >  1:117238‑M1/1‑45 (MQ=255)                                                                                                      gCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgc                                    >  1:226041‑M1/1‑39 (MQ=255)                                                                                                           ggCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggcaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcctgagcgtgagctgcccgct                               >  2:510507‑M1/1‑34 (MQ=255)                                                                                                                   tCTTGGGAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaggtctcgccggtaacacaactgtgcccttttccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtct                       <  2:10034‑M1/143‑118 (MQ=255)                                                                                                                       agaAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgtcggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtat                   <  1:77251‑M1/142‑122 (MQ=255)                                                                                                                           gtcgaccacgccCTGCATctccttcttctgcttttgagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgta               <  2:37404‑M1/141‑126 (MQ=255)                                                                                                                                        tGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtac  >  1:707918‑M1/1‑5 (MQ=255)                                                                                                                                            |                                                                                                                                           TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCC  >  NC_000913/4001505‑4001782

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑