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	47 [0]	[0] 47	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)	46 (0.910)	39/282	0.2	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

CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001501‑4001782                                                                                                                                               |                                                                                                                                            cAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCat                                                                                                                                             <  2:810753/143‑1 (MQ=255)    gatGGATTGATCTTGTCGAATCGGAGGACGACGAGCGACTGCGCGTACAATTTGAGTTTGGCCAGAGCCTGGCAACCCGCCCGGAACGGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctc                                                                                                                                          <  1:674111‑M1/142‑4 (MQ=255)     attgattGATCTCGTCGAACGGGACGACGCGGAGCGACTGTGCGGACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                         <  1:524534‑M1/140‑5 (MQ=255)     aTGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                         <  2:258844‑M1/143‑5 (MQ=255)     aTGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                         <  1:590992‑M1/143‑5 (MQ=255)           tGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGACGCCGACGACGGCCTGCATctccttcttc                                                                                                                                   >  2:738868‑M1/1‑133 (MQ=255)               cTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGCCGCCGGCCTGCATctccttcttcttcg                                                                                                                               >  2:876553‑M1/1‑129 (MQ=255)                ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCACCCCGCCCGGAACTGGAACCCATCGAAGCATCTGGACGTTTCTTTGAAGACGACGACGGCAAGCACctgcttcttcctcgt                                                                                                                              >  2:444150‑M1/1‑128 (MQ=255)                  gTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                            >  2:586303‑M1/1‑126 (MQ=255)                  gTCGAACCGCACGCCGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                            <  2:765966‑M1/143‑18 (MQ=255)                         cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGTCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                     <  1:238931‑M1/143‑25 (MQ=255)                             cgacgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGCAGACATCGAAGCATCGGCACGTCCCTTTGAAGACGACGACGGCCTGCATctacttcttcttctccgccgaggcgggg                                                                                                                 >  1:478378‑M1/1‑115 (MQ=255)                               acgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGACAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                     >  2:569840‑M1/1‑113 (MQ=255)                               acgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGACAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                     <  1:569840‑M1/137‑25 (MQ=255)                                 gacgaGCGACGGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatc                                                                                                             <  2:348622‑M1/143‑33 (MQ=255)                                  acgaGCGACTGCGCGTACAATCTGGACTTGGCCAGAGCCTGGCAACTCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatca                                                                                                            <  1:212434‑M1/143‑34 (MQ=255)                                          cTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggta                                                                                                    <  1:311497‑M1/143‑42 (MQ=255)                                              gcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactc                                                                                                <  1:349574‑M1/143‑46 (MQ=255)                                                    aaTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTATGAAGACGACGACGGCCTGCATctcctgcttcttctttgaagatgcgcaagatcacgccggtaactccactgt                                                                                          >  1:753765‑M1/1‑92 (MQ=255)                                                       cggaaCTGGACGAGGGCCTGGCGACCCGCCCGGAACTGGAAGACATCGAAGCAGCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggc                                                                                       <  1:807437‑M1/141‑55 (MQ=255)                                                        tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgccgaagatcacgccggtaactccactgtggca                                                                                      >  2:819508‑M1/1‑88 (MQ=255)                                                        tGAACTGGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggca                                                                                      <  1:505712‑M1/143‑56 (MQ=255)                                                             ttGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttac                                                                                 >  2:522847‑M1/1‑83 (MQ=255)                                                                   agagCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGGAGACGACGACGGCCTGCATctccttcttcttctttgaagctgcggaagatcacgccggtaactccactgtggcatttcccatcgg                                                                           >  2:44506‑M1/1‑77 (MQ=255)                                                                        ccggcaACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatg                                                                      <  2:83286‑M1/141‑72 (MQ=255)                                                                        cTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatg                                                                      <  2:381215‑M1/143‑72 (MQ=255)                                                                          ggCAACCCGCCGGGACCTGCAAGACATCGAAGCAGCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctttgaagatgcggaacatcacgccggtaactccactgtggcatttaccatccgtgatggtcgt                                                                    <  2:848240‑M1/143‑74 (MQ=255)                                                                           gCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgcgacggcc                                                                   >  2:847053‑M1/1‑69 (MQ=255)                                                                                    ccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtcggttta                                                          >  2:359191‑M1/1‑60 (MQ=255)                                                                                     cgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGCCGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtgccagttccaacccgtgctggtcgcaggtttac                                                         >  1:204422‑M1/1‑59 (MQ=255)                                                                                       gAACTGGAAGACATCGGAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactc                                                       <  2:508025‑M1/143‑87 (MQ=255)                                                                                        aaCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactct                                                      >  2:294693‑M1/1‑56 (MQ=255)                                                                                          cTGGAAAACATCGAAGCATCAGCACGCTGCTTGGAAGACGACGACGGCCTGCAGctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgc                                                    <  2:63939‑M1/143‑90 (MQ=255)                                                                                           tGGAAGACATCCAAGCATCGGCACGCTTCTTTGAAGACGACGACGGCCTGCATctcctgcttcttctttgaagatgcggaagatcacgccggtaactccactgcggcatttaccatccgtgatggtcgtctgtttactctgcg                                                   <  2:271794‑M1/143‑91 (MQ=255)                                                                                            ggAAGACATCGAAGCGTCGGCACGTTTCTGTGAAGACGACGACGGCCTGCATctcctgcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                  <  2:356857‑M1/143‑92 (MQ=255)                                                                                               aGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgggcgtctgtttactctgtgtggg                                               >  1:494843‑M1/1‑49 (MQ=255)                                                                                                  cATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccaccgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgt                                            <  2:826578‑M1/143‑98 (MQ=255)                                                                                                                    tttCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcgcgcttttggg                          >  2:135803‑M1/1‑28 (MQ=255)                                                                                                                      tCTTTGAAGCCGACGGCGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtct                        <  2:433190‑M1/143‑118 (MQ=255)                                                                                                                      tCTGTGAAGACGACGACGGCCTGCATctccttctccttctttgaagatgcggaagatcacgccggtaactcctctgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtct                        <  2:190590‑M1/143‑118 (MQ=255)                                                                                                                             aGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaacatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgt                 >  2:848618‑M1/1‑19 (MQ=255)                                                                                                                              gacgacgacGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgta                >  2:424303‑M1/1‑18 (MQ=255)                                                                                                                               acggcgacgGCCTGCATcgccttcgacggctcggaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtat               <  2:65570‑M1/143‑127 (MQ=255)                                                                                                                                       gggcggcagctccgtcttcttcttggaagatgcggacgatcatgccggtaactctacggtggcattaaccagccgtgatggtcgtctgttgactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgccc       <  1:770012‑M1/138‑135 (MQ=255)                                                                                                                                        gCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccg      >  2:822132‑M1/1‑8 (MQ=255)                                                                                                                                            gCATctccttcttgttctttgaagatgcggaagatttcgccgttaattccactgtggcatttaccatccgtgttggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagc  <  1:418510‑M1/143‑140 (MQ=255)                                                                                                                                            gCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgacctgcccgcttttcttctgtatcgtctgcgtgccctcagc  >  1:105098‑M1/1‑4 (MQ=255)                                                                                                                                               |                                                                                                                                            CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001501‑4001782

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