Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 4,001,644 Δ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 = 40016431 (0.020)47 (0.890) 39/284 0.3 97.9% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 1 (0.020)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

GATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001508‑4001782
                                                                                                                                             |                                                                                                                                     
gtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacggttctttgaagacgacgacggcctgcctct                                                                                                                                      >  2:439767‑M2/142‑143 (MQ=255)
  atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                    <  2:100500‑M2/4‑1 (MQ=255)
    ggattgatcttggcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCctt                                                                                                                                  <  2:274519‑M2/6‑1 (MQ=255)
      attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttct                                                                                                                                <  1:871038‑M2/8‑1 (MQ=255)
        cgatctggtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtgtctttgaagacgacgacggcctgcatCTCcttcttc                                                                                                                              <  1:295614‑M2/10‑1 (MQ=255)
          gTCGAAGCGGAAGACGAGGCGCGTCTGCGCGGTCCGTCTGAAGTTGACTGGGGCCTGGCACTACGCCCGGACCTGGCAGACCTGGAAGCATCGGCACGTTTCTTTGAAGCCGACATCGGCCTGCAT‑‑‑‑‑CTCCTTCTTCTTCTTTg                                                                                                                       <  1:324328/143‑1 (MQ=255)
          atcttgtcgaaccggacgactacgagcgactgcgcgcacaatctgaccttggtcagagcctggcaacccgcccggtactgcaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttctt                                                                                                                            <  2:126599‑M2/12‑1 (MQ=255)
          atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcctCTCcttcttctt                                                                                                                            >  1:375489‑M2/132‑143 (MQ=255)
             ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                         >  2:845575‑M2/129‑143 (MQ=255)
                      cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                >  1:88710‑M2/120‑143 (MQ=255)
                      cggaccacgacgagcgaccgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                <  2:100912‑M2/24‑1 (MQ=255)
                         acgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGa                                                                                                             <  2:90966‑M2/27‑1 (MQ=255)
                           gacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaaggcatcgaagcatgggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGTAg                                                                                                           >  2:395422‑M2/115‑143 (MQ=255)
                                cgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagtcgtcgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCAc                                                                                                      <  2:496371‑M2/34‑1 (MQ=255)
                                           gcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacggttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCcggtaccgc                                                                                           >  2:504686‑M2/99‑139 (MQ=255)
                                            cgtacaatccgaacttggccacagcctggcaacccgcccggaaccggaagacatcgacgcatcggcacgttgctctgaagactacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTcc                                                                                          <  2:795749‑M2/46‑1 (MQ=255)
                                                 actctgaacctcgccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacgtcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgt                                                                                     <  1:724880‑M2/51‑1 (MQ=255)
                                                  cgcggaacgtgggcagagcctggcaacccttccggaactggaagacatcgaagcatcggcacggttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCCCtgtg                                                                                    <  1:663112‑M2/52‑1 (MQ=255)
                                                  atcggaactctgccagagcctggcaacccgcccggaactggaagccatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                    <  2:541962‑M2/52‑1 (MQ=255)
                                                    gtgagcttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGc                                                                                  <  1:607012‑M2/54‑1 (MQ=255)
                                                     tgaacttggcccgagcctggcaaggtgtatgatagttgaaaacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACGCCACTGTGGCa                                                                                 >  1:653684‑M2/89‑143 (MQ=255)
                                                     tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCa                                                                                 >  2:677794‑M2/89‑143 (MQ=255)
                                                          ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacggtgctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAc                                                                            >  2:637322‑M2/84‑143 (MQ=255)
                                                                 gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcacCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCcgc                                                                     >  1:743438‑M2/77‑142 (MQ=255)
                                                                               gcccggaactggaagacatcgcagcatcggcacgtttctttgacgatgacgaggttcttcatCTTCTTCTTCTTCTTTGAAGAGGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                       <  2:663842‑M2/81‑1 (MQ=255)
                                                                                          gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                            >  2:874556‑M2/52‑143 (MQ=255)
                                                                                           aagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCGGTGATGGTCGTCTGTTTACTCtgcgtgg                                           >  2:842248‑M2/51‑142 (MQ=255)
                                                                                            agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                          >  2:431710‑M2/50‑143 (MQ=255)
                                                                                            agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                          >  1:467181‑M2/50‑143 (MQ=255)
                                                                                                atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgcgcgtg                                      >  1:623549‑M2/46‑143 (MQ=255)
                                                                                                atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtg                                      <  2:600246‑M2/98‑1 (MQ=255)
                                                                                                atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTGACTCTgcgtgagcgtg                                      >  2:26508‑M2/46‑143 (MQ=255)
                                                                                                atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACGCCACCGTGGCATTTACCATCCGTGATggtcgtctgtttactacgagtccgcggg                                      >  1:97931‑M2/46‑130 (MQ=255)
                                                                                                  cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCCTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGGGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                    <  2:45346‑M2/100‑1 (MQ=255)
                                                                                                          cggcacgttgctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGc                            <  2:54799‑M2/108‑1 (MQ=255)
                                                                                                            gcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtt                          >  2:435080‑M2/34‑143 (MQ=255)
                                                                                                             cacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGCTCACGCCGGTCACTCCACTGTGGCATTTACCATCCGTGATGCCCGTCTGTTTACTCCGCGTGACCGTCACCCGCCCGCttt                         >  1:417155‑M2/33‑143 (MQ=255)
                                                                                                               cgtttctttgaagactacgccggtctgcgtCGCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAGTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTc                       <  1:170814‑M2/113‑1 (MQ=255)
                                                                                                               cgtttctgttaagtcgctgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTc                       <  1:867386‑M2/113‑1 (MQ=255)
                                                                                                                 cttctttgatgacgacgactgtctgcatCTTCTTCTTCTTCTTTGAGGTTGCGGAAGATCACGTCGGTAACTTCACTGTGTCATTTGCCATCCGTGGTGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGt                     <  1:695443‑M2/115‑1 (MQ=255)
                                                                                                                   tctttggagacgacgacggcctgcctCGCCTTGTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                   <  2:153374‑M2/117‑1 (MQ=255)
                                                                                                                   tcgttgaggacggcggctgcctgcatCCCCTTCTTCTTCTTTGAATATGCGGAAGATCACGCCGGTAACTACACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                   <  1:696363‑M2/117‑1 (MQ=255)
                                                                                                                      ttgaagacgacgacggcctgcctCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTa                <  1:287557‑M2/120‑1 (MQ=255)
                                                                                                                         aagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCg             <  2:323781‑M2/123‑1 (MQ=255)
                                                                                                                            gcgacggcggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAt          <  1:526938‑M2/126‑1 (MQ=255)
                                                                                                                              gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACGATAGGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATAGTATGc        >  2:783410‑M2/16‑143 (MQ=255)
                                                                                                                                  acggcctgcatCTCCTTCGTCTTCGTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGc    <  2:375489‑M2/132‑1 (MQ=255)
                                                                                                                                    ggcctgcatCGCCTCCCTCTTCGGTGGAGATGCGGAAGATCACGCCGGTAACTCCGCTGTGGCATGGACCATCCGTGATGGTCGTCTGTTTTCTCTGCGTGAGCGTTAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGccc  <  2:175300‑M2/134‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                     
GATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCC  >  NC_000913/4001508‑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: 
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.