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 = 40016441 (0.020)37 (0.730) 33/282 0.6 97.4% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 1 (0.020)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

ATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGT  >  NC_000913/4001514‑4001785
                                                                                                                                        |                                                                                                                                       
gattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgcagcatcggcccgtgcctttcgagcccacgccggcctgctcCTCcttc                                                                                                                                   >  1:149446‑M2/137‑143 (MQ=255)
          gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggcactggaagacatcgaagcatcggcgcgggtctttgaagacgacgacggccggcacCTCCTTCGTCTTCTTTg                                                                                                                         >  2:286041‑M2/127‑143 (MQ=255)
          gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                         <  1:116115‑M2/17‑1 (MQ=255)
          gtcgaaccggacgacgacgagcgactgcgcggacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                         <  1:486220‑M2/17‑1 (MQ=255)
                ccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagatttcgacggcctgcatCTCCTTCTTCTTCTTTGAAGATg                                                                                                                   >  1:880747‑M2/121‑143 (MQ=255)
                ccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATg                                                                                                                   <  1:20747‑M2/23‑1 (MQ=255)
                 cggacggcgaccagcgactgcgcgtagcatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                  <  1:346754‑M2/24‑1 (MQ=255)
                         gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaatacatcgaagaaacggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCCGCAGATc                                                                                                          >  1:304080‑M2/112‑143 (MQ=255)
                             agccactgcgcgtacaatctgaactcggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGc                                                                                                      <  1:441671‑M2/36‑1 (MQ=255)
                              gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGcc                                                                                                     >  1:915252‑M2/107‑143 (MQ=255)
                               cgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccagatgtggttgacaacggatcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                    <  1:927147‑M2/38‑1 (MQ=255)
                                 actgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGCTCACGCCGGt                                                                                                  >  2:311839‑M2/104‑143 (MQ=255)
                                  ctgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTa                                                                                                 >  1:882223‑M2/103‑143 (MQ=255)
                                                    cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                               <  1:267168‑M2/59‑1 (MQ=255)
                                                         ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATc                                                                          >  1:61912‑M2/80‑143 (MQ=255)
                                                         ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATc                                                                          >  1:772519‑M2/80‑143 (MQ=255)
                                                         ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCATAATGGAAGATGCGGAAGAGCACGCCGGTGACTCCACTGTGGCAGTTACCATc                                                                          >  1:426751‑M2/80‑143 (MQ=255)
                                                          cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacggcgacggcctgcatCCCCTTCGTCTTCTTTGAAGATGCGGAAGAGAACGCCGGTACCGACGCTGTGGCCTGTaccacac                                                                         >  2:779359‑M2/79‑140 (MQ=255)
                                                                cgggcaacccccccggaacaggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                   <  1:229371‑M2/71‑1 (MQ=255)
                                                                  ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTATTCTTTGAAGATGCGGAAGAGCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGCtggg                                                                 >  1:82782‑M2/71‑142 (MQ=255)
                                                                     aacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCAGTGATGgtcgt                                                              >  1:240969‑M2/68‑143 (MQ=255)
                                                                       cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacgccccgaatATTCTTCTTCTGCTTTGAAGATGCGGCAGATCACGCCGGTACCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                            >  1:349843‑M2/66‑143 (MQ=255)
                                                                         cgcccggaactggaagacatcgacgcatcggcacgtttcttggaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                          <  1:286041‑M2/80‑1 (MQ=255)
                                                                                aactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActct                                                   <  1:437208‑M2/87‑1 (MQ=255)
                                                                                   tggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcg                                                <  1:532637‑M2/90‑1 (MQ=255)
                                                                                    ggaatacatcgaagcatcggcacgtttctttgaagtcgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                               <  1:465102‑M2/91‑1 (MQ=255)
                                                                                      aagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                             >  2:927019‑M2/51‑143 (MQ=255)
                                                                                      aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCAT‑‑‑‑‑CTCCTCCTTCTTCTTTGAAGGTGCGGAAGCTCACGCCGGTAACTCCGCTGTGGCATTTCCCACACGGCATGGCCGCGTGGTTACGCTgcgtgaccgcg                                        >  2:441940/1‑141 (MQ=255)
                                                                                           ggcgaagcatcggcgcgcttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCATGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtg                                        <  2:874030‑M2/98‑1 (MQ=255)
                                                                                                  catcggcacgtttctttgaagactacgacggcctgcatCTTCTTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                 <  2:304080‑M2/105‑1 (MQ=255)
                                                                                                   atcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGAGGGAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGccc                                >  2:155658‑M2/38‑143 (MQ=255)
                                                                                                               ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTg                    <  1:761022‑M2/118‑1 (MQ=255)
                                                                                                                tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGCGCGTGAACTCCCCGCTTTTCGTCTGt                   >  1:151201‑M2/25‑143 (MQ=255)
                                                                                                                  tgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                 >  1:327724‑M2/23‑143 (MQ=255)
                                                                                                                  tgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACGGCCCGCTGGTCGTCGGTAt                 >  2:312541‑M2/23‑143 (MQ=255)
                                                                                                                    aagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCg               >  2:536414‑M2/21‑143 (MQ=255)
                                                                                                                       acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTCCGTGCGCGTGCACTGCCCGCTTTTCGTCTGTATcgtgt            >  1:203053‑M2/18‑141 (MQ=255)
                                                                                                                                 cctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCCTCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATAGTATGCGTGCTCGt  >  1:705603‑M2/8‑143 (MQ=255)
                                                                                                                                        |                                                                                                                                       
ATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGT  >  NC_000913/4001514‑4001785

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