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)55 (1.000) 46/284 0.1 98.2% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 1 (0.020)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGC  >  NC_000913/4001507‑4001780
                                                                                                                                              |                                                                                                                                   
agtgtggattgagcgtgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatc                                                                                                                                     <  2:620684‑M2/1‑1 (MQ=255)
 gtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatct                                                                                                                                    >  1:772645‑M2/142‑143 (MQ=255)
   atggagtggtcttggcgaacgggacgacgacgagcgactgcgcgtacaatctgaactgggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                  <  2:447447‑M2/4‑1 (MQ=255)
     ggtctgatcgtgtcgcatcggacggcgacgcgcgactgcgcgtacaatcggaacttggccggagcctggcaacccgcccggcactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCctt                                                                                                                                <  1:129849‑M2/6‑1 (MQ=255)
           atcttgtcgaaccggacgacgacgagcggctgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttctt                                                                                                                          >  2:408648‑M2/132‑143 (MQ=255)
              ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgcagcatcggcacgtttctttgaagacgacgccggcctgcatCCCcttcttcttctt                                                                                                                       >  2:792063‑M2/129‑143 (MQ=255)
                           cgacgacgggcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                          <  1:400402‑M2/28‑1 (MQ=255)
                               gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctctgaagacgacgacggcctgcatCTCCTTCTTCTCCTTTGAAGATGCGGAagacc                                                                                                      >  1:153976‑M2/112‑141 (MQ=255)
                                    gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGcc                                                                                                 <  1:174188‑M2/37‑1 (MQ=255)
                                       actgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagtatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGt                                                                                              <  2:72341‑M2/40‑1 (MQ=255)
                                        ctgcgcgtacaatctgaacttggccagagcctggcaaccagcacggaactggaagacatcgaaacatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGCTGCGGAAGATCTCTCCGgtc                                                                                             >  2:205523‑M2/103‑142 (MQ=255)
                                         ttcgcgtacaatctgaactccgccagcgcctggcaacccgcccggaactggaagccctcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                            <  2:108681‑M2/42‑1 (MQ=255)
                                           cgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTCACt                                                                                          >  1:146585‑M2/100‑143 (MQ=255)
                                                   atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                  <  1:937897‑M2/52‑1 (MQ=255)
                                                   atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                  <  1:743883‑M2/52‑1 (MQ=255)
                                                      tgaagtgggccagagcctggcaacccgcccggcactggaggtcgtggtcgcagcggcacgtttctttgaactcgacgacggcctgcatCTCCTTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCa                                                                               <  2:904600‑M2/55‑1 (MQ=255)
                                                         acttggccagagcctggcaacccgcccggaactggaagagaccgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTACCTCCACTGTGGCAttt                                                                            <  1:936746‑M2/58‑1 (MQ=255)
                                                         acttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCAttt                                                                            <  2:47634‑M2/58‑1 (MQ=255)
                                                          cttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                           <  2:729728‑M2/59‑1 (MQ=255)
                                                               ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATc                                                                      >  2:863203‑M2/80‑143 (MQ=255)
                                                                    gcctggcagctcgcccggcactggaagacatcgaagcatcgtcacgtttctttgaagacgacgacggcctgcatGTCCTGCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACCCCACTGTGGCATTTACCATCCGTGa                                                                 <  1:242026‑M2/69‑1 (MQ=255)
                                                                     cctggcaacccgcccggaactggaagacatcgaagcagcggcacggttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGAt                                                                <  1:401963‑M2/70‑1 (MQ=255)
                                                                      ctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCCTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTATGGCATTTACCATCCGtgacg                                                               >  2:760183‑M2/73‑141 (MQ=255)
                                                                        ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgt                                                             >  1:791833‑M2/71‑143 (MQ=255)
                                                                         gcaaccctcccggaactggaagacatcgaagcgtcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtc                                                            <  1:301731‑M2/74‑1 (MQ=255)
                                                                             tcgGAACTTGAATACATCGATGCATAGGCACGTTTCTTTGAAGACGATGACTGCTTGCAT‑‑‑‑‑CTCCTTGTTATTTTTTTAAGATGCGGAAGATCACGCCGGTAACTCTACTTTGGCATTTACTATCCGTGATGGTCGTCTGTTTa                                                   <  1:463122/142‑1 (MQ=255)
                                                                              ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATTACAATGTGATCTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                       >  2:944199‑M2/65‑143 (MQ=255)
                                                                              ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                       >  1:652391‑M2/65‑143 (MQ=255)
                                                                              ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                       <  2:208378‑M2/79‑1 (MQ=255)
                                                                               cgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                      <  1:841058‑M2/80‑1 (MQ=255)
                                                                                gcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                     <  1:449305‑M2/81‑1 (MQ=255)
                                                                                 cccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTgttc                                                    >  1:558803‑M2/62‑142 (MQ=255)
                                                                                    ggaacgtgaagacatcgaagcatcggcgcgtttctttgaagacggcgacgggctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCCCGCCGGTAACTCCACTGTGGCATTTACCATCCGTGGTGGTCGTCTTTTTAct                                                 <  1:388988‑M2/85‑1 (MQ=255)
                                                                                       actggaagacatcgaagcatcggcacgattctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTg                                              >  1:527295‑M2/56‑143 (MQ=255)
                                                                                        ctggaagacagcgaaaaaccggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                             <  2:253003‑M2/89‑1 (MQ=255)
                                                                                           gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                          >  2:453153‑M2/52‑143 (MQ=255)
                                                                                             agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                        >  2:356596‑M2/50‑143 (MQ=255)
                                                                                              gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAGCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                       <  2:203611‑M2/95‑1 (MQ=255)
                                                                                                 atcgaagcatcggcacgtttctttgaagacgacgacgtcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtg                                    <  2:839628‑M2/98‑1 (MQ=255)
                                                                                                  tcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTCACTCCACTGTGGCATTTACCATCCGTGATGGTCGCCTGTTTACTCTgcgtgagcgtga                                   >  1:696004‑M2/45‑143 (MQ=255)
                                                                                                   cgacgcaccggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTTTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                  <  2:229819‑M2/100‑1 (MQ=255)
                                                                                                   cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                  <  2:746941‑M2/100‑1 (MQ=255)
                                                                                                   cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                  <  2:772645‑M2/100‑1 (MQ=255)
                                                                                                    gaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAAc                                 <  2:676378‑M2/101‑1 (MQ=255)
                                                                                                             gcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGGCGTCTGTTTACTCTGCGTGAGCGTGAACTACCCGCtt                        >  1:738332‑M2/34‑143 (MQ=255)
                                                                                                                cgttacgatgaagacgacgacggcctgcatCTCCTTCGTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCCCTGTGGCATTTACTATCCGTGATGGTCGTCTGGTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTc                     <  2:136292‑M2/113‑1 (MQ=255)
                                                                                                                 gtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACGGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGGGCGTGGACTGCCCGCTTTTCg                    >  2:938486‑M2/30‑143 (MQ=255)
                                                                                                                 gtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAGCACGCCGGTAACTCCCCTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCg                    >  2:803242‑M2/30‑143 (MQ=255)
                                                                                                                   ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTc                  >  2:49394‑M2/28‑143 (MQ=255)
                                                                                                                   ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTc                  >  2:47221‑M2/28‑143 (MQ=255)
                                                                                                                    tctttgaaggcgacgacggcctgcatCTCCTTCTGCTTCTTTGATGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                 <  1:234260‑M2/117‑1 (MQ=255)
                                                                                                                         gaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATc            >  2:786253‑M2/22‑143 (MQ=255)
                                                                                                                               gacgactgcttgcatCGCCTTCTTCTTCGGTGGGGATGCGGAAGATCAGGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCGGCGTGAGCGTGAACTGTCCGCTTTTCGTCTGTATCGTATGc      <  1:377453‑M2/128‑1 (MQ=255)
                                                                                                                               gacgacggcctgcgaCTCCTGCTTCTCCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGc      <  2:146585‑M2/128‑1 (MQ=255)
                                                                                                                               gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACTATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGc      <  2:746965‑M2/128‑1 (MQ=255)
                                                                                                                                   acggcccgcagCTCCCTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGCCATTCAACCTCCGTTAGGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGc  <  2:854908‑M2/132‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                   
GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGC  >  NC_000913/4001507‑4001780

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.