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

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

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

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