Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 4,001,645 Δ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 = 40016441 (0.020)46 (1.030) 41/284 0.1 97.9% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 1 (0.020)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCC  >  NC_000913/4001511‑4001789
                                                                                                                                           |                                                                                                                                           
atggattgatcttgtcgcaccggacgacgccgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                          <  2:658365‑M2/4‑1 (MQ=255)
   gcttgatcttgtcgaaccggaccacgacgagccacttcgcgtacaatctgaacttggccagcgcatggcaacccgcccggaactggaagacatcgacgcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                       <  2:647652‑M2/7‑1 (MQ=255)
       gatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttct                                                                                                                                   >  1:398885‑M2/133‑143 (MQ=255)
             gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTg                                                                                                                             <  1:357216‑M2/17‑1 (MQ=255)
                       tcgacgacgagctgctgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcgtcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGa                                                                                                                   <  2:454644‑M2/27‑1 (MQ=255)
                        cgacgacgagcgactgcgcgtacaatctgaacttggccagtgcctggcaacccgcccggaactcgaagacatcgaagcatcggcacgtttctttgacgacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                  <  1:35582‑M2/28‑1 (MQ=255)
                        cgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                  <  2:237818‑M2/28‑1 (MQ=255)
                        cgacgacgaccgactgcccgtaccatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                  <  1:387439‑M2/28‑1 (MQ=255)
                          acgacgagctactgcgcgtacgagctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                                <  1:24158‑M2/30‑1 (MQ=255)
                            gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcctcggcacgtgtctttgaagacgacgccggcctgcatCCCCTTCTTCCTCCTTGACGATGCGGAAGATc                                                                                                              >  2:430770‑M2/112‑143 (MQ=255)
                             acgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcccgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCa                                                                                                             <  2:582390‑M2/33‑1 (MQ=255)
                               gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggcagacatcgaagcaccggcacgtttctttgaagccgacgacggcctgcatCTCCTCCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                           >  2:558683‑M2/109‑143 (MQ=255)
                                agcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGc                                                                                                          <  1:751932‑M2/36‑1 (MQ=255)
                                 gcgcgtgcgcgtacgatctgaacgtggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCATGcc                                                                                                         <  2:733492‑M2/37‑1 (MQ=255)
                                 gcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGcc                                                                                                         >  1:131968‑M2/107‑143 (MQ=255)
                                  cgactgcgcgtacaagctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                        <  1:554471‑M2/38‑1 (MQ=255)
                                      tgcgcgtacaatctgagcttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTaa                                                                                                    <  2:736215‑M2/42‑1 (MQ=255)
                                          caTCTGAACTCGGCCAGAGCCTGGCAACCCGCCCGGAACGGGGAGACATCGAAGCCTCGGTCCGTTTCATTGTATTG‑‑GATGGCATCATTTGTCTCCTTCTTCTTCTTTGAAGACGCGGTATATCACGCCGGTAACTCCACtgt                                                                                           <  2:784378/142‑1 (MQ=255)
                                             tggatctgatcttcgacagggcctggcaacgcgcgcggaactggaagacatggaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGTGGAAGATCACGCCGGTAACTCCACt                                                                                             <  2:4913‑M2/49‑1 (MQ=255)
                                                 tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCGACTGTgg                                                                                         >  1:487185‑M2/91‑143 (MQ=255)
                                                  cgaaactcggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCGTCTGTGAAAATCATCAACATCACGCCGGTAACTCCACTGTGGc                                                                                        <  2:99277‑M2/54‑1 (MQ=255)
                                                           gccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGAGGCGGAAGCTCACGCCGGCAACTCCACTGTGGCAGTTACCAt                                                                               >  2:614384‑M2/81‑143 (MQ=255)
                                                               gagccgtgctacccgccgggaactggaatgcgtctgagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGt                                                                           <  2:127738‑M2/67‑1 (MQ=255)
                                                                 gcctggcaagccggtctgaattggaagacatcgaagcatgggcacgtttctttggagacgacgacggcctgcatCTCCGTCTTCTTCTTTGAAGATGCGGAAGATGACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                         <  2:81320‑M2/69‑1 (MQ=255)
                                                                   gtggcacgccgcccggagctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTGCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACGGTGGCATTTACCATCCGTGATg                                                                       <  2:359291‑M2/71‑1 (MQ=255)
                                                                   ctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                       >  1:680134‑M2/73‑143 (MQ=255)
                                                                   ctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGACGCGGAATATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                       <  1:765287‑M2/71‑1 (MQ=255)
                                                                      gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtc                                                                    <  1:598323‑M2/74‑1 (MQ=255)
                                                                            cgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                              <  1:337657‑M2/80‑1 (MQ=255)
                                                                                   aactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTActct                                                       <  1:17456‑M2/87‑1 (MQ=255)
                                                                                      tggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcg                                                    <  2:435981‑M2/90‑1 (MQ=255)
                                                                                       ggaggagatcgaagcgtcggtatgcttgttcgaagacgacgagggcctgcacCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCATGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                   <  2:343187‑M2/91‑1 (MQ=255)
                                                                                       ggaagacatcgaagcatctgcacgtttctttgaagacgacgccggcctgcatCTCCTTCTGCTTCTTGGAAGATGCGGAATATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                   <  2:550685‑M2/91‑1 (MQ=255)
                                                                                       ggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTCACCATCCGTGATGGTCGTGGATTTACTCTgcgt                                                   >  2:378686‑M2/53‑143 (MQ=255)
                                                                                        gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                                  >  2:292404‑M2/52‑143 (MQ=255)
                                                                                        gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                                  >  1:103849‑M2/52‑143 (MQ=255)
                                                                                         cagacactgaagcatcggcacttttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                                 <  2:752020‑M2/93‑1 (MQ=255)
                                                                                         aagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                                 >  1:138083‑M2/51‑143 (MQ=255)
                                                                                                 gaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTCACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAAc                                         >  2:582987‑M2/43‑143 (MQ=255)
                                                                                                    gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCCTGCGCGTGAGCTGc                                      >  2:415033‑M2/40‑143 (MQ=255)
                                                                                                        cggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGATGTTCCCGCCCGTCCCTCCACTGTGGCATTTCCCATTCTTGATGGTAGACTGTTTACTCTGCGTGAGTGTGAACTGACcga                                  >  1:64814‑M2/36‑142 (MQ=255)
                                                                                                            acgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtttt                              <  1:262808‑M2/112‑1 (MQ=255)
                                                                                                               tttctttgaagacgacgccgccctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGt                           <  1:154534‑M2/115‑1 (MQ=255)
                                                                                                               tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTGACTCCCCTGTGGCATTTCCCCTCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGACCTGCCCGCTTTTCGt                           >  1:203558‑M2/29‑143 (MQ=255)
                                                                                                                   tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCCTTTACCATCCGTGATGGTCGTCTGTCTACTCTGCCCGAGCGTGACCTGCCCGCTTTTCGTCCGt                       >  1:742275‑M2/25‑143 (MQ=255)
                                                                                                                    ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGCCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCCTCTgtg                      >  2:653858‑M2/24‑142 (MQ=255)
                                                                                                                                        catCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGAGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGCAGcc  >  2:494159‑M2/4‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCC  >  NC_000913/4001511‑4001789

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.