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.030)33 (0.890) 30/264 0.3 97.1% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 1 (0.030)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAG  >  NC_000913/4001508‑4001791
                                                                                                                                              |                                                                                                                                             
agtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                         <  6:135552‑M2/7‑1 (MQ=255)
agtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttc                                                                                                                                         <  4:38823‑M2/7‑1 (MQ=255)
           atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttc                                                                                                                                      <  1:117891‑M2/10‑1 (MQ=255)
           atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttc                                                                                                                                      >  2:117891‑M2/132‑141 (MQ=255)
                gtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                             <  8:136582‑M2/19‑1 (MQ=255)
                    aaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCgg                                                                                                                      >  6:66685‑M2/123‑148 (MQ=255)
                        ggacgacgacgcgcgactgcgcgtacaatctgaacttggccagagcctggccacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgccgccggcctgcatCTCCGTCTGCTTCTTTGATGATGCGGaagtt                                                                                                                 >  8:34159‑M2/119‑147 (MQ=255)
                        ggacgacgacgagcgactgcgcgtacaatctgatcttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGaagtt                                                                                                                 >  8:32188‑M2/119‑147 (MQ=255)
                               gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCg                                                                                                          >  5:72034‑M2/112‑149 (MQ=255)
                                           cgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtg                                                                                              >  4:38049‑M2/100‑149 (MQ=255)
                                                acaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCa                                                                                         >  6:27642‑M2/95‑149 (MQ=255)
                                                      tgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCGTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAc                                                                                    >  5:149303‑M2/89‑148 (MQ=255)
                                                        aacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAc                                                                                    >  7:48873‑M2/87‑146 (MQ=255)
                                                        aacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTAc                                                                                    <  8:48873‑M2/60‑1 (MQ=255)
                                                            agagCCAGGAAACCCGCCCGTAACTCTAAGAAATCTAAGAATCTGCACGTTTCTTAGAAGCCGACGAAGGACTTCAT‑‑‑‑‑CTCCTTCTTATTCATTTAAGATGCGGAAGATCACGCCGGTAACTCAACTGTGGCATTTACCATCCGTGATgg                                                                        <  7:32188/149‑1 (MQ=255)
                                                                 agagcctggcaacccgcccggacctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTGCTTCTTCTTTGAAGATGCGGAATATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGCGATgg                                                                        <  7:34159‑M2/72‑1 (MQ=255)
                                                                  gagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATgg                                                                        >  1:121930‑M2/77‑148 (MQ=255)
                                                                    gcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtc                                                                      <  2:121930‑M2/74‑1 (MQ=255)
                                                                          caacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                               <  5:27642‑M2/81‑1 (MQ=255)
                                                                              ccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTAct                                                           >  5:1533‑M2/65‑149 (MQ=255)
                                                                                    ggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                     >  2:110906‑M2/59‑149 (MQ=255)
                                                                                               acatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACt                                          <  5:66685‑M2/102‑1 (MQ=255)
                                                                                                 atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                                 <  8:28791‑M2/95‑1 (MQ=255)
                                                                                                 atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgagc                                                 >  7:28791‑M2/46‑140 (MQ=255)
                                                                                                   cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGccc                                      >  4:44833‑M2/44‑149 (MQ=255)
                                                                                                     aagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                    <  6:72034‑M2/108‑1 (MQ=255)
                                                                                                       gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                    >  1:143407‑M2/40‑147 (MQ=255)
                                                                                                       gcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGc                                    <  2:143407‑M2/108‑1 (MQ=255)
                                                                                                             ggacgtttctttgaagacgacgactgcctgcatCTCCTTCTTCTTATTTGATGATGCGGTAGATCACGCCGGTAACTCCACTTTTGCATTTACCATCCGTGATGTTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCACGCTTTTCGTc                            <  6:149303‑M2/116‑1 (MQ=255)
                                                                                                                  tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                         <  6:64986‑M2/71‑1 (MQ=255)
                                                                                                                  tttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATg                                                                         >  5:64986‑M2/29‑99 (MQ=255)
                                                                                                                       ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTAt                  >  1:95982‑M2/24‑149 (MQ=255)
                                                                                                                             acgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGc            <  3:44833‑M2/132‑1 (MQ=255)
                                                                                                                                       cctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAg  >  5:39345‑M2/8‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAG  >  NC_000913/4001508‑4001791

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ 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.