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.020)	46 (0.900)	39/284	0.3	97.9%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001650 =	1 (0.020)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

GTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4001504‑4001644 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001650‑4001783                                                                                                                                                                                                                                                                                     GTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTTCTTTGAAGACGACGACGGCCTGCCTCT                                                                                                                                      >  1:439767/1‑143   ATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCC                                                                                                                                    <  1:100500/143‑1     GGATTGATCTTGGCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTT                                                                                                                                  <  1:274519/143‑1       ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCT                                                                                                                                <  2:871038/143‑1         cGATCTGGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTGTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTC                                                                                                                              <  2:295614/142‑1           ATCTTGTCGAACCGGACGACTACGAGCGACTGCGCGCACAATCTGACCTTGGTCAGAGCCTGGCAACCCGCCCGGTACTGCAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTT                                                                                                                            <  1:126599/143‑1           ATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCCTCTCCTTCTTCTT                                                                                                                            >  2:375489/1‑143              TTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTT                                                                                                                         >  1:845575/1‑143                       CGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGC                                                                                                                >  2:88710/1‑143                       CGGACCACGACGAGCGACCGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGC                                                                                                                <  1:100912/143‑1                          ACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGA                                                                                                             <  1:90966/143‑1                            GACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGGCATCGAAGCATGGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGTAG                                                                                                           >  1:395422/1‑143                                 CGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGTCGTCGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCAC                                                                                                      <  1:496371/143‑1                                            GCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAccgc                                                                                           >  1:504686/1‑139                                             CGTACAATCCGAACTTGGCCACAGCCTGGCAACCCGCCCGGAACCGGAAGACATCGACGCATCGGCACGTTGCTCTGAAGACTACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCC                                                                                          <  1:795749/143‑1                                                  acTCTGAACCTCGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGTCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGT                                                                                     <  2:724880/141‑1                                                   cgcgGAACGTGGGCAGAGCCTGGCAACCCTTCCGGAACTGGAAGACATCGAAGCATCGGCACGGTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCCCTGTG                                                                                    <  2:663112/139‑1                                                   ATCGGAACTCTGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGCCATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTG                                                                                    <  1:541962/143‑1                                                     gTGAGCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGC                                                                                  <  2:607012/142‑1                                                      TGAACTTGGCCCGAGCCTGGCAAGGTGTATGATAGTTGAAAACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACGCCACTGTGGCA                                                                                 >  2:653684/1‑143                                                      TGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCA                                                                                 >  1:677794/1‑143                                                           TTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTGCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAC                                                                            >  1:637322/1‑143                                                                  GAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCACCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGc                                                                     >  2:743438/1‑142                                                                                GCCCGGAACTGGAAGACATCGCAGCATCGGCACGTTTCTTTGACGATGACGAGGTTCTTCATCTTCTTCTTCTTCTTTGAAGAGGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTT                                                       <  1:663842/143‑1                                                                                           GAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTG                                            >  1:874556/1‑143                                                                                            AAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCGGTGATGGTCGTCTGTTTACTCTGCGTGg                                           >  1:842248/1‑142                                                                                             AGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG                                          >  1:431710/1‑143                                                                                             AGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAG                                          >  2:467181/1‑143                                                                                                 ATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGCGCGTG                                      >  2:623549/1‑143                                                                                                 ATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTG                                      <  1:600246/143‑1                                                                                                 ATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTGACTCTGCGTGAGCGTG                                      >  1:26508/1‑143                                                                                                   CGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCCTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGGGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAA                                    <  1:45346/143‑1                                                                                                           CGGCACGTTGCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGC                            <  1:54799/143‑1                                                                                                             GCACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTT                          >  1:435080/1‑143                                                                                                              CACGTTTCTTTGAAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGCTCACGCCGGTCACTCCACTGTGGCATTTACCATCCGTGATGCCCGTCTGTTTACTCCGCGTGACCGTCACCCGCCCGCTTT                         >  2:417155/1‑143                                                                                                                CGTTTCTTTGAAGACTACGCCGGTCTGCGTCGCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAGTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTC                       <  2:170814/143‑1                                                                                                                CGTTTCTGTTAAGTCGCTGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTC                       <  2:867386/143‑1                                                                                                                  cTTCTTTGATGACGACGACTGTCTGCATCTTCTTCTTCTTCTTTGAGGTTGCGGAAGATCACGTCGGTAACTTCACTGTGTCATTTGCCATCCGTGGTGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGT                     <  2:695443/142‑1                                                                                                                    tcgTTGAGGACGGCGGCTGCCTGCATCCCCTTCTTCTTCTTTGAATATGCGGAAGATCACGCCGGTAACTACACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCT                   <  2:696363/140‑1                                                                                                                    TCTTTGGAGACGACGACGGCCTGCCTCGCCTTGTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCT                   <  1:153374/143‑1                                                                                                                       TTGAAGACGACGACGGCCTGCCTCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTA                <  2:287557/143‑1                                                                                                                          AAGACGACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCG             <  1:323781/143‑1                                                                                                                             gCGACGGCGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAT          <  2:526938/142‑1                                                                                                                               GACGACGGCCTGCATCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACGATAGGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATAGTATGC        >  1:783410/1‑143                                                                                                                                   ACGGCCTGCATCTCCTTCGTCTTCGTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGC    <  1:375489/143‑1                                                                                                                                     GGCCTGCATCGCCTCCCTCTTCGGTGGAGATGCGGAAGATCACGCCGGTAACTCCGCTGTGGCATGGACCATCCGTGATGGTCGTCTGTTTTCTCTGCGTGAGCGTTAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  <  1:175300/143‑1                                                                                                                                                                                                                                                                                     GTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4001504‑4001644 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001650‑4001783

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.