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

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

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.