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.040)49 (0.990) 43/284 0.1 96.1% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 2 (0.040)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGC  >  NC_000913/4001504‑4001780
                                                                                                                                           |                                                                                                                                         
tATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGCAGACATCGAAGCATCGGCACGTTTCTTTGAAGGCGACGACGGCCTGCATctc                                                                                                                                        >  2:730289‑M1/1‑140 (MQ=255)
         aTCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGGCGGCCTGCATctccttcttctt                                                                                                                               >  2:430834‑M1/1‑131 (MQ=255)
               tCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                         <  1:315552‑M1/143‑19 (MQ=255)
                 gAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCCCGTTTCTTTGAAGACGCCGACGGCCTGCATctccttcttcttctttgcag                                                                                                                       >  1:108728‑M1/1‑123 (MQ=255)
                 gAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCAGCGCGGAACTGGAAGACATCGAAGCATCGGCGCGGCTCTTGGATGACGACGACGGCCTGCCTctcctccgtctgctttgctg                                                                                                                       >  1:819511‑M1/1‑123 (MQ=255)
                  aaCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTCCTTTGACGACGACGACGGCCTGCGTccccttcttcttctttgacga                                                                                                                      >  2:209042‑M1/1‑122 (MQ=255)
                     cGGACGACGACGAGCGACTGCGCGTACAATCTGAGCTTGGCCAGAGCCTGGCAACCCGCGCGGAACTGGAAGACATCGAAGCATCGGCACGTTTGTTTGAAGACGACGGCGGCCTGCATctgcgtcttccacttttaaggtgc                                                                                                                   >  2:200075‑M1/1‑119 (MQ=255)
                     cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctcctttgaagatgc                                                                                                                   >  1:529247‑M1/1‑119 (MQ=255)
                     cGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACCACGGCCTGCATctccttcttcttctttgaagatgc                                                                                                                   >  1:71733‑M1/1‑119 (MQ=255)
                        acggcgacgAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgga                                                                                                                <  2:708647‑M1/143‑28 (MQ=255)
                          gacgacgaGCGACTGCGCGTACAATCTGACCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaag                                                                                                              <  1:259630‑M1/143‑30 (MQ=255)
                           acgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaaga                                                                                                             >  1:520158‑M1/1‑113 (MQ=255)
                             gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCGTctccttcttcgtctttgaagatgcgggaggtc                                                                                                           >  1:64339‑M1/1‑111 (MQ=255)
                             gacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagctgcggaagctc                                                                                                           >  1:670922‑M1/1‑111 (MQ=255)
                                 aGCGACTGCGCGTACAATCTGAACCTGGCCTGAGCCGCGCAACCCGCCCGCAACTGGAAGACATCGAAGCATCGACACTTTTCTTTGAAGACGACGACTTCCTGCAT‑‑‑‑‑CTCCTTCGTCTTCTTTGAAGATGCGGAAGATCACGc                                                                                                  <  1:651706/143‑1 (MQ=255)
                                  gCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgacgatgcggaagctcccccc                                                                                                      >  1:200020‑M1/1‑106 (MQ=255)
                                  gCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCAActccttcttcttctttgaagatgcggaagctcacgcc                                                                                                      >  1:652907‑M1/1‑106 (MQ=255)
                                   gggcggcgcgggcaACCTGAACTTGGCGAGGGCCTGGCAACCCGCCCGGAGCTGGAAGACATCGAAGCGTCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccg                                                                                                     <  1:689157‑M1/138‑39 (MQ=255)
                                      cTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagaggcggaagatcacgccggca                                                                                                  >  1:144615‑M1/1‑102 (MQ=255)
                                         cgcgTACCATCTGAGCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaact                                                                                               <  2:238726‑M1/143‑45 (MQ=255)
                                         cgcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaact                                                                                               >  1:739343‑M1/1‑99 (MQ=255)
                                             tACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGGCGGCCTGCATctccttcttcttctttgaagatgcggacgagcacgccggtaactccac                                                                                           >  1:153725‑M1/1‑95 (MQ=255)
                                              aCAATCGGGCCTTGGCCCGAGCCTGCCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGTCACGTTTGTTTGAAGACGACGACGGCCTGCATctccttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                          <  2:710667‑M1/143‑50 (MQ=255)
                                                  tCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccagtgtgg                                                                                      >  2:171613‑M1/1‑90 (MQ=255)
                                                    tGAACTTCGCCAGAGCCTTGCAACCCCCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggca                                                                                    <  2:714750‑M1/143‑56 (MQ=255)
                                                     gAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcaagccggtaactccactgtgtcat                                                                                   >  2:485783‑M1/1‑87 (MQ=255)
                                                       aCTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACAGACCGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcattt                                                                                 <  1:670049‑M1/143‑59 (MQ=255)
                                                              ctgagtcTGGCAATCCGCCCGGAACAGGAAGACATCGAAGCATCGGCACGTTTCGTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcc                                                                          <  1:594420‑M1/141‑66 (MQ=255)
                                                              cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgtggaagatcacgccggtaactccactgtggcatttaccatcc                                                                          >  1:359953‑M1/1‑78 (MQ=255)
                                                                   ccTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgat                                                                     >  1:839104‑M1/1‑73 (MQ=255)
                                                                     tGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaaggtgcggaagctcactccggtaactccactgcggcatttagcatccgagagca                                                                   >  2:91978‑M1/1‑71 (MQ=255)
                                                                          gcccgtccgGAGCCGTAAGCCCTCGATGCATCGGCACGTTTCTTGGAAGACGACGGCTGCCTGCAT‑‑‑‑‑CTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATcacggatggtaaat                                                         <  2:849486/142‑14 (MQ=255)
                                                                            ccgcccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctc                                                            >  1:451169‑M1/1‑64 (MQ=255)
                                                                                    aaCTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactct                                                    >  1:188558‑M1/1‑56 (MQ=255)
                                                                                       tGGAAGACATCGACGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcg                                                 <  2:261885‑M1/143‑91 (MQ=255)
                                                                                        ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatgggcggctgtttactctgcgt                                                >  2:623559‑M1/1‑52 (MQ=255)
                                                                                          aaGCCATGGAAGCATCGGCAGGTTTCTTTGAAGACGACGACGGCCTGCATttccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtga                                              <  2:561029‑M1/143‑94 (MQ=255)
                                                                                               aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttcctctgcgtgagcgtg                                         >  1:15888‑M1/1‑45 (MQ=255)
                                                                                                  gAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgcagatgcggaaggtcacggcggtaactccactgtgccatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaac                                      >  2:349366‑M1/1‑42 (MQ=255)
                                                                                                     gCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgggcgtgaactgc                                   >  1:858380‑M1/1‑39 (MQ=255)
                                                                                                            cACGTTCCTTTGCAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcgcaacatcacgccggtaacgccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttt                            <  1:196037‑M1/143‑112 (MQ=255)
                                                                                                               gTTTCTTTGAAGACGACGACGGCCTGCATctccttcttattctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcggcggtttactctgcgtgagcggggactgcccgcttttcg                         >  2:354626‑M1/1‑29 (MQ=255)
                                                                                                                tttCTTTGAAGACGACGACGTCCTGCATcttcttcttcttcttagaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgt                        <  1:42162‑M1/143‑116 (MQ=255)
                                                                                                                  tCTTTGGAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtct                      <  2:385417‑M1/143‑118 (MQ=255)
                                                                                                                     ttGAAGACGACGACGGCCTGCAGctccgtcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgta                   <  2:262196‑M1/143‑121 (MQ=255)
                                                                                                                       gAAGACGACGACTGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttatcatccgtaatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatc                 <  1:596961‑M1/143‑123 (MQ=255)
                                                                                                                                gacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgggcgtgacctagctgcgtttcgtctgttgaggggcccac        >  1:619694‑M1/1‑12 (MQ=255)
                                                                                                                                gacgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtgtcggatgcgtg        >  2:538927‑M1/1‑12 (MQ=255)
                                                                                                                                 acgGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgc       <  1:294107‑M1/143‑133 (MQ=255)
                                                                                                                                     ccTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccgctgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctcttcgtctgtatcgtatgcgtgcccgt   >  1:604739‑M1/1‑7 (MQ=255)
                                                                                                                                      cTGCATctccttcttcttcattgaagatgaggaagatcacgccggtaactccactgtggcagttaccatccgtgatgctagagtgtttactctgcgtgagcgtgaactgaccgcttgtcgtccgcatcgaatgcgtgagggtg  >  2:99570‑M1/1‑6 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGC  >  NC_000913/4001504‑4001780

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.