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	45 [0]	[0] 45	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)	44 (0.900)	38/282	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

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

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