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)35 (0.840) 32/284 0.3 97.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

CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCG  >  NC_000913/4001501‑4001783
                                                                                                                                              |                                                                                                                                            
ctggctGGATTTATCTTTTCGATCCGGACGACGAGGGGCGATTGCGCGTACAATCTGAACTTGGCTAGAGCCTGGCAAGCCTCTCGTAACTGGCAGACATCGAAGCATCGGCACGTTTCTTTGAAGACTACAACGGCCTGCat                                                                                                                                              <  2:151686/138‑1 (MQ=255)
    gcagtttgctctCGTCGAAGCTGACAACGCCGAGCGACTGCTCGTACAATCTGAACTTGGCCAGAGCCTGGCACCCCGCCCGGAACTGGAATACATCGAAGCATCGGTACGTTTCTTTGAAGACGACGACGGCCTGCATctcc                                                                                                                                          <  2:212600‑M1/138‑5 (MQ=255)
        attactcttgtCTACCCGGACGACGACGAGCGACTGCGCGTCCAACCCGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttct                                                                                                                                      <  1:345550‑M1/138‑9 (MQ=255)
             tCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttc                                                                                                                                 >  1:424931‑M1/1‑130 (MQ=255)
                  ccgAACCGGACGACGGCGAGCGACTGCGCGTACAATCTGAACCTGGCCAGAGCCTGGCAGCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttga                                                                                                                            <  1:464990‑M1/142‑19 (MQ=255)
                   cGACCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                           <  1:81461‑M1/143‑20 (MQ=255)
                   cGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaa                                                                                                                           <  1:581910‑M1/143‑20 (MQ=255)
                     aaCCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaaga                                                                                                                         <  1:238631‑M1/143‑22 (MQ=255)
                            cgacgacgaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaa                                                                                                                  <  1:352261‑M1/143‑29 (MQ=255)
                                   gaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctcctttgaagatgcggaagatcacg                                                                                                           >  2:466608‑M1/1‑108 (MQ=255)
                                   gaGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttctccttctctgaacatgcggaacgtcacg                                                                                                           >  2:584716‑M1/1‑108 (MQ=255)
                                       gACTGCGCGGACAAGCTGAAGTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccgg                                                                                                       <  1:316637‑M1/143‑40 (MQ=255)
                                             gcgTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgacgctgcggacggtcacgccggtagctg                                                                                                 >  1:296069‑M1/1‑98 (MQ=255)
                                                    aTCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtg                                                                                          <  1:237793‑M1/143‑53 (MQ=255)
                                                      cTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgtggaagatcacgccggtaactccactgtggc                                                                                        >  1:395295‑M1/1‑89 (MQ=255)
                                                       tGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGCAGCATCGGCACGTTTCTTGGAAGACGACGACGGCCTGCATctccttcttcttctttgaagaggcggaagatcacgccggtaaagccactgagcgc                                                                                       >  1:670761‑M1/1‑88 (MQ=255)
                                                                cccgagcctGGCAACCCGCCCGGAACTGGACGCCATCGACGCATCGGCACGTTTCTTTGAAGACGACGACCGCCTGCATctccttcttcttctttgcagatgcggaagatcacgccggtaactccactgtggcatttaccatc                                                                              <  1:176229‑M1/140‑65 (MQ=255)
                                                                 cAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatcc                                                                             >  2:517973‑M1/1‑78 (MQ=255)
                                                                   gagCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgt                                                                           <  1:176889‑M1/143‑68 (MQ=255)
                                                                            aaCCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGGTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttcttggaagatgcggaagatcacgccggtaactccactctggcatttaccatccgtgatggtcgt                                                                  <  2:694786‑M1/143‑77 (MQ=255)
                                                                              cccgcccgGAGCGGGAAGACATCGAAGCATCGGCACGTTTCTTTTAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtct                                                                <  2:402298‑M1/143‑79 (MQ=255)
                                                                                   ccgGAACTGGAAGACATCGAAGCATCGGCACGTTTCTCTACCGACGACGACCGCCTG‑‑‑A‑‑CTCTCCTTCTTCTTCGTTGAAGATGCGGCCCATCCCGCCGGCCACTACACGGCGGCATTCGCCACCCGCGACGgccgccacacta                                                      >  2:614812/1‑135 (MQ=255)
                                                                                      gAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgcgatggtcgtctgtttactc                                                        >  2:531430‑M1/1‑57 (MQ=255)
                                                                                       aaCTGGAAGACAGCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactct                                                       <  1:111158‑M1/143‑88 (MQ=255)
                                                                                           ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgt                                                   <  1:13540‑M1/143‑92 (MQ=255)
                                                                                           ggAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcgtctttgaagatgcggaagatcacgccggtacctccgctgtggcatttgccatccgtgatggtcggctgtttactctgcgt                                                   >  1:549282‑M1/1‑52 (MQ=255)
                                                                                               gACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagc                                               >  1:252733‑M1/1‑48 (MQ=255)
                                                                                                  aTCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatacgtgatggtcgtctgtttactctgcgtgagcgtgaac                                            >  1:177761‑M1/1‑45 (MQ=255)
                                                                                                              gCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgtggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                <  2:395295‑M1/143‑111 (MQ=255)
                                                                                                              gCACGTTTCTTTGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgctt                                <  2:675406‑M1/143‑111 (MQ=255)
                                                                                                                      cTTTGAAGACGACGCCGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtacctccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctg                        <  1:315235‑M1/143‑119 (MQ=255)
                                                                                                                       tttGAAGACGACGACGGCCTGCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttctcggctgt                       >  1:177128‑M1/1‑24 (MQ=255)
                                                                                                                              acgacgtcgGCCTGCATctcctttttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctatttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtat                <  1:696586‑M1/143‑127 (MQ=255)
                                                                                                                                tgcgtcggcctgCATctccttcttctgccttgacgttgcggcagctctcgccggtacctcctctgtggcacttaccctccgtgatggtcgtctgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgc              <  2:722643‑M1/138‑129 (MQ=255)
                                                                                                                                    gcggCATGCATctccttcttcttctctgaagatgcggaagatcacgccggtaactcccctgtggcatttaccatccgtgatggtcgtttgtttactctgcgtgagcgtgaactgcccgcttttcgtctgtatcgtatgcgtgc          <  1:145406‑M1/142‑133 (MQ=255)
                                                                                                                                           gCATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgcgcgtgaactgcccgcttttcgtctgtatcgtatgcgcgcccgtagc   >  2:615002‑M1/1‑4 (MQ=255)
                                                                                                                                            cATctccttcttcttctttgaagatgcggaagatcacgccggtaactccactgtggcatttaccatccgtgatggtcgtctgtttactctgcgtgagcgtgaaatattcgctgcccgtctgtatcgtatgcgtgcccgttgcc  >  1:302224‑M1/1‑3 (MQ=255)
                                                                                                                                              |                                                                                                                                            
CAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCG  >  NC_000913/4001501‑4001783

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.