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)48 (0.850) 41/284 0.4 98.0% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 1 (0.020)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGC  >  NC_000913/4001508‑4001788
                                                                                                                                              |                                                                                                                                          
agtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagagatcgacgcatcggcccgtttctttgaagacgacgccggcctgcatc                                                                                                                                            >  1:843940‑M2/143‑143 (MQ=255)
agtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaacgggaagacatcgaagcatcggcacgtttctttgaaggcgacgacggcctgcatc                                                                                                                                            >  1:362314‑M2/143‑143 (MQ=255)
    tggattgatcttggcgaaccggacgacgaccagcgactgcgcgtacaagctgaactgggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCct                                                                                                                                        <  1:376610‑M2/5‑1 (MQ=255)
      ggttgatcgtgccgaaccggacgactacgagcgacggcgcgtacagtctgaacttggccagagcctggctttccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                      <  1:517725‑M2/7‑1 (MQ=255)
             cttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                               <  2:417062‑M2/14‑1 (MQ=255)
                 tcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                           <  1:492014‑M2/18‑1 (MQ=255)
                     accggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGAt                                                                                                                       >  1:939095‑M2/122‑143 (MQ=255)
                      caccgacgaGGTACTGCGCCTACAAGCTGAACCTGGCCGGAGCCGGGCATGCCGCCCGGAACTGGAAGACATCGGAGGATCGGCACGTTTCTTCGAAGACGCCGACGGCCTGCAT‑‑‑‑‑CTCCTTCTTCTTCTTTGAAGATGCTGaa                                                                                                                 <  2:887298/140‑1 (MQ=255)
                         gacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctgtgcagacgacgacggcctgcatCTGATTCGTCGTCTTTGGAGACGCgg                                                                                                                   >  2:158209‑M2/118‑143 (MQ=255)
                            cgcgacgggcgactgcgcctacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgcagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAg                                                                                                                <  2:398542‑M2/29‑1 (MQ=255)
                               gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGGAGAtgcggacgctc                                                                                                             >  1:595345‑M2/112‑138 (MQ=255)
                                   agcgactgcgcgtacaatctgaactgggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGc                                                                                                         <  1:495756‑M2/36‑1 (MQ=255)
                                    gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGcc                                                                                                        <  1:695324‑M2/37‑1 (MQ=255)
                                    gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGAGGCGGGAGATCACGcc                                                                                                        >  1:768736‑M2/107‑143 (MQ=255)
                                        ctgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACgccgcga                                                                                                    >  1:727386‑M2/103‑140 (MQ=255)
                                         tgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                   >  1:807247‑M2/102‑143 (MQ=255)
                                         tgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                   >  1:1015653‑M2/102‑143 (MQ=255)
                                           cgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                 <  2:694521‑M2/44‑1 (MQ=255)
                                            gcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatagcagcatcggcacgggtctttggagacgacgacggccggcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAGCTCGCCGGGAACTc                                                                                                >  1:49940‑M2/99‑143 (MQ=255)
                                              gtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagaggacgacggcctgcatCTCCTTCTTCTTCTTTGAAGCTGCGGAAGATCACGCCGGTAACTCCa                                                                                              >  1:508066‑M2/97‑143 (MQ=255)
                                                    tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGACGATGCGGAAGATCACGCCGGTAACTCCACTGTgg                                                                                        >  1:919780‑M2/91‑143 (MQ=255)
                                                      tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCa                                                                                      <  2:908132‑M2/55‑1 (MQ=255)
                                                      tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCa                                                                                      >  1:495573‑M2/89‑143 (MQ=255)
                                                        aacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCAtt                                                                                    <  1:454770‑M2/57‑1 (MQ=255)
                                                          cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                  <  1:247386‑M2/59‑1 (MQ=255)
                                                           ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCGTCTTCTTCTTTGAAGATGCGGGAGATGACGCCGGTAACTCCACTTTGGCATCTAc                                                                                 >  1:262865‑M2/84‑143 (MQ=255)
                                                                           aaccggcccggaactggaagacatcgaagcatcggcacgtttctttgaagactacgacggcctgcacCTCTTAATACTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTTGCATTTACCATCCGTGATGgtcgt                                                                 <  2:254182‑M2/76‑1 (MQ=255)
                                                                              ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                              >  2:787177‑M2/65‑143 (MQ=255)
                                                                                  ccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                          >  2:304766‑M2/61‑143 (MQ=255)
                                                                                          ggaagacatcgaagcctcggcacgtttccttgcagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTCGCATTTACCCTCCGTGATGGTCGTCTGTTTACTCTgcgt                                                  <  1:238486‑M2/91‑1 (MQ=255)
                                                                                           gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                                 >  1:541048‑M2/52‑143 (MQ=255)
                                                                                           gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCAGTGATGGTCGTCTGTTTACTCTgcgtg                                                 >  1:49092‑M2/52‑143 (MQ=255)
                                                                                            aagacatcgaagcatcgtcacgtttctttgaagaccacgacggcctgcatCTCCTCCCTCTCCTGTGAAAATGCAAAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                                <  1:231550‑M2/93‑1 (MQ=255)
                                                                                            aagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                                <  1:304766‑M2/93‑1 (MQ=255)
                                                                                                   cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                         <  1:688862‑M2/100‑1 (MQ=255)
                                                                                                   cgaagcatcggcacgtttcgttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                         <  2:1015653‑M2/100‑1 (MQ=255)
                                                                                                         attggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGACCACGCCGGTAACTCCACTGCGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGccc                                   >  1:908397‑M2/38‑143 (MQ=255)
                                                                                                             gcacgtttcttggaagacgtggccggcctgcatCTCCTTCGTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtt                               <  2:76473‑M2/110‑1 (MQ=255)
                                                                                                                  tttctttgaagacgacggcggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGt                          <  2:299686‑M2/115‑1 (MQ=255)
                                                                                                                  tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCAAGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTCCTGTGGGTTGGGGTGGACTGCCCGCTTTTAGt                          >  1:420244‑M2/29‑143 (MQ=255)
                                                                                                                    tctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTGCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                        <  2:867417‑M2/117‑1 (MQ=255)
                                                                                                                    tctttgaagacgacgacggcctgcatCTCCTTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCGACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                        <  2:200070‑M2/117‑1 (MQ=255)
                                                                                                                     ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAGTCAACTGTGGCATTTAACATCCGTGATGGTCGTCGGTTTACTCTGCGTGAGCGTGGACTGCCGTCTTTTCGGCTg                       >  1:297591‑M2/26‑143 (MQ=255)
                                                                                                                        tgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTAt                    >  1:838003‑M2/23‑143 (MQ=255)
                                                                                                                         gcagacgacgacggcctgcatCTCCGTCTTCGTCTTTGAAGATGCGGAAGAGCACGCCGGTAACTCCACTGTGGCAGTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATc                   <  1:918813‑M2/122‑1 (MQ=255)
                                                                                                                           agacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTGTCGTCTGTATCGt                 >  1:108220‑M2/20‑143 (MQ=255)
                                                                                                                            gacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTa                >  2:695487‑M2/19‑143 (MQ=255)
                                                                                                                                  gacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTg          <  1:787177‑M2/131‑1 (MQ=255)
                                                                                                                                          gcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTCTCGGCTGTATCGTCTGCGTGCCCGTAGc  >  1:5794‑M2/5‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                          
GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGC  >  NC_000913/4001508‑4001788

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ 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.