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	36 [0]	[0] 35	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)	35 (0.870)	32/284	0.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

CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCG  >  NC_000913/4001502‑4001784                                                                                                                                               |                                                                                                                                             ctggctGGATTTATCTTTTCGATCCGGACGACGAGGGGCGATTGCGCGTACAATCTGAACTTGGCTAGAGCCTGGCAAGCCTCTCGTAACTGGCAGACATCGAAGCATCGGCACGTTTCTTTGAAGACTACAACGGCCTGCat                                                                                                                                              <  2:151686/138‑1 (MQ=255)     gcagtttgctctCGTCGAAGCTGACAACGCCGAGCGACTGCTCGTACAATCTGAACTTGGCCAGAGCCTGGCACCCCGCCCGGAACTGGAATACATCGAAGCATCGGTACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                          <  2:212600‑M1/138‑5 (MQ=255)         attactcttgtCTACCCGGACGACGACGAGCGACTGCGCGTCCAACCCGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                                      <  1:345550‑M1/138‑9 (MQ=255)              tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttc                                                                                                                                 >  1:424931‑M1/1‑130 (MQ=255)                   ccgAACCGGACGACGGCGAGCGACTGCGCGTACAATCTGAACCTGGCCAGAGCCTGGCAGCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                            <  1:464990‑M1/142‑19 (MQ=255)                    cGACCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                           <  1:81461‑M1/143‑20 (MQ=255)                    cGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                           <  1:581910‑M1/143‑20 (MQ=255)                      aaCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaaga                                                                                                                         <  1:238631‑M1/143‑22 (MQ=255)                             cgacgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                  <  1:352261‑M1/143‑29 (MQ=255)                                    gaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctcctttgaagatgcggaagatcacg                                                                                                           >  2:466608‑M1/1‑108 (MQ=255)                                    gaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctccttctctgaacatgcggaacgtcacg                                                                                                           >  2:584716‑M1/1‑108 (MQ=255)                                        gACTGCGCGGACAAGCTGAAGTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgg                                                                                                       <  1:316637‑M1/143‑40 (MQ=255)                                              gcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgacgctgcggacggtcacgccggtagctg                                                                                                 >  1:296069‑M1/1‑98 (MQ=255)                                                     aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                          <  1:237793‑M1/143‑53 (MQ=255)                                                       cTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgtggaagatcacgccggtaactccactgtggc                                                                                        >  1:395295‑M1/1‑89 (MQ=255)                                                        tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGCAGCATCGGCACGTTTCTTGGAAGACGACGACGGCCTGCATctccttcttcttctttgaagaggcggaagatcacgccggtaaagccactgagcgc                                                                                       >  1:670761‑M1/1‑88 (MQ=255)                                                                 cccgagcctGGCAACCCGCCCGGAACTGGACGCCATCGACGCATCGGCACGTTTCTTTGAAGACGACGACCGCCTGCATctccttcttcttctttgcagatgcggaagatcacgccggtaactccactgtggcatttaccatc                                                                              <  1:176229‑M1/140‑65 (MQ=255)                                                                  cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcc                                                                             >  2:517973‑M1/1‑78 (MQ=255)                                                                    gagCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgt                                                                           <  1:176889‑M1/143‑68 (MQ=255)                                                                             aaCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttggaagatgcggaagatcacgccggtaactccactctggcatttaccatccgtgatggtcgt                                                                  <  2:694786‑M1/143‑77 (MQ=255)                                                                               cccgcccgGAGCGGGAAGACATCGAAGCATCGGCACGTTTCTTTTAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                                <  2:402298‑M1/143‑79 (MQ=255)                                                                                       gAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgcgatggtcgtctgtttactc                                                        >  2:531430‑M1/1‑57 (MQ=255)                                                                                        aaCTGGAAGACAGCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactct                                                       <  1:111158‑M1/143‑88 (MQ=255)                                                                                            ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                   <  1:13540‑M1/143‑92 (MQ=255)                                                                                            ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcgtctttgaagatgcggaagatcacgccggtacctccgctgtggcatttgccatccgtgatggtcggctgtttactctgcgt                                                   >  1:549282‑M1/1‑52 (MQ=255)                                                                                                gACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagc                                               >  1:252733‑M1/1‑48 (MQ=255)                                                                                                   aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatacgtgatggtcgtctgtttactctgcgtgagcgtgaac                                            >  1:177761‑M1/1‑45 (MQ=255)                                                                                                               gCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgtggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                <  2:395295‑M1/143‑111 (MQ=255)                                                                                                               gCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                <  2:675406‑M1/143‑111 (MQ=255)                                                                                                                       cTTTGAAGACGACGCCGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtacctccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctg                        <  1:315235‑M1/143‑119 (MQ=255)                                                                                                                        tttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttctcggctgt                       >  1:177128‑M1/1‑24 (MQ=255)                                                                                                                               acgacgtcgGCCTGCATctcctttttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctatttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtat                <  1:696586‑M1/143‑127 (MQ=255)                                                                                                                                 tgcgtcggcctgCATctccttcttctgccttgacgttgcggcagctctcgccggtacctcctctgtggcacttaccctccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc              <  2:722643‑M1/138‑129 (MQ=255)                                                                                                                                     gcggCATGCATctccttcttcttctctgaagatgcggaagatcacgccggtaactcccctgtggcatttaccatccgtgatggtcgtttgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgc          <  1:145406‑M1/142‑133 (MQ=255)                                                                                                                                            gCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgcgcgtgaactgcccgcttttcgtctgtatcgtatgcgcgcccgtagc   >  2:615002‑M1/1‑4 (MQ=255)                                                                                                                                             cATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaaatattcgctgcccgtctgtatcgtatgcgtgcccgttgcc  >  1:302224‑M1/1‑3 (MQ=255)                                                                                                                                               |                                                                                                                                             CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCG  >  NC_000913/4001502‑4001784

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