Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	4,001,645	Δ5 bp	coding (220‑224/951 nt)	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.030)	31 (0.830)	25/284	0.5	96.9%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001650 =	1 (0.030)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTA  >  NC_000913/4001510‑4001786                                                                                                                                       |                                                                                                                                               aTTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGCAGACGACGACGGCCTCCATctccttct                                                                                                                                        >  1:142262‑M1/1‑135 (MQ=255)      tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGACCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCCGCATctccttcttctcc                                                                                                                                   >  2:394009‑M1/1‑130 (MQ=255)        ttGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGGCGGCGGACCGCATctccttctccttctt                                                                                                                                 >  2:387452‑M1/1‑128 (MQ=255)        ttGTCAAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGTCATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACCGCCTGCATctccctcctcttcct                                                                                                                                 >  2:560040‑M1/1‑128 (MQ=255)          gTCGATCCTGACGACGACGAGCGACTGCGCGGACAAGCTGAACGTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttg                                                                                                                               <  1:168513‑M1/143‑18 (MQ=255)            cggaccgGACGACGACGAGCGACTGCGCGTCGAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                             <  1:413605‑M1/140‑20 (MQ=255)            cGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                             <  1:266309‑M1/143‑20 (MQ=255)                            gaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctcctccttcttctctacagctacggaagagcaca                                                                                                             >  2:215079‑M1/1‑108 (MQ=255)                              gCGACGGCGCGTAGAATCTGAACTTGGCAAGAGCCTGGCAACCCGCCCGGAACTGGAAGACAGCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgcc                                                                                                           <  1:429964‑M1/143‑38 (MQ=255)                                   cgcccgtagaCTCCGAGCTTGGCCGGAGCCTGGCCACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaa                                                                                                      <  2:578820‑M1/139‑43 (MQ=255)                                     cgcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagctgcggaagatcacgccggtaact                                                                                                    >  1:622957‑M1/1‑99 (MQ=255)                                             aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                            <  2:239298‑M1/143‑53 (MQ=255)                                             aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                            <  2:238137‑M1/143‑53 (MQ=255)                                                               ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgat                                                                          >  1:432714‑M1/1‑73 (MQ=255)                                                               ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgat                                                                          >  1:655040‑M1/1‑73 (MQ=255)                                                                 tGGCAACCCGCCCGGAACTGGAAGACATTGAAGCATCGGTACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcgtttaccatccgtcatgg                                                                        >  1:197509‑M1/1‑71 (MQ=255)                                                                 tGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgg                                                                        <  1:187584‑M1/143‑73 (MQ=255)                                                                   gCAATCCGCCCGGAACTGGAAGACATCGAAGCATCGGCAGGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttgtttgaagatgcgaaagtgcacgccggtatctccactgtggcatttaccatccgtgatggtc                                                                      <  2:22066‑M1/143‑75 (MQ=255)                                                                            ccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTCGAAGACGACGCCGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgttta                                                             <  1:489502‑M1/143‑84 (MQ=255)                                                                              ggAACTGGAAGACATCGACGCATCGGCACTTTTCTTTGAAGCC‑‑C‑‑CTGCCTTGCGACAC‑CTCCTTCTTGGGCGAGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAct                                                      <  2:297632/143‑1 (MQ=255)                                                                              ggAACGGGAAGACAGCGAAGCATCGGAACGTTTCTGGGAAGACGACGACGGCCTGCAGctccttcttcttcttttaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttact                                                           <  1:8618‑M1/143‑86 (MQ=255)                                                                                   tGGAAGACATCGAAGCTGCGGTTCGTTTGTTGGAAGGCGACGACGGCCTGCATctccgtctgcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                                      <  1:1242‑M1/143‑91 (MQ=255)                                                                                   tGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                                      <  1:653295‑M1/143‑91 (MQ=255)                                                                                    ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                     <  2:451344‑M1/143‑92 (MQ=255)                                                                                        gACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtgactccactgtggcatttatcatccgtgatggtcgtctgtttactctgagtgagc                                                 >  1:608054‑M1/1‑48 (MQ=255)                                                                                         aCATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggacgatcacgccggtaactccactgtggcatttaccttccgtgatggtcgtctgtttactctgcctgagcc                                                >  1:107106‑M1/1‑47 (MQ=255)                                                                                           aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtg                                              >  1:45374‑M1/1‑45 (MQ=255)                                                                                           aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATatcctacttattatttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgggagcgtg                                              >  1:37876‑M1/1‑45 (MQ=255)                                                                                                      ggCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgct                                   <  1:365331‑M1/143‑110 (MQ=255)                                                                                                                 ttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgta                        >  1:80818‑M1/1‑23 (MQ=255)                                                                                                                         gacgccgGCCTGGAGctcctgcttcttcattgaggatgcggaagctcacgccggtgacggcactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc                <  2:39919‑M1/143‑129 (MQ=255)                                                                                                                                       tctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgcccgtagcgag  >  1:379793‑M1/1‑1 (MQ=255)                                                                                                                                       |                                                                                                                                               ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTA  >  NC_000913/4001510‑4001786

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: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.