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 = 40016432 (0.030)62 (0.800) 57/280 0.2 96.9% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 2 (0.030)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGC  >  NC_000913/4001510‑4001787
                                                                                                                                           |                                                                                                                                          
atggattgatcttgtcgaaccggacgacggcgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactgtaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                         <  2:1008735‑M2/4‑1 (MQ=255)
   gattgatcttgtggagccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacagcgaagcatcgcctcgtttctgttaagacgacgacggcctgcatCTCcttc                                                                                                                                      <  2:1287574‑M2/7‑1 (MQ=255)
        atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatgtgaacttggccagagcctggcaacccgcccggaactggaagacatggaagcatcggcacgttgcgttgaagaggacgacggcctgcagatccttcttggt                                                                                                                                 >  2:285643‑M2/132‑140 (MQ=255)
        atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttctt                                                                                                                                 <  2:859890‑M2/12‑1 (MQ=255)
         tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacggcggcctgcatCTacttcttcttt                                                                                                                                >  2:949684‑M2/131‑142 (MQ=255)
          ctggtcgaaccggacgacgacgagcgacggcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttct                                                                                                                               <  2:262505‑M2/14‑1 (MQ=255)
            tgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                             >  2:902791‑M2/128‑143 (MQ=255)
              tcgagcaggccgccgacggccgactgcgcgtacaatctgaccttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                           <  1:1010207‑M2/18‑1 (MQ=255)
               cgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                          <  2:1064605‑M2/19‑1 (MQ=255)
               cgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                          <  2:1047921‑M2/19‑1 (MQ=255)
                gaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggcactggaagacatcgaagcatcggcacgattctttgaagacggcgacggcctgcatCTCCTTCTTcttcctttacg                                                                                                                         >  1:717158‑M2/124‑139 (MQ=255)
                    cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTACTACTTTGCAGGTGc                                                                                                                     >  2:92295‑M2/120‑143 (MQ=255)
                        cgacgacgagcgcctgcgcgtacaatctgaacttgtctagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                 <  1:845990‑M2/28‑1 (MQ=255)
                         gacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggcactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctggatCTCCTTCTTCTTCTTTGCAGATGGCGAAg                                                                                                                >  1:805765‑M2/115‑143 (MQ=255)
                           cgacgagcgactgcgcgtacactctcaacctggccagagcctggcaacccgcccggaactggaagacatcggggcatcggcacttttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAt                                                                                                              <  1:816136‑M2/31‑1 (MQ=255)
                              cgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCAc                                                                                                           <  2:1150581‑M2/34‑1 (MQ=255)
                               gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggcactggaagacatcgaagcatcggcacgtttctttgatgacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGAGGCGGGAGAGCACg                                                                                                          >  2:687246‑M2/109‑143 (MQ=255)
                                      tgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                   >  2:1251336‑M2/102‑143 (MQ=255)
                                       gcgcgtaccatctgaacttggcaagagcctggcaacccgcccggaactggaagacatcgaagcatcggcccgtttctttgaacacgacgacgccctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGCTCACGCCGGTAAc                                                                                                  <  2:367360‑M2/43‑1 (MQ=255)
                                             acaatctgcacttggccagagccttgcgacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                            <  1:240176‑M2/49‑1 (MQ=255)
                                             acaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                            >  2:878908‑M2/95‑143 (MQ=255)
                                               aatctgaacttggccagagcccggcaacccgccgggaactgtaagacatcgaagcatcggcacgtttctttgaagacgacgacggccttcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgt                                                                                          <  2:87391‑M2/51‑1 (MQ=255)
                                                  ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgggtctttgaagacgacgacggcctgcagCTCCTCCTTCTTCTTTGAAGATGCGGAAGAGCACGCCGGTAACTCCACAGTGGc                                                                                       >  2:162833‑M2/90‑143 (MQ=255)
                                                       cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGCCATTTa                                                                                  >  2:132868‑M2/85‑143 (MQ=255)
                                                        ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAc                                                                                 >  2:589706‑M2/84‑143 (MQ=255)
                                                            ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATc                                                                             <  1:700173‑M2/64‑1 (MQ=255)
                                                             cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTGCCATcc                                                                            >  2:875365‑M2/79‑143 (MQ=255)
                                                               gagcctggcaccccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGt                                                                          <  2:234392‑M2/67‑1 (MQ=255)
                                                                ggcctggcaacccgcccggaagtggaggacatcgaggcatcggcacttttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                         <  1:1175554‑M2/68‑1 (MQ=255)
                                                                agccgggcaccccgcccggatctggaagacatcgaagcatcgtcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                         <  1:845082‑M2/68‑1 (MQ=255)
                                                                 gcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                        <  2:646561‑M2/69‑1 (MQ=255)
                                                                   ctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgacgacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                      >  2:240746‑M2/73‑143 (MQ=255)
                                                                     ctcgcccgGAGCGGGTGGTCATCGGGGCAGCGGCAAGTTTCTTTGAAGCCGTCGACGG‑‑TGCGCAT‑‑‑CTCAGTCTCCTCCTTTGATGATGCGGAAGCTCGCGCCGGTAACTCGACTGGGGCATTTACCATCCGTGATGGGCGTCt                                                               <  1:1291798/141‑1 (MQ=255)
                                                                       caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                        >  2:876968‑M2/69‑137 (MQ=255)
                                                                       caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                        <  1:876968‑M2/69‑1 (MQ=255)
                                                                        aacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcgt                                                                 <  2:1136290‑M2/76‑1 (MQ=255)
                                                                          tccgcccggagcgggaagacatcgaagcatcggcacctttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                               <  1:186817‑M2/78‑1 (MQ=255)
                                                                              cccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGttt                                                           <  2:940858‑M2/82‑1 (MQ=255)
                                                                                     ctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGCTGGTCGTCTGTTTACTCTgc                                                    >  2:1094749‑M2/55‑143 (MQ=255)
                                                                                     cggggggagatcgaggcatcggcacgtttctgtgaagacgacggccgcctgcatCCCCTCCGTCTGCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                    <  1:282973‑M2/89‑1 (MQ=255)
                                                                                         aagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                                <  2:598481‑M2/93‑1 (MQ=255)
                                                                                           gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                              <  2:1315621‑M2/95‑1 (MQ=255)
                                                                                           gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                              <  2:550776‑M2/95‑1 (MQ=255)
                                                                                             catcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTCCTTTGCAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgt                                            <  2:474597‑M2/97‑1 (MQ=255)
                                                                                              atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAATCCACTGTGGCATTTACCATCCGCGATGGCCGTCTGTTTACTCTgcgcgagcgtg                                           >  1:532248‑M2/46‑143 (MQ=255)
                                                                                               tcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcggga                                          >  2:432283‑M2/45‑143 (MQ=255)
                                                                                                cgaagcatcgggacgtttctttggagccggcgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTACCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                         <  1:495505‑M2/100‑1 (MQ=255)
                                                                                                 gaagcatcggcacgtttctttgaagacgacgacggcctgcctCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTactctgagtgagcggggac                                        >  1:991466‑M2/43‑138 (MQ=255)
                                                                                                    gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTGTTCTGCTTTGAAGATGCGGAAGAGCACGCCGGTAACCCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                                     <  2:580667‑M2/104‑1 (MQ=255)
                                                                                                     catcggcacgtttctgtgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACCCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                    <  2:1044827‑M2/105‑1 (MQ=255)
                                                                                                          gcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCCCTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtt                               >  1:425985‑M2/34‑143 (MQ=255)
                                                                                                          gcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGCCGGTTTTCTCTGCGTGAGCGTGAACTGCCCGCtt                               >  1:950545‑M2/34‑143 (MQ=255)
                                                                                                              tttGAAGACGACGACGGCCTGCAT‑‑‑‑‑CTCCTTCTTCACCTACAAAAATGCGGAAGATCACGCCGGTAACTCCACTGTAACACACACCAACCATAATCAACGTCTGTTTACTCTCCGTGACCATGAACTGCCCGCTTTTCGTCTGt                      >  1:896354/1‑143 (MQ=255)
                                                                                                               tttctttggagacgtcgtggccctgcacCTCCGTCTGCTTCTTGGGAGGTGCGGAAGATCACGCCGGTAACTCCGCTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGt                          <  1:1218263‑M2/115‑1 (MQ=255)
                                                                                                                 tctttgaagacgacgccggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                        <  1:949684‑M2/117‑1 (MQ=255)
                                                                                                                  ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                  >  1:437324‑M2/26‑116 (MQ=255)
                                                                                                                  ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                  <  2:437324‑M2/91‑1 (MQ=255)
                                                                                                                          acgacgacggcctgcatCTCCTTCTTCTTATTTGCAGATACGGAGGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAGCTGCCCACTTTTCGCCTGTCGCGtag               >  2:258515‑M2/18‑142 (MQ=255)
                                                                                                                            gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTCTGc             >  2:1354021‑M2/16‑143 (MQ=255)
                                                                                                                            gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGc             <  2:280845‑M2/128‑1 (MQ=255)
                                                                                                                               gacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGGACTGCCCGCTTTTCGTCTGTCTCGTTTGCGTg          >  2:644166‑M2/13‑143 (MQ=255)
                                                                                                                                  ggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTAAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGccc       <  1:240746‑M2/134‑1 (MQ=255)
                                                                                                                                    cctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAGCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGt     >  2:576099‑M2/8‑143 (MQ=255)
                                                                                                                                       gcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGc  >  2:186910‑M2/5‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGC  >  NC_000913/4001510‑4001787

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.