Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,001,644 Δ5 bp coding (220‑224/951 nt) corA → magnesium/nickel/cobalt transporter

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4001644 4001648 5 52 [0] [0] 53 corA magnesium/nickel/cobalt transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40016430 (0.000)52 (0.950) 43/284 0.2 100% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 0 (0.000)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001507‑4001786
                                                                                                                                              |                                                                                                                                         
agtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatc                                                                                                                                           <  2:577561‑M2/1‑1 (MQ=255)
 gtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatct                                                                                                                                          >  2:421556‑M2/142‑143 (MQ=255)
   atggtttgatcttgtcgaaccggacgacggcgatcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                        <  1:439438‑M2/4‑1 (MQ=255)
   atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                        <  2:362567‑M2/4‑1 (MQ=255)
   atggattgagcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                        <  1:421833‑M2/4‑1 (MQ=255)
        ttgatcttgtcgaaccggacgacgacgagcgactgcgcgtaccatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcgtcggcacgtctctttgacgacgacgccgacctgcaaCCccttctc                                                                                                                                   >  1:659797‑M2/135‑142 (MQ=255)
        ttgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggcctgcacCTCcttctt                                                                                                                                   >  1:696700‑M2/135‑143 (MQ=255)
        ctgatctggtcgaaccggaggacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttctt                                                                                                                                   <  1:866897‑M2/9‑1 (MQ=255)
            tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcctCCCcctcttcctc                                                                                                                               >  1:847666‑M2/131‑143 (MQ=255)
                 ccgacccggacgacgacgagcgagtgcgggtacgatcggaacttggccagagcctggcaacccgcccggaactggaagacctcgaagcatcggcacgtttctttgaagacgacgccggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                          <  2:283181‑M2/18‑1 (MQ=255)
                  cgaaccggacgacgacgagcgacggcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                         <  2:785246‑M2/19‑1 (MQ=255)
                          acgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaccccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGa                                                                                                                 <  1:608748‑M2/27‑1 (MQ=255)
                          acgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGa                                                                                                                 <  1:281475‑M2/27‑1 (MQ=255)
                           cgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                <  2:545737‑M2/28‑1 (MQ=255)
                             gctgcgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgtagcatcggcacgtttctttgaagtcgacgacggcttgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                              <  1:508820‑M2/30‑1 (MQ=255)
                             acgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                              <  2:419460‑M2/30‑1 (MQ=255)
                                tgcagcggctgcgcgtgcaagctgaagttggccagagtctggcaacccgcccggaactggaagacatcgcagcagcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCa                                                                                                           <  1:466247‑M2/33‑1 (MQ=255)
                                  gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGACGATGCGGTCGCTCACg                                                                                                         >  2:43957‑M2/109‑143 (MQ=255)
                                         tgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                  <  2:120544‑M2/42‑1 (MQ=255)
                                               tacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCAc                                                                                            <  2:44136‑M2/48‑1 (MQ=255)
                                               tacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtctctttgaagacgacgccggcgtgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCGCGTAACTCCAc                                                                                            >  1:878360‑M2/96‑143 (MQ=255)
                                                  aatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGGAACTCCACtgt                                                                                         >  1:791071‑M2/93‑143 (MQ=255)
                                                      tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCAGTGTGGCa                                                                                     >  2:824884‑M2/89‑143 (MQ=255)
                                                      tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACAGAGGCa                                                                                     >  2:816106‑M2/89‑143 (MQ=255)
                                                      tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTGCTTTGAAGATGCGGAGGATCACGCCGGTAACTCCACTGTGGCa                                                                                     >  2:592629‑M2/89‑143 (MQ=255)
                                                               ccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATc                                                                            >  2:668998‑M2/80‑143 (MQ=255)
                                                                  gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGt                                                                         <  2:930821‑M2/67‑1 (MQ=255)
                                                                    gcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                       <  2:223358‑M2/69‑1 (MQ=255)
                                                                     cctggcaccccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGGGAt                                                                      >  1:388333‑M2/74‑143 (MQ=255)
                                                                      ttggcaacgcgcccgtaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                     <  2:475142‑M2/71‑1 (MQ=255)
                                                                       tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTGCTTGTTCTTTGAAGCTGCGGAAGATCACGCCGGTAACTCCACtacgtcatttaccgcgcgagagag                                                                    >  1:60179‑M2/72‑132 (MQ=255)
                                                                              ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCCTTTACCCTCCGTGATGGTCGTCTg                                                             >  2:769541‑M2/65‑143 (MQ=255)
                                                                                gcccggaactggaagacatcgaagcatcggcgcgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                           <  1:729745‑M2/81‑1 (MQ=255)
                                                                                  ccggaactggaagacatcgaagcatcggcacgaaaaacttaagtctacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                         <  2:396896‑M2/83‑1 (MQ=255)
                                                                                    ggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAct                                                       <  2:918950‑M2/85‑1 (MQ=255)
                                                                                      aactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActct                                                     <  1:421556‑M2/87‑1 (MQ=255)
                                                                                        ctggaagacatcgaagcagcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                   <  2:866025‑M2/89‑1 (MQ=255)
                                                                                           gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTGCTCTgcgtg                                                >  2:13635‑M2/52‑143 (MQ=255)
                                                                                            acgacatcgaatcctcggcacgtttctttgcagacgaccacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                               <  1:473989‑M2/93‑1 (MQ=255)
                                                                                              gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                             >  2:6144‑M2/49‑143 (MQ=255)
                                                                                                catcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgt                                           <  1:611346‑M2/97‑1 (MQ=255)
                                                                                                 atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgcgagcgtg                                          >  1:71622‑M2/46‑143 (MQ=255)
                                                                                                   cgacgcatcggcacgtttctttgaagacgacgacggcctgcatCTCCGTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                        <  1:358785‑M2/100‑1 (MQ=255)
                                                                                                          tcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGGGTgagagtgaagtcccac                                 >  2:347980‑M2/37‑136 (MQ=255)
                                                                                                           cggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATTTTCGGCGGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                <  2:907784‑M2/108‑1 (MQ=255)
                                                                                                            ggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGCGATGGTCGTCTGTTTGCTCTGCGTGCGCGTGCACTGTCCGCt                               >  2:326563‑M2/35‑143 (MQ=255)
                                                                                                                         gaagaggacgacggcctgcatCTCCTTCTTCTTCTTTTAAGATGCGGAAGATCACGCCGGTAACTCCTCTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATc                  <  1:392207‑M2/122‑1 (MQ=255)
                                                                                                                          atgacgacggcggcgtgctgCTCGTTCTTCTTCTTTGAGGAGGCGGATGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAAGTGCCCGCTTTTCGTCTGTATCg                 <  1:916955‑M2/123‑1 (MQ=255)
                                                                                                                            gacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATAACGCCGGTAGATCCACTGTGGCATTTACAATCAGTGATGGTCGTCAGTTTGCTCTGCGGGAGCGTGAACTGCCTGCTTTTGGTCTGTAGCGTa               >  2:193577‑M2/19‑143 (MQ=255)
                                                                                                                                 cgacggcctgcatCTCCTTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATGTACCATCCGTGTTGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGt          <  2:378340‑M2/130‑1 (MQ=255)
                                                                                                                                     ggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACCGTGGCATTTACCATCCGTGATGGTCGTCTGTTGCCTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTGTCGTACGCGTGgcc      >  2:63511‑M2/10‑143 (MQ=255)
                                                                                                                                        cttcaccTCCCTCTTCATCTTTGGACAGGCGGAATATCACGCCTGTACCTCCACGGTGGCATTTACCATCCGTGATTGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTa   <  1:824884/137‑1 (MQ=255)
                                                                                                                                         tgcatTTTCTTCTTCTTCTTTGAAGATGCGGAAGATCTCGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCTTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAg  <  1:389846‑M2/138‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                         
GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAG  >  NC_000913/4001507‑4001786

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.