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 48 [0] [0] 48 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)48 (0.770) 42/282 0.5 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

ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCA  >  NC_000913/4001509‑4001789
                                                                                                                                            |                                                                                                                                            
tatggattgattctgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatctc                                                                                                                                            <  2:960867‑M2/3‑1 (MQ=255)
    gattgatcttctcgatccggacgacgacgatcgactgcgcgtacaatctgaacttggccagcgcctggcaatccgcccgtaacttgatgacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                        <  2:111211‑M2/7‑1 (MQ=255)
         atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttctt                                                                                                                                   >  1:979153‑M2/132‑143 (MQ=255)
          tcttgtcgaaccggacgacgacgagcgactgcgcgcacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                  <  2:425154‑M2/13‑1 (MQ=255)
            ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtatctttgaagacgacgacagcctgcatCTCcttcttcatctt                                                                                                                                >  1:738501‑M2/129‑143 (MQ=255)
                     cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                       <  2:670979‑M2/24‑1 (MQ=255)
                      ggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCg                                                                                                                      >  1:59087‑M2/119‑143 (MQ=255)
                           acgacgagcgacggcgcgtacaatgtgaacttggacagagccgggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctggatCTCCTTCTTCTTCTTTGAAGATGCGGAAGa                                                                                                                 <  2:573559‑M2/30‑1 (MQ=255)
                                   cgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                         <  1:1012742‑M2/38‑1 (MQ=255)
                                   cgactgcgcgtacaatctcaccttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                         <  1:468226‑M2/38‑1 (MQ=255)
                                   cgactgcgcgtacaagctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                         <  1:935904‑M2/38‑1 (MQ=255)
                                    gactgcgcggacaatcggaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCgg                                                                                                        <  2:887999‑M2/39‑1 (MQ=255)
                                            gtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTGCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCa                                                                                                >  1:955407‑M2/97‑143 (MQ=255)
                                              acaatctgaacttggctagtgtctttcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                              <  1:16318‑M2/49‑1 (MQ=255)
                                                aatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGCGGACGATCACGCCGGTCACTCCACtgt                                                                                            >  2:794789‑M2/93‑143 (MQ=255)
                                                 atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                           >  1:591138‑M2/92‑143 (MQ=255)
                                                   ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATGACGCGGGTAACGAGACTGTGGc                                                                                         >  1:710642‑M2/90‑143 (MQ=255)
                                                         ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGGAGATCACGCCGGTCACTCCACTGTGGCATTtaa                                                                                   >  1:106110‑M2/84‑142 (MQ=255)
                                                         ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCGCTGTGGcattttc                                                                                   >  1:988347‑M2/84‑141 (MQ=255)
                                                         gtggccagagccgggcaacccgcccggaactggaagacatcgaagcatcggcacgtctctttgaagacgacgacggcccgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGTATTTAc                                                                                   <  1:432955‑M2/60‑1 (MQ=255)
                                                              cagagccttgcaatcctctcgtatcttgaatatatcgaagcatcggcacggttcttggaagacgacgacggcctgcatCTCCTTCTGCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATcc                                                                              <  1:742498‑M2/65‑1 (MQ=255)
                                                                 agcctggcggtccgcccggacctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGTGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                           <  2:874535‑M2/68‑1 (MQ=255)
                                                                  gcccggcaacccgcccggaactggaagacagcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGa                                                                          <  1:175253‑M2/69‑1 (MQ=255)
                                                                    ctcgtggcccgccctgaactggaagacgtcgaagcatcggcacgtttctttgaagacgacgactgcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                        <  2:524171‑M2/71‑1 (MQ=255)
                                                                     tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                               >  1:182172‑M2/72‑103 (MQ=255)
                                                                     tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                               <  2:182172‑M2/32‑1 (MQ=255)
                                                                           cccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                 <  1:343205‑M2/78‑1 (MQ=255)
                                                                            ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                                >  1:240803‑M2/65‑143 (MQ=255)
                                                                            ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTg                                                                >  2:33511‑M2/65‑143 (MQ=255)
                                                                                ccggaactggaagacatcgacgcatcggcacgtttgtttgaagacgacgacggccggcatCGCCTTCTTCTTCTTTGAAGAGCCTGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                            <  1:20495‑M2/83‑1 (MQ=255)
                                                                                    aactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTGCCATCCGCGATGGTCGCCTGTTTActct                                                        >  2:1071571‑M2/57‑143 (MQ=255)
                                                                                      ctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgc                                                      <  1:751399‑M2/89‑1 (MQ=255)
                                                                                       tggcagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTCACTCTgcg                                                     >  1:87847‑M2/54‑143 (MQ=255)
                                                                                       tggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGGTGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGCTGGTCGTCTGTTTACTCTgcg                                                     >  2:855074‑M2/54‑143 (MQ=255)
                                                                                         gaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtg                                                   <  2:979153‑M2/92‑1 (MQ=255)
                                                                                           agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                                 >  2:729724‑M2/50‑143 (MQ=255)
                                                                                                     gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGc                                       <  2:41148‑M2/104‑1 (MQ=255)
                                                                                                         cggcacgtttctttgacggcgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGACCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTCCTCTGCGTGAGCGTGAACTGCCcgg                                   >  2:957202‑M2/36‑142 (MQ=255)
                                                                                                         cggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGCAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCGTCCGTGCTGGTCGTCTGTTTACTCTGCGTGAGCGTGAACtgcccca                                   >  2:792552‑M2/36‑141 (MQ=255)
                                                                                                           gcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGGGCGTGAACTGCCGGctg                                 >  1:55950‑M2/34‑142 (MQ=255)
                                                                                                                 ttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAGGATGCGGAAGGACACGCCGGTAACTCCACTGTGGCATTTACCATCCTTGTTGGTCGGCTGTTTACTCTGAGTGAGCGTGAACGGCCCGCTTTtcgcg                           >  2:176971‑M2/28‑141 (MQ=255)
                                                                                                                    tttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGGCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTCCATCTGt                        >  1:88608‑M2/25‑143 (MQ=255)
                                                                                                                        aagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGCTGGTCGTCTGTTTACTCTCCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCg                    >  2:735850‑M2/21‑143 (MQ=255)
                                                                                                                              acgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGACCGTCAAAAGTCCGCTTTTCGTCTGTATCGTATGCg              <  2:955407‑M2/129‑1 (MQ=255)
                                                                                                                              acgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCGCTTTACTCTGCGTGAGCGTAACCTGCCCCCTGCTCGTCTGTCTCGTCTGCg              >  2:35236‑M2/15‑143 (MQ=255)
                                                                                                                               cgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACGATCCGTGATGGGCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTGTCGTATGCGt             >  1:133907‑M2/14‑143 (MQ=255)
                                                                                                                                       tgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGGCTGTTTACTCTGCGTGTGCGTGAACTGCCTGCTTTTCGACTATATCGAATGAGTGCCGGTAg     >  1:333412‑M2/6‑143 (MQ=255)
                                                                                                                                          atCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTGTGCGTGCCCGTGgccg  >  1:97388‑M2/3‑142 (MQ=255)
                                                                                                                                            |                                                                                                                                            
ATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCA  >  NC_000913/4001509‑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.