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)44 (0.980) 36/284 0.2 95.7% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001649 = 2 (0.040)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAG  >  NC_000913/4001510‑4001790
                                                                                                                                           |                                                                                                                                             
atggattgagctggtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                            <  1:193594‑M2/4‑1 (MQ=255)
atcgattggtcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttgggcatagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTcc                                                                                                                                            <  2:471077‑M2/4‑1 (MQ=255)
 tggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacaacggcctgcatCTCct                                                                                                                                           >  2:182972‑M2/139‑143 (MQ=255)
   gattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccgggactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                         >  1:157846‑M2/137‑143 (MQ=255)
         tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                   >  1:716312‑M2/131‑143 (MQ=255)
         tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggaacgtttctttgaagacgacgacggcctgcagCTCcttcttcttc                                                                                                                                   >  2:131631‑M2/131‑143 (MQ=255)
           ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                                 <  1:424927‑M2/15‑1 (MQ=255)
            tgtcgaaccggacgacgacgagcgactgcgtgtacaaccttaacttttccatatcctggcaacccgcccggaattggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                <  2:270007‑M2/16‑1 (MQ=255)
            tgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                <  1:78312‑M2/16‑1 (MQ=255)
            tgccaacccggacgactacgagcggctgcgcgtacaagctgaacttggccagagcctgggaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                <  2:300131‑M2/16‑1 (MQ=255)
            gttcgctccgtgcgacgacgcgcgactgcgcctacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                <  2:754059‑M2/16‑1 (MQ=255)
             gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacggttctttgaagacgacgacggcatgcatCTCCTTCTTGTTGTTTg                                                                                                                               >  1:696530‑M2/127‑143 (MQ=255)
                    cggacgacgacgagcgaccgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGc                                                                                                                        <  1:558893‑M2/24‑1 (MQ=255)
                        cgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaa                                                                                                                    <  1:665974‑M2/28‑1 (MQ=255)
                         gacgacgagcgactgcgcgtacaatccgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTTTTTGAAGATGCGGAAg                                                                                                                   <  1:75368‑M2/29‑1 (MQ=255)
                           cgacgagcgcctgggcggacaatctgaacttggccagagcctggcaagccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAt                                                                                                                 <  2:148732‑M2/31‑1 (MQ=255)
                           cgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGAt                                                                                                                 <  2:638947‑M2/31‑1 (MQ=255)
                            gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcatgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                                >  2:174531‑M2/112‑143 (MQ=255)
                            gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgcttctttgaagacgacgccggcctgcatCTCCTTCTTCGTCTTCGAAGATGCGGAAGATc                                                                                                                >  2:453683‑M2/112‑143 (MQ=255)
                               gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgacggccgccatCACCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                             >  2:162006‑M2/109‑143 (MQ=255)
                                             cggagcgtGGCCAGAGCCTGGCATCCCGCCGGGAGCTGGACGTCGTCGGGGCTTCGGCACGTTTCTGTGAGGGCGACGACGGCCTGCAT‑‑‑‑‑CTCCTTCTTCTTCTTGGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGc                                                                                          <  2:384075/136‑1 (MQ=255)
                                                atctgaacttggtcagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                            >  1:663214‑M2/92‑143 (MQ=255)
                                                     aacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCAtt                                                                                       <  1:673866‑M2/57‑1 (MQ=255)
                                                       cttggccagagcctggcaacccgcccggaactggaagtcatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTa                                                                                     <  2:662642‑M2/59‑1 (MQ=255)
                                                           gccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCAt                                                                                 >  2:62761‑M2/81‑143 (MQ=255)
                                                             cagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCATCTTCTTTGAAGATGCGGAAGATGACGCGGGTAGCGCCACTGTGGCGTTTGCcatgc                                                                               >  1:303682‑M2/79‑141 (MQ=255)
                                                                  cctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGAt                                                                          >  1:17896‑M2/74‑143 (MQ=255)
                                                                   aaCCCGCCCGGAACTGGAAGCCATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGCACTG‑‑‑‑‑TGAGGGGGCTTTCTTCTTTGAAGTTGCGGAAGGTCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcgt                                                                    <  1:181056/143‑1 (MQ=255)
                                                                      gcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCCTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtc                                                                      <  1:704762‑M2/74‑1 (MQ=255)
                                                                          cccgcccggaactgcaagacctcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCt                                                                  <  2:529238‑M2/78‑1 (MQ=255)
                                                                            cgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGt                                                                <  2:112513‑M2/80‑1 (MQ=255)
                                                                                cggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCtgcttcc                                                            >  2:40230‑M2/60‑141 (MQ=255)
                                                                                           gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                                 <  1:464241‑M2/95‑1 (MQ=255)
                                                                                           gacagcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                                 <  1:19383‑M2/95‑1 (MQ=255)
                                                                                             catctaagcatcggcacgtttctctgaagacgagtacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgt                                               <  1:396111‑M2/97‑1 (MQ=255)
                                                                                              atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagcgtg                                              <  1:707696‑M2/98‑1 (MQ=255)
                                                                                                ggaagcatcggcacgttgctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTGGAAGATGCGGAAGAGCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGaa                                            <  1:247704‑M2/100‑1 (MQ=255)
                                                                                                         ggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCt                                   >  1:404365‑M2/35‑143 (MQ=255)
                                                                                                           cacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCttt                                 >  1:389261‑M2/33‑143 (MQ=255)
                                                                                                               tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGt                             >  1:135791‑M2/29‑143 (MQ=255)
                                                                                                                    ttgaagacgacgtcggcctgcatCTCCTTCTCCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGGTGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTa                        <  2:359560‑M2/120‑1 (MQ=255)
                                                                                                                           cgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGGACTGCCCGCTTTTCGTCTGTACCGtgtg                 >  1:635978‑M2/17‑143 (MQ=255)
                                                                                                                               gacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTg             <  1:368891‑M2/131‑1 (MQ=255)
                                                                                                                                acggcctgcatCTCCTTCCTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGc            <  1:131631‑M2/132‑1 (MQ=255)
                                                                                                                                  tccctgcatCTCCTTCTTCTTCTGTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCACCCGTGATGGTCGTCTGTTCACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGccc          <  1:62047‑M2/134‑1 (MQ=255)
                                                                                                                                          tCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCGGTGATGGTCGTCTGTTTACTCTGCGTGCGCGTGACCTGCCCGCTTTTGGTCTGTATCGTATCCGTCCCCGTagccta  >  2:349866‑M2/2‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                             
TTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAG  >  NC_000913/4001510‑4001790

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.